zjietiao.m

用MATLAB仿真雷达的函数

%产生线形调频信号
t=100e-6;   %线形调频信号时间长度，单位：s
fs=40e6;   %采样频率，单位：Hz
fc=9e6;   %线形调频的起始频率，单位：Hz
f0=10e6;   %本振信号的中心频率，单位：Hz 
B=2e6;   %线形调频信号的调频带宽，单位：Hz
ft=0:1/fs:t-1/fs;   %调制信号，对于线形调频来说，调制信号就是时间序列
N=length(ft);
k=B/fs*2*pi/max(ft);   %计算调制灵敏度，也就是理解为线形调频的步进系数
y=modulate(ft,fc,fs,'fm',k);
y_fft_result=fft(y);
figure,subplot(2,1,1),plot(ft,y),
xlabel('t（单位：秒）'),ylabel('y（单位：伏）'),title('线形调频信号y(t)');
subplot(2,1,2),plot((0:fs/N:fs/2-fs/N),abs(y_fft_result(1:N/2))),
xlabel('频率f(单位：Hz）');title('线形调频信号y(t)的频谱');
%%%%%%%%%%%%正交解调%%%%%%%%%
n=0:N-1;
local_oscillator_i=cos(n*f0/fs*2*pi);  %I路本振信号
local_oscillator_q=sin(n*f0/fs*2*pi);  %q路本振信号
fbb_i=local_oscillator_i.*y;   %I路解调
fbb_q=local_oscillator_q.*y;   %q路解调
window=chebwin(51,40);   %采用50阶cheby窗的FIR低通滤波器
[b,a]=fir1(50,2*B/fs,window);
fbb_i=filter(b,a,fbb_i);   %I路信号经过低通滤波器
fbb_q=filter(b,a,fbb_q);   %Q路信号经过低通滤波器
figure,subplot(2,1,1),plot(ft,fbb_i),xlabel('t（单位：秒）');title('解调后的I路信号');
subplot(2,1,2),plot(ft,fbb_q),xlabel('t（单位：秒）');title('解调后的Q路信号');
fbb=fbb_i+j*fbb_q;
fbb_fft_result=fft(fbb);
figure,plot((0:fs/N:fs/2-fs/N),abs(fbb_fft_result(1:N/2))),
xlabel('频率f（单位：Hz）');title('解调后的信号的频谱');