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

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

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

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

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

K=400; %調(diào)頻斜率

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

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

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

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

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

for i=2:1:M

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

end

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

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

for i=2:1:M

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

end

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

for i=1:1:M

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

end

%下面一段畫出相應(yīng)的圖形

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

figure(2);freqz(y)

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

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

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

結(jié)束語

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