數(shù)字通信系統(tǒng)課程設(shè)計

1、三峽大學(xué)數(shù)字通信系統(tǒng)課程設(shè)計基于matlab的2psk調(diào)制與解調(diào)系統(tǒng)學(xué) 院： 計算機與信息學(xué)院 學(xué) 號： 2010156203 班級序號： 72 姓 名： 蔡紅亮 指導(dǎo)老師： 夏平 日 期： 2013年1月2日 數(shù)字通信系統(tǒng)課程設(shè)計任務(wù)書題 目基于matlab的2psk調(diào)制與解調(diào)系統(tǒng)設(shè) 計目 的1、綜合應(yīng)用信號與系統(tǒng)、隨機信號分析、通信原理等課程知識，使學(xué)生建立通信系統(tǒng)的整體概念；2、培養(yǎng)學(xué)生系統(tǒng)設(shè)計與系統(tǒng)開發(fā)的思想；3、培養(yǎng)學(xué)生利用matlab軟件平臺進行通信系統(tǒng)仿真的能力。設(shè) 計要 求1、每人一組，完成指定設(shè)計問題的系統(tǒng)詳細功能；2、對通信系統(tǒng)有整體的較深入的理解，深入理解設(shè)計系統(tǒng)的基本原

2、理，畫出組成的每個通信子系統(tǒng)的原理框圖；3、提出設(shè)計方案；4、完成仿真軟件的編制；5、在matlab仿真平臺上，實現(xiàn)仿真系統(tǒng)正常運行；6、提交詳細的設(shè)計報告。主 要設(shè) 備計算機1臺，matlab仿真軟件一套。主要參考文獻1 張輝，曹麗娜等. 現(xiàn)代通信原理與技術(shù)（第三版）.西安電子科技大學(xué)出版社，2009.2 周炯槃等. 通信原理(第三版). 北京郵電大學(xué)出版社，2005. 3 樊昌信. 通信原理（第五版）.國防工業(yè)出版社，2006.4 曹志剛等.現(xiàn)代通信原理. 清華大學(xué)出版社，1992.5 leon w. couch. digital and analog communication syst

3、ems(影印第5版). 清華大學(xué)出版社，1999 .6 stephen g. wilson. digital modulation and coding（影印版）.電子工業(yè)出版社，1998.課程設(shè)計進度安排（起止時間、工作內(nèi)容）目 錄一、課程設(shè)計主要內(nèi)容二、系統(tǒng)設(shè)計的原理三、課程設(shè)計思路及過程四、課程設(shè)計實驗結(jié)果五、課程設(shè)計結(jié)果分析六、參考文獻 七、附錄1、 課程設(shè)計的主要內(nèi)容了解并掌握數(shù)字頻帶傳輸系統(tǒng)2psk的調(diào)制與解調(diào)系統(tǒng)，通過matlab來實現(xiàn)整個通信系統(tǒng)的設(shè)計，進一步加深對2psk系統(tǒng)的理解，并熟練掌握matlab語言。2、 課程設(shè)計原理數(shù)字信號的傳輸方式分為基帶傳輸和帶通傳輸，在實際

