AM及SSB調(diào)制與解調(diào)

1、通信原理課程設(shè)計(jì)設(shè)計(jì)題目：AM及SSB調(diào)制與解調(diào)及抗噪聲性能分析班級(jí):學(xué)生姓名:學(xué)生學(xué)號(hào):指導(dǎo)老師:一、引言31.1 概述31.2 課程設(shè)計(jì)的目的31.3 課程設(shè)計(jì)的要求3二、AMM制與解調(diào)及抗噪聲性能分析42.1 AM調(diào)制與解調(diào)42.1.1 AM調(diào)制與解調(diào)原理42.1.2 調(diào)試過程62.2 相干解調(diào)的抗噪聲性能分析92.2.1 抗噪聲性能分析原理92.2.2 調(diào)試過程10三、SSBM制與解調(diào)及抗噪聲性能分析123.1 SSB調(diào)制與解調(diào)原理123.2 SSB調(diào)制解調(diào)系統(tǒng)抗噪聲性能分析133.3 調(diào)試過程15四、心得體會(huì)19五、參考文獻(xiàn)19一、引言1.1 概述通信原理是通信工程專業(yè)的一門極為重要

2、的專業(yè)基礎(chǔ)課，但內(nèi)容抽象，基本概念較多，是一門難度較大的課程，通過MATLA昉真能讓我們更清晰地理解它的原理，因此信號(hào)的調(diào)制與解調(diào)在通信系統(tǒng)中具有重要的作用。本課程設(shè)計(jì)是AM及SSB調(diào)制解調(diào)系統(tǒng)的設(shè)計(jì)與仿真，用于實(shí)現(xiàn)AM及SSB信號(hào)的調(diào)制解調(diào)過程，并顯示仿真結(jié)果，根據(jù)仿真顯示結(jié)果分析所設(shè)計(jì)的系統(tǒng)性能。在課程設(shè)計(jì)中，幅度調(diào)制是用調(diào)制信號(hào)去控制高頻載波的振幅，使其按調(diào)制信號(hào)的規(guī)律變化，其他參數(shù)不變。同時(shí)也是使高頻載波的振幅載有傳輸信息的調(diào)制方式。1.2 課程設(shè)計(jì)的目的在此次課程設(shè)計(jì)中，我需要通過多方搜集資料與分析：(1)掌握模擬系統(tǒng)AM和SSB調(diào)制與解調(diào)的原理；(2)來理解并掌握AM和SSB調(diào)制解

3、調(diào)的具體過程和它在MATLA訃的實(shí)現(xiàn)方法；(3)掌握應(yīng)用MATLA的析系統(tǒng)時(shí)域、頻域特性的方法，進(jìn)一步鍛煉應(yīng)用MATLAB1行編程仿真的能力。通過這個(gè)課程設(shè)計(jì)，我將更清晰地了解AMBSSB的調(diào)制解調(diào)原理，同時(shí)加深對MATLAB款通信原理輔助教學(xué)操作的熟練度。1.3 課程設(shè)計(jì)的要求(1)熟悉MATLAB勺使用方法，掌握AM言號(hào)的調(diào)制解調(diào)原理，以此為基礎(chǔ)用MATLA褊程實(shí)現(xiàn)信號(hào)的調(diào)制解調(diào)；(2)設(shè)計(jì)實(shí)現(xiàn)AM調(diào)制與解調(diào)的模擬系統(tǒng)，給出系統(tǒng)的原理框圖，對系統(tǒng)的主要參數(shù)進(jìn)行設(shè)計(jì)說明；(3)采用MATLA的言設(shè)計(jì)相關(guān)程序，實(shí)現(xiàn)系統(tǒng)的功能，要求采用一種方式進(jìn)行仿真，即直接采用MATLAB言編程的靜態(tài)方式。要

