通信原理課程設(shè)計(jì)---基于matlab的模擬調(diào)制系統(tǒng)仿真與測(cè)試

1、<p>　　課 程 設(shè) 計(jì) 報(bào) 告</p><p>　　基于MATLAB的模擬調(diào)制系統(tǒng)仿真與測(cè)試</p><p>　　2012年12 月 27 日</p><p><b>　　目 錄</b></p><p>　　1、設(shè)計(jì)題目·······

2、··········3</p><p>　　設(shè)計(jì)原理·················3</p><p>　　設(shè)計(jì)主程序&

3、#183;···············4</p><p>　　設(shè)計(jì)調(diào)用子函數(shù)··············8</p>

4、;<p>　　Matlab程序運(yùn)行結(jié)果 ···········10</p><p>　　參考文獻(xiàn)················&

5、#183;13</p><p><b>　　一.設(shè)計(jì)題目</b></p><p>　　基于MATLAB的模擬調(diào)制系統(tǒng)仿真與測(cè)試</p><p><b>　　二.設(shè)計(jì)原理</b></p><p><b>　　1.1 模擬調(diào)制</b></p><p>　　?

6、模擬調(diào)制包括幅度調(diào)制（DSB，SSB，AM）和相角調(diào)制（頻率和相位調(diào)制）。</p><p>　　?幅度調(diào)制（線性調(diào)制）是正弦載波的幅度隨著調(diào)制信號(hào)而改變的調(diào)制方案。</p><p>　　?若調(diào)制信號(hào) ,頻譜為 ，帶寬為 ，SSB調(diào)制的帶寬為 ，DSB和AM調(diào)制的帶寬為2 。VSB-AM的帶寬在 -2 區(qū)間內(nèi)。</p><p>　　?AM調(diào)制與DSB調(diào)制在許多

7、方面十分相似，唯一區(qū)別在于AM調(diào)制用 代替DSB的 。</p><p><b>　　? </b></p><p><b>　　? </b></p><p><b>　　調(diào)制器模型如圖所示</b></p><p>　　?假設(shè)被調(diào)信號(hào) 是零均值信號(hào)，則調(diào)制信號(hào)的功率為 。&l

8、t;/p><p>　　? 為調(diào)制信號(hào)平均功率， 為邊帶功率， 為載波功率。</p><p>　　?調(diào)制信號(hào)中用于發(fā)送信息的功率和總功率的比值稱為調(diào)制效率， </p><p>　　?解調(diào)器輸入信噪比定義為, 為理想帶通濾波器帶寬,在理想信道中，當(dāng)帶通濾波器幅頻特性為常數(shù)1時(shí)， ， ，即可由解調(diào)器輸入信噪比計(jì)算出信道噪聲的單邊帶功率譜密度，而當(dāng)系統(tǒng)抽樣速率為時(shí)，產(chǎn)生的

9、高斯白噪聲帶寬為 ，由此可計(jì)算出信道中高斯白噪聲的平均功率，即方差 ，從而利用有關(guān)知識(shí)可以產(chǎn)生信道中所疊加的高斯白噪聲。</p><p>　　1.2信道加性高斯白噪聲</p><p>　　信道中加性高斯白噪聲功率由于其均值為0，故其方差 即為其平均功率，其中 為高斯白噪聲單邊帶功率譜密度；B為信道帶寬。在圖1中，信道噪聲的功率譜密度圖可以看出，當(dāng) =16Hz時(shí)，疊加于信道的高斯噪聲帶寬為8

10、Hz，信道中的加性高斯白噪聲通過(guò)帶寬為2Hz，幅度為1的理想帶通濾波器后，輸出的窄帶噪聲的平均功率即為相干解調(diào)器輸入噪聲的平均功率 ，其功率譜密度不變，仍為 ?？梢钥闯?， 與信道加性高斯白噪聲功率 之間有一定的關(guān)系，其共同點(diǎn)是其功率譜密度相同，從圖1也可以觀察出來(lái)，為-20dB。因此，在理想通信系統(tǒng)中，利用已給解調(diào)器輸入信噪比及已調(diào)信號(hào)功率和帶寬，可以計(jì)算出 ，從而算出信道加性高斯白噪聲的方差，由于其均值為0，故該方差為其平均功率，利用

