【正文】 end cp=[cp A]。%每個值 200 個點 else g(n)==1。 if g(n)==0。 t=0:40000*pi/199:40000*pi。 mod1=[]。 g=randint(1,max)。 clf。 設(shè)接收端帶通濾波器輸出波形 y((t)為： （公式 310） 經(jīng)過相干解調(diào)后，送入抽樣判決器的輸入波形為： （公式 311） 由于 ()t是均值為 0，方差為 2n? 的高斯噪聲，所以 x(t)的一維概率密度函數(shù)為 （公式 312） （公式 313） 由最佳判決門限分析可知，在發(fā)送“ 1”和發(fā)送“ 0”的概率相等時，最佳判決門限 b*=，發(fā)“ 1”而錯判為“ 0”的概率為 （公式 314） 0 1( 0 / 1) ( 0 ) ( )P P x f x d x??? ? ? ?? ?12erfc r?發(fā) 送 端 信 道 帶 通 濾 波 器 相 乘 器 低 通 濾 波 器 抽 樣 判 決 器 定時 脈沖 )(tsT )(tni )(tyi ct?cos2 y(t) x(t) 輸出 [ ( ) ] c o s ( ) s in , 1() [ ( ) ] c o s ( ) s in 0c c s cc c s ca n t t n t tyt a n t t n t t?????? ? ? ? ??發(fā) 送 “ ” 時， 發(fā) 送 “ ” 時( ) , 1() ( ) , 0cca n txt a n t??? ? ???發(fā) 送 “ ” 符 號發(fā) 送 “ ” 符 號時發(fā)送“ 12 )(e xp2 1)( 2 21 ?????? ???nnaxxf ???”時發(fā)送“ 02 )(e xp2 1)( 2 20 ?????? ???nnaxxf ???武漢理工大學 MATLAB 課程設(shè)計 15 同理，發(fā)送“ 0”而錯判為“ 1”的概率為 （公式 315） 故 2PSK 信號相干解調(diào)時系統(tǒng)的總誤碼率為 （公式 316） 在大信噪比條件下，上式可近似為 （公式 317） 用 matlab 實現(xiàn) PSK 系統(tǒng)的調(diào)制解調(diào) clear all。根據(jù)設(shè)計要求考慮不同信噪比的高斯 白噪聲對 PSK 系統(tǒng)的影響。 帶通 濾波器 相乘器 低通 濾波器 抽樣 判決器 a c d e b cosWctwcos epsk(t) 定時 脈沖 輸出 武漢理工大學 MATLAB 課程設(shè)計 14 圖 35PSK 信號相干解調(diào)系統(tǒng)性能分析模型 當信號經(jīng)過信道傳輸時會受到噪聲的影響，這是不可避免的。另外，在隨機信號碼元序列中，信號波形有可能出現(xiàn)長時間連續(xù)的正弦波形，致使在接收端無法辨認信號碼元的起止時刻。這種現(xiàn)象稱為 2PSK 方式的 “ 倒π ” 現(xiàn)象或 “ 反相工 作 ” 。所以得到下式 2 2 22 11( ) ( ) c o s ( ) c o s ( ) ( ) c o s 222P K S c c cc t e t t s t t s t s t t? ? ?? ? ? ? (公式 38) 通過低通濾波器后 21( ) ( )2d t s t? （ 公式 39） 最后通過抽樣判決器恢復出數(shù)字信號。 （ 2） PSK 解調(diào)過程分析 解調(diào)原理框圖如圖 34 所示。因此， 2PSK 信號的時域表達式為 )c o s ()(2 ?? ncP S K tAte ?? （公式 33） 其中， ?n表示第 n 個符號的絕對相位： ?n=0（發(fā)送“ 0”時，概率為 P），?n = ? （發(fā)送“ 1 ”時 , 概率為 1 — P ） 因此，公 式 33 可以改寫為??????? tAtAtccP S Ke ??c o sc o s)(2 （公式 34） 由于表示信號的兩種碼元的波形相同極性相反，故 2PSK 信號一般可以表述為一個雙極性全占空比矩形脈沖序列與一個正弦載波的相乘 ，即ttS ctP S Ke ?c o s)()(2 ? （公式 35） 其中， )()( Ta sn n ntgtS ?? ? （公式 36） 這里， g(t)是脈寬為 Ts 的單個矩形脈沖，而 an 的統(tǒng)計特性為 an=1(概率為P),an=0(概率為 1P).即 發(fā)送二進制符號 ” 0” 時 , )(2 tePSK取 0 相位 ,發(fā)送二進制符號 ” 1” 時 , )(2 tePSK取 ? 相位 .這種以載波的不同相位直接去表示相應(yīng)二進制數(shù)字信號的調(diào)制方式 ,稱為絕對相移方式。 end %符 號位的判斷 end 武漢理工大學 MATLAB 課程設(shè)計 12 調(diào)制解調(diào)原理 （ 1） PSK 調(diào)制過程分析 相移鍵控（ PSK）是利用載波的相位變化來傳遞數(shù)字信息，而振幅和頻率保持不變。 end if z(i)0 code(i,1)=1。 code(i,7)=str2num(k(3))。 %編碼段內(nèi) 碼為二進制 code(i,5)=str2num(k(1))。 %求出在段內(nèi)的位置 if q==0 code(i,(5:8))=[0,0,0,0]。 %量化間隔 b=[1,1,2,4,8,16,32,64]。 for i=1:length(S) N(i)=bin2dec(num2str(code(i,2:4)))+1。(Q(i)=2048) code(i,4)=1。(Q(i)128)||(Q(i)=256)amp。 %在第三四七八段內(nèi)，段位碼第二位為“ 1” end if (Q(i)=16)amp。(Q(i)128)||(Q(i)=512)amp。(Q(i)=2048) code(i,2)=1。 %量化 code=zeros(length(S),8)。 %求 S 的最大值 S=abs(S/MaxS)。 %量化 z=sign(S)。3000HZ 抽樣波形 39。 subplot(2,1,2) %抽樣 stem(ts,Sign_m)。 sin(2*pi*300*t)39。 figure subplot(2,1,1) plot(t,Signal)。 武漢理工大學 MATLAB 課程設(shè)計 10 Signal = sin(2*pi*300*t)。 fs = 3000。在表 33 和表 34 中給出了段落碼和段內(nèi)碼的編碼規(guī)則。其中第 2到 4 位（ c2c3c4）是段落碼，共計 3 位，可以表示 8 種斜率的段落；其他 5 位（ c5~c9）為段內(nèi)碼，可以表示每一段落內(nèi)的 32 種量化電平。其中第一位 c1 表示量化值的極性正負。它們能夠表示 8 個十進制數(shù)，從 0 到 7，如表 32 所列。 Iw 的個數(shù)決定于編碼的位數(shù)，現(xiàn)在共有三個不同的 Iw 值。每比較一次，得出 1 位二進制碼。編碼器的輸入信號抽樣脈沖值在 0和 之間。在用電路實現(xiàn)時，最常用的是一種方案稱為逐次比較法編碼，其基本原理方框圖如量化器 m(kT) mq(kT) 武漢理工大學 MATLAB 課程設(shè)計 8 圖 所示。對應(yīng)于各轉(zhuǎn)折點的橫坐標 x 值可 以按（公式 32）來計算，結(jié)果列于表 中。由于本課程設(shè)計要求用 ? 律來編碼，故對 A 律不作介紹。通常使用的壓縮器中，大多采用對數(shù)式壓縮。信號抽樣值小時，量化間隔也小，信號抽樣值大時，量化間隔也變大。均勻量化時，量化間隔是確定的，量化器的平均輸出信號量噪比隨量化電平數(shù) M 的增大而提高，故這種均勻量化器對于小輸入信號很不利。在實際中，量化過程常是和后續(xù)的編碼過程結(jié)合在一起完成的，不一定存在量化器。用這 M 個量化電平表示連續(xù)抽樣值的方法稱為量化。若僅用 N 個二進制數(shù)字碼元來代表此抽樣值的大小，則 N 個二進制碼元只能代表 M= 2N 個不同的抽樣值。 （ 2）量化 設(shè)模擬信號的抽樣值為 m(KT)，其中 T 是抽 樣周期， k 為整數(shù)。抽樣結(jié)果得到的是一系列周期性的沖激脈沖，其面積和模擬信號的取值成正比。 圖 31 PCM編碼原理框圖 （ 1）抽樣 模擬信號通常是在時 間上連續(xù)的信號，一在系列離散的點上，對這種信號抽取樣值稱為抽樣。 帶通濾波器 相乘器 低通濾波器 抽樣判決器 碼反變換器 tc?cos 定時脈 沖 輸出 a c d e f b )(2 teDPSK 帶通濾波器 相乘器 低通濾波器 抽樣判決器 輸出 定時脈 沖 )(2 teDPSK a b c d e 延遲 Ts 武漢理工大學