4、應(yīng)用中，大多數(shù)信道具有帶通特性而不能直接傳輸基帶信號。為了使數(shù)字信號在帶通信道中傳輸，必須使用數(shù)字基帶信號對載波進行調(diào)制，以使信號與信道的特性相匹配。這種用數(shù)字基帶信號控制載波，把數(shù)字基帶信號變換為數(shù)字帶通信號的過程稱為數(shù)字調(diào)制。數(shù)字調(diào)制技術(shù)的兩種方法：利用模擬調(diào)制的方法去實現(xiàn)數(shù)字式調(diào)制，即把數(shù)字調(diào)制看成是模擬調(diào)制的一個特例，把數(shù)字基帶信號當(dāng)作模擬信號的特殊情況處理；利用數(shù)字信號的離散取值特點通過開關(guān)鍵控載波，從而實現(xiàn)數(shù)字調(diào)制。這種方法通常稱為鍵控法，比如對載波的相位進行鍵控，便可獲得相移鍵控（psk）基本的調(diào)制方式。 乘法器)(2tspsk雙極性不歸零t0wcos)(tf碼型變換單/雙 圖-

5、1模擬調(diào)制的方法tcwcos)(tf)(2tspsk開關(guān)電路移相01800p 圖-2相移鍵控方法調(diào)制過程：數(shù)字調(diào)相：如果兩個頻率相同的載波同時開始振蕩，這兩個頻率同時達到正最大值，同時達到零值，同時達到負最大值，它們應(yīng)處于同相狀態(tài)；如果其中一個開始得遲了一點，就可能不相同了。如果一個達到正最大值時，另一個達到負最大值，則稱為反相。一般把信號振蕩一次（一周）作為360度。如果一個波比另一個波相差半個周期，我們說兩個波的相位差180度，也就是反相。當(dāng)傳輸數(shù)字信號時，1碼控制發(fā)0度相位，0碼控制發(fā)180度相位。載波的初始相位就有了移動，也就帶上了信息。相移鍵控：利用載波的相位變化來傳遞數(shù)字信息，而振

6、幅和頻率保持不變。在2psk中，通常用初始相位0和分別表示二進制“1”和“0”。因此，2psk信號的時域表達式為 (t)=acos其中，表示第n個符號的絕對相位：= 因此，上式可以改寫為 當(dāng)碼元寬度和載波周期相等時，2psk調(diào)制波形如下： 圖-3 2psk調(diào)制波形解調(diào)過程2psk信號的解調(diào)方法是相干解調(diào)法。由于psk信號本身就是利用相位傳遞信息的，所以在接收端必須利用信號的相位信息來解調(diào)信號。下圖2-3中給出了一種2psk信號相干接收設(shè)備的原理框圖。圖中經(jīng)過帶通濾波的信號在相乘器中與本地載波相乘，然后用低通濾波器濾除高頻分量，在進行抽樣判決。判決器是按極性來判決的。即正抽樣值判為1，負抽樣值判

7、為0.圖-4 2psk相干解調(diào)系統(tǒng)框圖0101sttabcd0tttte00110 圖-5 2psk信號相干解調(diào)各點時間波形可見，2psk信號相干解調(diào)的過程實際上是輸入已調(diào)信號與本地載波信號進行極性比較的過程，故常稱為極性比較法解調(diào)。 由于2psk信號實際上是以一個固定初相的末調(diào)載波為參考的，因此，解調(diào)時必須有與此同頻同相的同步載波。如果同步載波的相位發(fā)生變化，如0相位變?yōu)橄辔换蛳辔蛔優(yōu)?相位，則恢復(fù)的數(shù)字信息就會發(fā)生“0”變“1”或“1”變“0”，從而造成錯誤的恢復(fù)。這種因為本地參考載波倒相，而在接收端發(fā)生錯誤恢復(fù)的現(xiàn)象稱為“倒”現(xiàn)象或“反向工作”現(xiàn)象。絕對移相的主要缺點是容易產(chǎn)生相位模糊，

8、造成反向工作。這也是它實際應(yīng)用較少的主要原因。 三、課程設(shè)計思路及過程本課題主要是通過matlab來實現(xiàn)整個2psk系統(tǒng)，在用matlab語言編程的時候主要分為以下幾個模塊。模塊一、通過碼型變換產(chǎn)生雙極性數(shù)字基帶信號，首先調(diào)用matlab函數(shù)產(chǎn)生一組隨機序列，然后將離散的隨機序列變換成單極性的二進制數(shù)字基帶信號，將單極性信號取反并與上面信號相減即可得到我們所需要的二進制雙極性數(shù)字基帶信號，生成的二進制信號帶寬為1。模塊二、載波信號的產(chǎn)生，采用的載頻為10hz，所以基帶信號周期是載波周期的10倍，載波振幅為1。模塊三、通過模擬調(diào)制的方法得到2psk信號：（ 注：在這設(shè)置載波信號的初相為0）模塊四

9、、已調(diào)信號通過信道傳輸，通過matlab函數(shù)產(chǎn)生一個隨機信號來模擬高斯白噪聲，并與已調(diào)信號相加，在這里我們只考慮加性噪聲，不考慮乘型噪聲。模塊五、相干解調(diào)過程（注：這里沒有載波提取過程，而是直接用本地載波去與已調(diào)信號相乘）模塊六、用matlab來設(shè)計理想的巴特沃斯低通濾波器模塊七、根據(jù)位同步信號進行抽樣判決模塊八、2psk信號的頻譜分析模塊九、采用窗函數(shù)方法對信號進行功率譜密度分析4、 課程設(shè)計的試驗結(jié)果1、2psk調(diào)制、解調(diào)系統(tǒng)信號波形圖圖-6二進制數(shù)字基帶信號與載波信號 圖-7 2psk信號與接收機接收到的信號波形圖-8通過解調(diào)、低通輸出波形圖-9判決輸出波形將抽樣輸出的基帶信號波形與原信

