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

1、精選優(yōu)質(zhì)文檔-傾情為你奉上通信原理課程設(shè)計(jì)設(shè)計(jì)題目：AM及SSB調(diào)制與解調(diào)及抗噪聲性能分析班 級(jí)： 學(xué)生： 學(xué)生學(xué)號(hào)：指導(dǎo)老師： 目 錄一、引言1.1 概述通信原理是通信工程專(zhuān)業(yè)的一門(mén)極為重要的專(zhuān)業(yè)基礎(chǔ)課，但容抽象，基本概念較多，是一門(mén)難度較大的課程，通過(guò)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)過(guò)程，并顯示仿真結(jié)果，根據(jù)仿真顯示結(jié)果分析所設(shè)計(jì)的系統(tǒng)性能。在課程設(shè)計(jì)中，幅度調(diào)制是用調(diào)制信號(hào)去控制高頻載波的振幅，使其按調(diào)制信號(hào)的規(guī)律變化，其他參數(shù)不變。同時(shí)也是使

2、高頻載波的振幅載有傳輸信息的調(diào)制方式。1.2 課程設(shè)計(jì)的目的在此次課程設(shè)計(jì)中，我需要通過(guò)多方搜集資料與分析：(1) 掌握模擬系統(tǒng)AM和SSB調(diào)制與解調(diào)的原理；(2) 來(lái)理解并掌握AM和SSB調(diào)制解調(diào)的具體過(guò)程和它在MATLAB中的實(shí)現(xiàn)方法；(3) 掌握應(yīng)用MATLAB分析系統(tǒng)時(shí)域、頻域特性的方法，進(jìn)一步鍛煉應(yīng)用MATLAB進(jìn)行編程仿真的能力。通過(guò)這個(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)；(

3、2) 設(shè)計(jì)實(shí)現(xiàn)AM調(diào)制與解調(diào)的模擬系統(tǒng)，給出系統(tǒng)的原理框圖，對(duì)系統(tǒng)的主要參數(shù)進(jìn)行設(shè)計(jì)說(shuō)明；(3) 采用MATLAB語(yǔ)言設(shè)計(jì)相關(guān)程序，實(shí)現(xiàn)系統(tǒng)的功能，要求采用一種方式進(jìn)行仿真，即直接采用MATLAB語(yǔ)言編程的靜態(tài)方式。要求采用兩種以上調(diào)制信號(hào)源進(jìn)行仿真，并記錄各個(gè)輸出點(diǎn)的波形和頻譜圖；(4) 對(duì)系統(tǒng)功能進(jìn)行綜合測(cè)試，整理數(shù)據(jù)，撰寫(xiě)課程設(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í)域表示式：頻譜： 調(diào)制器模型如圖所示：圖1-1 調(diào)制器模

4、型AM的時(shí)域波形和頻譜如圖所示：時(shí)域 頻域 圖1-2 調(diào)制時(shí)、頻域波形AM信號(hào)的頻譜由載頻分量、上邊帶、下邊帶三部分組成。它的帶寬是基帶信號(hào)帶寬的2倍。在波形上，調(diào)幅信號(hào)的幅度隨基帶信號(hào)的規(guī)律而呈正比地變化，在頻譜結(jié)構(gòu)上，它的頻譜完全是基帶信號(hào)頻譜在頻域的簡(jiǎn)單搬移。所謂相干解調(diào)是為了從接受的已調(diào)信號(hào)中，不失真地恢復(fù)原調(diào)制信號(hào)，要求本地載波和接收信號(hào)的載波保證同頻同相。相干載波的一般模型如下： 將已調(diào)信號(hào)乘上一個(gè)與調(diào)制器同頻同相的載波，得 由上式可知，只要用一個(gè)低通濾波器，就可以將第1項(xiàng)與第2項(xiàng)分離，無(wú)失真的恢復(fù)出原始的調(diào)制信號(hào) 相干解調(diào)的關(guān)鍵是必須產(chǎn)生一個(gè)與調(diào)制器同頻同相位的載波。如果同頻同相

