通信對抗原理大作業(yè)-資料下載頁
2025-01-07 03:20本頁面
【正文】 )*sum(y.^2) subplot(2,2,3), plot(t1,y),title(39。噪聲調(diào)頻干擾時域波形 39。) axis([0,]) subplot(2,2,4), J=fft(y)。 plot(f,(abs(J(n+1)))) title(39。噪聲調(diào)頻干擾已調(diào)波功率譜 39。) 仿真結(jié)果如下： 0 1 2 3 4 5x 1 06505調(diào)制噪聲波形0 2 4 6 8x 1 07 8 0 6 0 4 0 2 00調(diào)制噪聲功率譜0 1 2 3 4 5x 1 06101噪聲調(diào)頻干擾時域波形0 2 4 6 8x 1 0702022400060008000噪聲調(diào)頻干擾已調(diào)波功率譜 3 梳狀干擾 如果在某個頻帶內(nèi)有多個離散的窄帶干擾，形成多個窄帶譜峰，則它們被稱為梳狀譜干擾。其信號表達式為 ? ?? ? ??? Ln Ln nnnn tttAtJtJ 1 1 ))(cos()()()( ?? 式中， Jn（ t）為第 n 個窄帶調(diào)幅干擾信號； An（ t）為窄帶干擾信號的包絡(luò)，φ n（ t）為干擾信號的相位；ω n 為干擾角頻率。 仿真程 序： n=randn(1,10e3)； J=0。 fp=*10e5。 fs=*10e5。 Fs=3*10e5。 fj=100*10e5。 wp=2*pi*fp/Fs。 ws=2*pi*fs/Fs。 rs=40。 B=wswp。 beta=*(rs21)^+*(rs21)。 N=ceil((rs8)/)。 wc=(wp+ws)/2/pi。 h=fir1(N,wc,kaiser(N+1,beta))。 H=fft(h)。 w=0:2*pi/length(H):2*pi2*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=*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(39。時域 39。) xlabel(39。t39。),ylabel(39。J39。) figure(2) plot(f,abs(fft(J))) title(39。頻域 39。) xlabel(39。f39。),ylabel(39。H39。) 仿真結(jié)果如下： 0 1 2 3 4 5 6x 1 0542024時域tJ0 1 2 3 4 5 6x 1 08050100150頻域fH 4 錐形音頻干擾 前面介紹的幾種信號屬于壓制式干擾，而錐形音頻干擾屬于欺騙式干擾，它由高斯噪聲與錐形脈沖相乘產(chǎn)生，該信號的表達式為)()()( tutdUtJ nj? 仿真程序 ： n=randn(1,10e3)。 fp=*10e5。 fs=*10e5。 Fs=3*10e5。 A=1。 Ta=2。 fj=100*10e5。 wp=2*pi*fp/Fs。 ws=2*pi*fs/Fs。 rs=40。 B=wswp。 beta=*(rs21)^+*(rs21)。 N=ceil((rs8)/)。 wc=(wp+ws)/2/pi。 h=fir1(N,wc,kaiser(N+1,beta))。 H=fft(h)。 w=0:2*pi/length(H):2*pi2*pi/length(H)。 s=conv(n,h)。 s=s/max(abs(s))。 t=0:1/(2*fj):length(s)/(2*fj)1/(2*fj)。 k=2。 d=2*pi*k*t。 J=A*d.*s。 f=[0:length(J)1]*2*fj/length(J)。 figure(1) plot(t,J) title(39。時域 39。) xlabel(39。t39。),ylabel(39。J39。) figure(2) plot(f,abs(fft(J))) title(39。頻域 39。) xlabel(39。f39。),ylabel(39。H39。) 仿真結(jié)果如下： 0 1 2 3 4 5 6x 1 05505x 1 04時域tJ0 0 . 5 1 1 . 5 2x 1 0800 . 0 10 . 0 20 . 0 30 . 0 4頻域fH
公司管理相關(guān)推薦
鄂ICP備17016276號-1