【導(dǎo)讀】調(diào)試和仿真波形分析。調(diào)制信號(hào)：300Hz正弦信號(hào)，經(jīng)過(guò)μ律PCM編碼；載波頻率：100kHz；全文用A4紙打印，圖紙應(yīng)符合繪圖規(guī)范。進(jìn)也是通信系統(tǒng)性能提高的重要途徑。本文是基于PCM編碼的PSK和DPSK信號(hào)的。擬信號(hào)轉(zhuǎn)換為離散的數(shù)字信號(hào)，然后運(yùn)用MATLAB來(lái)仿真得到相應(yīng)波形。調(diào)制方法的性能，并通過(guò)比較仿真與理論計(jì)算的性能，證明了仿真的可靠性。

　　

【正文】 t(3,2,4)。 plot(jiet)。grid on axis([0 15*length(g) 2 2])。 title(39。相乘后信號(hào)波形 39。) figure(2)。 subplot(3,2,4)。 plot(abs(fft(jiet)))。 axis([0 15*length(g) 0 1000])。 title(39。相乘后信號(hào)頻譜 39。)。 %低通濾波器 fp=500。fs=700。rp=3。rs=20。fn=11025。 ws=fs/(fn/2)。 wp=fp/(fn/2)。%計(jì)算歸一化角頻率 [n,wn]=buttord(wp,ws,rp,rs)。%計(jì)算階數(shù)和截止頻率 [b,a]=butter(n,wn)。%計(jì)算 H(z） jt=filter(b,a,jiet)。 figure(1)。 subplot(3,2,5)。 plot(jt)。 grid on axis([0 15*length(g) 2 2])。 title(39。經(jīng)低通濾波器后信號(hào)波形 39。) figure(2)。 subplot(3,2,5)。 plot(abs(fft(jt)))。 axis([0 15*length(g) 0 1000])。 title(39。經(jīng)低通濾波器后信號(hào)頻譜 39。)。 %抽樣判決 for m=1:200*length(g)。 if jt(m)0。 jt(m)=1。 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計(jì)說(shuō)明書(shū) 22 else jt(m)=0。 jt(m)=0。 end end figure(1)。 subplot(3,2,6)。 plot(jt)。 grid on axis([0 15*length(g) 2 2])。 title(39。經(jīng)抽樣判決后信號(hào)波形 39。) figure(2)。 subplot(3,2,6)。 plot(abs(fft(jt)))。 axis([0 15*length(g) 0 1000])。 title(39。經(jīng)抽樣判決后信號(hào)頻譜 39。)。 %誤碼率隨信噪比的變化曲線(xiàn) snr = 0 : 2: 20。%信噪比范圍 len_snr = length(snr)。 for i = 1:len_snr SNR_eb = exp(snr(i)*log(10)/10)。 theo_err_prb(i) = (1/2)*erfc(sqrt(SNR_eb))。%PSK 系統(tǒng)的誤碼率與信噪比的關(guān)系 end figure(3)。 semilogy(snr,theo_err_prb,39。ko39。)。 grid on。 title(39。誤碼率隨信噪比的變化曲線(xiàn) 39。)。 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計(jì)說(shuō)明書(shū) 23 附錄 3 %dpsk 調(diào)制解調(diào) fs = 300。 Time_Hold_On =。 Num_Unit = fs * Time_Hold_On。 High_Level = ones ( 1, Num_Unit )。 Low_Level = zeros ( 1, Num_Unit )。 w =100000。 A = 1。 max=50 %定義 max 長(zhǎng)度 t=0:1/3000:。 Sign_Set =pcm(sin(2*pi*300*t))。 Lenth_Of_Sign = length ( Sign_Set )。 st = zeros ( 1, Num_Unit * Lenth_Of_Sign )。 sign_orign = zeros ( 1, Num_Unit * Lenth_Of_Sign )。 sign_result = zeros ( 1, Num_Unit * Lenth_Of_Sign )。 t = 0 : 10/fs : 10*Time_Hold_On * Lenth_Of_Sign 10/fs。 %產(chǎn)生基帶信號(hào) for I = 1 : Lenth_Of_Sign if Sign_Set(I) == 1 sign_orign( (I1)*Num_Unit + 1 : I*Num_Unit) = High_Level。 else sign_orign( (I1)*Num_Unit + 1 : I*Num_Unit) = Low_Level。 end end %調(diào)制部分 for I = 1 : Lenth_Of_Sign if Sign_Set(I) == 1 st( (I1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w * t( (I1)*Num_Unit + 1 : I*Num_Unit ) + pi )。 else st( (I1)*Num_Unit + 1 : I*Num_Unit) = A * cos ( 2 * pi * w * t( (I1)*Num_Unit + 1 : I*Num_Unit ) )。 end end figure(1)。 subplot ( 3,2,1 )。 plot(t, sign_orign)。 axis( [ 0 , *Time_Hold_On *( Lenth_Of_Sign + 1), (A / 2), A + (A / 2) ] )。 title ( 39。PCM 編碼信號(hào)序列 39。 )。 figure(2)。 subplot (3,2,1 )。 plot(abs(fft(sign_orign)))。 axis([0,200*Time_Hold_On*(Lenth_Of_Sign+1),0,400])。 title(39。PCM 編碼信號(hào)頻譜 39。)。 grid figure(1)。 subplot ( 3, 2, 2 )。 plot ( t, st )。 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計(jì)說(shuō)明書(shū) 24 axis( [ 0 ,*Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 title ( 39。調(diào)制后的信號(hào) 39。 )。 figure(2)。 subplot (3,2,2)。 plot(abs(fft(st)))。 axis([0,200*Time_Hold_On*(Lenth_Of_Sign+1),0,400])。 title(39。調(diào)制后信號(hào)的頻譜 39。)。 grid tz=awgn(st,10)。%信號(hào)調(diào)制中加入白噪聲，信噪比為 10 %相乘 dt = st .* cos ( 2 * pi * w * t )。 figure(1)。 subplot(3,2,3)。 plot ( t, dt )。 axis( [ 0 , *Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 title ( 39。相乘后的波形 39。 )。 grid figure(2)。 subplot(3,2,3)。 plot(abs(fft(dt)))。 axis([0,200*Time_Hold_On*(Lenth_Of_Sign+1),0,400])。 title(39。相乘后的波形頻譜 39。)。 %低通濾波部分 [N,Wn] = buttord( 2*pi*50, 2*pi*150,3,25,39。s39。)。 %臨界頻率采用角頻率表示 [b,a]=butter(N,Wn,39。s39。)。 [bz,az]=impinvar(b,a,fs)。 %映射為數(shù)字的 dt = filter(bz,az,dt)。 figure(1)。 subplot(3,2,4)。 plot ( t, dt )。 axis( [ 0 ,* Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 title ( 39。低通濾波后的波形 39。 )。 figure(2)。 subplot(3,2,4)。 plot(abs(fft(dt)))。 axis([0,200*Time_Hold_On*(Lenth_Of_Sign+1),0,400])。 title(39。低通濾波后信號(hào)的頻譜 39。) grid %抽樣判決 amp。 逆碼變換部分 for I = 1 : Lenth_Of_Sign if dt((2*I1)*Num_Unit/2) 0 sign_result( (I1)*Num_Unit + 1 : I*Num_Unit) = High_Level。 else sign_result( (I1)*Num_Unit + 1 : I*Num_Unit) = Low_Level。 end end。 figure(1) subplot(3,2,5)。 plot ( t, sign_result )。 武漢理工大學(xué)《 Matlab 應(yīng)用》課程設(shè)計(jì)說(shuō)明書(shū) 25 axis( [ 0 , *Time_Hold_On *( Lenth_Of_Sign + 1), 3*(A / 2), 3*(A / 2) ] )。 title ( 39。逆碼變換后的波形 39。 )。 figure(2)。 subplot(3,2,5)。 plot(abs(fft(sign_result)))。 axis([0,200*Time_Hold_On*(Lenth_Of_Sign+1),0,400])。 title(39。經(jīng)逆碼變換后信號(hào)的頻譜 39。) grid %誤碼率隨信噪比變化曲線(xiàn) snr=0:2:20。%信噪比范圍 len_snr=length(snr)。 for i=1:len_snr SNR_eb=exp(snr(i)*log(10)/10)。 theo_err_prb(i)=erfc(sqrt(SNR_eb))。%DPSK 系統(tǒng)的誤碼率與信噪比的關(guān)系 end figure(3)。 semilogy(snr,theo_err_prb,39。ko39。)。 grid on。 title(39。誤碼率隨信噪比的變化曲線(xiàn) 39。)。