【導(dǎo)讀】通信更是成為通信的主流，因此我們決定對(duì)數(shù)字信號(hào)的調(diào)制方式進(jìn)行簡(jiǎn)單的研究。制系統(tǒng)有數(shù)字基帶傳輸系統(tǒng)與數(shù)字帶通傳輸系統(tǒng)兩種。數(shù)字基帶傳輸系統(tǒng)多用于低頻傳。因?yàn)楹芏嘈诺啦痪?。?shù)字調(diào)制有兩種方法：一種是用模擬調(diào)制方法來(lái)實(shí)現(xiàn)數(shù)字式調(diào)。本設(shè)計(jì)主要對(duì)鍵控法進(jìn)行仿真研究。們的相關(guān)參數(shù)，例如誤碼率和頻帶寬度等特性，然后對(duì)它們的性能進(jìn)行比較。

　　

【正文】 為 a,b sg1=filter(b,a,signal)。 %信號(hào)通過(guò)該 BPF sg2=abs(sg1)。 %信號(hào)通過(guò)全波整流器 Fs=1e3。 %采樣頻率 [b,a]=butter(2,10*2/Fs) %設(shè)計(jì)巴特沃斯低通濾波器 sg3=filter(b,a,sg2)。 %信號(hào)通過(guò)該 LPF b=。 %判決門(mén)限 LL=tzd/2。 for i=1:n if sg3((i1)*tzd+LL)。 %取 sg2的中間的點(diǎn)作為判決點(diǎn) sg4(i)=1。 else sg4(i)=0。 end end %得到判決信號(hào) [numbers,pe] =symerr(g,sg4)。 %利用函數(shù)得到誤碼率和誤碼數(shù) ratio=ratio+numbers。 end r1=ratio/(n*10)。 %誤碼數(shù)除以總點(diǎn)數(shù)為誤碼率 pel(1,snrb+1)=r1。 %11個(gè)信噪比對(duì)應(yīng)的 11個(gè)誤碼率存入數(shù)組 pel end figure。 %畫(huà)圖 x=0:1:10。 x1=10.^((x+7)./10)。 %分貝值轉(zhuǎn)化為真值 y=*(erfc(sqrt(x1/4)))。 %2ASK信號(hào)非相關(guān)解調(diào)理論誤碼率計(jì)算 semilogy(x,pel,39。r39。,x,y,39。b39。)。 legend(39。simulation39。,39。theoritical case39。)。 xlabel(39。信噪比 39。) ylabel(39。誤碼率 39。)。grid on。 title(39。2ASK信號(hào)非相干解調(diào)時(shí)信噪比與誤碼率的關(guān)系 39。)。 2FSK 相干解調(diào) %構(gòu)造載波 得 已調(diào)信號(hào) % a=randsrc(1,8,[0:1])。%產(chǎn)生 8個(gè)隨機(jī)的數(shù) l=linspace(0,2*pi,50)。%利用 linspace函數(shù)創(chuàng)建數(shù)組 f1=sin(4*l)。 %生成載波 1 f2=sin(2*l)。%生成載波 2 t=linspace(0,10*pi,400)。 %定義時(shí)間軸 length為 10pi out=1:400。%規(guī)定已調(diào)信號(hào) length b=1:400。% 規(guī)定基帶信號(hào) length w1=1:400。% 規(guī)定載波 1length w2=1:400。% 規(guī)定載波 2length %生成 FSK信號(hào) % for i=1:8 if a(i)==0 38 for j=1:50 out(j+50*(i1))=f2(j)。 %碼元為 0,將載波 2輸出 end else for j=1:50 out(j+50*(i1))=f1(j)。 %碼元為 1,將載波 1輸出 end end end for i=1:8 for j=1:50 b(j+50*(i1))=a(i)。 %b作為基帶信號(hào)輸出 w1(j+50*(i1))=f1(j)。 %w1作為載波 1輸出 w2(j+50*(i1))=f2(j)。 %w2作為載波 2輸出 end end subplot(3,2,1),plot(t,b),axis([0 10*pi ]), xlabel(39。t39。),ylabel(39。幅度 39。),title(39?；鶐盘?hào) 39。)。 grid on。 subplot(3,2,2),plot(t,w1),axis([0 10*pi ]), xlabel(39。t39。),ylabel(39。幅度 39。),title(39。載波 139。)。 grid on。 subplot(3,2,3),plot(t,w2),axis([0 10*pi ]), xlabel(39。t39。),ylabel(39。幅度 39。),title(39。載波 239。)。 grid on。 subplot(3,2,4),plot(t,out),axis([0 10*pi ]), xlabel(39。t39。),ylabel(39。幅度 39。),title(39。FSK波形 39。)。 grid on。 %axis設(shè)定 XY軸起始值和終值 noise=awgn(out,60,39。measured39。) 。%產(chǎn)生噪音 ,信噪比 60 subplot(325)。plot(t,noise)。 ylabel(39。噪音 +信號(hào) 39。)。 xlabel(39。t39。)。 axis([0 10*pi 1 1])。grid on。 %信號(hào)通過(guò) BPF% Fs=400。%抽樣頻率 400HZ t=(1:400)*10*pi/Fs。%時(shí)軸步進(jìn) [b1,a1]=ellip(4,40,[25,35]*2/Fs)。%設(shè)計(jì) IIR帶通濾波器 1 [b2,a2]=ellip(4,40,[10,20]*2/Fs)。%設(shè)計(jì) IIR帶通濾波器 2 sf1=filter(b1,a1,noise)。%信號(hào)通過(guò)該濾波器 1 sf2=filter(b2,a2,noise)。%信號(hào)通過(guò)該濾波器 2 figure(2)。%兩路信號(hào)的解調(diào)過(guò)程 subplot(321)。 plot(t,sf1)。%通過(guò)該 BPF1的波形 xlabel(39。t39。)。 ylabel(39。幅度 39。)。title(39。f1通過(guò) BPF后的波形 39。)。 39 axis([0 10*pi ])。grid on。 subplot(322)。 plot(t,sf2)。%通過(guò) BPF2的波形 xlabel(39。t39。)。 ylabel(39。幅度 39。)。title(39。f2通過(guò) BPF后的波形 39。)。 axis([0 10*pi ])。grid on。 f1=[f1 f1 f1 f1 f1 f1 f1 f1]。 f2=[f2 f2 f2 f2 f2 f2 f2 f2]。 s1=sf1.*f1。 s2=sf2.*f2。 subplot(323)。 plot(t,s1)。%通過(guò)該帶通濾波器的波形 xlabel(39。t39。)。 ylabel(39。通過(guò)相乘器后波形 39。)。 axis([0 10*pi 1 1])。grid on。 subplot(324)。 plot(t,s2)。%通過(guò)該帶通濾波器的波形 xlabel(39。t39。)。 ylabel(39。通過(guò)相乘器后 波形 39。)。 axis([0 10*pi 1 1])。grid on。 %信號(hào)通過(guò) LPF% Fs=400。%抽樣頻率 400HZ t=(1:400)*10*pi/Fs。%時(shí)軸步進(jìn) [b1,a1]=ellip(4,40,[25]*2/Fs)。%設(shè)計(jì) IIRLPF [b2,a2]=ellip(4,40,[10]*2/Fs)。%設(shè)計(jì) IIRLPF sf1=filter(b1,a1,s1)。%信號(hào)通過(guò)該濾波器 sf2=filter(b2,a2,s2)。%信號(hào)通過(guò)該濾波器 figure(2)。 subplot(325)。 plot(t,sf1)。%通過(guò)該低通濾波器的波形 xlabel(39。t39。)。 ylabel(39。幅度 39。)。 title(39。sf1通過(guò) LPF后的波形 39。)。 axis([0 10*pi 1 1])。grid on。 figure(2)。 subplot(326)。 plot(t,sf2)。%通過(guò)該低通濾波器的波形 xlabel(39。t39。)。 ylabel(39。幅度 39。)。title(39。sf2通過(guò) LPF后的波形 39。)。 axis([0 10*pi 1 1])。grid on。 %抽樣判決 % fin=1:400。 for i=1:8 for j=1:50 if sf1(j+50*(i1)) sf2(j+50*(i1)) fin(j+50*(i1))=1。 %若 sf2sf1,碼元為 1 else fin(j+50*(i1))=0。 %若 sf2sf1,碼元為 0 end end 40 end figure(1)。 subplot(326)。 plot(t,fin)。%通過(guò)抽樣判決器的波形 xlabel(39。t39。)。 ylabel(39。幅度 39。)。title(39。抽樣判決后波形 39。)。 axis([2 10*pi ])。 誤碼率： rs=1e3。%時(shí)間軸頻率步進(jìn) fc=1e2。%載波頻率 100HZ tzd=1e2。 t=0:1/rs:(tzd1/rs)。 for snrb=0:1:10 %不同信噪比 ratio=0。%初始誤碼數(shù)設(shè)為 0 for k=1:10 %十次循環(huán) 得 10000碼元 n=1e3。%一次產(chǎn)生碼元數(shù) g=randint(1,n)。%產(chǎn)生 1000個(gè)碼元 tz1=g(ceil(10*t+(1/rs))).*cos(2*pi*fc*t)。%載波頻率為 100HZ tz2=~g(ceil(10*t+(1/rs))).*cos(2*pi*2*fc*t)。%載波頻率為 200HZ tz=tz1+tz2。%得到調(diào)制信號(hào) tz， 100個(gè)點(diǎn)表示 1個(gè)碼元 signal=awgn(tz,snrb)。 %解調(diào)過(guò)程 % Fs=1e3。%采樣頻率 [b1,a1]=butter(2,[80,120]*2/Fs)。%設(shè)計(jì)巴特沃斯帶通濾波器 [b2,a2]=butter(2,[180,220]*2/Fs)。%設(shè)計(jì)巴特 沃斯帶通濾波器 sg1=filter(b1,a1,signal)。%100HZ信號(hào)通過(guò)該 BPF sag1=filter(b2,a2,signal)。%200HZ信號(hào)通過(guò)該 BPF sg2=2*sg1.*cos(2*pi*fc*t)。%100HZ信號(hào)通過(guò)相乘器 sag2=2*sag1.*cos(2*pi*2*fc*t)。%200HZ信號(hào)通過(guò)相乘器 Fs=1e3。%采樣頻率 [b3,a3]=butter(2,10*2/Fs)%設(shè)計(jì)巴特沃斯低通濾波器 [b4,a4]=butter(2,10*2/Fs)。%設(shè)計(jì)巴特 沃斯低通濾波器 sg3= filter(b3,a3,sg2)。%100HZ信號(hào)通過(guò)該 LPF sag3= filter(b4, a4, sag2)。%200HZ信號(hào)通過(guò)該 LPF %判決過(guò)程 % LL=tzd/2。 for i=1:1000 if sg3((i1)*tzd+LL)=sag3((i1)*tzd+LL)。%比較兩路信號(hào)輸出值 u(i)=1。 else u(i)=0。 end end %計(jì)算誤碼率過(guò)程 % [numbers,pe] =symerr(g,u)。 %誤碼率和誤碼數(shù) 41 ratio=ratio+numbers。 end r1=ratio/(n*10)。 %誤碼率 pel(1,snrb+1)=r1。 % 11個(gè)誤碼率存入數(shù)組 pel end figure。 %畫(huà)圖 x=0:1:10。 x1=10.^((x+7)./10)。 %分貝值轉(zhuǎn)化為真值 y=*(erfc(sqrt(x1/2)))。 %2FSK誤碼率計(jì)算 semilogy(x,pel,39。r39。,x,y,39。b39。)。 legend(39。simulation39。,39。theoritical case39。)。 xlabel(39。信噪比 39。)。ylabel(39。誤碼率 39。)。grid on。 title(39。2FSK信號(hào)相干解調(diào)時(shí)信噪比與誤碼率的關(guān)系 39。)。