【導(dǎo)讀】數(shù)字調(diào)制是通信系統(tǒng)中最為重要的環(huán)節(jié)之一，幾種數(shù)字調(diào)制方法的仿真模型。最后，在仿真的基礎(chǔ)上分析比較了。理論計算的性能，證明了仿真模型的可行性。Keywords:Digital，modulation，analysis，simulation，MATLAB，根據(jù)控制的載波參量的不同，數(shù)字調(diào)制有。調(diào)幅、調(diào)相和調(diào)頻三種基本形式，并可以派生出多種其他形式。為了進(jìn)行長途傳輸，必須對數(shù)字信號進(jìn)行載波調(diào)制，將。信號頻譜搬移到高頻處才能在信道中傳輸。使用數(shù)字調(diào)制技術(shù)。另外，由于數(shù)字通信具有建網(wǎng)靈活，容易采用數(shù)字差。，所以通信系統(tǒng)都有由模擬方式向數(shù)字方式過渡的趨勢。統(tǒng)的重要部分之一，對它的研究也是有必要的。于改進(jìn)系統(tǒng)，獲得更佳的傳輸性能。Matlab的編程功能簡單,并且很容易擴(kuò)展和創(chuàng)造新的。應(yīng)用Matlab可方便地解決復(fù)雜數(shù)值計算問題?；ビ煤蛿?shù)據(jù)交換。Simulink支持連續(xù)、離散及兩者混合的線性和非線性。仿真軟件包用差分方程和微分方程建模相比,更直觀、方便和靈活。另外，Matlab的圖形界面功能GUI

　　

