基于matlab的dsb調(diào)制與解調(diào)分析(編輯修改稿)
2024-09-06 17:10 本頁面
【文章內(nèi)容簡介】 聲于X中，SNR為信噪比，單位為dB，其值在假設(shè)X的功率為0dBM的情況下確定。信號的信噪比為信號中有用的信號功率與噪聲功率的比值，根據(jù)信號功率定義，采用MATLAB函數(shù)var實現(xiàn)，其函數(shù)常使用格式為：Y =VAR(X)，返回向量的方差，則信噪比為:SNR=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)制信號、已調(diào)信號、解調(diào)信號的波形、頻譜以及解調(diào)器輸入輸出信噪比的關(guān)系。MATLAB源程序如下：Fs=100000。 %抽樣頻率 Fc=30000。 %載波頻率N=1000。 %FFT長度n=0:N1。t=n/Fs。 %截止時間和步長x=a*sin(2*pi*300*t)。 %基帶調(diào)制信號 y=modulate(x,Fc,Fs,39。am39。)。 %抑制雙邊帶振幅調(diào)制yn=awgn(y,4)。 %加入高斯白噪聲yn1=awgn(y,10)。yn2=awgn(y,15)。yn3=awgn(y,20)。yn4=awgn(y,25)。y1=demod(y,Fc,Fs,39。am39。)。 %無噪聲已調(diào)信號解調(diào)yyn=demod(yn,30000,Fs,39。am39。)。 %加噪聲已調(diào)信號解調(diào)yyn1=demod(yn1,30000,Fs,39。am39。)。yyn2=demod(yn2,30000,Fs,39。am39。)。yyn3=demod(yn3,30000,Fs,39。am39。)。yyn4=demod(yn4,30000,Fs,39。am39。)。dy1=yny。 %高斯白噪聲snr1=var(y)/var(dy1)。 %輸入信噪比dy2=yyny1。 %解調(diào)后噪聲snr2=var(y1)/var(dy2)。 %輸出信噪比 dy11=yn1y。 snr11=var(y)/var(dy11)。 dy21=yyn1y1。 snr21=var(y1)/var(dy21)。 dy12=yn2y。 snr12=var(y)/var(dy12)。 dy22=yyn2y1。 snr22=var(y1)/var(dy22)。 dy13=yn3y。 snr13=var(y)/var(dy13)。 dy23=yyn3y1。 snr23=var(y1)/var(dy23)。 dy14=yn4y。 snr14=var(y)/var(dy14)。 dy24=yyn4y1。 snr24=var(y1)/var(dy24)。 in=[snr1,snr11,snr12,snr13,snr14]。 out=[snr2,snr21,snr22,snr23,snr24]。ff1=fft(x,N)。 %傅里葉變換mag1=abs(ff1)。 %取模f1=(0:length(ff1)1)39。*Fs/length(ff1)。 %頻率轉(zhuǎn)換ff2=fft(y,N)。mag2=abs(ff2)。f2=(0:length(ff2)1)39。*Fs/length(ff2)。ff3=fft(y1,N)。mag3=abs(ff3)。f3=(0:length(ff3)1)39。*Fs/length(ff3)。figure(1)。subplot(221) %繪制曲線plot(t,x)xlabel(39。調(diào)制信號波形39。)subplot(222)plot(f1,mag1)axis([0 1000 0 1000])xlabel(39。調(diào)制信號頻譜39。)subplot(223)plot(t,y)xlabel(39。已調(diào)信號波形39。)subplot(224)plot(f2,mag2)axis([0 40000 0 500])xlabel(39。已調(diào)信號頻譜39。)figure(2)。subplot(311)plot
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1