matlab仿真噪声,基于MATLAB的噪声仿真程序.docx

clear all;close all;clc;%清除变量

T=100e-6;%采样时间

fs=300e6;%采样频率

N=T*fs;%采样点数

detlf=20e6;%滤波器截止频率

f1=100e6;%调制信号中心频率

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));

end

s=s/M;

j=s;

figure(1)

t=0:1/fs:T-1/fs;

plot(t*1e6,xn);

xlabel('us');

title('视频噪声时域波形');

figure(2)

plot(f*1e-6,20*log10(j.^2/max(j.^2)));%视频噪声功率谱

axis([-1 22 -8 0]);

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');

s=0;

for i=1:100

zn=(1+m*cos(2*pi*xn)).*cos(2*pi*f1/fs*n);

j=abs(fft(zn));

s=s+j;

end

s=s/M;

j=s;

figure(4)

plot(f*1e-6,20*log10(j.^2/max(j.^2)));%噪声调幅干扰功率谱

title('噪声调幅干扰功率谱');

xlabel('MHZ');

axis([90 110 -200 0]);

sum(1)=0;

for i=1:N-1

sum(i+1)=xn(i)+sum(i);

end

wn=cos((2*pi*f1*t+2*pi*kfm*xn));%噪声调频干扰表达式

figure(5)

plot(t*1e6,wn);

title('噪声调频干扰时域波形');

xlabel('us');

s=0;

for i=1:100

xn=ifft(fft(random('Normal',0,1,1,N)).*p);

sum(1)=0;

for i=1:N-1;

sum(i+1)=xn(i)+sum(i);

wn=cos((2*pi*f1*t+2*pi*kfm*xn));

j=abs(fft(wn));

s=s+j;

s=s/M;

j=s;

figure(6)

plot(f*1e-6,20*log10(j.^2/max(j.^2)));%噪声调频干扰功率谱

axis([50 150 -150 0])

xlabel('MHZ');

title('噪声调频干扰功率谱');