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、AM調(diào)制與解調(diào)及抗噪聲性能分析4.2.1 AM調(diào)制與解調(diào)2.1.1 AM調(diào)制與解調(diào)原理調(diào)試過程62.2相干解調(diào)的抗噪聲性能分析 10抗噪聲性能分析原理 10調(diào)試過程11三、SSB調(diào)制與解調(diào)及抗噪聲性能分析 133.1 SSB調(diào)制與解調(diào)原理 133.2 SSB調(diào)制解調(diào)系統(tǒng)抗噪聲性能分析 143.3調(diào)試過程16四、心得體會(huì)20.五、參考文獻(xiàn)21引言1.1概述通信原理是通信工程專業(yè)的一門極為重要的專業(yè)基礎(chǔ)課，但容抽象，基本概念較多，是一

2、 門難度較大的課程，通過MATLAB仿真能讓我們更清晰地理解它的原理，因此信號(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)制解調(diào)的具體過程和它在MATLAB

3、中的實(shí)現(xiàn)方法；(3) 掌握應(yīng)用MATLAB分析系統(tǒng)時(shí)域、頻域特性的方法，進(jìn)一步鍛煉應(yīng)用MATLAB進(jìn)行編程仿真的能力。通過這個(gè)課程設(shè)計(jì)，我將更清晰地了解AM和SSB的調(diào)制解調(diào)原理，同時(shí)加深對(duì)MATLAB這款通信原理輔助教學(xué)操作的熟練度。1.3課程設(shè)計(jì)的要求(1) 熟悉MATLAB的使用方法，掌握AM信號(hào)的調(diào)制解調(diào)原理， 以此為基礎(chǔ)用 MATLAB編 程實(shí)現(xiàn)信號(hào)的調(diào)制解調(diào)；(2) 設(shè)計(jì)實(shí)現(xiàn)AM調(diào)制與解調(diào)的模擬系統(tǒng)， 給出系統(tǒng)的原理框圖， 對(duì)系統(tǒng)的主要參數(shù)進(jìn)行設(shè) 計(jì)說明；(3) 采用MATLAB語言設(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) 對(duì)系統(tǒng)功能進(jìn)行綜合測(cè)試，整理數(shù)據(jù)，撰寫課程設(shè)計(jì)論文。二、AM調(diào)制與解調(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) = Jo + ff?(/)cos(z?t Z =cos to f + 用(f) cos氓/頻譜：2調(diào)制器模型如圖所示：口 t >>賂t評(píng)f tAcos ct圖1-1調(diào)制器模型AM的時(shí)域波形和頻譜如圖所示：圖1-2調(diào)制時(shí)、頻域波形A

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

6、出原始的調(diào)M°(T)如。M(T)相干解調(diào)的關(guān)鍵是必須產(chǎn)生一個(gè)與調(diào)制器同頻同相位的載波。如果同頻同相位的條件得不到滿足，則會(huì)破壞原始信號(hào)的恢復(fù)。調(diào)試過程t=-1:0.00001:1; % 定義時(shí)長A1=6; %調(diào)制信號(hào)振幅A2=10; %外加直流分量f=3000;%載波頻率w0=2*f*pi;%角頻率Uc=cos(wO*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);%畫出載波

7、信號(hào)頻譜title('載波信號(hào)頻譜');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ū)間Uam 仁A2*(1+mes/A2).*cos(w0).*t);%AM 已調(diào)