11、它可以生成信道加性高斯白噪聲。轉(zhuǎn)換關(guān)系為：</p><p><b>　　， ，因此 ，</b></p><p>　　或 ，E是接收信號(hào)平均能量。</p><p>　　圖1 系統(tǒng)采樣頻率為16Hz時(shí)的噪聲</p><p><b>　　1.3 實(shí)驗(yàn)內(nèi)容</b></p><p>　

12、　調(diào)制信號(hào)為，利用AM調(diào)制方式調(diào)制載波 ，假設(shè) ，直流分量為3，采樣頻率1000Hz，解調(diào)器輸入信噪比為25dB，采用相干方式解調(diào)編寫matlab程序?qū)崿F(xiàn)AM信號(hào)的調(diào)制解調(diào)。</p><p>　　三、Matlab實(shí)現(xiàn)：</p><p><b>　　主程序</b></p><p><b>　　% AM調(diào)制解調(diào)</b><

13、/p><p>　　clear all;close all;echo on</p><p>　　%----------------系統(tǒng)仿真參數(shù)</p><p>　　A=3; %直流分量</p><p>　　fc=250; %載波頻率（Hz）</p><p>　　t0=0.15;%信號(hào)時(shí)長(zhǎng)</p>&l

14、t;p>　　snr=25; %解調(diào)器輸入信噪比dB</p><p>　　dt=0.001% 系統(tǒng)時(shí)域采樣間隔</p><p>　　fs=1/dt;%系統(tǒng)采樣頻率 </p><p>　　df = 0.2; %所需的頻率分辨率</p><p>　　t=0:dt:t0;</p><p>　　Lt=length(t

15、);%仿真過(guò)程中，信號(hào)長(zhǎng)度</p><p>　　snr_lin = 10^(snr/10); %解調(diào)器輸入信噪比</p><p>　　%-------------畫出調(diào)制信號(hào)波形及頻譜</p><p>　　% 產(chǎn)生模擬調(diào)制信號(hào)</p><p>　　m=[ones(1,t0/(3*dt)),-2*ones(1,t0/(3*dt)),zeros

16、(1,t0/(3*dt)+1)];</p><p>　　L=2*min(m);</p><p>　　R=2*max(abs(m))+A;</p><p>　　pause%畫出調(diào)制信號(hào)波形及頻譜</p><p><b>　　clf</b></p><p><b>　　figure(1)&l

