DDS不再孤單與Matlab“手拉手”一起創(chuàng )建仿真事業(yè)

clear all;clc %清除所有變量，清屏

M=2^11; %采樣的點(diǎn)數

a(1:1:M)=0; %頻率累加器的數值

b(1:1:M)=0; %相位累加器的數值

y(1:1:M)=0; %輸出的波形數據

K=400; %調頻斜率

Kc=10; %初始頻率控制字

N=10; %幅度量化位數

L=24; %相位累加器位數

%下面一段實(shí)現頻率累加器

a(1)=0.5*K; %初始頻率步進(jìn)量

for i=2:1:M

a(i)=a(i-1)+K;

end

%下面一段實(shí)現相位累加器

b(1)=Kc+a(1); %相位初始值

for i=2:1:M

b(i)=b(i-1)+(Kc+a(i));

end

%下面一段實(shí)現了查找表ROM以及進(jìn)行幅度量化

for i=1:1:M

y(i)=floor(2^N*cos(2*pi/(2^L)*b(i)));

end

%下面一段畫(huà)出相應的圖形

figure(1);plot(y);axis([0 M,-2^N-100 2^N+100]);

figure(2);freqz(y)

上面這段程序中，有很多的變量，包括調頻斜率K、頻率控制字Kc、幅度量化位數N以及相位累加器位數L等等，修改不同變量值可以得到各，不同的線(xiàn)性調頻信號。具體的實(shí)現要根據實(shí)際的需要來(lái)設置。4 實(shí)驗結果

根據上面的程序,取調頻斜率為400，頻率控制字為10，幅度量化為10位(和所使用的D/A配合)，相位累加器為24位，用Matlab仿真得到的線(xiàn)性調頻信號的波形和相應的幅頻響應如圖3和圖4所示。

圖4 針對圖3的幅頻特性曲線(xiàn)

結束語(yǔ)

本文所講的線(xiàn)性調頻信號的產(chǎn)生原理和方法有很好的可擴展性，修改頻率累加器的內容即可以實(shí)現其他的各種調頻信號。作為線(xiàn)性調頻信號，由于廣泛應用于高分辨率的雷達系統中，因此正確理解線(xiàn)性調頻信號的產(chǎn)生原理和掌握其產(chǎn)生的方法是很有現實(shí)意義的。在實(shí)際應用中主要采用專(zhuān)用DDS芯片或者FPGA來(lái)實(shí)現線(xiàn)性調頻信號(各有優(yōu)缺點(diǎn))，而本文則基于Matlab軟件的良好編程性來(lái)驗證這一方案，