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

《基于MATLAB的噪声仿真程序.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、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-1,detlf/fs*2));%滤波器频谱s=0;fori=1:100xn=ifft(fft(random('Normal',0,1,1,N)).*p);%高斯白噪声通过滤波器j=abs(fft(xn));ends=s/M;j=s;figu

2、re(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([-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*1e6,zn);title('噪声调幅干扰时域波形');xlabel('us');

3、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-6,20*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));%噪声调频干扰表

4、达式figure(5)plot(t*1e6,wn);title('噪声调频干扰时域波形');xlabel('us');s=0;fori=1:100xn=ifft(fft(random('Normal',0,1,1,N)).*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-6,20*log10(j.^2/max(j.^2)));%噪声调频干扰功率谱axis(

5、[50150-1500])xlabel('MHZ');title('噪声调频干扰功率谱');