17、t;/b></p><p>　　subplot(321); </p><p>　　plot(t,m(1:length(t)));% 畫出調(diào)制信號(hào)波形</p><p>　　axis([0 t0 -R/2 R/2]);</p><p>　　xlabel('t');</p><p>　　ylabel(&

18、#39;調(diào)制信號(hào)');</p><p>　　subplot(322); </p><p>　　[M,m,df1,f]=T2F(m,dt,df,fs);%求出調(diào)制信號(hào)頻譜</p><p>　　[Bw_eq]=signalband(M,df,t0);%求出信號(hào)等效帶寬</p><p>　　f_start=fc-Bw_eq;</p&g

19、t;<p>　　f_cutoff=fc+Bw_eq;</p><p>　　plot(f,fftshift(abs(M)))% 畫出調(diào)制信號(hào)頻譜</p><p>　　xlabel('f');ylabel('調(diào)制信號(hào)頻譜');</p><p>　　pause %畫出載波及頻譜</p><p>　　su

20、bplot(323); </p><p>　　c=cos(2*pi*fc*t);%載波</p><p>　　plot(t,c);</p><p>　　axis([0 t0 -1.2 1.2]);</p><p>　　xlabel('t');ylabel('載波');</p><p>　　

21、subplot(324)% 載波頻譜</p><p>　　[C,c,df1,f]=T2F(c,dt,df,fs);</p><p>　　plot(f,fftshift(abs(C)))% 畫出載波頻譜</p><p>　　xlabel('f');ylabel('載波頻譜');</p><p>　　pause%

22、畫已調(diào)信號(hào)及其頻譜</p><p>　　subplot(325)% 畫已調(diào)信號(hào)</p><p>　　u=(A+m(1:Lt)).*c(1:Lt);%已調(diào)信號(hào)</p><p>　　plot(t,u);% 畫出已調(diào)信號(hào)波形</p><p>　　axis([0 t0 -R R]);</p><p>　　xlabel('

23、;t');ylabel('已調(diào)信號(hào)');</p><p>　　subplot(326); </p><p>　　[U,u,df1,f]=T2F(u,dt,df,fs);</p><p>　　plot(f,fftshift(abs(U)))% 畫出已調(diào)信號(hào)頻譜</p><p>　　xlabel('f');

24、ylabel('已調(diào)信號(hào)頻譜');</p><p>　　%先根據(jù)所給信噪比產(chǎn)生高斯白噪聲</p><p>　　signal_power = power_x(u(1:Lt)); %已調(diào)信號(hào)的平均功率</p><p>　　noise_power=(signal_power*fs)/(snr_lin*4* Bw_eq);%求出噪聲方差（噪聲均值為0）&l

25、t;/p><p>　　noise_std = sqrt(noise_power); %噪聲標(biāo)準(zhǔn)偏差</p><p>　　noise = noise_std*randn(1,Lt); %產(chǎn)生噪聲</p><p>　　pause%畫出信道高斯白噪聲波形及頻譜，此時(shí)，噪聲已實(shí)現(xiàn)，為確知信號(hào)，可求其頻譜</p><p><b

26、>　　figure(2)</b></p><p>　　subplot(321); </p><p>　　plot(t,noise);% 畫出噪聲波形</p><p>　　axis([0 t0 -R R]);</p><p>　　xlabel('t');ylabel('噪聲信號(hào)');</p

27、><p>　　subplot(322); </p><p>　　[noisef,noise,df1,f]=T2F(noise,dt,df,fs);%噪聲頻譜</p><p>　　plot(f,fftshift(abs(noisef)))% 畫出噪聲頻譜</p><p>　　xlabel('f');ylabel('噪聲頻譜&

28、#39;);</p><p>　　pause%畫出疊加了噪聲的已調(diào)信號(hào)波形及頻譜</p><p>　　sam=u(1:Lt)+noise(1:Lt);%疊加了噪聲的已調(diào)信號(hào)</p><p>　　subplot(323); %畫出疊加了噪聲的已調(diào)信號(hào)波形</p><p>　　plot(t,sam);</p><p>　　

29、axis([0 t0 -R R]);</p><p>　　xlabel('t');</p><p>　　ylabel('信道中的信號(hào)');</p><p>　　subplot(324); </p><p>　　[samf,sam,df1,f]=T2F(sam,dt,df,fs);%求出疊加了噪聲的已調(diào)信號(hào)頻譜&

30、lt;/p><p>　　plot(f,fftshift(abs(samf)))% 畫出疊加了噪聲的已調(diào)信號(hào)頻譜</p><p>　　xlabel('f');</p><p>　　ylabel('信道中信號(hào)頻譜');</p><p>　　[H,f]=bp_f(length(sam),f_start,f_cutoff,

31、df1,fs,1);%求帶通濾波器</p><p>　　subplot(326); </p><p>　　plot(f,fftshift(abs(H)))% 畫出帶通濾波器</p><p>　　xlabel('f');ylabel('帶通濾波器');</p><p>　　pause%經(jīng)過(guò)理想帶通濾波器后的信號(hào)及

32、其頻譜</p><p>　　DEM = H.*samf; %濾波器輸出的頻譜</p><p>　　[dem]=F2T(DEM,fs);%濾波器的輸出波形 </p><p><b>　　figure(3)</b></p><p>　　subplot(321)%經(jīng)過(guò)理想帶通濾波器后的信號(hào)波形</p>

33、<p>　　plot(t,dem(1:Lt))%畫出經(jīng)過(guò)理想帶通濾波器后的信號(hào)波形</p><p>　　axis([0 t0 -R R]);</p><p>　　xlabel('t');</p><p>　　ylabel('理想BPF輸出信號(hào)');</p><p>　　[demf,dem,df1,f

34、]=T2F(dem(1:Lt),dt,df,fs);%求經(jīng)過(guò)理想帶通濾波器后信號(hào)頻譜</p><p>　　subplot(322)</p><p>　　plot(f,fftshift(abs(demf)));% 畫出經(jīng)過(guò)理想帶通濾波器后信號(hào)頻譜</p><p>　　xlabel('f');</p><p>　　ylabel(&

35、#39;理想BPF輸出信號(hào)頻譜');</p><p>　　%--------------和本地載波相乘，即混頻</p><p>　　pause%混頻后的信號(hào)，先畫本地載波及其頻譜</p><p>　　subplot(323)</p><p>　　plot(t,c(1:Lt));</p><p>　　axis([

36、0 t0 -1.2 1.2]);</p><p>　　xlabel('t');</p><p>　　ylabel('本地載波');</p><p>　　subplot(324)% 載波頻譜</p><p>　　[C,c,df1,f]=T2F(c(1:Lt),dt,df,fs);</p><p

37、>　　plot(f,fftshift(abs(C)))% 畫出載波頻譜</p><p>　　xlabel('f');</p><p>　　ylabel('本地載波頻譜');</p><p>　　pause%再畫混頻后信號(hào)及其頻譜</p><p>　　der=dem(1:Lt).*c(1:Lt);%混頻&l

38、t;/p><p>　　subplot(325)%畫出混頻后的信號(hào)</p><p>　　plot(t,der);</p><p>　　axis([0 t0 -R R]);</p><p>　　xlabel('t');</p><p>　　ylabel('混頻后的信號(hào)');</p>

39、<p>　　subplot(326)</p><p>　　[derf,der,df1,f]=T2F(der,dt,df,fs);%求混頻后的信號(hào)頻譜</p><p>　　plot(f,fftshift(abs(derf)))%畫出混頻后的信號(hào)的頻譜</p><p>　　xlabel('f');</p><p>　　

40、ylabel('混頻后信號(hào)頻譜');</p><p>　　pause%畫出理想低通濾波器</p><p><b>　　figure(4)</b></p><p>　　[LPF,f]=lp_f(length(der),Bw_eq,df1,fs,2);%求低通濾波器</p><p>　　subplot(322

41、)</p><p>　　plot(f,fftshift(abs(LPF)));% 畫出理想低通濾波器</p><p>　　xlabel('f');</p><p>　　ylabel('理想LPF');</p><p>　　pause%混頻信號(hào)經(jīng)理想低通濾波器后的頻譜及波形</p><p>

42、;　　DM = LPF.*derf; %理想低通濾波器輸出的頻譜</p><p>　　[dm]=F2T(DM,fs);%濾波器的輸出波形 </p><p>　　subplot(323)</p><p>　　plot(t,dm(1:Lt));%畫出經(jīng)過(guò)低通濾波器后的解調(diào)出的波形</p><p>　　axis([0 t0 -R R]);<

43、;/p><p>　　xlabel('t');</p><p>　　ylabel('LPF輸出信號(hào)');</p><p>　　subplot(324)</p><p>　　[dmf,dm,df1,f]=T2F(dm(1:Lt),dt,df,fs);%求LPF輸出信號(hào)的頻譜</p><p>　　

44、plot(f,fftshift(dmf));%畫出LPF輸出信號(hào)的頻譜</p><p>　　xlabel('f');</p><p>　　ylabel('LPF輸出信號(hào)頻譜');</p><p>　　axis([-fs/2 fs/2 0 0.5]);</p><p>　　%--------去除解調(diào)信號(hào)中的直流分

45、量</p><p>　　pause%去除解調(diào)信號(hào)中的直流分量</p><p>　　dmd=dm(1:Lt)-mean(dm(1:Lt));</p><p>　　subplot(325)</p><p>　　plot(t,dmd);%畫出恢復(fù)信號(hào)(去除直流分量)</p><p>　　axis([0 t0 -R/2 R/2

46、]);</p><p>　　xlabel('t');</p><p>　　ylabel('恢復(fù)信號(hào)');</p><p>　　[dmdf,dmd,df1,f]=T2F(dmd,dt,df,fs);%求恢復(fù)信號(hào)的頻譜</p><p>　　subplot(326)</p><p>　　plo

47、t(f,fftshift(dmdf));%畫出恢復(fù)信號(hào)的頻譜</p><p>　　xlabel('f');</p><p>　　ylabel('恢復(fù)信號(hào)的頻譜');</p><p>　　axis([-fs/2 fs/2 0 0.2]);</p><p>　　subplot(321); </p>&

48、lt;p>　　plot(t,m(1:Lt));% 畫出調(diào)制信號(hào)波形</p><p>　　axis([0 t0 -R/2 R/2]);</p><p>　　xlabel('t');</p><p>　　ylabel('調(diào)制信號(hào)');</p><p><b>　　調(diào)用子函數(shù)</b><

49、;/p><p><b>　　序列的傅立葉變換</b></p><p>　　function [M,m,df]=fftseq(m,ts,df)</p><p>　　%各參數(shù)含義與子函數(shù)T2F中的完全相同，完成</p><p>　　fs = 1/ts;</p><p>　　if nargin ==2<

50、;/p><p><b>　　n1 =0;</b></p><p><b>　　else</b></p><p>　　n1 = fs/df;</p><p><b>　　end </b></p><p>　　n2 = length(m);</p>

51、<p>　　n = 2^(max(nextpow2(n1),nextpow2(n2)));</p><p>　　M = fft(m,n);</p><p>　　m = [m,zeros(1,n-n2)];</p><p>　　df = fs/n;</p><p><b>　　計(jì)算信號(hào)功率</b></p&

52、gt;<p>　　function p=power_x(x)</p><p><b>　　%x:輸入信號(hào)</b></p><p>　　%p:返回信號(hào)的x功率</p><p>　　p=(norm(x).^2)./length(x);</p><p>　　信號(hào)從頻域轉(zhuǎn)換到時(shí)域</p><p&

53、gt;　　function [m]=F2T(M,fs)</p><p>　　%-------------------------輸入?yún)?shù)</p><p><b>　　%M：信號(hào)的頻譜</b></p><p>　　%fs:系統(tǒng)采樣頻率</p><p>　　%--------------------輸出(返回)參數(shù)</

