《通信對抗原理》大作業(yè)(第四組)

1、通信對抗原理大作業(yè)第四組第8、9題組長：呂才兵 02113058 組員：林家泉 02113057 韓樹榮 02113001 賈文凱 02113003 陳 員 02113038第8題：干擾信號產(chǎn)生。仿真產(chǎn)生4種以上的通信干擾信號。畫出信號的時(shí)域、頻域波形。本文主要介紹通過matlab產(chǎn)生噪聲調(diào)幅、噪聲調(diào)頻、梳狀干擾及錐形音頻干擾等四種干擾信號。下面將分別介紹這四種干擾信號。1 噪聲調(diào)幅干擾信號噪聲調(diào)幅干擾信號的時(shí)域表達(dá)式為：其中，調(diào)制噪聲為零均值，方差為，在區(qū)間分布的廣義平穩(wěn)隨機(jī)過程，為均勻分布，且為與獨(dú)立的隨機(jī)變量，為常數(shù)。噪聲調(diào)幅定理：式中，為調(diào)制噪聲的相關(guān)函數(shù)。噪聲調(diào)幅信號的總功率為：它等

2、于載波功率與調(diào)制噪聲功率的一半的和。其又可改寫為：式中，為載波功率；，為有效調(diào)制系數(shù)。噪聲調(diào)幅信號的功率譜可由噪聲調(diào)幅定理經(jīng)傅立葉變換求得：式中，為調(diào)制噪聲的功率譜，第一項(xiàng)代表載波的功率譜，后兩項(xiàng)代表調(diào)制噪聲功率譜的對稱平移。仿真程序：%噪聲調(diào)幅干擾function y=ZSAM(u0,N,wpp); %函數(shù)定義行 if nargin=0 %判斷輸入變量為零 wpp=0;u0=1; end fj=35e6;fs=4*fj; Tr=520e-6; %載波頻率，采樣頻率，總時(shí)間 t1=0:1/fs:3*Tr-1/fs; N=length(t1);%采樣時(shí)間，采樣點(diǎn)數(shù) u=wgn(1,N,wpp);

3、 %產(chǎn)生高斯白噪聲，wpp為以dbw為單位指定輸出噪聲的強(qiáng)度df1=fs/N;n=0:N/2;f=n*df1; wp=10e6; %通帶截止頻率ws=14e6; %阻帶截止頻率rp=1; rs=60; %通帶最大衰減系數(shù)，阻帶最小衰減系數(shù)n1,wn1=buttord(wp/(fs/2),ws/(fs/2),rp,rs); %用于計(jì)算巴特沃斯數(shù)字濾波器的階數(shù)N和3dB截止頻率wc。b,a=butter(n1,wn1); % 計(jì)算N階巴特沃斯數(shù)字濾波器系統(tǒng)函數(shù)分子、分母多項(xiàng)式的系數(shù)向量b、a。u1=filter(b,a,u); %Y = filter(B,A,X) ，輸入X為濾波前序列，Y為濾波結(jié)

