iq.m

雷达数字接收机正交解调的matlab仿真

%%产生调频信号，中心频率f0=20MHz,带宽2MHz,采样率80MHz

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(s)'),ylabel('y(V)'),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('线性调频信号频谱');



%%%%%%%%%正交解调%%%%%%%%%%%%%%%

n=0:N-1;
s1=cos(n*f0/fs*2*pi);       % I路本振
s2=sin(n*f0/fs*2*pi);       % Q路本振
sI=s1.*y;
sQ=s2.*y;
window=chebwin(51,40);      % 50阶cheby 滤波器
[b,a]=fir1(50,2*B/fs,window);
sI=filter(b,a,sI);
sQ=filter(b,a,sQ);
%%%%%%%%%%%%%

figure(2);
subplot(3,1,1),plot(ft,sI),title('解调后I路信号');% xlabel('t');
subplot(3,1,2),plot(ft,sQ),title('解调后Q路信号');% xlabel('t');
sig=sI+j*sQ;
sig_fft=fft(sig);
subplot(3,1,3),plot((0:fs/N:fs/2-fs/N),abs(sig_fft(1:N/2))),xlabel('频率f(Hz)'),title('解调信号频谱')