54、p><p>　　%m:傅里葉逆變換后的信號(hào)，注意其長(zhǎng)度為2的整數(shù)次冪，利用其畫波形時(shí)，要注意選取m的一部分，選取長(zhǎng)度和所給時(shí)間序列t的長(zhǎng)度要一致，plot(t,m(1:length(t)))，否則會(huì)出錯(cuò)。</p><p>　　m = real(ifft(M))*fs; </p><p>　　信號(hào)從時(shí)域轉(zhuǎn)換到頻域</p><p>　　functio

55、n [M,m,df1,f]=T2F(m,ts,df,fs)</p><p>　　%------------------------輸入?yún)?shù)</p><p><b>　　%m:信號(hào)</b></p><p>　　%ts:系統(tǒng)時(shí)域采樣間隔</p><p>　　%df:所需的頻率分辨率</p><p>　

56、　%fs:系統(tǒng)采樣頻率</p><p>　　%---------------------輸出(返回)參數(shù)</p><p>　　%M:傅里葉變換后的頻譜序列</p><p>　　%m: 輸入信號(hào)參與過(guò)傅里葉變換后對(duì)應(yīng)的序列，需要注意的是，該序列與輸入信號(hào)m的區(qū)別，其長(zhǎng)度是不一樣的，輸入的m長(zhǎng)度不一定是2的整數(shù)次冪，而傅里葉變換要求輸入信號(hào)長(zhǎng)度為2的整數(shù)次冪，故傅里葉變

