matlab二進(jìn)制數(shù)字調(diào)制與解調(diào)系統(tǒng)課程設(shè)計(jì)報(bào)告

1、<p><b>　　一．設(shè)計(jì)題目：</b></p><p>　　二進(jìn)制數(shù)字調(diào)制與解調(diào)系統(tǒng)的設(shè)計(jì)</p><p><b>　　二．主要內(nèi)容：</b></p><p>　　二進(jìn)制數(shù)字調(diào)制與解調(diào)系統(tǒng)的設(shè)計(jì)</p><p>　　MATLAB及SIMULINK建模環(huán)境簡介 </p>&

2、lt;p>　　MATLAB 是美國MathWorks公司出品的商業(yè)數(shù)學(xué)軟件，用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計(jì)算的高級技術(shù)計(jì)算語言和交互式環(huán)境，主要包括MATLAB和SIMULINK兩大部分。</p><p>　　Simulink是MATLAB最重要的組件之一，它提供一個(gè)動(dòng)態(tài)系統(tǒng)建模、仿真和綜合分析的集成環(huán)境。在該環(huán)境中，無需大量書寫程序，而只需要通過簡單直觀的鼠標(biāo)操作，就可構(gòu)造出復(fù)雜的系統(tǒng)。Si

3、mulink具有適應(yīng)面廣、結(jié)構(gòu)和流程清晰及仿真精細(xì)、貼近實(shí)際、效率高、靈活等優(yōu)點(diǎn)，并基于以上優(yōu)點(diǎn)Simulink已被廣泛應(yīng)用于控制理論和數(shù)字信號處理的復(fù)雜仿真和設(shè)計(jì)。同時(shí)有大量的第三方軟件和硬件可應(yīng)用于或被要求應(yīng)用于Simulink。 </p><p>　　Simulink 是MATLAB中的一種可視化仿真工具， 是一種基于MATLAB的框圖設(shè)計(jì)環(huán)境，是實(shí)現(xiàn)動(dòng)態(tài)系統(tǒng)建模、仿真和分析的一個(gè)軟件包，被廣泛應(yīng)用于線性系

4、統(tǒng)、非線性系統(tǒng)、數(shù)字控制及數(shù)字信號處理的建模和仿真中。Simulink可以用連續(xù)采樣時(shí)間、離散采樣時(shí)間或兩種混合的采樣時(shí)間進(jìn)行建模，它也支持多速率系統(tǒng)，也就是系統(tǒng)中的不同部分具有不同的采樣速率。為了創(chuàng)建動(dòng)態(tài)系統(tǒng)模型，Simulink提供了一個(gè)建立模型方塊圖的圖形用戶接口(GUI) ，這個(gè)創(chuàng)建過程只需單擊和拖動(dòng)鼠標(biāo)操作就能完成，它提供了一種更快捷、直接明了的方式，而且用戶可以立即看到系統(tǒng)的仿真結(jié)果。</p><p>

5、;<b>　　三．具體要求</b></p><p>　　a．利用所學(xué)的《通信原理及應(yīng)用》的基礎(chǔ)知識，設(shè)計(jì)一個(gè)2ASK數(shù)字調(diào)制器。完成對2ASK的調(diào)制與解調(diào)仿真電路設(shè)計(jì)，并對其仿真結(jié)果進(jìn)行分析。要求理解2ASK信號的產(chǎn)生，掌握2ASK信號的調(diào)制原理和實(shí)現(xiàn)方法并畫出實(shí)現(xiàn)框圖。</p><p>　　b．設(shè)計(jì)一個(gè)2FSK數(shù)字調(diào)制器。要求給出2FSK的產(chǎn)生原理框圖（調(diào)頻法、鍵控

6、法）、Matlab仿真調(diào)制解調(diào)的原理框圖，給出信號的頻譜圖、調(diào)制前與解調(diào)后數(shù)據(jù)波形圖</p><p>　　c．設(shè)計(jì)一個(gè)2PSK數(shù)字調(diào)制器。給出信號的頻譜圖、調(diào)制前與解調(diào)后數(shù)據(jù)波形圖.</p><p>　　d. 盡可能做出加噪前后相關(guān)波形。(加分項(xiàng))</p><p>　　數(shù)字通信系統(tǒng)的基本模型 </p><p>　　從消息傳輸角度看,該系統(tǒng)包

