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

clearall;closeall;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-1detlf/fs*2));%滤波器频谱s=0;fori=1:100xn=ifft(fft(random('Normal'011N)).*p);%高斯白噪声通过滤波器j=abs(fft(xn));ends=s/M;j=s;figure(1)t=0:1/fs:T-1/fs;plot(t*1e6xn);xlabel('us');title('视频噪声时域波形');figure(2)plot(f*1e-620*log10(j.^2/max(j.^2)));%视频噪声功率谱axis([-122-80]);xlabel('MHZ');title('视频噪声功率谱')n=1:N;zn=(1+m*cos(2*pi*xn)).*cos(2*pi*f1/fs*n);%噪声调幅干扰表达式figure(3)plot(t*1e6zn);title('噪声调幅干扰时域波形');xlabel('us');s=0;fori=1:100zn=(1+m*cos(2*pi*xn)).*cos(2*pi*f1/fs*n);j=abs(fft(zn));s=s+j;ends=s/M;j=s;figure(4)plot(f*1e-620*log10(j.^2/max(j.^2)));%噪声调幅干扰功率谱title('噪声调幅干扰功率谱');xlabel('MHZ');axis([90110-2000]);sum(1)=0;fori=1:N-1sum(i+1)=xn(i)+sum(i);endwn=cos((2*pi*f1*t+2*pi*kfm*xn));%噪声调频干扰表达式figure(5)plot(t*1e6wn);title('噪声调频干扰时域波形');xlabel('us');s=0;fori=1:100xn=ifft(fft(random('Normal'011N)).*p);sum(1)=0;fori=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-620*log10(j.^2/max(j.^2)));%噪声调频干扰功率谱axis([50150-1500])xlabel('MHZ');title('噪声调频干扰功率谱');