4、求采用兩種以上調(diào)制信號(hào)源進(jìn)行仿真，并記錄各個(gè)輸出點(diǎn)的波形和頻譜圖；(4)對系統(tǒng)功能進(jìn)行綜合測試，整理數(shù)據(jù)，撰寫課程設(shè)計(jì)論文。二、AMM制與解調(diào)及抗噪聲性能分析2.1 AM調(diào)制與解調(diào)2.1.1 AM調(diào)制與解調(diào)原理幅度調(diào)制是由調(diào)制信號(hào)去控制高頻載波的幅度，使正弦載波的幅度隨著調(diào)制信號(hào)而改變的調(diào)制方案，屬于線性調(diào)制。AM信號(hào)的時(shí)域表示式：sAM（Z）=4）+陋（）cos3/=4COS3J+m（Z）COS35頻譜：SAM（=乃4回（0+牡）+5（一籃）+然）+M（刃一曲）調(diào)制器模型如圖所示：Acosctc圖1-1調(diào)制器模型AM的時(shí)域波形和頻譜如圖所示：(3頻域圖1-2調(diào)制時(shí)、頻域波形2倍。AM信號(hào)的頻

5、譜由載頻分量、上邊帶、下邊帶三部分組成。它的帶寬是基帶信號(hào)帶寬的在波形上，調(diào)幅信號(hào)的幅度隨基帶信號(hào)的規(guī)律而呈正比地變化，在頻譜結(jié)構(gòu)上，它的頻譜完全是基帶信號(hào)頻譜在頻域內(nèi)的簡單搬移。所謂相干解調(diào)是為了從接受的已調(diào)信號(hào)中，不失真地恢復(fù)原調(diào)制信號(hào)，要求本地載波和接收信號(hào)的載波保證同頻同相。相干載波的一般模型如下：sp(t)LPFct=cosrtc將已調(diào)信號(hào)乘上一個(gè)與調(diào)制器同頻同相的載波，得2.Sam(t)*COSWct=A0m(t)cosWct11二2隔m(t)-Aom(t)cos2Wct由上式可知，只要用一個(gè)低通濾波器，就可以將第1項(xiàng)與第2項(xiàng)分離，無失真的恢復(fù)出原始的調(diào)制信一1_Mo(T)=AoM

6、(T)2相干解調(diào)的關(guān)鍵是必須產(chǎn)生一個(gè)與調(diào)制器同頻同相位的載波。如果同頻同相位的條件得不到滿足，則會(huì)破壞原始信號(hào)的恢復(fù)。2.1.2 調(diào)試過程t=-i：0.0000i：i；%定義時(shí)長A1=6;%調(diào)制信號(hào)振幅A2=10;%外加直流分量f=3000;%載波頻率w0=2*f*pi;%角頻率Uc=cos(w0*t);%載波信號(hào)subplot(5,2,1);plot(t,Uc);%畫載波信號(hào)title('載波信號(hào));axis(0,0.01,-1,1);%坐標(biāo)區(qū)間T1=fft(Uc);%傅里葉變換subplot(5,2,2);plot(abs(T1);%畫出載波信號(hào)頻譜title('載波信號(hào)頻

7、譜);axis(5800,6200,0,200000);%坐標(biāo)區(qū)間mes=A1*cos(0.002*w0*t);%調(diào)制信號(hào)subplot(5,2,3);plot(t,mes);%畫出調(diào)制信號(hào)title('調(diào)制信號(hào));T2=fft(mes);%傅里葉變換subplot(5,2,4);plot(abs(T2);%畫出調(diào)制信號(hào)頻譜title('調(diào)制信號(hào)頻譜');axis(198000,202000,0,1000000);%坐標(biāo)區(qū)間Uam1=A2*(1+mes/A2).*cos(w0).*t);%AM已調(diào)信號(hào)subplot(5,2,5);plot(t,Uam1);%畫出已調(diào)信號(hào)

8、title('已調(diào)信號(hào));T3=fft(Uam1);%已調(diào)信號(hào)傅里葉變換subplot(5,2,6);plot(abs(T3);%畫出已調(diào)信號(hào)頻譜title('已調(diào)信號(hào)頻譜');axis(5950,6050,0,900000);%坐標(biāo)區(qū)間sn1=20;%信噪比db1=A1A2/(2*(10A(sn1/10);%計(jì)算對應(yīng)噪聲方差n1=sqrt(db1)*randn(size(t);%生成高斯白噪聲Uam=n1+Uam1;%疊加噪聲后的已調(diào)信號(hào)Dam=Uam.*cos(w0*t);%又AM已調(diào)信號(hào)進(jìn)行解調(diào)subplot(5,2,7);plot(t,Dam);%濾波前的AM解

9、調(diào)信號(hào)title(,濾波前的AM解調(diào)信號(hào)波形,);T4=fft(Dam);%求人”信號(hào)的頻譜subplot(5,2,8);plot(abs(T4);%濾波前的AM解調(diào)信號(hào)頻譜title(,濾波前的AM解調(diào)信號(hào)頻譜,);axis(187960,188040,0,600000);Ft=2000;%采樣頻率fpts=100120;%通帶邊界頻率fp=100Hz阻帶截止頻率fs=120Hzmag=10;dev=0.010.05;%通帶波動(dòng)1%,阻帶波動(dòng)5%n21,wn21,beta,ftype=kaiserord(fpts,mag,dev,Ft);%kaiserord估計(jì)采用凱塞窗設(shè)計(jì)的FIR濾波器的

10、參數(shù)b21=fir1(n21,wn21,Kaiser(n21+1,beta);%由fir1設(shè)計(jì)濾波器z21=fftfilt(b21,Dam);%FIRf氐通濾波subplot(5,2,9);plot(t,z21,'r');%濾波后的AM解調(diào)信號(hào)title(,濾波后的AM解調(diào)信號(hào)波形,);axis(0,1,-1,10);T5=fft(z21);%求人”信號(hào)的頻譜subplot(5,2,10);plot(abs(T5),'r');%畫出濾波后的AM解調(diào)信號(hào)頻譜title(,濾波后的AM解調(diào)信號(hào)頻譜,);axis(198000,202000,0,500000);運(yùn)行結(jié)

11、果：。mokmw哪哪m嘛m則柏珊用0i''°WWAAWVW41111111.'-1«46也42Q1,114580.SIE卿YfrrrfrrrrjjAAAAAAAAAAAA-1488H4201.2H5.6。1I"1101111卿卿撇蜘U蒯制Q6150咖4雕哪11偉1,365例即2娜加挪成捕即能雕IrIIIT1III5,.AA0i.i_JL_Uu_|明躺蒯5«5«6聊m冊m唧冊q浦躺轆岫翎ilm儆m懈佛i.刪唧k翩M二1M.6m4201.2H5.5。1VIIIIIII$WWA/0'iiiiiiii'0111J

12、J10.5M1Il1UK1,M期22.005州MW12.2相干解調(diào)的抗噪聲性能分析2.2.1 抗噪聲性能分析原理AM線性調(diào)制系統(tǒng)的相干解調(diào)模型如下圖所示。施工BPFLPFCOM圖3.5.1線性調(diào)制系統(tǒng)的相干解調(diào)模型圖中Sm可以是AM調(diào)幅信號(hào)，帶通濾波器的帶寬等于已調(diào)信號(hào)帶寬10。下面討論AM調(diào)制系統(tǒng)的抗噪聲性能11。AM信號(hào)的時(shí)域表達(dá)式為Sam(t)=Aom(t)coswct通過分析可得AM信號(hào)的平均功率為.22.S_2.m(Si)AM一22又已知輸入功率Ni=n0B,其中b表示已調(diào)信號(hào)的帶寬。由此可得AM信號(hào)在解調(diào)器的輸入信噪比為(Si.Ni)AMA；m2(t)A；m2(t)2n0Bam4n

13、0fHAM信號(hào)經(jīng)相干解調(diào)器的輸出信號(hào)為因此解調(diào)后輸出信號(hào)功率為，一、-2,、12,、(So)am=m0(t)m4在上圖中輸入噪聲通過帶通濾波器之后，變成窄帶噪聲ni(t),經(jīng)乘法器相乘后的輸出噪聲為np(t)coswct=nc(t)coswct-ns(t)sinwctcoswct11=30)以門c(t)cos2wct-ns(t)sin2wct=-Ni4因此解調(diào)器的輸出噪聲功率為212N0=n0(t)nc(t)4可得AM信號(hào)經(jīng)過解調(diào)器后的輸出信噪比為(S0.No)AMm2(t)noBm2(t)2nofh由上面分析的解調(diào)器的輸入、輸出信噪比可得AM信號(hào)的信噪比增益為GAMSoNSiNi-22m(t

14、)22Aom2(t)2.2.2調(diào)試過程clf;t=o：o.oi:2;fc=5o；A=10;fa=5;mt=A*cos(2*pi*fa.*t);xzb=5;snr=1O.A(xzb/1O);db=AA2./(2*snr);nit=sqrt(db).*randn(size(mt);psmt=(A+mt).*cos(2*pi*fc.*t);psnt=psmt+nit;xzb1=3O;snr1=1O.A(xzb1/1O);db1=AA2./(2*snr1);nit1=sqrt(db1).*randn(size(mt);psnt1=psmt+nit1;subplot(2,2,1);plot(t,nit,

15、'g');title('小信噪比高斯白躁聲,);xlabel('t');ylabel('nit');subplot(2,2,2);plot(t,psnt,'b');title('疊加小信噪比已調(diào)信號(hào)波形');%清除窗口中的圖形%定義變量區(qū)間%給出相干載波的頻率%定義輸入信號(hào)幅度%定義調(diào)制信號(hào)頻率%輸入調(diào)制信號(hào)表達(dá)式%輸入小信躁比(dB)%由信躁比求方差%產(chǎn)生小信噪比高斯白躁聲%輸出調(diào)制信號(hào)表達(dá)式%輸出疊加小信噪比已調(diào)信號(hào)波形%輸入大信躁比(dB)%由信躁比求方差%產(chǎn)生大信噪比高斯白躁聲%輸出已調(diào)信號(hào)波形%劃

16、分畫圖區(qū)間%畫出輸入信號(hào)波形xlabel('時(shí)間');ylabel('輸出調(diào)制信號(hào),);subplot(2,2,3);plot(t,nit1,'r');%length用于長度匹配title(,大信噪比高斯白躁聲,);%畫出輸入信號(hào)與噪聲疊加波形xlabel('t');ylabel('nit');subplot(2,2,4);plot(t,psnt1,'k');title('疊加大信噪比已調(diào)信號(hào)波形');%畫出輸出信號(hào)波形xlabel('時(shí)間');ylabel('輸出調(diào)

17、制信號(hào)');運(yùn)行結(jié)果：小信噪比高斯白躁聲20-賽加小信噪比已調(diào)信號(hào)波形40-20111Q0.511.52t大信噪比高斯白躁聲1-4Q1100511.52時(shí)間嶷加大信喋比已調(diào)信號(hào)波形40-100.511.5211.52時(shí)間由上圖可見，當(dāng)輸入信號(hào)一定時(shí)，隨著噪聲的加強(qiáng)，接收端輸入信號(hào)被干擾得越嚴(yán)重。而相應(yīng)的輸出波形相對于發(fā)送端的波形誤差也越大。而當(dāng)噪聲過大時(shí)信號(hào)幾難分辨。這信噪比變小導(dǎo)致的在實(shí)際的信號(hào)傳輸過程中當(dāng)信道噪聲過大將會(huì)導(dǎo)致幅度相位等各失真當(dāng)然由于非線性元件如濾波器等的存在。非線性失真也會(huì)隨噪聲加大而變大。、SSB調(diào)制與解調(diào)及抗噪聲性能分析3.1SSB調(diào)制與解調(diào)原理單邊帶調(diào)制信號(hào)是

18、將雙邊帶信號(hào)中的一個(gè)邊帶濾掉而形成的。根據(jù)方法的不同，產(chǎn)生SSB信號(hào)的方法有：濾波法和相移法。由于濾波法在技術(shù)上比較難實(shí)現(xiàn)，所以在此我們將用相移法對SS則制與解調(diào)系統(tǒng)進(jìn)行討論與設(shè)計(jì)。相移法和SSB信號(hào)的時(shí)域表示設(shè)單頻調(diào)制信號(hào)為m(t)=Amcosmt載波為c(t)=cosct則其雙邊帶信號(hào)DSB信號(hào)的時(shí)域表示式為sdsb(t)=Amcosmtcosct11A2Amcos('cm)t2Amcos('c1m)t若保留上邊帶，則有1、，SuSB(t)=Amcos(Cm)t=2若保留下邊帶，則有一Amcosmcosct-Amsinmsinct_IIIIIIcIIIIIIc22小1A，、

19、，1slsb(t)=Aicos(C-I)t=-22，1，.，AIcosItcosctAIsinItsinct2將上兩式合并得:s=Amco$8/cosfij/isinsing由希爾伯特變換Amc?smt=Amsinmt故單邊帶信號(hào)經(jīng)過希爾伯特變換后得：1_1人.把上式推廣到一般情況，則得到“。)=jm(f)com8的傅金咻變換i而興偽式中囹是m的希爾伯特變換若M。)是mt)的傅里葉變換，則M?()=M(,)I-jsgn1上式中的-jsgn可以看作是希爾伯特濾波器傳遞函數(shù)，即Hh(')=M?()/M()=-jsgn移相法SSB調(diào)制器方框圖相移法是利用相移網(wǎng)絡(luò)，對載波和調(diào)制信號(hào)進(jìn)行適當(dāng)?shù)南?/p>

20、移，以便在合成過程中將其中的一個(gè)邊帶抵消而獲得SSB信號(hào)。相移法不需要濾波器具有陡峭的截止特性，不論載頻有多高，均可一次實(shí)現(xiàn)SSB調(diào)制。SSB信號(hào)的解調(diào)SSB信號(hào)的解調(diào)不能采用簡單的包絡(luò)檢波，因?yàn)镾SB信號(hào)是抑制載波的已調(diào)信號(hào)，它的包絡(luò)不能直接反映調(diào)制信號(hào)的變化，所以仍需采用相干解調(diào)。SSB信號(hào)的性能SSB信號(hào)的實(shí)現(xiàn)比AMDS期復(fù)雜，但SSB調(diào)制方式在傳輸信息時(shí)，不僅可節(jié)省發(fā)射功率，而且它所占用的頻帶寬度比AMDS聰少了一半。它目前已成為短波通信中一種重要的調(diào)制方式。3.2SSB調(diào)制解調(diào)系統(tǒng)抗噪聲性能分析噪聲功率1 1NoNin°B44這里，B=fH為SSB信號(hào)的帶通濾波器的帶寬。信

21、號(hào)功率SSB信號(hào)«)=gj不;而。)sin與相干載波相乘后，再經(jīng)低通濾波可得解調(diào)器輸出信號(hào)，、1，、m°(t)m(t)4因此，輸出信號(hào)平均功率八21-28b=n(t)=-m輸入信號(hào)平均功率為H=£：(，)二sin=；配,(O+g而“D因m(t)與m(t)的幅度相同，所以具有相同的平均功率，故上式C12Sim2(t)4單邊帶解調(diào)器的輸入信噪比為Si；m2(t)=4m2(t)Nin0B4n0B單邊帶解調(diào)器的輸出信噪比為c1m2(t)So16m(t)N。1noB4%B4制度增益GSSBSo/N。S/N號(hào)和噪聲中的正交分量因?yàn)樵赟SB系統(tǒng)中，信號(hào)和噪聲有相同表示形式，所以

22、相干解調(diào)過程中,均被抑制掉，故信噪比沒有改善。3.3調(diào)試過程先建立3個(gè)M文件1.afd_buttfunctionb,a=afd_butt(Wp,Ws,Rp,As);ifWp<=0error('Passbandedgemustbelargerthan0')endifWs<=Wperror('StopbandedgemustbelargerthanPassbandedge')endif(Rp<=0)|(As<0)error('PBrippleand/orSBattenuationustbelargerthan0')endN=c

23、eil(log10(10A(Rp/10)-1)/(10A(As/10)-1)/(2*log10(Wp/Ws);fprintf('n*ButterworthFilterOrder=%2.0fn',N)OmegaC=Wp/(10A(Rp/10)-1)A(1/(2*N);b,a=u_buttap(N,OmegaC);2. u_buttapfunctionb,a=u_buttap(N,Omegac);z,p,k=buttap(N);p=p*Omegac;k=k*OmegacAN;B=real(poly(z);b0=k;b=k*B;a=real(poly(p);3. imp_invr%脈

24、沖向應(yīng)不變法子程序functionb,a=imp_invr(c,d,T)R,p,k=residue(c,d);p=exp(p*T);b,a=residuez(R,p,k);b=real(b).*T;a=real(a);程序：> >t0=0.1;fs=12000;%t0采樣區(qū)間，fs采樣頻率> >fc=1000;Vm0=2.5;ma=0.25;%fc載波頻率，Vm0輸出載波電壓振幅，ma調(diào)幅度> >n=-t0/2:1/fs:t0/2;%定義變量區(qū)間> >N=length(n);> >A=4;%定義調(diào)制信號(hào)幅度> >x1=A*

25、cos(150*pi*n);%調(diào)制信號(hào)> >x2=hilbert(x1,N);%對x1做希爾伯特變換> >y=(Vm0*x1.*cos(2*pi*fc*n)-Vm0*x2.*sin(2*pi*fc*n)/2;%保留上邊帶的已調(diào)波信號(hào)> >xzb=2;%輸入小信噪比(dB)> >snr=10.A(xzb/10);> >h,l=size(x1);%求調(diào)制信號(hào)的維度> >fangcha=A*A./(2*snr);%由信噪比求方差> >nit=sqrt(fangcha).*randn(h,l);%產(chǎn)生高斯白噪聲>

26、 >yn=y+nit;%疊加小信噪比噪聲的已調(diào)波信號(hào)> >xzb=10;%輸入小信噪比(dB)> >snr1=10.A(xzb/10);> >h,l=size(x1);%求調(diào)制信號(hào)的維度> >fangcha=A*A./(2*snr1);%由信噪比求方差> >nit1=sqrt(fangcha).*randn(h,l);%產(chǎn)生高斯白噪聲> >yn1=y+nit1;%疊加小信噪比噪聲的已調(diào)波信號(hào)> >figure(1)> >subplot(2,2,1)%劃分畫圖區(qū)間> >plot(n,

27、x1)%畫出調(diào)制信號(hào)的波形> >title('調(diào)制信號(hào)');> >subplot(2,2,2)%劃分畫圖區(qū)間> >plot(n,y)%畫出已調(diào)波信號(hào)波形> >title(,已調(diào)波信號(hào),);> >subplot(2,2,3)%劃分畫圖區(qū)間> >plot(n,yn)%畫出疊加噪聲的已調(diào)波信號(hào)波形> >title(,疊加小信噪比噪聲的已調(diào)波信號(hào),);> >subplot(2,2,4)%劃分畫圖區(qū)間> >plot(n,yn1)%畫出疊加噪聲的已調(diào)波信號(hào)波形> >tit

28、le('疊加大信噪比噪聲的已調(diào)波信號(hào),);> >X=fft(x1);%調(diào)制信號(hào)x1的傅里葉變換> >Y=fft(y);%已調(diào)信號(hào)y的傅里葉變換> >Yn=fft(yn);%疊加小信噪比噪聲的已調(diào)信號(hào)yn的傅里葉變換> >Yn1=fft(yn1);%疊加大信噪比噪聲的已調(diào)信號(hào)yn的傅里葉變換> >w=0:2*pi/(N-1):2*pi;%定義變量區(qū)間> >figure(2)> >subplot(2,2,1)%劃分畫圖區(qū)間> >plot(w,abs(X)%畫出調(diào)制信號(hào)頻譜波形> >a

29、xis(0,pi/4,0,3000);%給出橫縱坐標(biāo)的范圍> >title('調(diào)制信號(hào)頻譜,);> >subplot(2,2,2)%劃分畫圖區(qū)間> >plot(w,abs(Y)%畫出已調(diào)波信號(hào)頻譜> >axis(pi/6,pi/4,0,2500);%給出橫縱坐標(biāo)的范圍> >title(,已調(diào)波信號(hào)頻譜,);> >subplot(2,2,3)%劃分畫圖區(qū)間> >plot(w,abs(Yn)%畫出疊加小信噪比噪聲的已調(diào)波信號(hào)頻譜> >axis(pi/6,pi/4,0,2500);%給出橫縱坐標(biāo)的

30、范圍> >title(,疊加小信噪比噪聲的已調(diào)波信號(hào)頻譜,);> >subplot(2,2,4)%劃分畫圖區(qū)間> >plot(w,abs(Yn1)%畫出疊加大信噪比噪聲的已調(diào)波信號(hào)頻譜> >axis(pi/6,pi/4,0,2500);%給出橫縱坐標(biāo)的范圍> >title(,疊加大信噪比噪聲的已調(diào)波信號(hào)頻譜,);> >y1=y-2*cos(1500*pi*n);> >y2=Vm0*y1.*cos(2*pi*fc*n);%將已調(diào)幅波信號(hào)的頻譜搬移到原調(diào)制信號(hào)的頻譜處> >wp=40/N*pi;ws=6

31、0/N*pi;Rp=1;As=15;T=1;%濾波器參數(shù)設(shè)計(jì)> >OmegaP=wp/T;OmegaS=ws/T;> >cs,ds=afd_butt(OmegaP,OmegaS,Rp,As);> >b,a=imp_invr(cs,ds,T);> >y=filter(b,a,y2);> >yn=y+nit;> >figure(3)> >subplot(2,1,1)%劃分畫圖區(qū)間> >plot(n,y)%畫出解調(diào)波波形> >title('解調(diào)波,);> >Y=fft(y)

32、;%解調(diào)波y的傅里葉變換> >subplot(2,1,2)%劃分畫圖區(qū)間> >plot(w,abs(Y)%畫出解調(diào)信號(hào)頻譜> >axis(0,pi/6,0,2500);%給出橫縱坐標(biāo)的范圍> >title('解調(diào)信號(hào)頻譜,);運(yùn)行結(jié)果：> *ButterworthFilterOrder=6420-2-4-0.0500.05調(diào)制信號(hào)已調(diào)波信號(hào)疊加小信噪比噪聲的已調(diào)波信號(hào)20100-10-20-0.0500.05疊加大信噪比噪聲的已調(diào)波信號(hào)調(diào)制信號(hào)頻譜300020001000000.20.40.625002000150010005000已調(diào)