【文章內(nèi)容簡介】 為 ? ?12(t)o nt? 故輸出噪聲功率為 201 1 1( t)4 4 4o i iN n N n B? ? ? 式中， B=2,為 DSB 的帶通濾波器的帶寬， 0n 為噪聲單邊功率譜密度。 解調(diào)器輸入信號平均功率為 ? ?2 2 21( t) [ ( t) c o s ] ( t)2i m cSs tmm?? ? ? 可得解調(diào)器的輸入信噪比 201 (t)2iimSN n B?，解調(diào)器的輸出信噪比 20(t)ooS mN n B? 因此制度增益為 2DSBG ? ，也就是說， DSB 信號的解調(diào)器使信噪比改善一倍。 武漢理工大學《 MATLAB 課程設計》報告 7 調(diào)制解調(diào)分析的 MATLAB 實現(xiàn) 信號 DSB 調(diào)制采用 MATLAB 函數(shù) modulate 實現(xiàn)，其函數(shù)格式為： Y = MODULATE(X,fc,fs,METHOD,OPT) X 為基帶調(diào)制信號， fc 為載波頻率， fs 為抽樣頻率， METHOD 為調(diào)制方式選擇， DSB調(diào)制時為 ’am’， OPT 在 DSB 調(diào)制時可不選， fs 需滿足 fs 2*fc + BW， BW為調(diào)制信號帶寬。 DSB 信號解調(diào)采用 MATLAB 函數(shù) demod 實現(xiàn)，其函數(shù)使用格式為： X = DEMOD(Y,fc,fs,METHOD,OPT) Y 為 DSB 已調(diào)信號， fc 為載波頻率， fs 為抽樣頻率， METHOD 為解調(diào)方式選擇， DSB解調(diào)時為 ’am’， OPT 在 DSB 調(diào)制時可不選。 觀察信號頻譜需對信號進行傅里葉變換，采用 MATLAB函數(shù) fft實現(xiàn)，其函數(shù)常使用格式為： Y=FFT(X,N)， X為時域函數(shù)， N為傅里葉變換點數(shù)選擇，一般取值。頻域變換后，對頻域函數(shù)取模，格式： Y1=ABS(Y)，再進行頻率轉換，轉換方法：f=(0:length(Y)1)’*fs/length(Y) 分析解調(diào)器的抗噪性能時，在輸入端加入高斯白噪聲，采用 MATLAB函數(shù) awgn實現(xiàn)，其函數(shù)使用格式為： Y =AWGN(X,S_N)，加高斯白噪聲于 X中， S_N為信噪比，單位為 dB，其值在假設 X的功率為 0dBM的情況下確定。 信號的信噪比為信號中有用的信號功率與噪聲功率的比值，根據(jù)信號功率定義，采用MATLAB函數(shù) var實現(xiàn)，其函數(shù)常使用格式為： Y =VAR(X)，返回向量的方差，則信噪比為 :S_N=VAR(X1)/VAR(X2)。 繪制曲線采用 MATLAB函數(shù) plot實現(xiàn)，其函數(shù)常使用格式： PLOT（ X， Y） ， X為橫軸變量， Y為縱軸變量，坐標范圍限定 AXIS([x1 x2 y1 y2])，軸線說明 XLABEL(‘ ‘)和 YLABEL(‘ ‘)。 用頻率 300HZ 正弦波調(diào)制頻率 30KHZ 的正弦波，采用同步解調(diào)，觀察調(diào)制信號、已武漢理工大學《 MATLAB 課程設計》報告 8 調(diào)信號、解調(diào)信號的波形、頻譜以及解調(diào)器輸入輸出信噪比的關系。 MATLAB 源程序如下： fc=30000。 fs=100000。 N=1000。 n=0:N1。 t=n/fs。 x=A*sin(2*pi*300*t)。 y=modulate(x,fc,fs,39。am39。)。 %抑制雙邊帶振幅調(diào)制 fft1=fft(x,N)。 %傅里葉變換 mag1=abs(fft1)。 %取模 f1=(0:length(fft1)1)39。*fs/length(fft1)。 fft2=fft(y,N)。 mag2=abs(fft2)。 f2=(0:length(fft2)1)39。*fs/length(fft2)。 figure(1)。 subplot(2,2,1)。 plot(t,x)。 xlabel(39。調(diào)制信號波形 39。) subplot(2,2,2)。 plot(f1,mag1)。 axis([0 600 0 600])。 xlabel(39。調(diào)制信號頻譜 39。) subplot(2,2,3)。 plot(t,y)。 xlabel(39。已調(diào)信號波形 39。) subplot(2,2,4)。 plot(f2,mag2)。 axis([28000 32021 0 400])。 xlabel(39。已調(diào)信號頻譜 39。) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% yn=awgn(y,4)。 %加入高斯白噪聲 武漢理工大學《 MATLAB 課程設計》報告 9 znn=demod(y,fc,fs,39。am39。)。 %無噪聲已調(diào)信號解調(diào) zn=demod(yn,fc,fs,39。am39。)。 %加噪聲已調(diào)信號解調(diào) fft3=fft(znn,N)。 mag3=abs(fft3)。 f3=(0:length(fft3)1)39。*fs/length(fft3)。 figure(2)。 subplot(3,1,1)。 plot(t,zn)。 xlabel(39。加噪聲解調(diào)信號波形 39。) subplot(3,1,2)。 plot(t,znn)。 xlabel(39。無噪聲解調(diào)信號波形 39。) subplot(3,1,3)。 plot(f3,mag3)。 axis([0 500 0 500])。 xlabel(39。解調(diào)信號頻譜 39。) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% yn1=awgn(y,8)。 yn2=awgn(y,12)。 yn3=awgn(y,16)。 yn4=awgn(y,20)。 zn1=demod(yn1,fc,fs,39。am39。)。 zn2=demod(yn2,fc,fs,39。am39。)。 zn3=demod(yn3,fc,fs,39。am39。)。 zn4=demod(yn4,fc,fs,39。am39。)。 dyi=yny。 %高斯白噪聲 s_ni=var(y)/var(dyi)。 %輸入信噪比 dyo=znznn。 %解調(diào)后噪聲 s_no=var(znn)/var(dyo)。 %輸出信噪比 dyi1=yn1