干扰信号MATLAB仿真ppt课件

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