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