FPGA產(chǎn)生正弦波方法四分之一周期COE文件

1、-. z.FPGA產(chǎn)生正弦波方法1/4周期COE文件 本文主要介紹通過(guò)查表法產(chǎn)生正弦波。涉及利用matlab生成四分之一周期的正弦相位、幅值對(duì)應(yīng)表，生成COE文件的IP核方法，F(xiàn)PGA編程產(chǎn)生數(shù)字正弦信號(hào)。正弦信號(hào)，是一個(gè)模擬信號(hào)。FPGA只能產(chǎn)生數(shù)字信號(hào)。因此需要用DA將數(shù)字量轉(zhuǎn)化為模擬量。這里采用modelsim的模擬波形顯示，不表達(dá)DA轉(zhuǎn)換相關(guān)知識(shí)了。產(chǎn)生正弦信號(hào)的方法有很多，這里用的是查找rom的方法，產(chǎn)生正弦信號(hào)。正弦信號(hào)，是一個(gè)幅值介于-1和1之間周期變化的模擬量，因此只需要將1/4周期的值存進(jìn)rom里，其波形就可以根據(jù)這1/4周期幅值的變化得到。matlab產(chǎn)生正弦信號(hào)的值。以0

2、.01為步長(zhǎng)，從0采集到pi/2，共158個(gè)點(diǎn)。* = linspace(0, pi/2,158); 因sin的幅值值是小數(shù)，而FPGA不能表示小數(shù)，我們把幅值乘以12位二進(jìn)制后取整實(shí)現(xiàn)幅值整數(shù)化。即幅值放大4096倍。例如幅值為0.3，而0.3*212=1228.8，取整為1229。當(dāng)然這樣表示會(huì)有一定誤差，但可以根據(jù)精度需要調(diào)整參數(shù)，使誤差可承受。 生成ROM。這里rom用的是ISE的IP rom。將matlab生成的數(shù)據(jù)存入rom里面，讀取rom的值，即可生成正弦信號(hào)。首先是生成rom的初始化文件。*ilin*的rom的初始化文件的后綴為.coe。這個(gè)特殊文件有固定的格式如下： 文件前面

3、兩行的容是固定的。MEMORY_INITIALIZATION_RADI*=10;表示下面的數(shù)字是10進(jìn)制數(shù)，后面的數(shù)據(jù)是依次存入rom的值，以逗號(hào)分開(kāi)，最后以一個(gè)分號(hào)完畢。生成初始化文件的matlab程序如下所示：fid = fopen(sin_rom.t*t,w);fprintf(fid,MEMORY_INITIALIZATION_RADI* = 10;n);fprintf(fid,MEMORY_INITIALIZATION_VECTOR =n);for i = 0:1:pi/2*100y = sin(i/100);rom =floor( y * 212);if i = 157fprintf

4、(fid,%d;,rom);elsefprintf(fid,%d,rom);endif mod(i,10)=0 & i = 0fprintf(fid,n);endendfclose(fid);matlab生成的文件是.t*t文件。將后綴直接改為.coe即可。然后復(fù)制到ISE分工程目錄下。初始化文件生成后，剩下就是編寫(xiě)verilog代碼。首先建一個(gè)工程，然后新建一個(gè)IP。找到rom，翻開(kāi)。設(shè)置按如下設(shè)置：這里是設(shè)置位寬，采用的是13位第一位為符號(hào)位，后面12位為數(shù)據(jù)位來(lái)顯示sin的值。因此這里是設(shè)置為13.深度是因?yàn)橐?58個(gè)值，所以這里設(shè)置為158.將剛剛生成的.coe文件載入，生成即可。接

5、著就是編寫(xiě)verilog代碼了。首次查看用ip生成的rom的例化代碼.sin_rom your_instance_name (.clka(clka), / input clka.addra(addra), / input 7 : 0 addra.douta(douta) / output 12 : 0 douta);從上面的程序可看出，只需要給時(shí)鐘信號(hào)，和輸入地址，就可以了。輸出的就是正弦的數(shù)字信號(hào)了。分析正弦信號(hào)，前1/4個(gè)周期，地址從0自加，一直加到1571/4個(gè)周期的點(diǎn)數(shù)。然后再自減，減到0。然后進(jìn)入到負(fù)半周，負(fù)數(shù)的表示是以二進(jìn)制的補(bǔ)碼來(lái)表示的，即絕對(duì)值數(shù)的二進(jìn)制取反在加1。編寫(xiě)地址自加

6、自減的代碼，然后再根據(jù)地址的值，判斷輸出值的正負(fù)，如為負(fù)數(shù)，則取反加1后再輸出。代碼，如下所示：module sin_top(input clk, /輸入時(shí)鐘信號(hào)input 9:0 address, /輸入地址信號(hào)output reg 12:0 data_out /輸出sin的數(shù)字值);reg 7:0 add;wire 12:0 douta;/ 以下是判斷地址的值always * beginif( address = 157 )add = address;else if( address = 314 )add = 10d314 - address;else if( address = 471

7、)add = address - 10d314;else if( address = 628 )add = 10d628 - address;elseadd = 0;end/例化之前生成的sin_romsin_rom u1_sin_rom (.clka(clk), / input clka.addra(add), / input 7 : 0 addra.douta(douta) / output 12 : 0 douta);/ 判斷輸出值是正數(shù)還是負(fù)數(shù)。always * beginif( address = 314)data_out = douta;else if( address = 628

8、 )data_out = douta + 1b1;elsedata_out = 0;endendmodule程序?qū)懞昧?，剩下就是要仿真了。仿真的testbench，只需要輸入時(shí)鐘和地址信號(hào)就可以了。地址信號(hào)一直加一，直到不小于628，剛好一個(gè)周期完畢。就返回0值，在繼續(xù)自加。測(cè)試代碼如下所示：module sin_top_test;/ Inputsreg clk;reg 9:0 address;/ Outputswire 12:0 data_out;/ Instantiate the Unit Under Test (UUT)sin_top uut (.clk(clk),.address(address),.data_out(data_out);always#5 clk = clk;initial begin/ Initialize Inputsclk = 0;address = 0;/ Wait 100 ns for global reset to finishwhile(1)begin(negedge(clk);if( address 628 )address