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

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

M=2^11; %采樣的點數(shù)

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

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

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

K=400; %調頻斜率

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

N=10; %幅度量化位數(shù)

L=24; %相位累加器位數(shù)

%下面一段實現(xiàn)頻率累加器

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

for i=2:1:M

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

end

%下面一段實現(xiàn)相位累加器

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

for i=2:1:M

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

end

%下面一段實現(xiàn)了查找表ROM以及進行幅度量化

for i=1:1:M

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

end

%下面一段畫出相應的圖形

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

figure(2);freqz(y)

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

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

圖4 針對圖3的幅頻特性曲線

結束語

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