梳状谱干扰matlab仿真,干扰信号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('噪声调相干扰频域波形');

谢 谢！

放映结束 感谢各位批评指导！

让我们共同进步