一定要理解带通采样定理的周期延拓

clear all;
close  all;

%采样参数
T=9.005e-6;%采样时间
fs=200e+6;%采样率
n=round(T*fs);%采样点个数
t=linspace(0,T,n);
f=linspace(-fs/2,fs/2,n);
A_lfm=1;%幅度
B=20e+6;%带宽20hz

k=B/T;%调频斜率

f_lfm =120e+6;%初始频率102M

%NCO 数字混频本振频率
f2=120e+6;

y_lfm= A_lfm*(cos(2*pi*f_lfm*t + pi*k*t.^2));%LFM信号

%  Hd=low_pass_95M;
%  y_lfm=filter(Hd,y_lfm);

%%产生NCO IQ两路信号
 y_cos=cos(2*pi*f2*t);
 y_sin=sin(2*pi*f2*t);

%混频生成I路信号
for i=1:n    
 y_I(i)=y_lfm(i)*y_cos(i);
end

%混频生成Q路信号
for i=1:n    
 y_Q(i)=y_lfm(i)*y_sin(i);
end

%频谱
fft_lfm=fftshift(fft(y_lfm));%
fft_cos=fftshift(fft(y_cos));
fft_sin=fftshift(fft(y_sin));
fft_yI=fftshift(fft(y_I));
fft_yQ=fftshift(fft(y_Q));

%滤出I路
Hd=low_pass_35M;
y_after_I=filter(Hd,y_I);
y_after_Q=filter(Hd,y_Q);
fft_after_yI=fftshift(fft(y_after_I));
fft_after_yQ=fftshift(fft(y_after_Q));

%图1
% figure(1);
% plot(t,y_lfm);
% title('输入LFM信号时域');
% xlabel('t/s');
% ylabel('幅度');
% 
% %图2
% figure(2);
% plot(t,y_cos);
% title('NCO产生cos信号时域');
% xlabel('t/s');
% ylabel('幅度');
% 
% %图3
% figure(3);
% plot(t,y_sin);
% title('NCO产生sin信号时域');
% xlabel('t/s');
% ylabel('幅度');
% 
% %图4
% figure(4);
% plot(t,y_I);
% title('混频后I路的时域');
% xlabel('t/s');
% ylabel('幅度');
% 
% %图5
% figure(5);
% plot(t,y_Q);
% title('混频后Q路的时域');
% xlabel('t/s');
% ylabel('幅度');
% 
% %图6
% figure(6);
% plot(t,y_after_I);
% title('I路滤波后信号时域');
% xlabel('t/s');
% ylabel('幅度');
% 
% %图7
% figure(7);
% plot(t,y_after_Q);
% title('Q路滤波后信号时域');
% xlabel('t/s');
% ylabel('幅度');

%图8
figure(8);
plot(f,abs(fft_lfm));
title('输入信号LFM频谱');
xlabel('f/Hz');
ylabel('幅度');

%图9
figure(9);
plot(f,abs(fft_cos));
title('NCO信号cos频谱');
xlabel('f/Hz');
ylabel('幅度');

%图11
figure(10);
plot(f,abs(fft_sin));
title('NCO信号sin频谱');
xlabel('f/Hz');
ylabel('幅度');

%图12
figure(11);
plot(f,abs(fft_yI));
title('混频信号I路频谱');
xlabel('f/Hz');
ylabel('幅度');

%图13
figure(12);
plot(f,abs(fft_yQ));
title('混频信号Q路频谱');
xlabel('f/Hz');
ylabel('幅度');

%图14
figure(13);
plot(f,abs(fft_after_yI));
title('滤波后I路信号频谱');
xlabel('f/Hz');
ylabel('幅度');

%图15
figure(14);
plot(f,abs(fft_after_yQ));
title('滤波后Q路信号频谱');
xlabel('f/Hz');
ylabel('幅度');

% s=200e6;
% Tr=9.005e-6;             % 设脉宽 
% n=round(Tr*fs);      % 采样点个数
% t=linspace(0,Tr,n);  
% 混频信号
% f1=100e6;
% y_I=cos(2*pi*f1*t);
% y_Q=sin(2*pi*f1*t);
% 
%     for i=1:1
%     S_1=y_Q;
%     ff=linspace(-fs/2,fs/2,n);
%     对测试中频信号进行混频
%     fft_y=fftshift(fft(y_I));
%     figure;
%     plot(ff,abs(fft_y));
%     title('测试中频信号的频谱');
%     ylim([0,1000]);%对X轴设定显示范围 
%     xlim([-1.2e8,1.2e8]);%对X轴设定显示范围 
%     end