10、號相比較可以發(fā)現(xiàn)兩個信號波形完全相同，說明整個2psk的調(diào)制、解調(diào)系統(tǒng)是正確的。五、課程設(shè)計結(jié)果分析1、2psk調(diào)制解調(diào)系統(tǒng)頻譜分析數(shù)字基帶信號為二進制非歸零碼，所以它的傅里葉變換應(yīng)該是門函數(shù)，信號帶寬為1/ts，其頻譜圖-9。載波信號為正弦信號，載頻為10hz，其頻譜如下。圖-9 數(shù)字基帶信號與載波信號頻譜2psk信號的實質(zhì)為基帶信號與載波相乘的結(jié)果，所以2psk的信號頻譜是基帶信號頻譜的線性搬移，中心頻率為fc.，帶寬為基帶信號帶寬的兩倍，如圖-10，2psk信號通過信道后，加上了高斯白噪聲，其頻譜如下。圖-10 2psk信號與接收機接收的信號頻譜相干解調(diào)實質(zhì)是將已調(diào)信號與提取的載波信號相

11、乘，所以解調(diào)輸出的信號的頻率成分包括基帶信號頻率和2倍的載頻信號和噪聲如圖-11，通過低通濾波器后，高頻成分被慮掉，剩下基帶信號和部分噪聲。圖-11 解調(diào)輸出與低通輸出的信號頻譜圖-12 判決輸出信號頻譜與基帶信號頻譜通過比較判決輸出信號的頻譜和基帶信號頻譜，得出結(jié)論，2psk的調(diào)制解調(diào)系統(tǒng)設(shè)計是正確的。2、 各階段信號功率譜密度分析 圖-12 2psk功率譜密度3、“倒”現(xiàn)象分析 在相干解調(diào)之前，需要根據(jù)已調(diào)信號來提取載波信號，如果同步載波的相位發(fā)生變化，則恢復(fù)的數(shù)字信息就會發(fā)生“0”變“1”或“1”變“0”，從而造成錯誤的恢復(fù)。這種因為本地參考載波倒相。所以在實際應(yīng)用中多采用2dpsk調(diào)制

12、解調(diào)系統(tǒng)，可以防止“倒”現(xiàn)象的發(fā)生。在matlab程序中，只用將提取的載波加上90度的相位偏移即可看到“倒”現(xiàn)象。 六、參考文獻1 張輝，曹麗娜等. 現(xiàn)代通信原理與技術(shù)（第三版）.西安電子科技大學(xué)出版社，2009.2 周炯槃等. 通信原理(第三版). 北京郵電大學(xué)出版社，2005. 3 樊昌信. 通信原理（第五版）.國防工業(yè)出版社，2006.4 曹志剛等.現(xiàn)代通信原理. 清華大學(xué)出版社，1992.5 leon w. couch. digital and analog communication systems(影印第5版). 清華大學(xué)出版社，1999 .6 stephen g. wilson.