57、換前需對(duì)輸入的m信號(hào)進(jìn)行補(bǔ)零，其長(zhǎng)度有所增加，故輸出參數(shù)中的m為補(bǔ)零后的輸入信號(hào)，其長(zhǎng)度與輸入?yún)?shù)m不一樣，但與M,f長(zhǎng)度是一樣的，并且，其與時(shí)間序列t所對(duì)應(yīng)的序列m(1:length(t))與輸入?yún)?shù)中的m是一致的。</p><p>　　%df1:返回的頻率分辨率</p><p>　　%f:與M相對(duì)應(yīng)的頻率序列</p><p>　　[M,m,df1]=fftseq

58、(m,ts,df);</p><p>　　f = [0:df1:df1*(length(m)-1)] -fs/2; %頻率向量</p><p><b>　　M=M/fs;</b></p><p><b>　　帶通濾波器</b></p><p>　　Function[H,f]=bp_f(n,f_

59、start,f_cutoff,df1,fs,p)</p><p>　　%帶通濾波器函數(shù) 輸入設(shè)計(jì)的濾波器參數(shù)，產(chǎn)生帶通濾波器頻率特性函數(shù)H和頻率向量f</p><p>　　%------------------------輸入?yún)?shù)</p><p>　　%n 帶通濾波器的輸入信號(hào)長(zhǎng)度</p><p>　　%f_start 通帶起始

