【正文】 通濾波器后信號波形 39。%計(jì)算歸一化角頻率 [n,wn]=buttord(wp,ws,rp,rs)。)。 title(39。 title(39。 grid on axis([0 15*length(g) 2 2])。 axis([0 15*length(g) 0 1000])。 plot(abs(fft(cp)))。 plot(tiaoz)。 %載波信號 else g(n)==1。PCM 編碼信號序列 39。 %s(t),碼元寬度 200 c=cos(f*t)。 f=2*100000*pi。 end %符號位的判斷 end for i=1:length(S) for j=1:9 code2(9*(i1)+j)=code(i,j)。 %如果輸入為零則輸出 “0” 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計(jì)說明書 19 else k=num2str(dec2bin(q1,5))。amp。(Q(i)=8160) code(i,3)=1。 %量 化 code=zeros(length(S),9)。 subplot(2,1,2) %抽樣 stem(ts,Sign_m)。 fs = 3000。但是由于我的理論知識沒有很好地掌握，對于仿真軟件也用的不太熟練，因此這次的課設(shè)做起來也沒有想象中輕松容易。 第二個(gè)為 DPSK調(diào)制信號的波形，將其與前面的二進(jìn)制序列作比較可以看出其符合要求。 PSK 調(diào)制解調(diào)仿真分析 圖 33 PSK調(diào)制解調(diào)仿真結(jié)果 為了使 每個(gè)圖形相對應(yīng)，方便看出波形調(diào)制解調(diào)前后及通過高斯白噪聲信道前后的波形變化， 圖 33 中將每個(gè)圖像的比例設(shè)置為相同值 ， 此時(shí) 2PSK 調(diào)制信號波形過于緊密。 （ 3） DPSK 相干解調(diào)系統(tǒng)性能分析 2DPSK 原相干解調(diào)法，又稱極性比較碼反變換法，其解調(diào)原理是：對 2DPSK 信號進(jìn)行相干解調(diào)，恢復(fù)出相對碼序列 {bn},再通過碼反變換器變換為絕對碼序列 {an}，從而恢復(fù)也發(fā)送的二進(jìn)制數(shù)字信息。也就是說， DPSK 信號的相位并不直接代表基帶信號，而前后碼元相對相位的差才唯一決定信息符號。 帶通 濾波器 相乘器 低通 濾波器 抽樣 判決器 a c d e b cosWctwcos epsk(t) 定時(shí) 脈沖 輸出 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計(jì)說明書 8 圖 24 PSK信號相干解調(diào)系統(tǒng)性能分析模型 當(dāng)信號經(jīng)過信道傳輸時(shí)會受到噪聲的影響，這是不可避免的。 （ 2） PSK 解調(diào)過程分析 解調(diào)原理框圖如圖 23 所示。其中第 2 到 4 位（ c2c3c4）是段落碼，共計(jì) 3 位，可以表示 8 種斜率的段落；其他 5 位（ c5~c9）為段內(nèi)碼，可以表示每一段落內(nèi)的 32 種量化電平。每比較一次，得出 1 位二進(jìn)制碼。 ? 律特性： )1ln( )1ln( ???? x (公式 21) 由于 ? 律不易用電子線路準(zhǔn)確實(shí)現(xiàn)，所以目前實(shí)用中采用特性近似的 15 折線代替 ?律，把縱坐標(biāo) y 從 0 到 1 之間劃分為 8 等份。均勻量化時(shí)，量化間隔是確定的，量化器的平均輸出信號量噪比隨量化電平數(shù) M 的增大而提高，故這種均勻量化器對于小輸入信號很不利。 （ 2）量化 設(shè)模擬信號的抽樣值為 m(KT)，其中 T 是抽樣周期， k 為整數(shù)。 圖 14 DPSK差分相干解調(diào)法原理框圖 由原理框圖可以看出方案一就是在 PSK 解調(diào)的基礎(chǔ)上增加了一個(gè)碼反變換器，而方案二是沒有用碼反變換器，在前面加了一個(gè) 延遲。 調(diào)制部分 方案選擇 PSK， DPSK 調(diào)制有兩種 方案 。 20xx 年 12 月 27 日至 20xx 年 12 月 30 日， 仿真 調(diào)試和設(shè)計(jì)說明書撰寫。武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計(jì)說明書 課程設(shè)計(jì)任務(wù)書 學(xué)生姓名： 專業(yè)班級： 指導(dǎo)教師： 工作單位： 題 目 : PSK 和 DPSK 信 號的仿真分析 初始條件： 在 MATLAB 仿真軟件系統(tǒng)中設(shè)計(jì)并實(shí)現(xiàn) PSK 和 DPSK 信號系統(tǒng)的仿真 ，并進(jìn)行調(diào)試和仿真 波形 分析。 20xx 年 1 月 8 日 ， 上交課程設(shè)計(jì)成果及報(bào)告，同時(shí)進(jìn)行答辯。 一是模擬調(diào)制法，其原理框圖如圖 11 所示。因?yàn)榍懊嬉呀?jīng)做了 PSK 的調(diào)制解調(diào)，為了方便，故選擇方案一。此抽樣值仍然是一個(gè)取值連續(xù)的變量。在非均勻量化時(shí)，量化間隔是隨信號抽樣值的不同而變化的。對應(yīng)于各轉(zhuǎn)折點(diǎn)的橫坐標(biāo) x 值可以按（公式22）來計(jì)算，結(jié)果列于表 21 中。權(quán)值電流 Iw 是在電路中預(yù)先產(chǎn)生的。段 內(nèi)碼代表的 32 個(gè)量化電平是均勻劃分的，所以，這 8 位共能表示 25628? 種量化值。 圖 23 PSK解調(diào) 原理框圖 帶通濾波器的意義是讓有用信號（已調(diào)信號）通過，濾除一部分噪聲，所以有用信號在 a 處得到信號為 ： ? ?( ) cos ca t s t t?? (公式 27) 假設(shè)相干載波的基準(zhǔn)相位與 2PSK 信號的調(diào)制載波的基準(zhǔn)相位一致（通常默認(rèn)為 0 相位）。而通信系統(tǒng)中常見的熱噪聲近似為高斯白噪聲，且符合加性。 2DPSK 信號的產(chǎn)生方法，先對二進(jìn)制數(shù)字基帶信號進(jìn)行差分編碼，即把表示數(shù)字 信息序列的絕對碼變換成相對碼，然后再根據(jù)相對 碼進(jìn)行絕對調(diào)相，從而 產(chǎn)生二進(jìn)制差分相(1 ) ( 0 / 1 ) ( 0) ( 0 / 1 )eP P P P P?? ? ?12erfc r?00(1 / 0 ) ( 0 ) ( )P P x f x d x?? ? ? ?? ?12erfc r?0 1( 0 / 1) ( 0 ) ( )P P x f x dx??? ? ? ?? ?12erfc r?12 rePer? ??武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計(jì)說明書 10 稱鍵控信號。因此，碼反變換器輸入端的誤碼率 Pe 可由 2PSK 信號采用相干解調(diào)時(shí)的誤碼率公式來確定。調(diào)制信號波形的相位變化可以見圖 34。第三個(gè)是同步解調(diào)相乘后的波形，與理論知識聯(lián)系比較后，其符合題目要求。 由于對 simulink 電路模塊的仿真不熟悉，所以我選擇用程序仿真。 ts = 0:1/fs:。 title(39。 %代碼存儲矩陣（全零） %段落碼判斷程序 for i=1:length(S) if (Q(i)=480)amp。 %在第三四七八段內(nèi)，段 落 碼第二位為 “1” end if (Q(i)=32)amp。(Q(i)=8160) code(i,4)=1。 %編碼段內(nèi)碼為二進(jìn)制 code(i,5)=str2double(k(1))。 end end 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計(jì)說明書 20 附錄 2 %psk 調(diào)制解調(diào) clear all。 t=0:40000*pi/199:40000*pi。%載波信號 mod1=[mod1 c]。)。 B=ones(1,200)。 grid on。 axis([0 15*length(g) 0 1000])。 title(39。 title(39。加入白噪聲的 2PSK 信號頻譜 39。相乘后信號波形 39。 %低通濾波器 fp=500。%計(jì)算階數(shù)和截止頻率 [b,a]=butter(n,wn)。) figure(2)。 if jt(m)0。 title(39。)。)。 High_Level = ones ( 1, Num_Unit )。 sign_orign = zeros ( 1, Num_Unit * Lenth_Of_Sign )。 subplot ( 3,2,1 )。 plot(abs(fft(sign_orign)))。 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計(jì)說明書 24 axis( [ 0 ,*Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 title(39。 axis( [ 0 , *Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 title(39。)。低通濾波后的波形 39。) grid %抽樣判決 amp。逆碼變換后的波形 39。) grid %誤碼率隨信噪比變化曲線 snr=0:2:20。 grid on。誤碼率隨信噪比的變化曲線 39。 for i=1:len_snr SNR_eb=exp(snr(i)*log(10)/10)。 figure(2)。 else sign_result( (I1)*Num_Unit + 1 : I*Num_Unit) = Low_Level。 figure(2)。 %映射為數(shù)字的 dt = filter(bz,az,dt)。)。相乘后的波形 39。)。調(diào)制后的信號 39。 title(39。 axis( [ 0 , *Time_Hold_On *( Lenth_Of_Sign + 1), (A / 2), A + (A / 2) ] )。 t = 0 : 10/fs : 10*Time_Hold_On * Lenth_Of_Sign 10/fs。 w =100000。 title(39。%信噪比范圍 len_snr = length(snr)。) figure(2)。 武漢理工大學(xué)《 Matl