7、括了兩個(gè)重要交換,即消息與數(shù)字基帶信號之間的交換,數(shù)字基帶信號與信道信號之間的交換.通常前一種交換由發(fā)收端設(shè)備完成.而后一種交換則由調(diào)制和解調(diào)完成.</p><p><b>　　數(shù)字通信系統(tǒng)模型</b></p><p>　　一、2ASK調(diào)制解調(diào) </p><p><b>　　基本原理</b></p><

8、;p>　　2ASK是利用載波的幅度變化來傳遞數(shù)字信息，而其頻率和初始相位保持不變。</p><p>　　其信號表達(dá)式為： ，S (t)為單極性數(shù)字基帶信號。</p><p><b>　　2ASK幅移鍵控</b></p><p>　　幅移鍵控（ASK）相當(dāng)于模擬信號中的調(diào)幅，只不過與載頻信號相乘的

9、是二進(jìn)數(shù)碼而已。幅移就是把頻率、相位作為常量，而把振幅作為變量，信息比特是通過載波的幅度來傳遞的。由于調(diào)制信號只有0或1兩個(gè)電平，相乘的結(jié)果相當(dāng)于將載頻或者關(guān)斷，或者接通，它的實(shí)際意義是當(dāng)調(diào)制的數(shù)字信號"1時(shí)，傳輸載波；當(dāng)調(diào)制的數(shù)字信號為"0"時(shí)，不傳輸載波。</p><p>　　由圖可以看出2ASK信號的時(shí)間波形e2ASK（t）隨二進(jìn)制基帶信號s（t）通斷變化。所以又被稱為通斷鍵控

10、信號</p><p>　　2ASK信號的產(chǎn)生方法通常有兩種：模擬調(diào)制法和鍵控法。</p><p>　　模擬調(diào)制法使用乘法器實(shí)現(xiàn)</p><p>　　鍵控法使用開關(guān)電路實(shí)現(xiàn)</p><p><b>　　2ASK的調(diào)制方法</b></p><p>　　2ASK有兩種基本解調(diào)方法：相干解調(diào)法（同步檢測法

11、）和非相干解調(diào)法（包絡(luò)檢波法）。相干解調(diào)需要將載頻位置的已調(diào)信號頻譜重新搬回原始基帶位置，因此用相乘器與載波相乘來實(shí)現(xiàn)。為確保無失真還原信號，必須在接收端提供一個(gè)與調(diào)制載波嚴(yán)格同步的本地載波，這是整個(gè)解調(diào)過程能否順利完好進(jìn)行的關(guān)鍵。</p><p><b>　　相干解調(diào)</b></p><p><b>　　非相干解調(diào)</b></p>

12、<p>　　2ASK信號非相干解調(diào)過程的時(shí)間波形</p><p>　　振幅鍵控是利用載波的幅度變化來傳遞數(shù)字信息，而頻率和初始相位保持不變。</p><p><b>　　在2ASK中：</b></p><p>　　S2ask=m(t)*cos(2*pi*f*t)，</p><p>　　其中m(t)為數(shù)字信號，

13、后者為載波。</p><p>　　載波在二進(jìn)制基帶信號控制下通斷變化，所以又叫通-斷鍵控（OOK）。</p><p>　　2ASK的產(chǎn)生方法有兩種：模擬調(diào)制和鍵控法</p><p>　　而解調(diào)也有兩中基本方式：非相干解調(diào)（包絡(luò)檢波）和相干解調(diào)（同步檢測法）</p><p>　　DS2ask=s(t)*cos(2*pi*f*t)</p&g

14、t;<p>　　=0.5*m(t)+0.5*m(t)*cos(2*wc*t)</p><p>　　乘以相干載波后，只要濾去高頻部分就可以了</p><p>　　本次仿真使用相干解調(diào)方式：</p><p>　　2ask信號 帶通濾波器與與載波相乘低通濾波器 抽樣判決 輸出</p><p>　　以下就是matlab的仿真結(jié)果極其頻譜

