MATLAB通信原理仿真五章作業(yè)

1、通信原理仿真作業(yè)第五章 模擬調(diào)制1. AM、DSB調(diào)制及解調(diào)用matlab產(chǎn)生一個頻率為1Hz，功率為1的余弦信源，設(shè)載波頻率，,試畫出：l AM及DSB調(diào)制信號的時域波形；l 采用相干解調(diào)后的AM及DSB信號波形；l AM及DSB已調(diào)信號的功率譜；l 調(diào)整載波頻率及m0，觀察分的AM的過調(diào)與DSB反相點現(xiàn)象。l 在接收端帶通后加上窄帶高斯噪聲，單邊功率譜密度，重新解調(diào)。2. SSB調(diào)制及解調(diào)用matlab產(chǎn)生一個頻率為1Hz，功率為1的余弦信源，設(shè)載波頻率，,試畫出：l SSB調(diào)制信號的時域波形；l 采用相干解調(diào)后的SSB信號波形；l SSB已調(diào)信號的功率譜；l 在接收端帶通后加上窄帶高斯噪

2、聲，單邊功率譜密度，重新解調(diào)。3. FM調(diào)制及解調(diào)設(shè)輸入信號為，載波中心頻率為，VCO的壓控振蕩系數(shù)為，載波平均功率為1W。試畫出：l 已調(diào)信號的時域波形；l 已調(diào)信號的振幅譜；l 用鑒頻器解調(diào)該信號，并與輸入信號比較。 加入相同功率的高斯白噪聲信號，DSB信號解調(diào)較AM信號解調(diào)誤差較大。 二代碼 clc; close all; clear all; fm=1; fc=10; T=5; t=linspace(-5,5,1024);%采樣頻率 a=sqrt(2); m

3、t=a*cos(2*pi*fm*t);% 帶通濾波器 A=2; s_db=mt.*cos(2*pi*fc*t); s_am=(A+mt).*cos(2*pi*fc*t);  figure(1); subplot(2,1,1); 00.511.522.533.544.55-101t/sm(t)AM相干解調(diào)00.511.522.533.544.55-101t/sm(t)DSB相干解調(diào)plot(t,s_am); hold on;  plot(t,A+mt,'r-')

4、;  title('AM調(diào)制信號'); xlabel('t/s'); ylabel('幅度'); subplot(2,1,2); plot(t,s_db); title('DSB調(diào)制信號'); xlabel('t/s'); ylabel('幅度'); hold on;  plot(t,mt,'r-');   sp=s_am.*

