一维卷积噪声干扰

%%%%%%%%噪声卷积干扰%%%%%%%%%%%%%%%%%
clear all;close all;clc;

%%%%%%%初始化数据%%%%%%%%%
T=20e-6;
B=10e6;
f0=0;
Fs=5*B;
K=B/T;
Ts=1/Fs;
N=T*Fs;
%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%LFM信号%%%%%%%%%%%%%%%
t=linspace(0,T,N);
St=exp(j*pi*K*(t-T/2).^2);
figure,
subplot(121),plot(t*1e6,real(St));
xlabel('时间/us');
ylabel('幅度');
title('LFM时域信号');
grid on;

f=linspace(-Fs/2,Fs/2,N);
subplot(122),
plot(f*1e-6,fftshift(abs(fft(St))));
xlabel('频率/MHz');
ylabel('幅度');
grid on;
title('LFM信号频谱图');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%脉冲压缩%%%%%%%%%%%%%
St1=exp(j*pi*K*(t-10).^2);
y1=abs(xcorr(St,St));
M=length(y1);
axisX=(1:M)*Ts*1e6;
figure,
plot(axisX,abs(y1)/max(y1));
xlabel('时间/us');
ylabel('归一化幅度');
title('LFM脉冲压缩信号');
axis([10 30 0 1]);
grid on;
hold on;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%噪声卷积干扰%%%%%%%%%%%
Ls=2e-6;
M=fix(Ls*Fs);
Gnoise=wgn(1,M,0);      %均值为0的高斯视频噪声
J=conv(Gnoise,St);
so=J(1:length(St));
y3=abs(xcorr(so,St));
M=length(y3);
axisX=(1:M)*Ts*1e6;
figure,
subplot(121),
plot(axisX,real(y3)/max(y3));
xlabel('时间/us');
ylabel('归一化幅度');
title('(a)噪声时宽为2us');
axis([10 30 0 1]); 
%grid on;

%%%%%%%%%%%%%%%%%%%
Ls=4e-6;
M=fix(Ls*Fs);
Gnoise=wgn(1,M,0);      %均值为0的高斯视频噪声
J=conv(Gnoise,St);
so=J(1:length(St));
y3=abs(xcorr(so,St));
M=length(y3);
axisX=(1:M)*Ts*1e6;
subplot(122),
plot(axisX,real(y3)/max(y3));
xlabel('时间/us');
ylabel('归一化幅度');
title('（b)噪声时宽为4us');
axis([10 30 0 1]); 
%grid on;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%移频卷积噪声干扰%%%%%%%%
f1=1e6;       %移频量
f2=1.5e6;
St_j=exp(j*pi*K*(t-T/2).^2+j*2*pi*(0+f1)*(t-T/2));
St_j_1=exp(j*pi*K*(t-T/2).^2+j*2*pi*(0+f2)*(t-T/2));

J=conv(Gnoise,St_j);
so1=J(1:length(St))+St_j;
y3=abs(xcorr(so1,St));
M1=length(y3);
axisX=(1:M1)*Ts*1e6;
figure,
subplot(121),
plot(axisX,real(y3)/max(y3));
xlabel('时间/us');
ylabel('归一化幅度');
title('(a)移频量1MHz');
axis([10 30 0 1]); 

J=conv(Gnoise,St_j_1);
so1=J(1:length(St))+St_j_1;
y3=abs(xcorr(so1,St));
M1=length(y3);
axisX=(1:M1)*Ts*1e6;
subplot(122),
plot(axisX,real(y3)/max(y3));
xlabel('时间/us');
ylabel('归一化幅度');
title('(b)移频量1.5MHz');
axis([10 30 0 1]); 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%