干扰信号MATLAB仿真

• 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');
• • • • • • • • • •
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');
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('噪声调幅干扰频域波形');
Biblioteka Baidu 实验结果
• 视频噪声
噪声调幅干扰
• 噪声调频干扰
• 噪声调相干扰
程序代码
• • • • • • • • • • • 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;%积累次数
• 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('视频噪声频域波形');
干扰信号MATLAB仿真
实验目的
• 仿真产生4种以上的通信干扰信号。 • 画出信号的时域、频域波形。
实验思路与步骤
• 产生一个高斯白噪声， • 利用Matlab自带的fir1函数产生一个低通滤波器，限制高斯白噪声 的带宽，由此产生了视频噪声。 • 利用产生的视频噪声，分别代入噪声调幅干扰的时域表达式，并 且进行100次的积累后求平均值，对其进行快速傅里叶变换后， 由此画出噪声调幅干扰频域波形 • 重复上述步骤，分别代入噪声调频干扰和噪声调相干扰的时域表 达式，分别画出其时域波形和频域波形。
• figure(8) • f=(0:N-1)*fs/N; • y4=fft(on,N); • mag4=abs(y4); • plot(f*1e-6,mag4); • xlabel('MHZ'); • title('噪声调相干扰频域波形');
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');