【正文】 Matlab 應(yīng)用》課程設(shè)計說明書 24 axis( [ 0 ,*Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 figure(2)。 title(39。%信號調(diào)制中加入白噪聲，信噪比為 10 %相乘 dt = st .* cos ( 2 * pi * w * t )。 axis( [ 0 , *Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 grid figure(2)。 title(39。s39。)。 subplot(3,2,4)。低通濾波后的波形 39。 plot(abs(fft(dt)))。) grid %抽樣判決 amp。 figure(1) subplot(3,2,5)。逆碼變換后的波形 39。 plot(abs(fft(sign_result)))。) grid %誤碼率隨信噪比變化曲線 snr=0:2:20。%DPSK 系統(tǒng)的誤碼率與信噪比的關(guān)系 end figure(3)。 grid on。 。誤碼率隨信噪比的變化曲線 39。ko39。 for i=1:len_snr SNR_eb=exp(snr(i)*log(10)/10)。 title(39。 figure(2)。 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計說明書 25 axis( [ 0 , *Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 else sign_result( (I1)*Num_Unit + 1 : I*Num_Unit) = Low_Level。 title(39。 figure(2)。 axis( [ 0 ,* Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 %映射為數(shù)字的 dt = filter(bz,az,dt)。 %臨界頻率采用角頻率表示 [b,a]=butter(N,Wn,39。)。 plot(abs(fft(dt)))。相乘后的波形 39。 subplot(3,2,3)。)。 plot(abs(fft(st)))。調(diào)制后的信號 39。 subplot ( 3, 2, 2 )。 title(39。 figure(2)。 axis( [ 0 , *Time_Hold_On *( Lenth_Of_Sign + 1), (A / 2), A + (A / 2) ] )。 else st( (I1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w * t( (I1)*Num_Unit + 1 : I*Num_Unit ) )。 t = 0 : 10/fs : 10*Time_Hold_On * Lenth_Of_Sign 10/fs。 Lenth_Of_Sign = length ( Sign_Set )。 w =100000。 Time_Hold_On =。 title(39。 semilogy(snr,theo_err_prb,39。%信噪比范圍 len_snr = length(snr)。 title(39。) figure(2)。 plot(jt)。 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計說明書 22 else jt(m)=0。)。 plot(abs(fft(jt)))。 title(39。 figure(1)。 wp=fp/(fn/2)。rp=3。相乘后信號頻譜 39。 subplot(3,2,4)。grid on axis([0 15*length(g) 2 2])。 jiet=2*mod1.*tz。 axis([0 15*length(g) 0 1000])。)。 plot(tz)。)。 plot(abs(fft(tiaoz)))。PCM 編碼信號頻譜 39。 subplot(3,2,1)。title(39。 subplot(3,2,2)。 %載波信號 end cm=[cm B]。%每個值 200 個點 c=cos(f*t)。mod=[]。title(39。subplot(3,2,1)。 end cp=[cp A]。 if g(n)==0。 mod1=[]。 clf。 elseif z(i)0 code(i,1)=0。 code(i,7)=str2double(k(3))。 %求出 在段內(nèi)的位置 if q==0 code(i,(5:9))=[0,0,0,0,0]。 for i=1:length(S) N(i)=bin2dec(num2str(code(i,2:4)))+1。(Q(i)20xx)||(Q(i)=4064)amp。(Q(i)96)||(Q(i)=224)amp。amp。(Q(i)=8160) code(i,2)=1。 %歸一化 Q=8160*S。)。)。 Sign_m = sin(2*pi*300*ts)。在實驗課上學(xué)得的，我將發(fā)揮到其它 應(yīng)用 中去，也將在今后的學(xué)習(xí)和工作中不斷提高 和 完善；在 做課程設(shè)計時發(fā)現(xiàn)的不足，我將努力改善，通過學(xué)習(xí) 和 實踐等方式 ，不斷提高自己的綜合 能力！ 5 參考文獻(xiàn) [1]樊昌信，曹麗娜 .《通信原理》 .第六版 .國防工業(yè)出版社， 20xx [2]葛哲學(xué) .《精通 MATLAB》 .電子工業(yè)出版社 .20xx [3]韓利竹，王華 .《基于 MATLAB 的通信系統(tǒng)仿真》 .國防工業(yè)出版社 .20xx [4]Vinay .《數(shù)字信號處理及其 MATLAB 實現(xiàn)》 .電子工業(yè)出版社 .1998 [5]Sanjit K. Miltra.《 Digital Signal Processing Laboratory Using Matlab》 .McGrawHill 出版社 .20xx 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計說明書 18 附錄 1 %PCM 編碼程序 function code2=pcm(S) %S 為輸入信號 t = 0::。在 PSK 和 DPSK 的調(diào)制解調(diào)程序的編寫過程中，我在圖書館 和網(wǎng)上查了不少資料，借用了別人的程序，并進(jìn)行了部分修改和完善，最終仿真成功。本次課設(shè)要求我們運用 Matlab 這個專業(yè)的學(xué)習(xí)軟件將課本上的 知識進(jìn)行仿真，從而更透徹地理解課堂上學(xué)到的理論知識。第五個是經(jīng)抽樣判決及逆碼變換后的波形，也就是 DPSK 解調(diào)后的波形，由圖中可以看出它與第一個 PCM 編碼信號 序列 的波形是完全一致的，故是符合要求的。 圖 35 PSK各信號的頻譜仿真結(jié)果 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計說明書 14 圖 36 PSK信號誤碼率隨信噪比的變化曲線 DPSK 調(diào)制解調(diào)仿真分析 圖 37 DPSK調(diào)制解調(diào)仿真結(jié)果 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計說明書 15 如圖 37 所示，第一個為 PCM 編碼信號序列，與 PCM 編碼結(jié)果相符。第三個是通過高斯白噪聲信道后的PSK 信號波形，通過將其與調(diào)制后的 PSK 信號作比較可以看出，此波形較上個波形要模糊一些，故符合題目要求。由圖可以看出，編碼結(jié)果正確，與實際相符。由分析可知， Pe=erfc( r ) （公式 218） 用 Matlab 實現(xiàn) DPSK 系統(tǒng)的調(diào)制解調(diào) 用 Matlab 實現(xiàn) 2DPSK 信號 的調(diào)制解調(diào)的程序見附錄 3。其原理框圖如圖 13 所示。就是在 PSK 調(diào)制原理框圖前加了一個差分編碼器。 于是可以將一組二進(jìn)制數(shù)字信息與其對應(yīng)的 DPSK 信號的相位關(guān)系示例如下： 數(shù)字信息： 1 1 0 1 0 0 1 1 0 DPSK 相位：（ ? ） 0 ? ? 0 0 0 ? 0 0 或 （ 0） ? 0 0 ? ? ? 0 ? ? 由此示例可知，對于相同的基帶數(shù)字信息序列，由于序列初始碼元的參考相位不同，DPSK 信號的相位可以不同。在此過程中，我用函數(shù) randn 來添 加噪聲，此函數(shù)功能為向信號中添加噪聲功率為其方差的高斯白噪聲。 （ 3） PSK 相干解調(diào)系統(tǒng)性能分析 PSK 相干解調(diào)方式又稱為極性比較法，其性能分析模型如圖 24 所示。 但是，由于在 2PSK 信號的載波恢復(fù)過程中存在著的相位模糊，即恢復(fù)的本地載波與所需的相干載波可能同相，也可能反相，這種相位關(guān)系的不確定性將會造成解調(diào)出的數(shù)字基帶信號與發(fā)送的數(shù)字基帶信號正好相反，即 “ 1” 變?yōu)?“ 0” ， “ 0” 變?yōu)?“ 1” ，判決器輸出數(shù)字信號全部出錯。其調(diào)制原理框圖如圖 11 所示。 表 23段內(nèi)碼 量化間隔 段內(nèi)碼c5c6c7c8c9 量化間隔 段內(nèi)碼c5c6c7c8c9 31 11111 15 01111 30 11110 14 01110 29 11101 13 01101 28 11100 12 01100 27 11011 11 01011 26 11010 10 01010 25 11001 9 01001 24 11000 8 01000 23 10111 7 00111 22 10110 6 00110 21 10101 5 00101 20 10100 4 00100 19 10011 3 00011 18 10010 2 00010 17 10001 1 00001 16 10000 0 00000 武漢理工大學(xué)《 Matlab 應(yīng)用