干扰信号MATLAB仿真

干扰信号MATLAB仿真实验目的•仿真产生4种以上的通信干扰信号。•画出信号的时域、频域波形。实验思路与步骤•产生一个高斯白噪声，•利用Matlab自带的fir1函数产生一个低通滤波器，限制高斯白噪声的带宽，由此产生了视频噪声。•利用产生的视频噪声，分别代入噪声调幅干扰的时域表达式，并且进行100次的积累后求平均值，对其进行快速傅里叶变换后，由此画出噪声调幅干扰频域波形•重复上述步骤，分别代入噪声调频干扰和噪声调相干扰的时域表达式，分别画出其时域波形和频域波形。实验结果•视频噪声噪声调幅干扰•噪声调频干扰•噪声调相干扰程序代码•clearall;closeall;clc;%清除变量•T=50e-6;%采样时间•fs=150e6;%采样频率•N=T*fs;%采样点数•detlf=10e6;%滤波器截止频率•f1=50e6;%调制信号中心频率•m=0.5;%调制度•kfm=5e6;%调频斜率•kpm=5;%调相斜率•M=100;%积累次数•p=fft(fir1(N-1,detlf/fs*2));%滤波器频谱s=0;fori=1:100xn=ifft(fft(random('Normal',0,1,1,N)).*p);%高斯白噪声通过滤波器j=abs(fft(xn));s=s+j;endfigure(1)t=0:1/fs:T-1/fs;plot(t*1e6,xn);xlabel('us');•title('视频噪声时域波形');•figure(2)•f=(0:N-1)*fs/N;•y1=fft(xn,N);•mag1=abs(y1);•plot(f*1e-6,y1);%视频噪声频域波形•axis([-215-150150]);•xlabel('MHZ');•title('视频噪声频域波形');n=1:N;zn=(1+m*cos(2*pi*xn)).*cos(2*pi*f1/fs*n);%噪声调幅干扰表达式figure(3)plot(t*1e6,zn);title('噪声调幅干扰时域波形');xlabel('us');figure(4)f=(0:N-1)*fs/N;y2=fft(zn,N);mag2=abs(y2);plot(f*1e-6,mag2);xlabel('MHZ');title('噪声调幅干扰频域波形');•sum(1)=0;•fori=1:N-1;•sum(i+1)=xn(i)+sum(i);•end•xn=sum/fs;•wn=cos((2*pi*f1*t+2*pi*kfm*xn));%噪声调频干扰表达式•figure(5)•plot(t*1e6,wn);•title('噪声调频干扰时域波形');•xlabel('us');•figure(6)•f=(0:N-1)*fs/N;•y3=fft(wn,N);•mag3=abs(y3);•plot(f*1e-6,mag3);•xlabel('MHZ');•title('噪声调频干扰频域波形');•sum(1)=0;•fori=1:N-1;•sum(i+1)=xn(i)+sum(i);•end•xn=sum/fs;•on=cos(2*pi*f1*t+kpm*xn);%噪声调相干扰表达式•figure(7)•plot(t*1e6,on);•title('噪声调相干扰时域波形');•xlabel('us');••figure(8)•f=(0:N-1)*fs/N;•y4=fft(on,N);•mag4=abs(y4);•plot(f*1e-6,mag4);•xlabel('MHZ');•title('噪声调相干扰频域波形');