15、圖（省去了帶通filter）</p><p>　　可以看到解調(diào)后的信號與信源有一定的延時(shí)。</p><p>　　通過觀察頻譜圖，用放大鏡可以清楚的看到，2ask實(shí)現(xiàn)了頻譜的搬移，將基帶信號</p><p>　　搬移到了fc=150hz的頻率上，而且若只計(jì)頻譜的主瓣則有：</p><p>　　B2ask=2fs，fs=1/Ts</p>

16、;<p>　　其中Ts為一個(gè)碼元寬度</p><p>　　即：2ask信號的傳輸帶寬是碼元傳輸速率的2倍</p><p>　　Matlab程序?qū)崿F(xiàn)</p><p><b>　　clc;</b></p><p>　　clear all;</p><p>　　close all;<

17、/p><p><b>　　%信源</b></p><p>　　a=randint(1,15,2);</p><p>　　t=0:0.001:0.999;</p><p>　　m=a(ceil(15*t+0.01)); </p><p>　　subplot(511)</p><p&

18、gt;　　plot(t,m);</p><p>　　axis([0 1.2 -0.2 1.2]);</p><p>　　title('信源');</p><p><b>　　%載波</b></p><p><b>　　f=150;</b></p><p>　　

19、carry=cos(2*pi*f*t);</p><p><b>　　%2ASK調(diào)制</b></p><p>　　st=m.*carry;</p><p>　　subplot(512);</p><p>　　plot(t,st)</p><p>　　axis([0 1.2 -1.2 1.2])&l

20、t;/p><p>　　title('2ASK信號') </p><p><b>　　%加高斯噪聲</b></p><p>　　nst=awgn(st,70);</p><p><b>　　%解調(diào)部分</b></p><p>　　nst=nst.*carry; &l

21、t;/p><p>　　subplot(513)</p><p>　　plot(t,nst)</p><p>　　axis([0 1.2 -0.2 1.2]);</p><p>　　title('乘以相干載波后的信號')</p><p><b>　　%低通濾波器設(shè)計(jì)</b></p&