【正文】 id。 xlabel(39。相干解調(diào)后的波形 39。) n1=1024。 f1=(n1/2:1:n1/21)。 s=fft(y1,n1)。 subplot(4,2,2)。 plot(f1,abs(fftshift(s))),%axis([30 30 0 7]) xlabel(39。源信號頻譜 39。) n2=1024。 f2=(n2/2:1:n2/21)。 s=fft(y2,n2)。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 35 subplot(4,2,4)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。調(diào)制后的頻譜 39。) f2=(n2/2:1:n2/21)。 s=fft(y3,n2)。 subplot(4,2,6)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。經(jīng)過高斯信道后的頻譜 39。) f1=(n1/2:1:n1/21)。 s=fft(y4,n1)。 subplot(4,2,8)。 plot(f1,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。解調(diào)后的頻譜 39。) r=ErrorVec2(1)。 w=sprintf(39。符號誤差率 = %d\n39。,r)。 uiwait(msgbox(w,39。符號誤差率 39。,39。modal39。))。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 36 附 件 4 2MSK 仿真結(jié)果處理程序 附 件 [4] 2MSK仿真結(jié)果處理程序 function msk8() sim(39。msk939。)。 t1=ScopeData(:,1)。 y1=ScopeData(:,2)。 t3=ScopeData2(:,1)。 y3=ScopeData2(:,2)。 t4=ScopeData3(:,1)。 y4=ScopeData3(:,2)。 t5=ScopeData4(:,1)。 y5=ScopeData4(:,2)。 subplot(2,2,1)。 stairs(t1,y1),axis([0 20 1/2 3/2]),grid。 xlabel(39。源信號 39。) subplot(2,2,2)。 plot(t3,y3),grid。 xlabel(39。調(diào)制后的波形 39。) subplot(2,2,3)。 plot(t4,y4),grid。 xlabel(39。經(jīng)過高斯信道后的波形 39。) subplot(2,2,4)。 stairs(t5,y5),axis([0 20 1/2 3/2]),grid。 xlabel(39。相干解調(diào)后的波形 39。) figure n1=1024。 f1=(n1/2:1:n1/21)。 s=fft(y1,n1)。 subplot(2,2,1)。 plot(f1,abs(fftshift(s))),%axis([30 30 0 7]) xlabel(39。源信號頻譜 39。) n2=1024。 f2=(n2/2:1:n2/21)。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 37 s=fft(y3,n2)。 subplot(2,2,2)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。調(diào)制后的頻譜 39。) f2=(n2/2:1:n2/21)。 s=fft(y4,n2)。 subplot(2,2,3)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。經(jīng)過高斯信道后的頻譜 39。) f2=(n1/2:1:n1/21)。 s=fft(y5,n1)。 subplot(2,2,4)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。調(diào)制后的頻譜 39。) r=ErrorVec4(1)。 w=sprintf(39。符號誤差率 = %d\n39。,r)。 uiwait(msgbox(w,39。符號誤差率 39。,39。modal39。))。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 38 附 件 5 MDPSK 仿真結(jié)果處理程序 附 件 [5] MDPSK仿真結(jié)果處理程序 function mdpsk() sim(39。m_dpsk39。)。 t1=ScopeData(:,1)。 y1=ScopeData(:,2)。 t2=ScopeData1(:,1)。 y2=ScopeData1(:,2)。 t3=ScopeData2(:,1)。 y3=ScopeData2(:,2)。 t4=ScopeData3(:,1)。 y4=ScopeData3(:,2)。 subplot(4,2,1)。 stairs(t1,y1),axis([0 20 1/2 8]),grid。 xlabel(39。源信號 39。) subplot(4,2,3)。 plot(t2,y2),grid。 xlabel(39。調(diào)制后的波形 39。) subplot(4,2,5)。 plot(t3,y3),grid。 xlabel(39。經(jīng)過高斯信道后的波形 39。) subplot(4,2,7)。 stairs(t4,y4),axis([0 20 1/2 8]),grid。 xlabel(39。相干解調(diào)后 的波形 39。) n1=1024。 f1=(n1/2:1:n1/21)。 s=fft(y1,n1)。 subplot(4,2,2)。 plot(f1,abs(fftshift(s))),%axis([30 30 0 7]) xlabel(39。源信號頻譜 39。) n2=1024。 f2=(n2/2:1:n2/21)。 s=fft(y2,n2)。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 39 subplot(4,2,4)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。調(diào)制后的頻譜 39。) f2=(n2/2:1:n2/21)。 s=fft(y3,n2)。 subplot(4,2,6)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。經(jīng)過高斯信道后的頻譜 39。) f2=(n1/2:1:n1/21)。 s=fft(y4,n1)。 subplot(4,2,8)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。解調(diào)后的 頻譜 39。) r=ErrorVec4(1)。 w=sprintf(39。符號誤差率 = %d\n39。,r)。 uiwait(msgbox(w,39。符號誤差率 39。,39。modal39。))。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 40 附 件 6 各種 仿真 模型性能比較 程序 附 件 [6] 各仿真模型性能比較程序 x=1:11。 dpkk % 打開模型 for i=1:11 A=x(i)+26。 set_param(39。DPKK/AWGN Channel39。,39。EsNodB39。,num2str(A))。%給信道置數(shù) sim(39。DPKK39。)。 %運(yùn)行模型 y4(i)=ErrorVec(1)。 end save_system(39。DPKK39。)。 close_system(39。DPKK39。)。 %關(guān)閉模型 FKK for i=1:11 A=x(i)。 set_param(39。FKK/AWGN Channel39。,39。SNRdB39。,num2str(A))。 sim(39。FKK39。)。 y5(i)=ErrorVec(1)。 end save_system(39。FKK39。)。 close_system(39。FKK39。)。 mkk for i=1:11 A=x(i)+26。 set_param(39。MKK/AWGN Channel39。,39。EsNodB39。,num2str(A))。 sim(39。MKK39。)。 y6(i)=ErrorVec(1)。 end save_system(39。MKK39。)。 close_system(39。MKK39。)。 semilogy(x,y5,39。g39。,x,y4,39。b39。,x,y6,39。r39。) %畫圖 title(39。各調(diào)制模型性能比較 39。) xlabel(39。SNR(dB)39。) ylabel(39。ErrorVec39。) legend(39。2FSK39。,39。2DPSK39。,39。MSK39。)。 %圖示 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 41 附 件 7 2FSK 模型性能與理論值比較程序 附 件 [7] 2FSK模型性能與理論值比較程序 r=0::10。 for i=1:length(r), SNR=exp(r(i)*log(10)/10)。 % signaltonoise ratio theo_err_prb(i)=Qfunct(sqrt(SNR))。 %計算理論值 end x=1:8。 FKK for i=1:8 A=x(i)。 set_param(39。FKK/AWGN Channel39。,39。SNRdB39。,num2str(A))。 sim(39。FKK39。)。 y3(i)=ErrorVec(1)。 end save_system(39。FKK39。)。 close_system(39。FKK39。)。 semilogy(r,theo_err_prb,x,y3,39。*39。)。 title(39。FSK性能理論值與仿真值對比 39。) xlabel(39。SNR(dB)39。) ylabel(39。ErrorVec39。) legend(39。2FSK理論值 39。,39。2FSK仿真值 39。)。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 42 附 件 8 判決模塊 S 函數(shù)的 m 文件 附 件 [8] 判決模塊 S函數(shù)的 m文件： function [sys,x0,str,ts] = sfuntmpl(t,x,u,flag) switch flag, case 0, [sys,x0,str,ts]=mdlInitializeSizes。 case 3, sys=mdlOutputs(t,x,u)。 case{1,2,4, 9}, sys=[]。 otherwise error([39。Unhandled flag = 39。,num2str(flag)])。 end function [sys,x0,str,ts]=mdlInitializeSizes sizes = simsizes。 = 0。 = 1。 =1。 = 1。 = 1。 = 1。 sys = simsizes(sizes)。 x0 = [0]。 str = []。 ts = [0 0]。 function sys=mdlOutputs(t,x,u) if u=。 sys=1。 elseif u=。 sys=0。 end