8、信號(hào)subplot(5,2,5);plot(t,Uam1);%畫出已調(diào)信號(hào)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ū)間sn 1=20; %信噪比db仁A1A2/(2*(10A(s n1/10); % 計(jì)算對(duì)應(yīng)噪聲方差n1=sqrt(db1)*randn(size(t); % 生成高斯白噪聲Uam=n1+Uam1; %疊加噪聲后的已調(diào)信號(hào)Dam=Uam.*cos

9、(wO*t); %對(duì)AM已調(diào)信號(hào)進(jìn)行解調(diào)subplot(5,2,7);plot(t,Dam);%濾波前的AM解調(diào)信號(hào)title('濾波前的AM解調(diào)信號(hào)波形');T4=fft(Dam); %求AM信號(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=100 120; % 通帶邊界頻率fp=100Hz 阻帶截止頻率fs=120Hzmag=1 0;dev=0.01 0.05;%通帶波動(dòng)1%,阻

10、帶波動(dòng)5%n 21,w n21,beta,ftype=kaiserord(fpts,mag,dev,Ft);%kaiserord估計(jì)采用凱塞窗設(shè)計(jì)的 FIR濾波器的參數(shù)b21=fir1(n21,wn21,Kaiser(n21+1,beta);firl 設(shè)計(jì)濾波器z2 仁fftfilt(b21,Dam);%FIF低 通濾波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); %求AM信號(hào)的頻譜subplot(5,2,10);plot(

11、abs(T5),'r');%畫出濾波后的AM解調(diào)信號(hào)頻譜title('濾波后的AM解調(diào)信號(hào)頻譜');axis(198000,202000,0,500000);運(yùn)行結(jié)果：I I I I I 丨 I'I擁樁刖吧i n u 0.6 m 1D1-_1-1S 1.585 1J1 Iffi I 2伽加 2.(15 WH3U1 擁征 4J 4i U DJ D.E 0.) 1卅111 1111 一i卿泓迦訓(xùn)【M 汎EK孤 W 做.1111111111擁愛刖罷I徂ij(,f g.s 1)0.1 OJ W t.4 0.5 IE W 0.6 U 1訓(xùn)b I I I 2 I I

12、 S IWWTTWWT 一、一 |uik m ib m i.sa m w w ukj2.2相干解調(diào)的抗噪聲性能分析221抗噪聲性能分析原理AM線性調(diào)制系統(tǒng)的相干解調(diào)模型如下圖所示。圖線性調(diào)制系統(tǒng)的相干解調(diào)模型!j(C10。下面討論AM調(diào)制系圖中Sm(t)可以是am調(diào)幅信號(hào)，帶通濾波器的帶寬等于已調(diào)信號(hào)帶寬 統(tǒng)的抗噪聲性能11。AM信號(hào)的時(shí)域表達(dá)式為m(t)coswctSam (t)A0通過分析可得AM信號(hào)的平均功率為(Si )AMA0m2 (t)又已知輸入功率 Ni n0B ,其中B表示已調(diào)信號(hào)的帶寬。 由此可得AM信號(hào)在解調(diào)器的輸入信噪比為(Si : N i ) AMA m2(t)A2m2(

13、t)2n 0 B AM4n°fAM信號(hào)經(jīng)相干解調(diào)器的輸出信號(hào)為m°(t)2m(t)因此解調(diào)后輸出信號(hào)功率為(S0 ) AMn i(t)，經(jīng)乘法器相乘后的輸出噪聲為在上圖中輸入噪聲通過帶通濾波器之后，n p(t) ni (t)cosw ct n c(t)coswct- ns(t)si nw ctcoswct12n c(t)變成窄帶噪聲1n c(t)cos2wct -n s(t)s in 2w ct經(jīng)LPF后,no(t)1n c(t)因此解調(diào)器的輸出噪聲功率為Non2(t)丄而-Ni44可得AM信號(hào)經(jīng)過解調(diào)器后的輸出信噪比為AM信號(hào)的信噪比增益為22m,t)A m2(t)S0

14、N 0Si Ni% 青除窗口中的圖形%定義變量區(qū)間% 合出相干載波的頻率%定義輸入信號(hào)幅度%定義調(diào)制信號(hào)頻率% 俞入調(diào)制信號(hào)表達(dá)式%俞入小信躁比(dB)%h信躁比求方差滬生小信噪比高斯白躁聲% 俞出調(diào)制信號(hào)表達(dá)式%俞出疊加小信噪比已調(diào)信號(hào)波形%俞入大信躁比(dB)%由信躁比求方差滬生大信噪比高斯白躁聲%俞出已調(diào)信號(hào)波形%劃分畫圖區(qū)間%畫出輸入信號(hào)波形2 2(so Nomf mf由上面分析的解調(diào)器的輸入、輸出信噪比可得G AM調(diào)試過程elf;t=0:0.01:2;fe=50;A=10;fa=5;mt=A*cos(2*pi*fa.*t);xzb=5;sn r=10.A(xzb/10);db=AA2

15、./(2*s nr);n it=sqrt(db).*ra ndn( size(mt);psmt=(A+mt).*eos(2*pi*fe.*t);psn t=psmt+ nit;xzb1=30;sn r1=10.A(xzb1/10);db1=AA2./(2*s nr1);n it仁sqrt(db1).*ra ndn( size(mt);psn t1=psmt+ nitl;subplot(2,2,1);plot(t, nit,'g');title('小信噪比高斯白躁聲');xlabel(' t');ylabel(' nit');sub

16、plot(2,2,2);plot(t,ps nt,'b');title（'疊加小信噪比已調(diào)信號(hào)波形'）;xlabel（'時(shí)間'）;ylabel（'輸出調(diào)制信號(hào)'）;subplot（2,2,3）;plot（t ,n it1,'r'）;title（'大信噪比高斯白躁聲'）;xlabel（' t'）;ylabel（' ni t'）;subplot（2,2,4）;plot（t,ps nt1,'k'）;title（'疊加大信噪比已調(diào)信號(hào)波形'）;

17、xlabel（'時(shí)間'）;ylabel（'輸出調(diào)制信號(hào)'）;%le ngth用于長度匹配%畫出輸入信號(hào)與噪聲疊加波形%畫出輸出信號(hào)波形運(yùn)行結(jié)果:小信噪比高斯白躁聲400.511.622010苣 00.511.52時(shí)間20-2010 5壬 0-0 5犬信噪比專斯白蹦聲405120O-2Q由上圖可見，當(dāng)輸入信號(hào)一定時(shí)，隨著噪聲的加強(qiáng)，接收端輸入信號(hào)被干擾得越嚴(yán)重。而 相應(yīng)的輸出波形相對(duì)于發(fā)送端的波形誤差也越大。而當(dāng)噪聲過大時(shí)信號(hào)幾難分辨。這信噪比變小導(dǎo)致的 在實(shí)際的信號(hào)傳輸過程中當(dāng)信道噪聲過大將會(huì)導(dǎo)致幅度相位等各失真當(dāng)然由于非線性元件如濾波器等的存在。非線性失真也會(huì)

18、隨噪聲加大而變大。、SSB調(diào)制與解調(diào)及抗噪聲性能分析3.1 SSB調(diào)制與解調(diào)原理單邊帶調(diào)制信號(hào)是將雙邊帶信號(hào)中的一個(gè)邊帶濾掉而形成的。根據(jù)方法的不同，產(chǎn)生SSB信號(hào)的方法有：濾波法和相移法。由于濾波法在技術(shù)上比較難實(shí)現(xiàn)，所以在此我們將用相移法對(duì)SSB調(diào)制與解調(diào)系統(tǒng)進(jìn)行討論與設(shè)計(jì)。相移法和SSB信號(hào)的時(shí)域表示設(shè)單頻調(diào)制信號(hào)為 m(t) Am cos mt載波為 c(t) cos ct則其雙邊帶信號(hào)DSB信號(hào)的時(shí)域表示式為SDSB (t)Am cos mtcos ct1 12 Am cos( c m)t -Amcos( c m)t若保留上邊帶，則有1 1 1 sUSB (t)2 Am cos( c

19、 m)t 2 Am cos m COS 菽 2 Am sin m Sin 菽若保留下邊帶，則有Z 2 Am cos( Cm)t2Amcos mg 龍2Amsln mtsln 龍將上兩式合并得：£ 曲(f)二 * 蟲加 cos 血 J=F g *押 sin ojmt sin(oc t由希爾伯特變換Am cos mtAm sin mt故單邊帶信號(hào)經(jīng)過希爾伯特變換后得1 _ 1 *畑(0 = - cos cos 4 cos sm 叫把上式推廣到一般情況，則得到s則 =£血(r) com糾的傅里葉變換iM?億)為式中m(t)是m(t )的希爾伯特變換若M()是m(t)的傅里葉變換，

20、貝UM?( ) M( )jsgn上式中的-jsgn 可以看作是希爾伯特濾波器傳遞函數(shù)，即Hh( )M?( )/M( ) jsgn移相法SSB調(diào)制器方框圖相移法是利用相移網(wǎng)絡(luò)，對(duì)載波和調(diào)制信號(hào)進(jìn)行適當(dāng)?shù)南嘁?，以便在合成過程中將其中的一個(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)比AM、DSB要復(fù)雜，但SSB調(diào)制方式在傳輸信息時(shí)， 不僅可節(jié)省發(fā)射功率,

21、 而且它所占用的頻帶寬度比 AM、DSB減少了一半。它目前已成為短波通信中一種重要的調(diào)制方式。3.2 SSB調(diào)制解調(diào)系統(tǒng)抗噪聲性能分析噪聲功率11NoNinoB44這里，B = fH為SSB信號(hào)的帶通濾波器的帶寬。信號(hào)功率SSB信號(hào)5込(f) = m(f)cos cojt m (f) sin oj 2 2與相干載波相乘后，再經(jīng)低通濾波可得解調(diào)器輸出信號(hào)1m°(t)-m(t)4因此，輸出信號(hào)平均功率2 1 2S) mo(t) -m (t)16輸入信號(hào)平均功率為S嚴(yán) £ (f) = -cosa?£?干 A(f)sin 毋.于4因m(t)與m(t)的幅度相同，所以具有相

22、同的平均功率，故上式1Si -m (t)4單邊帶解調(diào)器的輸入信噪比為1 一 m (t)Si4m (t)Nin0B4n 0B單邊帶解調(diào)器的輸出信噪比為1 m?(t)S。16m (t)N。丄4noB4制度增益GssbSo/N。1S/叫因?yàn)樵赟SB系統(tǒng)中，信號(hào)和噪聲有相同表示形式，所以相干解調(diào)過程中，信號(hào)和噪聲中的正交分量 均被抑制掉，故信噪比沒有改善。3.3調(diào)試過程先建立3個(gè)M文件1. afd_buttfunction b,a = afd_butt(Wp,Ws,Rp,As);if Wp <= 0error( 'Passband edge must be larger than O&#

23、39; ) end if Ws <= Wperror( 'Stopband edge must be larger than Passband edge' ) endif (Rp <= 0) | (As < 0)error( 'PB ripple and/or SB attenuation ust be larger than 0')endN = ceil(log10(10A(Rp/10)-1)/(10A(As/10)-1)/(2*log10(Wp/Ws);fprintf( 'n* Butterworth Filter Order =

24、%2.0f n',N)OmegaC = Wp/(10A(Rp/10)-1)A(1/(2*N);b,a=u_butt ap(N,OmegaC);2. u_buttapfunction b,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. impnvr淋沖響應(yīng)不變法子程序function b,a=imp_invr(c,d,T) R,p,k=residue(c,d); p=exp(p*T);b,a=r

25、esiduez(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載波頻率，VmO輸出載波電壓振幅，ma調(diào)幅度>> n=-t0/2:1/fs:t0/2;%老義變量區(qū)間>> N=le ngth( n);>> A=4;淀義調(diào)制信號(hào)幅度>> x仁 A*cos(150*pi* n);%調(diào)制信號(hào)>> x2=hilbert(x1,N);%對(duì)x1做希爾伯特變換>> y=(

26、VmO*x1.*cos(2*pi*fc* n)-VmO*x2.*si n(2*pi*fc* n) )/2;%保留上邊帶的已調(diào)波信號(hào)>> xzb=2;%俞入小信噪比(dB)>> sn r=1O.A(xzb/1O);>> h,l=size(x1);%求調(diào)制信號(hào)的維度>> fan gcha=A*A./(2*s nr);%由信噪比求方差>> n it=sqrt(fa ngcha).*ra ndn( h,l);%產(chǎn)生高斯白噪聲>> yn=y+nit;%疊加小信噪比噪聲的已調(diào)波信號(hào)>> xzb=10;%輸入小信噪比(dB)

27、>> sn r1=10.A(xzb/10);>> h,l=size(x1);%求調(diào)制信號(hào)的維度>> fan gcha=A*A./(2*s nr1);%由信噪比求方差>> n it仁sqrt(fa ngcha).*ra ndn( h,l);%產(chǎn)生高斯白噪聲>> yn仁y+ni t1;%疊加小信噪比噪聲的已調(diào)波信號(hào)>> figure(1)>> subplot(2,2,1)%劃分畫圖區(qū)間>> plot( n,x1)%畫出調(diào)制信號(hào)的波形>> title('調(diào)制信號(hào)');>&

28、gt; 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)波形>> title('疊加大信噪比噪聲的已調(diào)波信號(hào))；>&

29、gt; X=fft(x1);%調(diào)制信號(hào)x1的傅里葉變換>> Y=fft(y);%已調(diào)信號(hào)y的傅里葉變換>> Yn=fft(yn);%疊加小信噪比噪聲的已調(diào)信號(hào)yn的傅里葉變換>> Yn仁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)頻譜波形>> axis(0,pi/4,0,3000);%合出橫縱坐標(biāo)的圍&g

30、t;> 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(Y n)%畫出疊加小信噪比噪聲的已調(diào)波信號(hào)頻譜>> axis(pi/6,pi/4,0,2500);%合出橫縱坐標(biāo)的圍>> title(

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

32、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 _in vr(cs,ds,T);>> y=filter(b,a,y2);>> yn=y+n it;>> figure(3)>> subplot(2,1,1)%劃分畫圖區(qū)間>> plot (n ,y)%畫出解調(diào)波波形>> title('解調(diào)波');>> Y=fft(y);%軍調(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é)果：* Butterworth Filter Order =6疊加小信噪比噪聲的已調(diào)波信號(hào)207已調(diào)波信號(hào)100-10-20 :-0.0500.05300020001000000.20.