13、 digital modulation and coding（影印版）.電子工業(yè)出版社，1998. 七、附錄 程序代碼實現(xiàn)：clear allclose alli=20;j=5000;t=linspace(0,20,5000);%在0到20之間等分成5000個點% % 功能實現(xiàn)模塊一：產(chǎn)生隨機二進制數(shù)字基帶信號 %rem=round(rand(1,i);%調(diào)用函數(shù)產(chǎn)生隨機序列signal1=t;%將離散的的隨機序列轉(zhuǎn)變成連續(xù)的單極性二進制碼元for n=1:20 if rem(n)=1; signal2(k)=0; else signal2(k)=1; endend;signal3=signa

14、l1-signal2;%將兩個單極性數(shù)字基帶信號相減即可獲得雙極性二進制數(shù)字基帶信號subplot(312);plot(250*t,signal3);title(雙極性二進制數(shù)字基帶信號);axis(0,5000,-2,2);grid on;% 功能實現(xiàn)模塊二：產(chǎn)生載波信號 % 載波周期與碼元寬度相同且均為1 %fc=10;%載波頻率為10hzs=sin(2*pi*fc*t); subplot(313);plot(s);title(載波信號波形);axis(0 5000 -1 1);grid on;% 功能實現(xiàn)模塊三：2psk調(diào)制部分:基帶信號與載波相乘 %psk=signal3.*s;fig

15、ure(2);subplot(211);plot(t*250,psk);title(已調(diào)信號波形);axis(0,5000,-2,2);% 功能實現(xiàn)模塊四：調(diào)制信號通過信道傳輸產(chǎn)生加性噪聲 %noise=rand(1,j);y1=psk+noise;%模擬加性噪聲subplot(212);plot(y1);axis(0,5000,-1,2);title(加噪后的波形);% 功能實現(xiàn)模塊五：相干解調(diào)實現(xiàn)過程 %y2=y1.*s;%接收機收到信號輸入乘法器與載波信號相乘figure(3)subplot(211);plot(y2);title(解調(diào)輸出波形);grid;% 功能實現(xiàn)模塊六：通過巴特沃

16、斯低通濾波器 %fp=100;fs=200;rp=3;rs=20;fn=10000;%設(shè)計巴特沃斯低通濾波器的參數(shù) ws=fs/(fn/2); wp=fp/(fn/2);%計算歸一化角頻率n,wn=buttord(wp,ws,rp,rs);%計算階數(shù)和截止頻率b,a=butter(n,wn);y3=2*filter(b,a,y2);subplot(212);plot(y3);title(經(jīng)低通濾波器后的輸出信號波形);grid on;% 功能實現(xiàn)模塊七：根據(jù)位同步信號進行抽樣判決 %y4=linspace(0,19);y4=0;for n=1:20; if y3(n*250)0; y4(n)=

17、1; else y4(n)=-1; endendfigure(4)subplot(211);stem(y4);axis(0 20 -2 2);title(抽樣后得到的離散采樣值);grid on;xlabel(時間/s);y5=t;for n=1:20 if y4(n)1; for m=j/i*(n-1)+1:j/i*n y5(m)=-1; end else for m=j/i*(n-1)+1:j/i*n y5(m)=1; end endendsubplot(212);plot(y5);axis(0 5000 -2 2);grid on;title(抽樣判決輸出基帶信號);xlabel(采樣點

18、數(shù));% 功能實現(xiàn)模塊八：各階段信號頻譜分析 %fs=250;n=1024; %采樣頻率和數(shù)據(jù)點數(shù)n=0:n-1;t=n/fs; %時間序列y=fft(signal3,n);%對信號進行快速fourier變換mag=abs(y);%求取fourier變換的振幅f=n*fs/n;figure(5);subplot(211)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(頻率/hz);ylabel(振幅);title(基帶信號頻譜（n=1024）);grid on; fs=250;n=1024;n=0:n-1;t=n/fs;y=fft(s,n);%

19、對信號進行快速fourier變換mag=abs(y);%求取fourier變換的振幅f=n*fs/n;subplot(212)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(頻率/hz);ylabel(振幅);title(載波信號頻譜（n=1024);grid on; fs=250;n=1024;n=0:n-1;t=n/fs;y=fft(psk,n);%對信號進行快速fourier變換mag=abs(y);%求取fourier變換的振幅f=n*fs/n;figure(6)subplot(211)plot(f(1:n/2),mag(1:n/2);

20、axis(0 50 0 500);xlabel(頻率/hz);ylabel(振幅);title(2psk信號頻譜（n=1024);grid on; fs=250;n=1024;n=0:n-1;t=n/fs;y=fft(y1,n);%對信號進行快速fourier變換mag=abs(y);%求取fourier變換的振幅f=n*fs/n;subplot(212)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(頻率/hz);ylabel(振幅);title(接收機接收到的信號頻譜（n=1024);grid on; fs=250;n=1024;n=0:n

21、-1;t=n/fs;y=fft(y2,n);%對信號進行快速fourier變換mag=abs(y);%求取fourier變換的振幅f=n*fs/n;figure(7)subplot(211)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(頻率/hz);ylabel(振幅);title(解調(diào)輸出信號頻譜（n=1024);grid on; fs=250;n=1024;n=0:n-1;t=n/fs;y=fft(y3,n);%對信號進行快速fourier變換mag=abs(y);%求取fourier變換的振幅f=n*fs/n;subplot(212)p

22、lot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(頻率/hz);ylabel(振幅);title(低通輸出信號頻譜（n=1024);grid on; fs=250;n=1024;n=0:n-1;t=n/fs;y=fft(y5,n);%對信號進行快速fourier變換mag=abs(y);%求取fourier變換的振幅f=n*fs/n;figure(8)subplot(211)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(頻率/hz);ylabel(振幅);title(抽樣判決輸出信號頻譜（n=10

23、24);grid on; fs=250;n=1024; %采樣頻率和數(shù)據(jù)點數(shù)n=0:n-1;t=n/fs; %時間序列y=fft(signal3,n);%對信號進行快速fourier變換mag=abs(y);%求取fourier變換的振幅f=n*fs/n;subplot(212)plot(f(1:n/2),mag(1:n/2);axis(0 50 0 500);xlabel(頻率/hz);ylabel(振幅);title(基帶信號頻譜（n=1024）);grid on;% 功能實現(xiàn)模塊九：各階段信號功率譜密度分析 %fs=1000; n=0:1/fs:1; nfft=1024; window=hamming(100); %海明窗noverlap=20; %數(shù)據(jù)無重疊 range=half;