1、 課 程 設 計 報 告基于matlab的模擬調制系統仿真與測試課程名稱 通信原理課程設計 專業班級 通信工程2班 學 號 14目 錄1、設計題目32、 設計原理33、 設計主程序44、 設計調用子函數85、 matlab程序運行結果 106、 參考文獻13一.設計題目基于matlab的模擬調制系統仿真與測試二.設計原理1.1 模擬調制l模擬調制包括幅度調制（dsb，ssb，am）和相角調制（頻率和相位調制）。l幅度調制（線性調制）是正弦載波的幅度隨著調制信號而改變的調制方案。l若調制信號 ,頻譜為 ，帶寬為 ，ssb調制的帶寬為 ，dsb和am調制的帶寬為2 。vsb-am的帶寬在 -2 區

2、間內。lam調制與dsb調制在許多方面十分相似，唯一區別在于am調制用 代替dsb的 。l l 調制器模型如圖所示 l假設被調信號 是零均值信號，則調制信號的功率為 。l 為調制信號平均功率， 為邊帶功率， 為載波功率。l調制信號中用于發送信息的功率和總功率的比值稱為調制效率， l解調器輸入信噪比定義為, 為理想帶通濾波器帶寬,在理想信道中，當帶通濾波器幅頻特性為常數1時， ， ，即可由解調器輸入信噪比計算出信道噪聲的單邊帶功率譜密度，而當系統抽樣速率為時，產生的高斯白噪聲帶寬為 ，由此可計算出信道中高斯白噪聲的平均功率，即方差 ，從而利用有關知識可以產生信道中所疊加的高斯白噪聲。1.2信道加

3、性高斯白噪聲信道中加性高斯白噪聲功率由于其均值為0，故其方差 即為其平均功率，其中 為高斯白噪聲單邊帶功率譜密度；b為信道帶寬。在圖1中，信道噪聲的功率譜密度圖可以看出，當 =16hz時，疊加于信道的高斯噪聲帶寬為8hz，信道中的加性高斯白噪聲通過帶寬為2hz，幅度為1的理想帶通濾波器后，輸出的窄帶噪聲的平均功率即為相干解調器輸入噪聲的平均功率 ，其功率譜密度不變，仍為 ?？梢钥闯觯?與信道加性高斯白噪聲功率 之間有一定的關系，其共同點是其功率譜密度相同，從圖1也可以觀察出來，為-20db。因此，在理想通信系統中，利用已給解調器輸入信噪比及已調信號功率和帶寬，可以計算出 ，從而算出信道加性高斯

4、白噪聲的方差，由于其均值為0，故該方差為其平均功率，利用它可以生成信道加性高斯白噪聲。轉換關系為： ， ，因此 ， 或 ，e是接收信號平均能量。 圖1 系統采樣頻率為16hz時的噪聲1.3 實驗內容調制信號為，利用am調制方式調制載波 ，假設 ，直流分量為3，采樣頻率1000hz，解調器輸入信噪比為25db，采用相干方式解調編寫matlab程序實現am信號的調制解調。三、matlab實現：主程序% am調制解調clear all;close all;echo on%-系統仿真參數 a=3; %直流分量fc=250; %載波頻率（hz）t0=0.15;%信號時長snr=25; %解調器輸入信噪比