5、cos(2*pi*fc*t); sp2=s_db.*cos(2*pi*fc*t); fs1=0.3;fp1=0.5;fpu=1.8;fsu=2.0;Fs=50;Rp=1;Rs=42; ws1=2*pi*fs1/Fs;wp1=2*pi*fp1/Fs;wpu=2*pi*fpu/Fs;wsu=2*pi*fsu/Fs; wd1=0.5*(wp1+ws1);wd2=0.5*(wpu+wsu); bt=min(abs(ws1-wp1),abs(wsu-wpu); n=ceil(6.6*pi/bt); wn=wd1/pi,wd2/pi;&

6、#160;hn=fir1(n-1,wn,'bandpass',hamming(n); y=conv(sp,hn); y2=conv(sp2,hn); y=2*y; y2=2*y2; figure(2); subplot(2,1,1); plot(t,y(1:1024); hold on; plot(t,mt,'r-'); axis(0 5 -1.5 1.5); xlabel('t/s'),ylabe

7、l('m(t)'),title('AM相干解調(diào)'); subplot(2,1,2); plot(t,y2(1:1024); hold on; plot(t,mt,'r-'); axis(0 5 -1.5 1.5); xlabel('t/s'),ylabel('m(t)'),title('DSB相干解調(diào)AM、DSB調(diào)制及解調(diào) 用matlab產(chǎn)生一個頻率為1Hz，功率為1的余弦信源()mt，設(shè)載波頻率1

8、0cHzw=，02m=,試畫出： l AM及DSB調(diào)制信號的時域波形； l 采用相干解調(diào)后的AM及DSB信號波形； l AM及DSB已調(diào)信號的功率譜； l 調(diào)整載波頻率及m0，觀察分的AM的過調(diào)與DSB反相點現(xiàn)象。 l 在接收端帶通后加上窄帶高斯噪聲，單邊功率譜密度00.1n=，重新解調(diào)。 一圖表輸出及結(jié)果分析  1AM及DSB調(diào)制信號的時域波形   AM信號的包絡(luò)波形與其調(diào)制信號的波形完全一樣。DSB調(diào)制信號在載波處反向。 

9、  2 采用相干解調(diào)后的AM及DSB信號波形  AM信號可用包絡(luò)檢波器檢波，濾除直流后可恢復出原信號。與AM相比，DSB信號由于不存在載波分量，全部功率都用來傳輸信息。 -5-4-3-2-1012345-505AM調(diào)制信號t/s幅度-5-4-3-2-1012345-202DSB調(diào)制信號t/s幅度00.511.522.533.544.55-101t/sm(t)AM相干解調(diào)00.511.522.533.544.55-101t/sm(t)DSB相干解調(diào)  3  AM及DSB已調(diào)信號的功率譜 

10、 AM信號所需的傳輸帶寬為調(diào)制信號帶寬的兩倍。DSB節(jié)省了載波功率，但所需的傳輸帶寬仍為調(diào)制信號帶寬的兩倍。 4 調(diào)整載波頻率及m0，觀察分的AM的過調(diào)與DSB反相點現(xiàn)象。   AM信號用于包絡(luò)檢波器檢波時，當調(diào)制信號的最大值大于直流信號時，會出現(xiàn)“過調(diào)制”現(xiàn)象，此時會發(fā)生失真。 DSB相位在0.2角頻率處頻率處明顯存在反相點現(xiàn)象。  5 在接收端帶通后加上窄帶高斯噪聲，單邊功率譜密度00.1n=，重新解調(diào)。 02040608010012005001000AM已調(diào)信號功率譜f/hzw0

11、204060801001200100200DSB已調(diào)信號功率譜f/hzw-5-4-3-2-1012345-505AM不過調(diào)t/s幅度-5-4-3-2-1012345-505AM過調(diào)t/s幅度00.10.20.30.40.50.60.70.80.910200400DSB調(diào)幅波w/(2*pi)幅度00.10.20.30.40.50.60.70.80.91-505w/(2*pi)相位 加入相同功率的高斯白噪聲信號，DSB信號解調(diào)較AM信號解調(diào)誤差較大。 二代碼 clc; close all; clear all; fm=

12、1; fc=10; T=5; t=linspace(-5,5,1024);%采樣頻率 a=sqrt(2); mt=a*cos(2*pi*fm*t);% 帶通濾波器 A=2; s_db=mt.*cos(2*pi*fc*t); s_am=(A+mt).*cos(2*pi*fc*t);  figure(1); subplot(2,1,1); 00.511.522.533.544.55-101t/sm(t)AM相干解調(diào)00.511.522.533.544.55-101t/s

13、m(t)DSB相干解調(diào)plot(t,s_am); hold on;  plot(t,A+mt,'r-');  title('AM調(diào)制信號'); xlabel('t/s'); ylabel('幅度'); subplot(2,1,2); plot(t,s_db); title('DSB調(diào)制信號'); xlabel('t/s'); ylabel('幅度');&#

14、160;hold on;  plot(t,mt,'r-');   sp=s_am.*cos(2*pi*fc*t); sp2=s_db.*cos(2*pi*fc*t); fs1=0.3;fp1=0.5;fpu=1.8;fsu=2.0;Fs=50;Rp=1;Rs=42; ws1=2*pi*fs1/Fs;wp1=2*pi*fp1/Fs;wpu=2*pi*fpu/Fs;wsu=2*pi*fsu/Fs; wd1=0.5*(wp1+ws1);wd2=0.5*(wpu+wsu); b

15、t=min(abs(ws1-wp1),abs(wsu-wpu); n=ceil(6.6*pi/bt); wn=wd1/pi,wd2/pi; hn=fir1(n-1,wn,'bandpass',hamming(n); y=conv(sp,hn); y2=conv(sp2,hn); y=2*y; y2=2*y2; figure(2); subplot(2,1,1); plot(t,y(1:1024); hold on; plot(t,mt,'r-'); axis(0 5 -1.5 1.5); xlabel('t/s'),ylabel('m(t)'),title('AM相干解調(diào)'); subplot(2,1,2); plot(t,y2(1:1024);