扰Рfunction y=noiseFM(uj,mf,wpp); Р Рif nargin==0 Р uj=1; mf=0.6;Р wpp=6; РendРfj=35e6;fs=4*fj;Tr=520e-6; Рbj=10e6;Рt1=0:1/fs:3*Tr-1/fs; N=length(t1);Рu=wgn(1,N,wpp); Рdf1=fs/N;n=0:N/2;f=n*df1; Рwp=10e6;Рws=13e6; Рrp=1; rs=60; Р[Nn,wn]=buttord(wp/(30e6/2),ws/(30e6/2),rp,rs); Р[b,a]=butter(Nn,wn); Рu1=filter(b,a,u); РfigureРsubplot(2,2,1),plot(t1,u1),title('调制噪声波形');axis([0,0.05e-4,-6,6])Рsubplot(2,2,2),j2=fft(u1); plot(f,10*log10(abs(j2(n+1)*2/N)))Рtitle('调制噪声功率谱'); Рi=1:N-1;ss=cumsum([0 u1(i)])Рss=ss*Tr/N; Рy=uj*cos(2*pi*fj*t1+2*pi*mf*bj*ss+100);Рp=(1/N)*sum(y.^2)Рsubplot(2,2,3), plot(t1,y),title('噪声调频干扰时域波形')Рaxis([0,0.05e-4,-1.5,1.5])Рsubplot(2,2,4), J=fft(y); plot(f,(abs(J(n+1))))Рtitle('噪声调频干扰已调波功率谱')РР结果:РР精品范文模板 可修改删除РР免责声明:图文来源于网络搜集,版权归原作者所以Р若侵犯了您的合法权益,请作者与本上传人联系,我们将及时更正删除。