5、dbdt=0.001% 系統時域采樣間隔fs=1/dt;%系統采樣頻率 df = 0.2; %所需的頻率分辨率t=0:dt:t0;lt=length(t);%仿真過程中，信號長度snr_lin = 10(snr/10); %解調器輸入信噪比%-畫出調制信號波形及頻譜% 產生模擬調制信號m=ones(1,t0/(3*dt),-2*ones(1,t0/(3*dt),zeros(1,t0/(3*dt)+1);l=2*min(m);r=2*max(abs(m)+a;pause%畫出調制信號波形及頻譜clffigure(1)subplot(321); plot(t,m(1:length(t);% 畫出調

6、制信號波形axis(0 t0 -r/2 r/2);xlabel(t);ylabel(調制信號);subplot(322); m,m,df1,f=t2f(m,dt,df,fs);%求出調制信號頻譜bw_eq=signalband(m,df,t0);%求出信號等效帶寬f_start=fc-bw_eq;f_cutoff=fc+bw_eq;plot(f,fftshift(abs(m)% 畫出調制信號頻譜xlabel(f);ylabel(調制信號頻譜);pause %畫出載波及頻譜subplot(323); c=cos(2*pi*fc*t);%載波plot(t,c);axis(0 t0 -1.2 1.2

7、);xlabel(t);ylabel(載波);subplot(324)% 載波頻譜c,c,df1,f=t2f(c,dt,df,fs);plot(f,fftshift(abs(c)% 畫出載波頻譜xlabel(f);ylabel(載波頻譜);pause% 畫已調信號及其頻譜subplot(325)% 畫已調信號u=(a+m(1:lt).*c(1:lt);%已調信號plot(t,u);% 畫出已調信號波形axis(0 t0 -r r);xlabel(t);ylabel(已調信號);subplot(326); u,u,df1,f=t2f(u,dt,df,fs);plot(f,fftshift(abs

8、(u)% 畫出已調信號頻譜xlabel(f);ylabel(已調信號頻譜);%先根據所給信噪比產生高斯白噪聲signal_power = power_x(u(1:lt); %已調信號的平均功率noise_power=(signal_power*fs)/(snr_lin*4* bw_eq);%求出噪聲方差（噪聲均值為0）noise_std = sqrt(noise_power); %噪聲標準偏差noise = noise_std*randn(1,lt); %產生噪聲pause%畫出信道高斯白噪聲波形及頻譜，此時，噪聲已實現，為確知信號，可求其頻譜figure(2)subplot(321); pl

9、ot(t,noise);% 畫出噪聲波形axis(0 t0 -r r);xlabel(t);ylabel(噪聲信號);subplot(322); noisef,noise,df1,f=t2f(noise,dt,df,fs);%噪聲頻譜plot(f,fftshift(abs(noisef)% 畫出噪聲頻譜xlabel(f);ylabel(噪聲頻譜);pause%畫出疊加了噪聲的已調信號波形及頻譜sam=u(1:lt)+noise(1:lt);%疊加了噪聲的已調信號subplot(323); %畫出疊加了噪聲的已調信號波形plot(t,sam);axis(0 t0 -r r);xlabel(t);

10、ylabel(信道中的信號);subplot(324); samf,sam,df1,f=t2f(sam,dt,df,fs);%求出疊加了噪聲的已調信號頻譜plot(f,fftshift(abs(samf)% 畫出疊加了噪聲的已調信號頻譜xlabel(f);ylabel(信道中信號頻譜);h,f=bp_f(length(sam),f_start,f_cutoff,df1,fs,1);%求帶通濾波器subplot(326); plot(f,fftshift(abs(h)% 畫出帶通濾波器xlabel(f);ylabel(帶通濾波器);pause%經過理想帶通濾波器后的信號及其頻譜dem = h.*

11、samf; %濾波器輸出的頻譜dem=f2t(dem,fs);%濾波器的輸出波形 figure(3)subplot(321)%經過理想帶通濾波器后的信號波形plot(t,dem(1:lt)%畫出經過理想帶通濾波器后的信號波形axis(0 t0 -r r);xlabel(t);ylabel(理想bpf輸出信號);demf,dem,df1,f=t2f(dem(1:lt),dt,df,fs);%求經過理想帶通濾波器后信號頻譜subplot(322)plot(f,fftshift(abs(demf);% 畫出經過理想帶通濾波器后信號頻譜xlabel(f);ylabel(理想bpf輸出信號頻譜);%-和

12、本地載波相乘，即混頻pause%混頻后的信號，先畫本地載波及其頻譜subplot(323)plot(t,c(1:lt);axis(0 t0 -1.2 1.2);xlabel(t);ylabel(本地載波);subplot(324)% 載波頻譜c,c,df1,f=t2f(c(1:lt),dt,df,fs);plot(f,fftshift(abs(c)% 畫出載波頻譜xlabel(f);ylabel(本地載波頻譜);pause%再畫混頻后信號及其頻譜der=dem(1:lt).*c(1:lt);%混頻subplot(325)%畫出混頻后的信號plot(t,der);axis(0 t0 -r r);

13、xlabel(t);ylabel(混頻后的信號);subplot(326)derf,der,df1,f=t2f(der,dt,df,fs);%求混頻后的信號頻譜plot(f,fftshift(abs(derf)%畫出混頻后的信號的頻譜xlabel(f);ylabel(混頻后信號頻譜);pause%畫出理想低通濾波器figure(4)lpf,f=lp_f(length(der),bw_eq,df1,fs,2);%求低通濾波器subplot(322)plot(f,fftshift(abs(lpf);% 畫出理想低通濾波器xlabel(f);ylabel(理想lpf);pause%混頻信號經理想低通

14、濾波器后的頻譜及波形dm = lpf.*derf; %理想低通濾波器輸出的頻譜dm=f2t(dm,fs);%濾波器的輸出波形 subplot(323)plot(t,dm(1:lt);%畫出經過低通濾波器后的解調出的波形axis(0 t0 -r r);xlabel(t);ylabel(lpf輸出信號);subplot(324)dmf,dm,df1,f=t2f(dm(1:lt),dt,df,fs);%求lpf輸出信號的頻譜plot(f,fftshift(dmf);%畫出lpf輸出信號的頻譜xlabel(f);ylabel(lpf輸出信號頻譜);axis(-fs/2 fs/2 0 0.5);%-去除

15、解調信號中的直流分量pause%去除解調信號中的直流分量dmd=dm(1:lt)-mean(dm(1:lt);subplot(325)plot(t,dmd);%畫出恢復信號(去除直流分量)axis(0 t0 -r/2 r/2);xlabel(t);ylabel(恢復信號);dmdf,dmd,df1,f=t2f(dmd,dt,df,fs);%求恢復信號的頻譜subplot(326)plot(f,fftshift(dmdf);%畫出恢復信號的頻譜xlabel(f);ylabel(恢復信號的頻譜);axis(-fs/2 fs/2 0 0.2);subplot(321); plot(t,m(1:lt)

16、;% 畫出調制信號波形axis(0 t0 -r/2 r/2);xlabel(t);ylabel(調制信號);調用子函數l 序列的傅立葉變換function m,m,df=fftseq(m,ts,df)%各參數含義與子函數t2f中的完全相同，完成fs = 1/ts;if nargin =2 n1 =0;else n1 = fs/df;end n2 = length(m);n = 2(max(nextpow2(n1),nextpow2(n2);m = fft(m,n);m = m,zeros(1,n-n2);df = fs/n;l 計算信號功率function p=power_x(x)%x:輸入信

17、號%p:返回信號的x功率p=(norm(x).2)./length(x);l 信號從頻域轉換到時域function m=f2t(m,fs)%-輸入參數%m：信號的頻譜%fs:系統采樣頻率%-輸出(返回)參數%m:傅里葉逆變換后的信號，注意其長度為2的整數次冪，利用其畫波形時，要注意選取m的一部分，選取長度和所給時間序列t的長度要一致，plot(t,m(1:length(t)，否則會出錯。m = real(ifft(m)*fs; l 信號從時域轉換到頻域function m,m,df1,f=t2f(m,ts,df,fs)%-輸入參數%m:信號%ts:系統時域采樣間隔%df:所需的頻率分辨率%fs

18、:系統采樣頻率%-輸出(返回)參數%m:傅里葉變換后的頻譜序列%m: 輸入信號參與過傅里葉變換后對應的序列，需要注意的是，該序列與輸入信號m的區別，其長度是不一樣的，輸入的m長度不一定是2的整數次冪，而傅里葉變換要求輸入信號長度為2的整數次冪，故傅里葉變換前需對輸入的m信號進行補零，其長度有所增加，故輸出參數中的m為補零后的輸入信號，其長度與輸入參數m不一樣，但與m,f長度是一樣的，并且，其與時間序列t所對應的序列m(1:length(t)與輸入參數中的m是一致的。%df1:返回的頻率分辨率%f:與m相對應的頻率序列m,m,df1=fftseq(m,ts,df);f = 0:df1:df1*(

19、length(m)-1) -fs/2; %頻率向量m=m/fs;l 帶通濾波器functionh,f=bp_f(n,f_start,f_cutoff,df1,fs,p)%帶通濾波器函數 輸入設計的濾波器參數，產生帶通濾波器頻率特性函數h和頻率向量f%-輸入參數%n 帶通濾波器的輸入信號長度%f_start 通帶起始頻率%f_cutoff 帶通濾波器的截止頻率%df1 頻率分辨率%fs 抽樣頻率%p 濾波器幅度%-輸出(返回)參數%h 帶通濾波器頻率響應%f 頻率向量%設計濾波器 n_cutoff = floor(f_cutoff/df1); n_start = floor(f_start/df1); f = 0:df1:df1*(n-1) -fs/2; %頻率向量h = zeros(size(f);h(n_start+1:n_cutoff) = p*ones(1,n_cutoff-n_start);h(length(f) - n_cutoff+1:length(f)-n_start) = p*ones(1,n_cutoff-n_start);l 低通濾波器function h,f=lp_f(n,f_cutoff,df1,fs,p)%低通濾波