【正文】 sys=0。 sys=1。 ts = [0 0]。 x0 = [0]。 = 1。 = 1。 = 1。 end function [sys,x0,str,ts]=mdlInitializeSizes sizes = simsizes。Unhandled flag = 39。 case{1,2,4, 9}, sys=[]。 基于 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。2FSK仿真值 39。2FSK理論值 39。ErrorVec39。SNR(dB)39。FSK性能理論值與仿真值對比 39。)。 semilogy(r,theo_err_prb,x,y3,39。FKK39。)。 end save_system(39。)。 sim(39。SNRdB39。FKK/AWGN Channel39。 FKK for i=1:8 A=x(i)。 % signaltonoise ratio theo_err_prb(i)=Qfunct(sqrt(SNR))。 %圖示 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 41 附 件 7 2FSK 模型性能與理論值比較程序 附 件 [7] 2FSK模型性能與理論值比較程序 r=0::10。MSK39。2DPSK39。2FSK39。ErrorVec39。SNR(dB)39。各調(diào)制模型性能比較 39。r39。b39。g39。)。 close_system(39。MKK39。 y6(i)=ErrorVec(1)。MKK39。,num2str(A))。,39。 set_param(39。)。 close_system(39。FKK39。 y5(i)=ErrorVec(1)。FKK39。,num2str(A))。,39。 set_param(39。)。 close_system(39。DPKK39。 %運行模型 y4(i)=ErrorVec(1)。DPKK39。,num2str(A))。,39。 set_param(39。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 40 附 件 6 各種 仿真 模型性能比較 程序 附 件 [6] 各仿真模型性能比較程序 x=1:11。modal39。符號誤差率 39。,r)。 w=sprintf(39。解調(diào)后的 頻譜 39。 subplot(4,2,8)。) f2=(n1/2:1:n1/21)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。 s=fft(y3,n2)。調(diào)制后的頻譜 39。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 39 subplot(4,2,4)。 f2=(n2/2:1:n2/21)。源信號頻譜 39。 subplot(4,2,2)。 f1=(n1/2:1:n1/21)。相干解調(diào)后 的波形 39。 stairs(t4,y4),axis([0 20 1/2 8]),grid。經(jīng)過高斯信道后的波形 39。 plot(t3,y3),grid。調(diào)制后的波形 39。 plot(t2,y2),grid。源信號 39。 stairs(t1,y1),axis([0 20 1/2 8]),grid。 y4=ScopeData3(:,2)。 y3=ScopeData2(:,2)。 y2=ScopeData1(:,2)。 y1=ScopeData(:,2)。)。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 38 附 件 5 MDPSK 仿真結(jié)果處理程序 附 件 [5] MDPSK仿真結(jié)果處理程序 function mdpsk() sim(39。modal39。符號誤差率 39。,r)。 w=sprintf(39。調(diào)制后的頻譜 39。 subplot(2,2,4)。) f2=(n1/2:1:n1/21)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。 s=fft(y4,n2)。調(diào)制后的頻譜 39。 subplot(2,2,2)。 f2=(n2/2:1:n2/21)。源信號頻譜 39。 subplot(2,2,1)。 f1=(n1/2:1:n1/21)。相干解調(diào)后的波形 39。 stairs(t5,y5),axis([0 20 1/2 3/2]),grid。經(jīng)過高斯信道后的波形 39。 plot(t4,y4),grid。調(diào)制后的波形 39。 plot(t3,y3),grid。源信號 39。 stairs(t1,y1),axis([0 20 1/2 3/2]),grid。 y5=ScopeData4(:,2)。 y4=ScopeData3(:,2)。 y3=ScopeData2(:,2)。 y1=ScopeData(:,2)。)。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 36 附 件 4 2MSK 仿真結(jié)果處理程序 附 件 [4] 2MSK仿真結(jié)果處理程序 function msk8() sim(39。modal39。符號誤差率 39。,r)。 w=sprintf(39。解調(diào)后的頻譜 39。 subplot(4,2,8)。) f1=(n1/2:1:n1/21)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。 s=fft(y3,n2)。調(diào)制后的頻譜 39。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 35 subplot(4,2,4)。 f2=(n2/2:1:n2/21)。源信號頻譜 39。 subplot(4,2,2)。 f1=(n1/2:1:n1/21)。相干解調(diào)后的波形 39。 stairs(t4,y4),axis([0 20 1/2 3/2]),grid。經(jīng)過高斯信道后的波形 39。 plot(t3,y3),grid。調(diào)制后的波形 39。 plot(t2,y2),grid。源信號 39。 stairs(t1,y1),axis([0 20 1/2 3/2]),grid。 y4=ScopeData3(:,2)。 y3=ScopeData2(:,2)。 y2=ScopeData1(:,2)。 y1=ScopeData(:,2)。)。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 34 附 件 3 2DPSK 仿真結(jié)果處理程序 附 件 [3] 2DPSK仿真結(jié)果處理程序 function psk() sim(39。modal39。符號誤差率 39。,r)。 w=sprintf(39。調(diào)制后的頻譜 39。 subplot(4,2,8)。) f2=(n1/2:1:n1/21)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。 s=fft(y3,n2)。調(diào)制后的頻譜 39。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 33 subplot(4,2,4)。 f2=(n2/2:1:n2/21)。源信號頻譜 39。 subplot(4,2,2)。 f1=(n1/2:1:n1/21)。相干解調(diào)后的波形 39。 stairs(t4,y4),axis([0 20 1/2 3/2]),grid。經(jīng)過高斯信道后的波形 39。 plot(t3,y3),grid。調(diào)制后的波形 39。 plot(t2,y2),grid。源信號 39。 stairs(t1,y1),axis([0 20 1/2 3/2]),grid。 y4=ScopeData3(:,2)。 y3=ScopeData2(:,2)。 y2=ScopeData1(:,2)。 y1=ScopeData(:,2)。)。) 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 32 附 件 2 2FSK 仿真結(jié)果處理程序 附 件 [2] 2FSK仿真結(jié)果處理程序 function bfsk2() sim(39。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。) s=fft(y7,n2)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。) s=fft(y1,n2)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。) s=fft(y2,n2)。 plot(f2,abs(fftshift(s))),%axis([400 400 0 300]) xlabel(39。 s=fft(y5,n2)。) n2=1024。 plot(f1,abs(fftshift(s))),%axis([30 30 0 7]) xlabel(39。 s=fft(y3,n1)。) figure(2) n1=1024。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 31 xlabel(39。) subplot(3,2,6)。 xlabel(39。) subplot(3,2,3)。 xlabel(39。) subplot(3,2,2)。) xlabel(39。 title(39。) subplot(3,2,1)。 xlabel(39。) subplot(3,2,4)。 xlabel(39。 %引入示波器的數(shù)據(jù) subplot(3,2,5)。 t7=ScopeData3(:,1)。 t6=ScopeData7(:,1)。 t5=ScopeData1(:,1)。 t3=ScopeData(:,1)。 t2=ScopeData2(:,1)。 t1=ScopeData5(:,1)。dm39。 基于 MATLAB 的數(shù)字調(diào)制系統(tǒng)仿真 29 參考文獻 [1] 樊昌信等 . 通信原理（第 5版） . 北京：國防工業(yè)出版社， 2020 [2] 曹志剛，錢亞生 . 現(xiàn)代通信原理 . 北京：清華大學 出版社， 2020 [3] 李建新等 . 現(xiàn)代通信系統(tǒng)分析與仿真 —— MATLAB通信工具箱 . 西安：西安電子科技大學出版社， 2020 [4] 王沫然 . Simulink4建模及動態(tài)仿真 . 北京：電子工業(yè)出版社， 2020 [5] ,，張航等譯 . 精通 MATLAB6. 北京：清華大學出版社， 2020 [6] 約翰 .G..普羅克斯，馬蘇德 .薩勒赫著 .劉樹棠譯 . 現(xiàn)在通信系統(tǒng) —— 使用MATLAB. 西安：西安交通大學出版社， 2020 [7] 戴虹，戴悟僧 . MATLAB在通信原理仿真中的應用 . 上海第二工業(yè)大學學報 . 2020年，第 1期 [8] 婁 莉 . GMSK 數(shù)字調(diào)制的仿真與分析 . 現(xiàn)代電子技術(shù) . 2020 年第 18 期總第185 期： P66— P68 [9] 李蕾，杜巖 . 基于 MATLAB 的數(shù)字基帶傳輸系統(tǒng)仿真實驗設計 . 山東工業(yè)大