60、頻率</p><p>　　%f_cutoff 帶通濾波器的截止頻率</p><p>　　%df1 頻率分辨率</p><p>　　%fs 抽樣頻率</p><p>　　%p 濾波器幅度</p><p>　　%----------------------輸出(返回)參數(shù)

61、</p><p>　　%H 帶通濾波器頻率響應(yīng)</p><p><b>　　%f 頻率向量</b></p><p>　　%設(shè)計(jì)濾波器 </p><p>　　n_cutoff =

62、 floor(f_cutoff/df1); </p><p>　　n_start = floor(f_start/df1); </p><p>　　f = [0:df1:df1*(n-1)] -fs/2; %頻率向量</p><p>　　H = zeros(size(f));</p><p>　　H(n_start+1:n_cu

63、toff) = p*ones(1,n_cutoff-n_start);</p><p>　　H(length(f) - n_cutoff+1:length(f)-n_start) = p*ones(1,n_cutoff-n_start);</p><p><b>　　低通濾波器</b></p><p>　　function [H,f]=lp_f

64、(n,f_cutoff,df1,fs,p)</p><p>　　%低通濾波器函數(shù) 輸入設(shè)計(jì)的濾波器參數(shù)，產(chǎn)生低通濾波器頻率特性函數(shù)H和頻率向量f</p><p>　　%------------------------輸入?yún)?shù)</p><p>　　%n 低通濾波器的輸入信號(hào)長(zhǎng)度</p><p>　　%f_cutoff 低通濾波器的截

65、止頻率</p><p>　　%df1 頻率分辨率</p><p>　　%fs 抽樣頻率</p><p>　　%p 濾波器幅度</p><p>　　%---------------------輸出(返回)參數(shù)</p><p>　　%H 低通濾波器頻率響應(yīng)</p>

66、;<p><b>　　%f 頻率向量</b></p><p>　　n_cutoff = floor(f_cutoff/df1); %設(shè)計(jì)濾波器</p><p>　　f = [0:df1:df1*(n-1)] -fs/2; %頻率向量</p><p>　　H = zeros(size(f));</p>

67、<p>　　H(1:n_cutoff) = p*ones(1,n_cutoff);</p><p>　　H(length(f) - n_cutoff+1:length(f)) = p*ones(1,n_cutoff);</p><p><b>　　計(jì)算信號(hào)有效帶寬</b></p><p>　　function [Bw_eq]=sig

68、nalband(sf,df,T)</p><p><b>　　%計(jì)算信號(hào)等效帶寬</b></p><p><b>　　%sf：信號(hào)頻譜</b></p><p><b>　　%df:頻率分辨率</b></p><p><b>　　%T：信號(hào)持續(xù)時(shí)間</b>&

69、lt;/p><p>　　sf_max=max(abs(sf));</p><p>　　Bw_eq=sum(abs(sf).^2)*df/T/sf_max.^2</p><p>　　四 Matlab程序運(yùn)行結(jié)果：</p><p>　　figure(1) ：</p><p><b>　　figure(2)</b

70、></p><p><b>　　、</b></p><p>　　figure(3)：</p><p>　　figure(4)：</p><p><b>　　五.參考文獻(xiàn)</b></p><p>　　[1]西蒙?赫金.通信系統(tǒng)(第4版)[M].北京:電子工業(yè)出版社,200