干扰信号MATLAB仿真.ppt

干扰信号MATLAB仿真 实验目的 仿真产生4种以上的通信干扰信号。 画出信号的时域、频域波形。 实验思路与步骤   产生一个高斯白噪声， 利用Matlab自带的fir1函数产生一个低通滤波器，限制高斯白噪声的带宽，由此产生了视频噪声。 利用产生的视频噪声，分别代入噪声调幅干扰的时域表达式，并且进行100次的积累后求平均值，对其进行快速傅里叶变换后，由此画出噪声调幅干扰频域波形 重复上述步骤，分别代入噪声调频干扰和噪声调相干扰的时域表达式，分别画出其时域波形和频域波形。 实验结果 视频噪声 噪声调幅干扰 噪声调频干扰 噪声调相干扰 程序代码 clear all;close all;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; for i=1:100 xn=ifft(fft(random('Normal',0,1,1,N)).*p);%高斯白噪声通过滤波器 j=abs(fft(xn)); s=s+j; end figure(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([-2 15 -150 150]); 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; for i=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; for i=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('噪声调相干扰频域波形'); 谢 谢！ 放映结束 感谢各位批评指导！ 让我们共同进步