AM與DSB信號(hào)的調(diào)制解調(diào)(含代碼)

./AM、DSB調(diào)制及解調(diào)用matlab產(chǎn)生一個(gè)頻率為1Hz,功率為1的余弦信源,設(shè)載波頻率,,試畫出：AM及DSB調(diào)制信號(hào)的時(shí)域波形；采用相干解調(diào)后的AM及DSB信號(hào)波形；AM及DSB已調(diào)信號(hào)的功率譜；調(diào)整載波頻率及m0,觀察分的AM的過調(diào)與DSB反相點(diǎn)現(xiàn)象。在接收端帶通后加上窄帶高斯噪聲,單邊功率譜密度,重新解調(diào)。一圖表輸出及結(jié)果分析1AM及DSB調(diào)制信號(hào)的時(shí)域波形AM信號(hào)的包絡(luò)波形與其調(diào)制信號(hào)的波形完全一樣。DSB調(diào)制信號(hào)在載波處反向。2采用相干解調(diào)后的AM及DSB信號(hào)波形AM信號(hào)可用包絡(luò)檢波器檢波,濾除直流后可恢復(fù)出原信號(hào)。與AM相比,DSB信號(hào)由于不存在載波分量,全部功率都用來傳輸信息。3AM及DSB已調(diào)信號(hào)的功率譜AM信號(hào)所需的傳輸帶寬為調(diào)制信號(hào)帶寬的兩倍。DSB節(jié)省了載波功率,但所需的傳輸帶寬仍為調(diào)制信號(hào)帶寬的兩倍。4調(diào)整載波頻率及m0,觀察分的AM的過調(diào)與DSB反相點(diǎn)現(xiàn)象。AM信號(hào)用于包絡(luò)檢波器檢波時(shí),當(dāng)調(diào)制信號(hào)的最大值大于直流信號(hào)時(shí),會(huì)出現(xiàn)"過調(diào)制"現(xiàn)象,此時(shí)會(huì)發(fā)生失真。DSB相位在0.2π角頻率處頻率處明顯存在反相點(diǎn)現(xiàn)象。5在接收端帶通后加上窄帶高斯噪聲,單邊功率譜密度,重新解調(diào)。加入相同功率的高斯白噪聲信號(hào),DSB信號(hào)解調(diào)較AM信號(hào)解調(diào)誤差較大。二代碼clc;closeall;clearall;fm=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>;plot<t,s_am>;holdon;plot<t,A+mt,'r--'>;title<'AM調(diào)制信號(hào)'>;xlabel<'t/s'>;ylabel<'幅度'>;subplot<2,1,2>;plot<t,s_db>;title<'DSB調(diào)制信號(hào)'>;xlabel<'t/s'>;ylabel<'幅度'>;holdon;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>;bt=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>>;holdon;plot<t,mt,'r--'>;axis<[05-1.51.5]>;xlabel<'t/s'>,ylabel<'m<t>'>,title<'AM相干解調(diào)'>;subplot<2,1,2>;plot<t,y2<1:1024>>;holdon;plot<t,mt,'r--'>;axis<[05-1.51.5]>;xlabel<'t/s'>,ylabel<'m<t>'>,title<'DSB相干解調(diào)'>;S=fft<s_am,1024>;p=S.*conj<S>/1024;S2=fft<s_db,1024>;p2=S2.*conj<S2>/1024;f=0:0.1:102.3;figure<3>;subplot<2,1,1>;plot<f,p>;title<'AM已調(diào)信號(hào)功率譜'>;xlabel<'f/hz'>;ylabel<'w'>;subplot<2,1,2>;plot<f,p2>;title<'DSB已調(diào)信號(hào)功率譜'>;xlabel<'f/hz'>;ylabel<'w'>;s_am1=<A-1+mt>.*cos<2*pi*fc*t>;figure<4>;subplot<2,1,1>;plot<t,s_am>;holdon;plot<t,A+mt,'r--'>;title<'AM不過調(diào)'>;xlabel<'t/s'>;ylabel<'幅度'>;subplot<2,1,2>;plot<t,s_am1>;title<'AM過調(diào)'>;xlabel<'t/s'>;ylabel<'幅度'>;f=<0:1023>/1024;S1=fft<s_db,1024>;xp=abs<S1>;gl=angle<S1>;figure<5>;subplot<2,1,1>;plot<f,xp>;title<'DSB調(diào)幅波'>;xlabel<'w/<2*pi>'>;ylabel<'幅度'>;subplot<2,1,2>;plot<f,gl>;xlabel<'w/<2*pi>'>;ylabel<'相位'>;holdon;f=0:0.1:102.3;s_am=s_am+wgn<1,length<t>,0.1>;s_db=s_db+wgn<1,length<t>,0.1>;sp1=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>;bt=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<sp1,hn>;y2=conv<sp2,hn>;y=2*y;y2=2*y2;figure<6>;subplot<2,1,1>;plot<t,y<1:1024>>;holdon;plot<t,mt,'r--'>;axis<[05-1.51.5]>;xlabel<'t/s'>,ylabel<'m<t>'