5、位的條件得不到滿(mǎn)足，則會(huì)破壞原始信號(hào)的恢復(fù)。2.1.2調(diào)試過(guò)程t=-1:0.00001:1; %定義時(shí)長(zhǎng) 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); %畫(huà)載波信號(hào)title('載波信號(hào)'); axis(0,0.01,-1,1); %坐標(biāo)區(qū)間T1=fft(Uc); %傅里葉變換subplot(5,2,2); plot(abs(T1);%畫(huà)出載波信號(hào)頻譜title('載波信號(hào)頻譜'); axis(5800,

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

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

8、#39;濾波前的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(,0,);Ft=2000; %采樣頻率fpts=100 120; %通帶邊界頻率fp=100Hz 阻帶截止頻率fs=120Hzmag=1 0; dev=0.01 0.05; %通帶波動(dòng)1%，阻帶波動(dòng)5%n21,wn21,beta,ftype=kaiserord(fpts,mag,dev,Ft);%kaiserord估計(jì)采用凱塞窗設(shè)計(jì)的FIR濾波器的

9、參數(shù)b21=fir1(n21,wn21,Kaiser(n21+1,beta); %由fir1設(shè)計(jì)濾波器z21=fftfilt(b21,Dam); %FIR低通濾波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(abs(T5),'r');%畫(huà)出濾波后的AM解調(diào)信號(hào)頻譜title('濾波后的AM解調(diào)信號(hào)頻譜'); axis(,0,)

10、; 運(yùn)行結(jié)果:2.2 相干解調(diào)的抗噪聲性能分析2.2.1抗噪聲性能分析原理AM線性調(diào)制系統(tǒng)的相干解調(diào)模型如下圖所示。圖3.5.1 線性調(diào)制系統(tǒng)的相干解調(diào)模型圖中可以是AM調(diào)幅信號(hào)，帶通濾波器的帶寬等于已調(diào)信號(hào)帶寬10。下面討論AM調(diào)制系統(tǒng)的抗噪聲性能11。AM信號(hào)的時(shí)域表達(dá)式為 通過(guò)分析可得AM信號(hào)的平均功率為 又已知輸入功率, 其中B表示已調(diào)信號(hào)的帶寬。由此可得AM信號(hào)在解調(diào)器的輸入信噪比為 AM信號(hào)經(jīng)相干解調(diào)器的輸出信號(hào)為 因此解調(diào)后輸出信號(hào)功率為 在上圖中輸入噪聲通過(guò)帶通濾波器之后，變成窄帶噪聲，經(jīng)乘法器相乘后的輸出噪聲為 經(jīng)LPF后， 因此解調(diào)器的輸出噪聲功率為 可得AM信號(hào)經(jīng)過(guò)解調(diào)器

11、后的輸出信噪比為 由上面分析的解調(diào)器的輸入、輸出信噪比可得AM信號(hào)的信噪比增益為 2.2.2 調(diào)試過(guò)程clf; %清除窗口中的圖形t=0:0.01:2; %定義變量區(qū)間fc=50; %給出相干載波的頻率A=10; %定義輸入信號(hào)幅度f(wàn)a=5; %定義調(diào)制信號(hào)頻率mt=A*cos(2*pi*fa.*t); %輸入調(diào)制信號(hào)表達(dá)式xzb=5; %輸入小信躁比(dB)snr=10.(xzb/10); db=A2./(2*snr); %由信躁比求方差nit=sqrt(db).*randn(size(mt); %產(chǎn)生小信噪比高斯白躁聲psmt=(A+mt).*cos(2*pi*fc.*t); %輸出調(diào)制信

12、號(hào)表達(dá)式psnt=psmt+nit; %輸出疊加小信噪比已調(diào)信號(hào)波形xzb1=30; %輸入大信躁比(dB)snr1=10.(xzb1/10); db1=A2./(2*snr1); %由信躁比求方差nit1=sqrt(db1).*randn(size(mt) ); %產(chǎn)生大信噪比高斯白躁聲psnt1=psmt+nit1; %輸出已調(diào)信號(hào)波形subplot(2,2,1); %劃分畫(huà)圖區(qū)間plot(t,nit,'g'); %畫(huà)出輸入信號(hào)波形title('小信噪比高斯白躁聲');xlabel(' t');ylabel(' nit');s

13、ubplot(2,2,2);plot(t,psnt,'b');title('疊加小信噪比已調(diào)信號(hào)波形');xlabel('時(shí)間');ylabel('輸出調(diào)制信號(hào)');subplot(2,2,3);plot(t,nit1,'r'); %length用于長(zhǎng)度匹配title('大信噪比高斯白躁聲'); %畫(huà)出輸入信號(hào)與噪聲疊加波形xlabel(' t');ylabel('nit');subplot(2,2,4);plot(t,psnt1,'k');title

14、('疊加大信噪比已調(diào)信號(hào)波形'); %畫(huà)出輸出信號(hào)波形xlabel('時(shí)間');ylabel('輸出調(diào)制信號(hào)');運(yùn)行結(jié)果：由上圖可見(jiàn)，當(dāng)輸入信號(hào)一定時(shí)，隨著噪聲的加強(qiáng)，接收端輸入信號(hào)被干擾得越嚴(yán)重。而相應(yīng)的輸出波形相對(duì)于發(fā)送端的波形誤差也越大。而當(dāng)噪聲過(guò)大時(shí)信號(hào)幾難分辨。這信噪比變小導(dǎo)致的在實(shí)際的信號(hào)傳輸過(guò)程中當(dāng)信道噪聲過(guò)大將會(huì)導(dǎo)致幅度相位等各失真當(dāng)然由于非線性元件如濾波器等的存在。非線性失真也會(huì)隨噪聲加大而變大。三、SSB調(diào)制與解調(diào)及抗噪聲性能分析 3.1 SSB調(diào)制與解調(diào)原理單邊帶調(diào)制信號(hào)是將雙邊帶信號(hào)中的一個(gè)邊帶濾掉而形成的。根據(jù)方法的不

15、同，產(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)為載波為則其雙邊帶信號(hào)DSB信號(hào)的時(shí)域表示式為若保留上邊帶，則有若保留下邊帶，則有將上兩式合并得：由希爾伯特變換故單邊帶信號(hào)經(jīng)過(guò)希爾伯特變換后得:把上式推廣到一般情況，則得到 式中若M(w)是m(t)的傅里葉變換，則上式中的-jsgnw可以看作是希爾伯特濾波器傳遞函數(shù)，即移相法SSB調(diào)制器方框圖相移法是利用相移網(wǎng)絡(luò)，對(duì)載波和調(diào)制信號(hào)進(jìn)行適當(dāng)?shù)南嘁疲员阍诤铣蛇^(guò)程中將其中的一個(gè)邊帶抵消而獲得SSB信號(hào)。相移法不需要濾波器具

16、有陡峭的截止特性，不論載頻有多高，均可一次實(shí)現(xiàn)SSB調(diào)制。SSB信號(hào)的解調(diào) SSB信號(hào)的解調(diào)不能采用簡(jiǎn)單的包絡(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ā)射功率，而且它所占用的頻帶寬度比AM、DSB減少了一半。它目前已成為短波通信中一種重要的調(diào)制方式。3.2 SSB調(diào)制解調(diào)系統(tǒng)抗噪聲性能分析 噪聲功率這里，B = fH 為SSB 信號(hào)的帶通濾波器的帶寬。信號(hào)功率SSB信號(hào)與相干載波相乘后，再經(jīng)低通濾波可得解調(diào)器輸出信號(hào)因此，輸出信號(hào)平均功率輸入

17、信號(hào)平均功率為單邊帶解調(diào)器的輸入信噪比為單邊帶解調(diào)器的輸出信噪比為制度增益因?yàn)樵赟SB系統(tǒng)中，信號(hào)和噪聲有相同表示形式，所以相干解調(diào)過(guò)程中，信號(hào)和噪聲中的正交分量均被抑制掉，故信噪比沒(méi)有改善。3.3 調(diào)試過(guò)程先建立3個(gè)M文件1. afd_buttfunction b,a = afd_butt(Wp,Ws,Rp,As);if Wp <= 0 error('Passband edge must be larger than 0')endif Ws <= Wp error('Stopband edge must be larger than Passband ed

18、ge')endif (Rp <= 0) | (As < 0) error('PB ripple and/or SB attenuation ust be larger than 0')endN = ceil(log10(10(Rp/10)-1)/(10(As/10)-1)/(2*log10(Wp/Ws);fprintf('n* Butterworth Filter Order = %2.0f n',N)OmegaC = Wp/(10(Rp/10)-1)(1/(2*N);b,a=u_buttap(N,OmegaC); 2. u_buttapfu

19、nction b,a = u_buttap(N,Omegac);z,p,k = buttap(N);p = p*Omegac;k = k*OmegacN;B = real(poly(z);b0 = k;b = k*B;a = real(poly(p);3. imp_invr%脈沖響應(yīng)不變法子程序function b,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采樣頻率>>

20、 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*cos(150*pi*n); %調(diào)制信號(hào)>> x2=hilbert(x1,N); %對(duì)x1做希爾伯特變換>> y=(Vm0*x1.*cos(2*pi*fc*n)-Vm0*x2.*sin(2*pi*fc*n)/2; %保留上邊帶的已調(diào)波信號(hào)>> xzb=2; %輸入小信噪比(d

21、B)>> snr=10.(xzb/10);>> h,l=size(x1); %求調(diào)制信號(hào)的維度>> fangcha=A*A./(2*snr); %由信噪比求方差>> nit=sqrt(fangcha).*randn(h,l); %產(chǎn)生高斯白噪聲>> yn=y+nit; %疊加小信噪比噪聲的已調(diào)波信號(hào)>> xzb=10; %輸入小信噪比(dB)>> snr1=10.(xzb/10);>> h,l=size(x1); %求調(diào)制信號(hào)的維度>> fangcha=A*A./(2*snr1); %由

22、信噪比求方差>> nit1=sqrt(fangcha).*randn(h,l); %產(chǎn)生高斯白噪聲>> yn1=y+nit1; %疊加小信噪比噪聲的已調(diào)波信號(hào)>> figure(1)>> subplot(2,2,1) %劃分畫(huà)圖區(qū)間>> plot(n,x1) %畫(huà)出調(diào)制信號(hào)的波形>> title('調(diào)制信號(hào)');>> subplot(2,2,2) %劃分畫(huà)圖區(qū)間>> plot(n,y) %畫(huà)出已調(diào)波信號(hào)波形>> title('已調(diào)波信號(hào)');>>

23、; subplot(2,2,3) %劃分畫(huà)圖區(qū)間>> plot(n,yn) %畫(huà)出疊加噪聲的已調(diào)波信號(hào)波形>> title('疊加小信噪比噪聲的已調(diào)波信號(hào)');>>subplot(2,2,4) %劃分畫(huà)圖區(qū)間>> plot(n,yn1) %畫(huà)出疊加噪聲的已調(diào)波信號(hào)波形>> title('疊加大信噪比噪聲的已調(diào)波信號(hào)');>> X=fft(x1); %調(diào)制信號(hào)x1的傅里葉變換>> Y=fft(y); %已調(diào)信號(hào)y的傅里葉變換>> Yn=fft(yn); %疊加小信噪比噪

24、聲的已調(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) %劃分畫(huà)圖區(qū)間>> plot(w,abs(X) %畫(huà)出調(diào)制信號(hào)頻譜波形>> axis(0,pi/4,0,3000); %給出橫縱坐標(biāo)的圍>> title('調(diào)制信號(hào)頻譜');>> subplot(2,2,2) %劃分畫(huà)圖區(qū)間>> plot(w,ab

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

26、(Yn1) %畫(huà)出疊加大信噪比噪聲的已調(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=60/N*pi;Rp=1;As=15;T=1; %濾波器參數(shù)設(shè)計(jì)>> OmegaP=wp/T;OmegaS=ws/T;>> cs,

27、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) %劃分畫(huà)圖區(qū)間>> plot(n,y) %畫(huà)出解調(diào)波波形>> title('解調(diào)波');>> Y=fft(y); %解調(diào)波y的傅里葉變換>> subplot(2,1,2) %劃分畫(huà)圖區(qū)間>> plot(w,abs(Y) %畫(huà)出解調(diào)信號(hào)頻譜>> axis(0,pi/6,0,2500);