22、gt;<p>　　wp=2*pi*2*f*0.5;</p><p>　　ws=2*pi*2*f*0.9; </p><p><b>　　Rp=2;</b></p><p><b>　　As=45;</b></p><p>　　[N,wc]=buttord(wp,ws,Rp,As,

23、9;s');</p><p>　　[B,A]=butter(N,wc,'s');</p><p><b>　　%低通濾波</b></p><p>　　h=tf(B,A); %轉(zhuǎn)換為傳輸函數(shù) </p><p>　　dst=lsim(h,nst,t);</p><p>　　s

24、ubplot(514)</p><p>　　plot(t,dst)</p><p>　　axis([0 1.2 -0.2 1.2]);</p><p>　　title('經(jīng)過低通濾波器后的信號');</p><p><b>　　%判決器</b></p><p><b>　

25、　k=0.25;</b></p><p>　　pdst=1*(dst>0.25);</p><p>　　subplot(515)</p><p>　　plot(t,pdst)</p><p>　　axis([0 1.2 -0.2 1.2]); </p><p>　　title('經(jīng)過抽樣判決后

26、的信號')</p><p><b>　　%頻譜觀察</b></p><p><b>　　%調(diào)制信號頻譜</b></p><p><b>　　T=t(end);</b></p><p><b>　　df=1/T;</b></p><

27、;p>　　N=length(st);</p><p>　　f=(-N/2:N/2-1)*df;</p><p>　　sf=fftshift(abs(fft(st)));</p><p><b>　　figure(2)</b></p><p>　　subplot(411)</p><p>　　

28、plot(f,sf)</p><p>　　title('調(diào)制信號頻譜')</p><p><b>　　%信源頻譜 </b></p><p>　　mf=fftshift(abs(fft(m)));</p><p>　　subplot(412)</p><p>　　plot(f,mf)

29、</p><p>　　title('信源頻譜')</p><p>　　% 乘以相干載波后的頻譜 </p><p>　　mmf=fftshift(abs(fft(nst)));</p><p>　　subplot(413)</p><p>　　plot(f,mmf)</p><p>

30、;　　title('乘以相干載波后的頻譜')</p><p>　　%經(jīng)過低通濾波后的頻譜</p><p>　　dmf=fftshift(abs(fft(dst)));</p><p>　　subplot(414)</p><p>　　plot(f,dmf)</p><p>　　title('經(jīng)過

31、低通濾波后的頻譜 );</p><p>　　二、2FSK調(diào)制解調(diào) </p><p>　　頻移鍵控是利用載波的頻率來傳遞數(shù)字信號，在2FSK中，載波的頻率隨著二進(jìn)制基帶信號在f1和f2兩個(gè)頻率點(diǎn)間變化，頻移鍵控是利用載波的頻移變化來傳遞數(shù)字信息的。在2FSK中，載波的頻率隨基帶信號在f1和f2兩個(gè)頻率點(diǎn)間變化。故其表達(dá)式為：</p><p>　　典型波形如下圖所示。由

32、圖可見。2FSK信號可以看作兩個(gè)不同載頻的ASK信號的疊加。因此2FSK信號的時(shí)域表達(dá)式又可以寫成：</p><p>　　2FSK數(shù)字系統(tǒng)的調(diào)制原理</p><p>　　2F S K調(diào)制就是使用兩個(gè)不同的頻率的載波信號來傳輸一個(gè)二進(jìn)制信息序列。可以用二進(jìn)制“1”來對應(yīng)于載頻f1，而“0”用來對應(yīng)于另一相載頻w2的已調(diào)波形，而這個(gè)可以用受矩形脈沖序列控制的開關(guān)電路對兩個(gè)不同的獨(dú)立的頻率源w1

33、、f2進(jìn)行選擇通。如下原理圖：</p><p><b>　　2FSK的解調(diào)方式</b></p><p>　　2FSK的解調(diào)方式有兩種:相干解調(diào)方式和非相干解調(diào)方式.下面我們將詳細(xì)的介紹:</p><p><b>　　1 非相干解調(diào)</b></p><p>　　經(jīng)過調(diào)制后的2FSK數(shù)字信號通過兩個(gè)頻率

34、不同的帶通濾波器f1、f2濾出不需要的信號，然后再將這兩種經(jīng)過濾波的信號分別通過包絡(luò)檢波器檢波，最后將兩種信號同時(shí)輸入到抽樣判決器同時(shí)外加抽樣脈沖，最后解調(diào)出來的信號就是調(diào)制前的輸入信號。其原理圖如下圖所示： </p><p><b>　　2 相干解調(diào) </b></p><p>　　根據(jù)已調(diào)信號由兩個(gè)載波f1、f2調(diào)制而成，則先用兩個(gè)分別對f1、f2帶通的濾波器對已調(diào)

35、信號進(jìn)行濾波，然后再分別將濾波后的信號與相應(yīng)的載波f1、f2相乘進(jìn)行相干解調(diào)，再分別低通濾波、用抽樣信號進(jìn)行抽樣判決器即可。原理圖如下：</p><p>　　Matlab程序?qū)崿F(xiàn)</p><p>　　Fc=150; %載頻</p><p>　　Fs=40; %系統(tǒng)采樣頻率</p><p>　　Fd=1; %碼速率</p&g

36、t;<p><b>　　N=Fs/Fd; </b></p><p><b>　　df=10;</b></p><p>　　numSymb=25;%進(jìn)行仿真的信息代碼個(gè)數(shù)</p><p>　　M=2; %進(jìn)制數(shù)</p><p>　　SNRpBit=60;%信噪比</p&

37、gt;<p>　　SNR=SNRpBit/log2(M);%60 </p><p>　　seed=[12345 54321];</p><p>　　numPlot=15; </p><p>　　x=randsrc(numSymb,1,[0:M-1]);%產(chǎn)生25個(gè)二進(jìn)制隨機(jī)碼</p><p><b>　　figur

38、e(1)</b></p><p>　　stem([0:numPlot-1],x(1:numPlot),'bx');%顯示15個(gè)碼元，桿圖，從x的前十五個(gè)隨機(jī)數(shù)中選取</p><p>　　title(' 二進(jìn)制隨機(jī)序列') </p><p>　　xlabel('Time');</p><p

39、>　　ylabel('Amplitude') ; </p><p><b>　　%調(diào)制</b></p><p>　　y=dmod(x,Fc,Fd,Fs,'fsk',M,df);%數(shù)字帶通調(diào)制</p><p>　　numModPlot=numPlot*Fs; %15*40</p><p&g

40、t;　　t=[0:numModPlot-1]./Fs;%數(shù)組除法（仿真時(shí)間） </p><p><b>　　figure(2)</b></p><p>　　plot(t,y(1:length(t)),'b-');</p><p>　　axis([min(t) max(t) -1.5 1.5]);</p><p

41、>　　title('調(diào)制后的信號')</p><p>　　xlabel('Time');</p><p>　　ylabel('Amplitude');</p><p>　　%在已調(diào)信號中加入高斯白噪聲</p><p>　　randn('state',seed(2));

42、 %生成-2到+2之間的隨機(jī)數(shù)矩陣</p><p>　　y=awgn(y,SNR-10*log10(0.5)-10*log10(N),'measured',[],'dB');%在已調(diào)信號中加入高斯白噪聲</p><p><b>　　figure(3)</b></p><p>　　plot(t,y(1:lengt

43、h(t)),'b-');%畫出經(jīng)過信道的實(shí)際信號</p><p>　　axis([min(t) max(t) -1.5 1.5]);</p><p>　　title('加入高斯白噪聲后的已調(diào)信號')</p><p>　　xlabel('Time');</p><p>　　ylabel('

44、Amplitude');</p><p><b>　　%相干解調(diào)</b></p><p><b>　　figure(4)</b></p><p>　　z1=ddemod(y,Fc,Fd,Fs,'fsk/eye',M,df);</p><p>　　title('相干解調(diào)

45、后的信號的眼圖')</p><p>　　%帶輸出波形的相干M元頻移鍵控解調(diào)</p><p><b>　　figure(5)</b></p><p>　　stem([0:numPlot-1],x(1:numPlot),'bx');</p><p><b>　　hold on;</b&

46、gt;</p><p>　　stem([0:numPlot-1],z1(1:numPlot),'ro');</p><p><b>　　hold off;</b></p><p>　　axis([0 numPlot -0.5 1.5]);</p><p>　　title('相干解調(diào)后的信號原序列比

47、較')</p><p>　　legend('原輸入二進(jìn)制隨機(jī)序列','相干解調(diào)后的信號')</p><p>　　xlabel('Time');</p><p>　　ylabel('Amplitude');</p><p><b>　　%非相干解調(diào)</b&g

48、t;</p><p><b>　　figure(6)</b></p><p>　　z2=ddemod(y,Fc,Fd,Fs,'fsk/eye/noncoh',M,df);</p><p>　　title('非相干解調(diào)后的信號的眼圖')</p><p>　　%帶輸出波形的非相干M元頻移鍵控解

49、調(diào)</p><p><b>　　figure(7)</b></p><p>　　stem([0:numPlot-1],x(1:numPlot),'bx');</p><p><b>　　hold on;</b></p><p>　　stem([0:numPlot-1],z2(1:nu

50、mPlot),'ro');</p><p><b>　　hold off;</b></p><p>　　axis([0 numPlot -0.5 1.5]);</p><p>　　title('非相干解調(diào)后的信號')</p><p>　　legend('原輸入二進(jìn)制隨機(jī)序列'

51、;,'非相干解調(diào)后的信號')</p><p>　　xlabel('Time');</p><p>　　ylabel('Amplitude');</p><p><b>　　%誤碼率統(tǒng)計(jì)</b></p><p>　　[errorSym ratioSym]=symerr(x,z

52、1);</p><p><b>　　figure(8)</b></p><p>　　simbasebandex([0:1:5]);</p><p>　　title('相干解調(diào)后誤碼率統(tǒng)計(jì)')</p><p>　　[errorSym ratioSym]=symerr(x,z2); </p&g

53、t;<p><b>　　figure(9)</b></p><p>　　simbasebandex([0:1:5]);</p><p>　　title('非相干解調(diào)后誤碼率統(tǒng)計(jì)')</p><p><b>　　%濾除高斯白噪聲</b></p><p>　　Delay=3

54、;R=0.5;PropD=0; %滯后3s</p><p>　　[yf,tf]=rcosine(Fd,Fs,'fir',R,Delay); %升余弦函數(shù)</p><p>　　[yo2,to2]=rcosflt(y,Fd,Fs,'filter',yf);</p><p>　　%加入高斯白噪聲后的已調(diào)信號和經(jīng)過升余弦濾波器后的

55、已調(diào)信號</p><p>　　t=[0:numModPlot-1]./Fs;</p><p>　　figure(10)</p><p>　　plot(t,y(1:length(t)),'r-');</p><p><b>　　hold on;</b></p><p>　　plot(

56、to2,yo2,'b-');%濾出帶外噪聲</p><p><b>　　hold off;</b></p><p>　　axis([0 30 -1.5 1.5]);</p><p>　　xlabel('Time');</p><p>　　ylabel('Amplitude'

57、);</p><p>　　legend('加入高斯白噪聲后的已調(diào)信號','經(jīng)過升余弦濾波器后的已調(diào)信號')</p><p>　　title('升余弦濾波前后波形比較')</p><p>　　eyediagram(yo2,N);%眼圖</p><p>　　title('加入高斯白噪聲后的已

58、調(diào)信號的眼圖') </p><p><b>　　仿真結(jié)果 </b></p><p>　　三、2psk信號調(diào)制解調(diào)</p><p>　　2psk信號的調(diào)制不能采用包絡(luò)檢測的方法，只能進(jìn)行相干解調(diào)，其原理框圖如下：</p><p>　　不考慮噪聲時(shí)，帶通濾波器輸出可以表示為</p><p>　

59、　y(t)=cos(wct+Φn)</p><p>　　式中Φn為2psk信號某一碼元的初相。Φn=0時(shí)，代表數(shù)字“0”， Φn=π時(shí)，代表數(shù)字“1”。與同步載波COSwct相乘后，輸出為</p><p>　　Z(t)=COS(wct+Φn) COSwct=1/2cosΦn+1/2cos(2wct+Φn)</p><p>　　經(jīng)過低通濾波器濾除高頻分量，得解調(diào)輸出為

60、</p><p>　　根據(jù)發(fā)送端產(chǎn)生2psk信號時(shí)Φn代表數(shù)字信息1或0的規(guī)定，以及接收端x（t）與Φn的關(guān)系特性，抽樣判決器的判決準(zhǔn)則為 </p><p>　　其中，x為x（t）在抽樣時(shí)刻的值。</p><p>　　2psk信號相干解調(diào)的過程實(shí)際上就是輸入已調(diào)信號與本地載波信號進(jìn)行極性比較的過程，故常稱為極性比較解調(diào)。 </p><p>　

61、　Matlab程序?qū)崿F(xiàn)</p><p><b>　　clear ;</b></p><p>　　close all; </p><p>　　fs=8e5; %抽樣頻率</p><p>　　fm=20e3; %基帶頻率</p><p>　　n=2

62、*(6*fs/fm);</p><p>　　final=(1/fs)*(n-1);</p><p>　　fc=2e5; % 載波頻率</p><p>　　t=0:1/fs:(final);</p><p>　　Fn=fs/2; %耐奎斯特頻率 </p><p>&l

63、t;b>　　%用正弦波產(chǎn)生方波</b></p><p>　　twopi_fc_t=2*pi*fm*t; </p><p><b>　　A=1;</b></p><p><b>　　phi=0;</b></p><p>　　x = A * cos(twopi_fc_t + phi);

64、</p><p><b>　　% 方波</b></p><p><b>　　am=1;</b></p><p>　　x(x>0)=am;</p><p>　　x(x<0)=-1;</p><p><b>　　figure(1)</b><

65、/p><p>　　subplot(321);</p><p>　　plot(t,x);</p><p>　　axis([0 2e-4 -2 2]);</p><p>　　title('基帶信號');</p><p><b>　　grid on</b></p><p&

66、gt;　　car=sin(2*pi*fc*t); %載波</p><p>　　ask=x.*car; %載波調(diào)制</p><p>　　subplot(322);</p><p>　　plot(t,ask);</p><p>　　axis([0 200e-6 -2 2]);</p><p>　　t

67、itle('PSK信號');</p><p><b>　　grid on;</b></p><p>　　%=====================================================</p><p><b>　　vn=0.1;</b></p><p>　

68、　noise=vn*(randn(size(t))); %產(chǎn)生噪音</p><p>　　subplot(323);</p><p>　　plot(t,noise);</p><p><b>　　grid on;</b></p><p>　　title('噪音信號');</p><p&g

69、t;　　axis([0 .2e-3 -1 1]);</p><p>　　askn=(ask+noise); %調(diào)制后加噪</p><p>　　subplot(324);</p><p>　　plot(t,askn);</p><p>　　axis([0 200e-6 -2 2]);</p><p>　　t

70、itle('加噪后信號');</p><p><b>　　grid on;</b></p><p><b>　　%帶通濾波</b></p><p><b>　　fBW=40e3;</b></p><p>　　f=[0:3e3:4e5];</p>&

71、lt;p>　　w=2*pi*f/fs;</p><p>　　z=exp(w*j);</p><p>　　BW=2*pi*fBW/fs;</p><p>　　a=.8547; %BW=2(1-a)/sqrt(a)</p><p>　　p=(j^2*a^2);</p><p>　　gai

72、n=.135;</p><p>　　Hz=gain*(z+1).*(z-1)./(z.^2-(p));</p><p>　　subplot(325);</p><p>　　plot(f,abs(Hz));</p><p>　　title('帶通濾波器');</p><p><b>　　grid

73、 on;</b></p><p>　　Hz(Hz==0)=10^(8); %avoid log(0)</p><p>　　subplot(326);</p><p>　　plot(f,20*log10(abs(Hz)));</p><p><b>　　grid on;</b></p>

74、<p>　　title('Receiver -3dB Filter Response');</p><p>　　axis([1e5 3e5 -3 1]);</p><p><b>　　%濾波器系數(shù)</b></p><p>　　a=[1 0 0.7305]; %[1 0 p]</p>&l

75、t;p>　　b=[0.135 0 -0.135]; %gain*[1 0 -1]</p><p>　　faskn=filter(b,a,askn);</p><p><b>　　figure(2)</b></p><p>　　subplot(321);</p><p>　　plot(t,faskn);<

76、;/p><p>　　axis([0 100e-6 -2 2]);</p><p>　　title('通過帶通濾波后輸出');</p><p><b>　　grid on;</b></p><p>　　cm=faskn.*car; %解調(diào)</p><p>　　subplot

77、(322);</p><p>　　plot(t,cm);</p><p>　　axis([0 100e-6 -2 2]);</p><p><b>　　grid on;</b></p><p>　　title('通過相乘器后輸出');</p><p><b>　　%低通濾

78、波器</b></p><p><b>　　p=0.72;</b></p><p>　　gain1=0.14;%gain=(1-p)/2</p><p>　　Hz1=gain1*(z+1)./(z-(p));</p><p>　　subplot(323);</p><p>　　Hz1(H

79、z1==0)=10^(-8);%avoid log(0)</p><p>　　plot(f,20*log10(abs(Hz1)));</p><p><b>　　grid on;</b></p><p>　　title('LPF -3dB response');</p><p>　　axis([0 5e4

80、 -3 1]);</p><p><b>　　%濾波器系數(shù)</b></p><p>　　a1=[1 -0.72]; %(z-(p))</p><p>　　b1=[0.14 0.14]; %gain*[1 1] </p><p>　　so=filter(b1,a1,cm);</p><p

81、>　　so=so*10; %add gain</p><p>　　so=so-mean(so); %removes DC component</p><p>　　subplot(324);</p><p>　　plot(t,so);</p><p>　　axis([0 8e-4 -3.5 3.5]);</p

82、><p>　　title('通過低通濾波器后輸出');</p><p><b>　　grid on;</b></p><p><b>　　%比較器</b></p><p><b>　　High=2.5;</b></p><p><b&g

83、t;　　Low=-2.5;</b></p><p>　　vt=0; %設(shè)立比較標(biāo)準(zhǔn)</p><p><b>　　error=0;</b></p><p>　　len1=length(so);</p><p>　　for ii=1:len1</p><p&g

84、t;　　if so(ii) >= vt</p><p>　　Vs(ii)=High;</p><p><b>　　else</b></p><p>　　Vs(ii)=Low;</p><p><b>　　end</b></p><p><b>　　end<

85、;/b></p><p><b>　　Vo=Vs;</b></p><p>　　subplot(325);</p><p>　　plot (t,Vo), title('解調(diào)后輸出信號'),</p><p>　　axis([0 2e-4 -5 5])</p><p><b