信抗大作业汇编.docxVIP

信息对抗大作业 班级： 题目7、使用 Matlab 仿真产生多音干扰、噪声调频、梳状谱干扰的通信干扰信号。画出信号的时域、频域波形。 1 多音干扰信 仿真程序： 仿真结果如下： 2噪声调频干扰信号 噪声调频干扰信号的时域表达式为： 其中，调制噪声为零均值、广义平稳的随机过程，为均匀分布，且与相互独立的随机变量，为噪声调频信号的幅度，为噪声调频信号的中心频率，为调频斜率。 式中，为调幅函数的方差，其为 式中为的自相关函数，它可由调制噪声的功率谱变换求得。设其具有带限均匀谱，如下式所示： 则的功率谱为 式中，为调制噪声的谱宽，为有效调频指数，其中为有效调频带宽。 噪声调频信号功率谱的表达式为： 当时，由噪声调频信号功率谱表达式可得： 当时，由噪声调频信号功率谱表达式可得： 仿真程序： %噪声调频干扰 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(噪声调频干扰已调波功率谱) 仿真结果如下： 3梳状干扰 如果在某个频带内有多个离散的窄带干扰，形成多个窄带谱峰，则它们被称为梳状谱干扰。其信号表达式为 式中，Jn（t）为第n个窄带调幅干扰信号；An（t）为窄带干扰信号的包络，φn（t）为干扰信号的相位；ωn为干扰角频率。 仿真程序： n=randn(1,10e3)； J=0; fp=0.98*10e5; fs=0.99*10e5; Fs=3*10e5; fj=100*10e5; wp=2*pi*fp/Fs; ws=2*pi*fs/Fs; rs=40; B=ws-wp; beta=0.5842*(rs-21)^0.4+0.07886*(rs-21); N=ceil((rs-8)/2.285/B); wc=(wp+ws)/2/pi; h=fir1(N,wc,kaiser(N+1,beta)); H=fft(h); w=0:2*pi/length(H):2*pi-2*pi/length(H); s=conv(n,h); s=s/max(abs(s)); t=0:1/(2*fj):length(s)/(2*fj)-1/(2*fj); A=1.3*zeros(1,length(s)); dealt=50*10e5; for i=1:4 S=(A+s).*cos(2*pi*fj*t); fj=fj+dealt; J=J+S; end f=[0:length(J)-1]*2*fj/length(J); figure(1) plot(t,J) title(时域) xlabel(t),ylabel(J) figure(2) plot(f,abs(fft(J))) title(频域) xlabel(f),ylabel(H) 仿真结果如下： 题目6、使用 Matlab 进行两种不同数字调制信号特征值的计算比较。 1、2ASK信号仿真 对于二进制情况2ASK信号可表示为： 其中m(t)是单极性数字基带信号，取值为0和1，码元宽度为，码元速率是，是载波频率。m(t)是复包络瞬时幅度a(t)和瞬时相位φ（t）分别为： 2ASK信号瞬时参数仿