4、果序列，B/A 提供濾波器系數(shù)，B為分子， A為分母 figure%畫圖subplot(4,1,1),plot(t1,u1),title(噪聲調(diào)制波形);xlabel(時(shí)間),ylabel(幅度), axis(0,0.06e-4,-2,2)subplot(4,1,3), j2=fft(u1);plot(f,10*log10(abs(j2(n+1)*2/N),xlabel(頻率),ylabel(功率)title(調(diào)制噪聲功率譜);y=(u0+u1).*cos(2*pi*fj*t1+2); %噪聲調(diào)幅信號表達(dá)式 p=(1/N)*sum(y.2);subplot(4,1,2), plot(t1,y)

5、,title(噪聲調(diào)幅干擾時(shí)域波形);xlabel(時(shí)間),ylabel(幅度), axis(0,0.05e-4,-2,2)subplot(4,1,4), J=fft(y);plot(f,10*log10(abs(J(n+1),xlabel(頻率),ylabel(功率)title(已調(diào)波功率譜);仿真結(jié)果如圖8.1所示：圖8.1 噪聲調(diào)幅信號仿真結(jié)果2 噪聲調(diào)頻干擾信號噪聲調(diào)頻干擾信號的時(shí)域表達(dá)式為：其中，調(diào)制噪聲為零均值、廣義平穩(wěn)的隨機(jī)過程，為均勻分布，且與相互獨(dú)立的隨機(jī)變量，為噪聲調(diào)頻信號的幅度，為噪聲調(diào)頻信號的中心頻率，為調(diào)頻斜率。式中，為調(diào)幅函數(shù)的方差，其為式中為的自相關(guān)函數(shù)，它可由調(diào)

6、制噪聲的功率譜變換求得。設(shè)其具有帶限均勻譜，如下式所示：則的功率譜式中，為調(diào)制噪聲的譜寬，為有效調(diào)頻指數(shù)，其中為有效調(diào)頻帶寬。噪聲調(diào)頻信號功率譜的表達(dá)式為：當(dāng)時(shí)，由噪聲調(diào)頻信號功率譜表達(dá)式可得：當(dāng)時(shí)，由噪聲調(diào)頻信號功率譜表達(dá)式可得：仿真程序：%噪聲調(diào)頻干擾function y=ZSFM(uj,mf,wpp); if nargin=0 %輸入?yún)?shù)個(gè)數(shù) uj=1; mf=0.6;%調(diào)頻系數(shù)mf=Uw比Uc wpp=6; %以dBW為單位的指定輸出噪聲的強(qiáng)度endfj=35e6;fs=4*fj;Tr=520e-6; %載波頻率和采樣頻率以及搜索時(shí)間，Tr為采樣的總時(shí)間bj=10e6;t1=0:1/f

7、s:3*Tr-1/fs; N=length(t1);%時(shí)間的取值范圍及步長，求其長度u=wgn(1,N,wpp);%高斯白噪聲 wpp輸出強(qiáng)度，N為采樣個(gè)數(shù)df1=fs/N;n=0:N/2;f=n*df1; wp=10e6;%通帶截止頻率ws=13e6;%阻帶截止頻率rp=2; rs=50;%通帶最大衰減系數(shù)為2，阻帶最小衰減系數(shù)為50 Nn,wn=buttord(wp/(30e6/2),ws/(30e6/2),rp,rs); %Nn為濾波器的階數(shù)，wn為3dB截止頻率即功率下降一半是對應(yīng)的頻率b,a=butter(Nn,wn);% 計(jì)算濾波器系統(tǒng)函數(shù)分子分母多項(xiàng)式系數(shù)B和Au1=filter

8、(b,a,u);%一階濾波，ul為濾波結(jié)果序列，b為輸入x的系數(shù)，a表示輸出y的系數(shù)，u為濾波前的輸入序列 figuresubplot(4,1,1),plot(t1,u1,b-),title(調(diào)制噪聲波形);axis(0,0.06e-4,-6,6),xlabel(時(shí)間),ylabel(幅度)subplot(4,1,3),j2=fft(u1);plot(f,10*log10(abs(j2(n+1)*2/N),b-),xlabel(頻率),ylabel(功率)% 一個(gè)傅里葉變換，由時(shí)域轉(zhuǎn)為頻域 title(調(diào)制噪聲功率譜); i=1:N-1;ss=cumsum(0 u1(i)%返回由標(biāo)量u（i）指

9、定的維數(shù)累家和，分別從u（1）u（N-1）ss=ss*Tr/N; y=uj*cos(2*pi*fj*t1+2*pi*mf*bj*ss+100);%函數(shù)時(shí)域表達(dá)式p=(1/N)*sum(y.2)%函數(shù)的功率表達(dá)式subplot(4,1,2), plot(t1,y),xlabel(時(shí)間),ylabel(幅度),title(噪聲調(diào)頻干擾時(shí)域波形)axis(0,0.06e-4,-1.5,1.5)%x的范圍06，y的范圍-1.51.5subplot(4,1,4), J=fft(y); plot(f,(abs(J(n+1),xlabel(頻率),ylabel(功率)title(噪聲調(diào)頻干擾已調(diào)波功率譜)仿

10、真結(jié)果如圖8.2所示：圖8.2 噪聲調(diào)頻干擾信號仿真結(jié)果3 音頻干擾信號音頻干擾是指使用單個(gè)正弦波或者多個(gè)正弦波的干擾信號。當(dāng)只使用單個(gè)正弦波時(shí)，稱為單音干擾，當(dāng)使用多個(gè)正弦信號時(shí)，稱為多音干擾。1） 單音干擾單音干擾就是干擾信號只發(fā)射一個(gè)正弦波，因此，它是一個(gè)單頻連續(xù)波。單音干擾也叫點(diǎn)頻干擾。單音干擾的時(shí)域表達(dá)式為式中，是干擾角頻率，通常它與被干擾信號的載波角頻率相同。單音信號的功率譜是在干擾頻率處的單根譜線。其時(shí)域波形和頻譜如圖8.3-1所示。圖8.3-1 單音干擾時(shí)域波形和頻譜仿真程序如下：%單音干擾t=0:0.001:0.6;y=sin(2*pi*50*t);%y=x+1.5*rand

11、n(1,length(t);subplot(2,1,1);plot(t,y);title(單音干擾的信號波形);ylabel(幅度);xlabel(時(shí)間);Y=fft(y,512);P=Y.*conj(Y)/512;%計(jì)算信號功率譜密度f=1000*(0:255)/512;%計(jì)算橫軸頻率值subplot(2,1,2);plot(f,P(1:256)title(單音干擾的功率譜);ylabel(功率譜);xlabel(頻率/Hz);2) 多音干擾干擾機(jī)可以發(fā)射L1個(gè)正弦信號，這些音頻可以隨機(jī)分布，或者位于特定的頻率上。多音干擾是由L個(gè)獨(dú)立的正弦波信號疊加而產(chǎn)生的，其時(shí)域表達(dá)式為多音干擾的功率譜為

12、多根等間隔的譜線。獨(dú)立多音干擾的時(shí)域和頻域波形如圖8.3-2所示。圖8.3-2 多音干擾時(shí)域波形和頻譜仿真程序如下：%多音干擾Ujn=1;fc=50;Phase=100;fs=1000;t=0:1/fs:1;y=Ujn*sin(2*pi*fc*t+Phase);%+0.1*randn(1,length(t);for i=1:10 Ujn=Ujn+0.1; fc=fc+50; Phase=Phase+100; y=y+Ujn*sin(2*pi*fc*t+Phase);endsubplot(2,1,1);plot(t,y);title(多音干擾的信號波形);ylabel(幅度);xlabel(時(shí)間

13、); Y=fft(y,512);P=Y.*conj(Y)/512;%計(jì)算信號功率譜密度f=10*(0:255)/512;%計(jì)算橫軸頻率值subplot(2,1,2);plot(f,P(1:256);title(多音干擾的功率譜);%subplot(3,1,3);%plot(f,20*log10(P(1:256);%title(多音干擾的功率譜);ylabel(功率譜);xlabel(頻率/Hz);4 梳狀譜干擾信號如果在某個(gè)頻帶內(nèi)有多個(gè)離散的窄帶干擾，形成多個(gè)窄帶譜峰，則它們被稱為梳狀譜干擾。其信號表達(dá)式為式中，為第個(gè)窄帶調(diào)幅干擾信號；為窄帶干擾信號的包絡(luò)，為干擾信號的相位；為干擾角頻率。仿真

14、程序：%梳狀干擾n=randn(1,10e3); %產(chǎn)生標(biāo)準(zhǔn)正態(tài)分布的隨機(jī)數(shù)J=0; %第一個(gè)干擾信號的功率為零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=45; %給定的一個(gè)參數(shù)B=ws-wp; %帶寬beta=0.5842*(rs-21)0.4+0.07886*(rs-21); %窗函數(shù)的beta參數(shù)N=ceil(rs-8)/2.285/B); %濾波器階數(shù)wc=(wp+ws)/2/pi; %截

15、止頻率h=fir1(N,wc,kaiser(N+1,beta); %截止頻率為wc的N階凱瑟窗函數(shù) H=fft(h); %對凱瑟窗函數(shù)進(jìn)行傅里葉變換w=0:2*pi/length(H):2*pi-2*pi/length(H); s=conv(n,h); %對凱瑟窗函數(shù)進(jìn)行卷積運(yùn)算s=s/max(abs(s); t=0:1/(2*fj):length(s)/(2*fj)-1/(2*fj); %產(chǎn)生信號的時(shí)間長度與時(shí)間間隔A=1.3*zeros(1,length(s); dealt=50*10e5; %頻率間隔for i=1:100 %循環(huán)100次S=(A+s).*cos(2*pi*fj*t);

16、%干擾信號fj=fj+dealt; %干擾信號的頻率遞增J=J+S; %干擾信號累加end f=0:length(J)-1*2*fj/length(J);%figure(1)subplot(2,1,1)plot(t,J) title(時(shí)域) xlabel(t),ylabel(幅度) subplot(2,1,2)plot(f,abs(fft(J) title(頻域) xlabel(f),ylabel(功率譜) 仿真結(jié)果如圖8.4所示：圖8.4 梳狀譜干擾信號仿真結(jié)果5 錐形音頻干擾信號前面介紹的幾種信號屬于壓制式干擾，而錐形音頻干擾屬于欺騙式干擾，它由高斯噪聲與錐形脈沖相乘產(chǎn)生，該信號的表達(dá)式為

17、式中，是干擾信號的幅度；是均值為0，方差為的窄帶高斯噪聲，其譜寬為0.33kHz;是錐形脈沖序列：為錐形脈沖重復(fù)周期，為0.125s。仿真程序：%-%錐形音頻干擾%信號各項(xiàng)參數(shù)n=randn(1,10e3); %返回正態(tài)分布值fp=0.98*10e5; fs=1.02*10e5; Fs=3*10e5; A=1 ; Ta=2; fj=100*10e5; %后期的采樣頻率wp=2*pi*fp/Fs; ws=2*pi*fs/Fs; %-%窗函數(shù)的各項(xiàng)參數(shù)rs=40; %損耗系數(shù)B=ws-wp; %帶寬beta=0.5842*(rs-21)0.4+0.07886*(rs-21); %波紋N=ceil(

18、rs-8)/2.285/B); %由rs和B確定凱撒窗N值wc=(wp+ws)/2/pi; %fir1是用窗函數(shù)法設(shè)計(jì)線性相位FIRDF的工具箱函數(shù)，以實(shí)現(xiàn)線性相位FIRDF的標(biāo)準(zhǔn)窗函數(shù)法設(shè)計(jì)。 %hn=fir1（M，wc），返回6dB截止頻率為wc的M階（單位脈沖響應(yīng)h(n)長度N=M+1）FIR低通（wc為標(biāo)量）濾波器系數(shù)向量hn，默認(rèn) %選用哈明窗。得到噪聲信號進(jìn)過濾波器的信號s。 %- %-h=fir1(N,wc,kaiser(N+1,beta); %單位脈沖，由圖像知道凱撒窗在這里有一個(gè)低通濾波的作用。H=fft(h); w=0:2*pi/length(H):2*pi-2*pi/le

19、ngth(H); s=conv(n,h); %卷積得到s，s是一個(gè)經(jīng)過濾波器后的信號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) subplot(2,1,1)plot(t,J) title(時(shí)域) xlabel(時(shí)間),ylabel(幅度) subplot(2,1,2) plot(f,abs(fft(J) title(頻域) xlabel(頻率

20、),ylabel(功率譜) 仿真結(jié)果如圖8.5所示：圖8.5 錐形音頻干擾信號仿真結(jié)果第9題：干擾效果仿真。仿真產(chǎn)生各種通信信號和通信干擾信號，送入通信接收機(jī)進(jìn)行解調(diào)。在不同干信比下統(tǒng)計(jì)干擾效果。1 高斯白噪聲對DSB信號的干擾以下是對不同干信比條件下高斯白噪聲對DSB信號的干擾仿真。仿真程序如下：%DSB信號加入噪聲后解調(diào)程序:通信信號為DSB信號，干擾信號為高斯白噪聲ts=0.01;%定義變量區(qū)間步長t0=2;%定義變量區(qū)間終止值t=-t0+0.0001:ts:t0;%定義變量區(qū)間fc=10;%給出相干載波的頻A=2;%定義輸入信號幅度fa=1; %定義調(diào)制信號頻率mt=A*cos(2*p

21、i*fa.*t);%輸入調(diào)制信號表達(dá)式xzb=2;%輸入小信躁比(dB)，干信比即信噪比的倒數(shù)snr=10.(xzb/10); %snr是信噪比真值h,l=size(mt); %求調(diào)制信號的維數(shù)fangcha=A*A./(2*snr); %由信躁比求方差nit=sqrt(fangcha).*randn(h,l); %產(chǎn)生小信噪比（大干信比）高斯白躁聲 psmt=mt.*cos(2*pi*fc.*t); %輸出調(diào)制信號表達(dá)式psnt=psmt+nit; %輸出疊加小信噪比已調(diào)信號波形figure(1)subplot(3,1,1),plot(t,mt),title(調(diào)制信號)subplot(3,1

22、,2),plot(t,cos(2*pi*fc.*t),title(載波信號)subplot(3,1,3),plot(t,psmt),title(已調(diào)波（DSB信號）)xzb=20; %輸入大信躁比(dB)snr1=10.(xzb/10); h,l=size(mt); %求調(diào)制信號的維數(shù)fangcha1=A*A./(2*snr1); %由信躁比求方差nit1=sqrt(fangcha1).*randn(h,l); %產(chǎn)生大信噪比（小干信比）高斯白躁聲psnt1=psmt+nit1; %輸出已調(diào)信號波形figure(2)subplot(2,2,1); %劃分畫圖區(qū)間plot(t,nit,g); %

23、畫出輸入信號波形 title(大干信比高斯白躁聲);xlabel(Variable t);ylabel(Variable nit);subplot(2,2,2);plot(t,psnt,b);title(疊加大干信比已調(diào)信號波形);xlabel(Variable t); ylabel(Variable psnt); subplot(2,2,3); plot(t,nit1,r);%length用于長度匹配title(小干信比高斯白躁聲);%畫出輸入信號與噪聲疊加波形 xlabel(Variable t);ylabel(Variable nit);subplot(2,2,4);plot(t,psn

24、t1,k); title(疊加小干信比已調(diào)信號波形); %畫出輸出信號波形 xlabel(Variable t);ylabel(Variable psmt);figure(3)subplot(3,1,1),plot(t,psmt),title(已調(diào)波（DSB信號）)subplot(3,1,2),plot(t,psnt,b);title(疊加大干信比已調(diào)信號波形);subplot(3,1,3),plot(t,psnt1,k); title(疊加小干信比已調(diào)信號波形); %畫出輸出信號波形 figure(4)psnt2=psnt1.*cos(2*pi*fc.*t);psnt3=psnt.*cos(

25、2*pi*fc.*t);subplot(2,1,1), plot(t,psnt3),title(大干信比情況下解調(diào))subplot(2,1,2), plot(t,psnt2),title(小干信比情況下解調(diào))figure(5)n1,wn1=buttord(0.1,1,1,60,s); %用于計(jì)算巴特沃斯數(shù)字濾波器的階數(shù)N和3dB截止頻率wc。b,a=butter(n1,wn1); % 計(jì)算N階巴特沃斯數(shù)字濾波器系統(tǒng)函數(shù)分子、分母多項(xiàng)式的系數(shù)向量b、a。y=filter(b,a,psnt3) %Y = filter(B,A,X) ，輸入X為濾波前序列，Y為濾波結(jié)果序列，B/A 提供濾波器系數(shù)，B

26、為分子， A為分母 subplot(3,1,1),plot(t,mt),title(調(diào)制信號)subplot(3,1,2),plot(y),title(大干信比下解調(diào)信號)y1=filter(b,a,psnt2) subplot(3,1,3),plot(y1),title(小干信比下解調(diào)信號)干擾效果如下： 圖9.1-1 調(diào)制信號+載波產(chǎn)生通信信號DSB信號圖9.1-2 干擾信號為噪聲，在不同干信比的情況下通信信號波形圖9.1-3在不同干信比的情況下通信信號波形與原始通信信號比較圖9.1-4 對受干擾的通信信號解調(diào)圖9.1-5 在不同干信比下解調(diào)結(jié)果與調(diào)制信號對比結(jié)論：由圖9.1-5可得干信比

27、較小的情況下，解調(diào)信號更加清晰，更加接近原調(diào)制信號。干信比較大的情況下，解調(diào)信號失真較嚴(yán)重。2 噪聲調(diào)頻干擾信號對2PSK信號的干擾以下是噪聲調(diào)頻干擾信號對2PSK信號的干擾仿真。仿真代碼如下：%噪聲調(diào)頻干擾2PSKclear all;close all;clc;N=150;fc=500;%載波頻率fs=150;%采樣頻率t=0:1/fs:(N-1)/fs; fb=10;n_dB=150;echo off;% -%信號調(diào)制n1=fs/fb; % The number of sample per symbolL=N/n1; % The number of symbolsPN=randsrc(1,

28、L);for k=0:L-1for i=1:n1a(n1*k+i)=PN(k+1);endendfor j=1:Nif a(j)=1b(j)=1;elseb(j)=0;endendy=20*b.*cos(2*pi*fc*t); % signal of 2ASK askzerodb=y; Y=fft(askzerodb);f=0:2*pi/150:2*pi*149/150;figure(1)subplot(211);plot(t,askzerodb);xlabel(t/s);ylabel(y)title(調(diào)制信號的時(shí)域波形) subplot(212);plot(f,Y);xlabel(w/rad

29、);ylabel(H)title(調(diào)制信號的頻域波形)%-%噪聲信號y1=noise();figure(2)subplot(211);plot(t,y1);xlabel(t),ylabel(y)title(噪聲信號的時(shí)域波形) subplot(212);plot(f,abs(fft(y1);xlabel(w/rad);ylabel(H)title(噪聲信號的頻域波形)%-%合成信號y=y+y1;figure(3)subplot(211);plot(t,y);title(合成信號的時(shí)域波形)xlabel(t/s);ylabel(y)subplot(212);plot(f,abs(fft(y);t

30、itle(合成信號的頻域波形) xlabel(w/rad);ylabel(H) %- y=y.*cos(2*pi*fc*t);figure(4)subplot(211);plot(t,y);xlabel(t/s);ylabel(y)title(解調(diào)后的時(shí)域波形)subplot(212);plot(f,abs(fft(y);xlabel(w/rad);ylabel(H)title(解調(diào)后的頻域波形) %-%低通濾波器wp=1.4;ws=1.6;rs=40; B=ws-wp; beta=0.5842*(rs-21)0.4+0.07886*(rs-21); N=ceil(rs-8)/2.285/B)

31、; 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(y,h); s=s/max(abs(s); t=0:1/150:length(s)/150-1/150;f=0:2*pi/221:2*pi*220/221;figure(5)subplot(211)plot(t,s)xlabel(t/s);ylabel(y)title(濾波后的時(shí)域波形)subplot(212)plot(f,abs(fft(s)xlabel(w/rad);ylabel(

32、H)title(濾波后的頻域波形) 仿真效果如下：圖9.2-1 調(diào)制信號的時(shí)域和頻域波形圖9.2-2 噪聲信號的時(shí)域和頻域波形圖9.2-3 合成信號的時(shí)域和頻域波形圖9.2-4 合成信號解調(diào)后的時(shí)域和頻域波形圖9.2-5 濾波后的時(shí)域和頻域波形3 錐形音頻干擾信號對2ASK信號的干擾以下是錐形音頻干擾信號對2ASK信號的干擾仿真。仿真代碼如下：%錐形干擾信號干擾2ASK信號close all;clear all;clc; N=300;fc=500;%載波頻率fs=150;%采樣頻率t=0:1/fs:(N-1)/fs; fb=10;% Modulated signals.調(diào)制信號n1=fs/fb

33、; % The number of sample per symbol每個(gè)樣品的數(shù)量L=ceil(N/n1); % The number of symbols樣品總數(shù)PN=randsrc(1,L);for k=0:L-1for i=1:n1a(n1*k+i)=PN(k+1);endendfor j=1:Nif a(j)=1b(j)=1;elseb(j)=0;endendy=20*b.*cos(2*pi*fc*t); % signal of 2ASK askzerodb=y; Y=abs(fft(askzerodb);pyy=Y.*conj(Y); %相關(guān)函數(shù)f=0:2*pi/300:2*pi*

34、299/300;figure(1)subplot(211);plot(t,askzerodb);xlabel(t/s);ylabel(y)title(2ASK信號的時(shí)域波形)subplot(212);plot(f,Y);xlabel(w/rad);ylabel(H)title(2ASK信號的頻域波形)%-%產(chǎn)生錐形干擾信號 n=0.00005*randn(1,132); fp=0.98*10e5; fs=1.02*10e5; Fs1=3*10e5; A=1 ; Ta=2; fj=150; wp=2*pi*fp/Fs1; ws=2*pi*fs/Fs1; 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); t=0:1/fj:leng