噪声调幅与调频干扰信号仿真分析.doc

噪声调幅与调频干扰信号仿真分析一、噪声调幅干扰信号时域表达式和功率谱仿真分析噪声调幅干扰信号的时域表达式为：其中，调制噪声为零均值，方差为，在区间分布的广义平稳随机过程，为均匀分布，且为与独立的随机变量，，为常数噪声调幅定理：式中，为调制噪声的相关函数噪声调幅信号的总功率为：它等于载波功率()与调制噪声功率()一半的和其又可改写为：式中，，为载波功率；，为有效调制系数噪声调幅信号的功率谱可由噪声调幅定理经傅立叶变换求得：式中，为调制噪声的功率谱，第一项代表载波的功率谱，后两项代表调制噪声功率谱的对称平移用MATLAB仿真分析：程序：%噪声调幅干扰function y=noiseAM(u0,N,wpp); if nargin==0 wpp=0;u0=1; end fj=35e6;fs=4*fj; Tr=520e-6; 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=14e6; rp=1; rs=60; [n1,wn1]=buttord(wp/(fs/2),ws/(fs/2),rp,rs); [b,a]=butter(n1,wn1); u1=filter(b,a,u); p=0.1503*mean((u1.^2)) ;figuresubplot(2,2,1),plot(t1,u1),title('噪声调制波形'); axis([0,0.05e-4,-2,2])subplot(2,2,2), j2=fft(u1);plot(f,10*log10(abs(j2(n+1)*2/N)))title('调制噪声功率谱');rand('state', 0); y=(u0+u1).*cos(2*pi*fj*t1+2); p=(1/N)*sum(y.^2);subplot(2,2,3), plot(t1,y),title('噪声调幅干扰时域波形'); axis([0,0.05e-4,-2,2])subplot(2,2,4), J=fft(y);plot(f,10*log10(abs(J(n+1))))title('已调波功率谱');结果：二、噪声调幅干扰信号时域表达式和功率谱仿真分析噪声调频干扰信号的时域表达式为：其中，调制噪声为零均值、广义平稳的随机过程，为均匀分布，且与相互独立的随机变量，为噪声调频信号的幅度，为噪声调频信号的中心频率，为调频斜率。

式中，为调幅函数的方差，其为式中为的自相关函数，它可由调制噪声的功率谱变换求得设其具有带限均匀谱，如下式所示：则的功率谱为 式中，为调制噪声的谱宽，为有效调频指数，其中为有效调频带宽噪声调频信号功率谱的表达式为：当时，由噪声调频信号功率谱表达式可得：当时，由噪声调频信号功率谱表达式可得：用MATLAB仿真分析：程序：%噪声调频干扰function y=noiseFM(uj,mf,wpp); if nargin==0 uj=1; mf=0.6; wpp=6; endfj=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); figuresubplot(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('噪声调频干扰已调波功率谱')结果：。