干扰信号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;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(-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);endxn=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);endxn=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(噪声调相干扰频域波形);,