【正文】 t(t,b),axis([0 10*pi ]), xlabel(39。),title(39。 subplot(3,3,2),plot(t,w),axis([0 10*pi ]), xlabel(39。),title(39。 subplot(3,3,3),plot(t,out),axis([0 10*pi ]), xlabel(39。),title(39。 %已調(diào)信號加入高斯白噪聲 % noise=awgn(out,100,39。 plot(t,noise)。 xlabel(39。噪音 +信號 39。 %信號通過 BPF% 36 Fs=400。%信號通過 該濾波器 subplot(335)。)。title(39。grid on。%設(shè)計 IIRLPF sf=filter(b,a,sf)。t39。)。 axis([0 10*pi 3 3])。 %sf判決門限，碼元為 1 else sf(j+50*(i1))=0。t39。)。 axis([3 10*pi ])。 %載波頻率 100HZ tzd=1e2。 %1次產(chǎn)生碼元數(shù) g=randint(1,n)。 %采樣頻率 37 [b,a]=butter(2,[80,120]*2/Fs)。 %采樣頻率 [b,a]=butter(2,10*2/Fs) %設(shè)計巴特沃斯低通濾波器 sg3=filter(b,a,sg2)。 %取 sg2的中間的點作為判決點 sg4(i)=1。 end r1=ratio/(n*10)。 x1=10.^((x+7)./10)。,x,y,39。simulation39。 xlabel(39。)。)。 %生成載波 1 f2=sin(2*l)。% 規(guī)定基帶信號 length w1=1:400。 %碼元為 1,將載波 1輸出 end end end for i=1:8 for j=1:50 b(j+50*(i1))=a(i)。t39。基帶信號 39。t39。載波 139。t39。載波 239。t39。FSK波形 39。measured39。 ylabel(39。t39。 %信號通過 BPF% Fs=400。%設(shè)計 IIR帶通濾波器 2 sf1=filter(b1,a1,noise)。 plot(t,sf1)。 ylabel(39。f1通過 BPF后的波形 39。 subplot(322)。)。title(39。grid on。 s2=sf2.*f2。t39。)。 plot(t,s2)。 ylabel(39。grid on。%設(shè)計 IIRLPF [b2,a2]=ellip(4,40,[10]*2/Fs)。 subplot(325)。)。 title(39。grid on。%通過該低通濾波器的波形 xlabel(39。幅度 39。)。 for i=1:8 for j=1:50 if sf1(j+50*(i1)) sf2(j+50*(i1)) fin(j+50*(i1))=1。 plot(t,fin)。 ylabel(39。抽樣判決后波形 39。%時間軸頻率步進 fc=1e2。%初始誤碼數(shù)設(shè)為 0 for k=1:10 %十次循環(huán) 得 10000碼元 n=1e3。%載波頻率為 200HZ tz=tz1+tz2。%設(shè)計巴特沃斯帶通濾波器 [b2,a2]=butter(2,[180,220]*2/Fs)。%100HZ信號通過相乘器 sag2=2*sag1.*cos(2*pi*2*fc*t)。%100HZ信號通過該 LPF sag3= filter(b4, a4, sag2)。 else u(i)=0。 %誤碼率 pel(1,snrb+1)=r1。 %分貝值轉(zhuǎn)化為真值 y=*(erfc(sqrt(x1/2)))。b39。,39。信噪比 39。)。)。 title(39。ylabel(39。)。 legend(39。r39。 %畫圖 x=0:1:10。 %誤碼率和誤碼數(shù) 41 ratio=ratio+numbers。 for i=1:1000 if sg3((i1)*tzd+LL)=sag3((i1)*tzd+LL)。%采樣頻率 [b3,a3]=butter(2,10*2/Fs)%設(shè)計巴特沃斯低通濾波器 [b4,a4]=butter(2,10*2/Fs)。%100HZ信號通過該 BPF sag1=filter(b2,a2,signal)。 %解調(diào)過程 % Fs=1e3。%產(chǎn)生 1000個碼元 tz1=g(ceil(10*t+(1/rs))).*cos(2*pi*fc*t)。 t=0:1/rs:(tzd1/rs)。 axis([2 10*pi ])。)。t39。 %若 sf2sf1,碼元為 0 end end 40 end figure(1)。grid on。title(39。)。 subplot(326)。)。幅度 39。%通過該低通濾波器的波形 xlabel(39。%信號通過該濾波器 sf2=filter(b2,a2,s2)。%抽樣頻率 400HZ t=(1:400)*10*pi/Fs。)。t39。grid on。 ylabel(39。 plot(t,s1)。 f2=[f2 f2 f2 f2 f2 f2 f2 f2]。)。幅度 39。%通過 BPF2的波形 xlabel(39。 39 axis([0 10*pi ])。)。t39。%信號通過該濾波器 2 figure(2)。%時軸步進 [b1,a1]=ellip(4,40,[25,35]*2/Fs)。 axis([0 10*pi 1 1])。)。%產(chǎn)生噪音 ,信噪比 60 subplot(325)。 grid on。幅度 39。 grid on。幅度 39。 grid on。幅度 39。 grid on。幅度 39。 %w1作為載波 1輸出 w2(j+50*(i1))=f2(j)。% 規(guī)定載波 2length %生成 FSK信號 % for i=1:8 if a(i)==0 38 for j=1:50 out(j+50*(i1))=f2(j)。 %定義時間軸 length為 10pi out=1:400。%產(chǎn)生 8個隨機的數(shù) l=linspace(0,2*pi,50)。 title(39。) ylabel(39。theoritical case39。)。 %2ASK信號非相關(guān)解調(diào)理論誤碼率計算 semilogy(x,pel,39。 %11個信噪比對應(yīng)的 11個誤碼率存入數(shù)組 pel end figure。 end end %得到判決信號 [numbers,pe] =symerr(g,sg4)。 %判決門限 LL=tzd/2。 %信號通過該 BPF sg2=abs(sg1)。 %得到已調(diào)信號 tz， 100個點表示 1個碼元 signal=awgn(tz,snrb)。 %碼元周期 ,即 碼元速率 10 for snrb=0:1:10 %不同信噪比 ratio=0。 誤碼率： rs=1e3。抽樣判決后波形 39。 ylabel(39。 plot(t,sf)。 %抽樣判決 % b=。通過 LPF 后的波形 39。 ylabel(39。 plot(t,sf)。%抽樣頻率 400HZ t=(1:400)*10*pi/Fs。)。幅度 39。%通過該 BPF 的波形 xlabel(39。%時軸步進 [b,a]=ellip(4,40,[10,25]*2/Fs)。 axis([0 10*pi ])。)。幅度 39。) 。)。),ylabel(39。)。),ylabel(39。)。),ylabel(39。 %碼元為 1 則載波輸出 end end end for i=1:8 for j=1:50 b(j+50*(i1))=a(i)。%規(guī)定已調(diào)信號 length b=1:400。%產(chǎn)生 8 個隨機的數(shù) l=linspace(0,2*pi,50)。 title(39。) ylabel(39。theoritical case39。)。%2ASK 誤碼率計算 semilogy(x,pel,39。% 11 個誤碼率存入數(shù)組 pel end figure。 end end%得到判決后信號 sg4 [numbers,pe] =symerr(g,sg4)。%判決門限 LL=tzd/2。%信號通過該 BPF sg2=2*sg1.*cos(2*pi*fc*t)。%得到調(diào)制信號 tz signal=awgn(tz,snrb)。 for snrb=0:1:10 %不同信噪比 ratio=0。 誤碼率： rs=1e3。抽樣判決后波形 39。 ylabel(39。 plot(t,sf)。 %抽樣判決 % b=。通過 LPF 后的波形 39。 ylabel(39。 plot(t,sf)。%抽樣頻率 400HZ t=(1:400)*10*pi/Fs。)。幅度 39。%通過該相乘器的波形 xlabel(39。 %信號經(jīng)過相乘器 % f=[f f f f f f f f]。通過 BPF 后的波形 39。 ylabel(39。 plot(t,sf)。%抽樣頻率 400HZ t=(1:400)*10*pi/Fs。)。幅度 39。 33 xlabel(39。measured39。 ASK 波形 39。t39。載波 39。t39。調(diào)制信號 39。t39。 %碼元 0 則將 0 輸出 end else for j=1:50 out(j+50*(i1))=5*f(j)。%定義時軸 out=1:400。 [7] 樊昌信 . 通信原理教程 [M]. 北京 : 電子工業(yè)出版社 , 2021. [8] 劉淑堂（譯） . 現(xiàn)代通信系統(tǒng)使用 MATLAB大全 [M]. 西安 : 西安交通大學出版社 , 2021. [9] 張志涌 . 精通 [M]. 北京 : 北京航空航天大學出版社 , 2021. [10] Math Works. MATLAB Compiler, 2021. [11] 王興亮 . 數(shù)字通信原理與技術(shù) [M]. 西安 : 西安電子科技大學出版社 , 2021. [12] 楊世海 , 周良柱 . 基于 Matlab的 CDMA無線通信可視化仿真 [J]. 系統(tǒng)仿真學報 , 2021, 13(1): 98100. [13] 李建新 , 劉乃安 . 現(xiàn)代通信系統(tǒng)分析與仿真 MATLAB 通信工具箱 [M]. 西安電子科技大學出版社 , 2021. [14] 郭海燕 , 畢紅軍 . Matlab 在偽隨機碼的生成及仿真中的應(yīng)用 [J]. 計算機仿真， 2021, 21(3): 149152. [15] 徐炳祥 . 數(shù)字動態(tài)濾波器 [J]. 電子科學技 術(shù) , 1981, 4(2): 102158. [16] Peter Communications(Second Edition).China Machine Press [17] Adam Lender. Correlative Digital Communication Techniques. IEEE Transactions on Communication Technology, 1996 [18] , [19] Emmanuel Barrie Sigital House Of Electronics Industry [20] Richard Digital Siginal Press [21] of Quantized ,1984,27(3) 32 附 錄 2ASK 相干解調(diào) %構(gòu)造載波得已調(diào)信號 % a