基于matlab的am調(diào)制解調(diào)及抗噪性分析(編輯修改稿)
2024-12-13 21:56 本頁面
【文章內(nèi)容簡介】 y2=y0.*cos(2*pi*fc*n/fs)。 y3=fft(y2,N)。% fft 變換 q1=(0:N/21)*fs/N。 mx1=abs(y3(1:N/2))。 figure(2) subplot(2,1,1)。 plot(t,y2)。 title(39。已調(diào)信號的時(shí)時(shí)域波 39。)。 subplot(2,1,2)。 plot(q1,mx1)。 title(39。f1 已調(diào)信號的頻譜 39。)。 %繪圖 yc=cos(2*pi*fc*t)。 figure(3) subplot(2,1,1),plot(t,yc),title(39。載波 fc 時(shí)域波形 39。) N=512。 n=0:N1。 yc1=Am*cos(2*pi*fc*n/fs)。 y3=fft(yc1,N)。 q=(0:N/21)*fs/N。 mx=abs(y3(1:N/2))。 figure(3) subplot(2,1,2),plot(q,mx),title(39。載波 fc 頻譜 39。) y4=*randn(1,length(t))。%用 RANDN 產(chǎn)生高斯分布序列 w=y4.^2。 %噪聲功率 figure(4) subplot(2,1,1)。 plot(t,y4)。 title(39。高斯白噪聲時(shí)域波形 39。) y5=fft(y4,N)。 q2=(0:N/21)*fs/N。 mx2=abs(y5(1:N/2))。 subplot(2,1,2),plot(q2,mx2),title(39。高斯白噪聲頻域波形 39。) y6=y2+y4。 figure(5) subplot(2,1,1),plot(t,y6),title(39。疊加后的調(diào)制信號時(shí)域波形 39。) q3=q1。 mx3=mx1+mx2。 subplot(2,1,2),plot(q3,mx3),title(39。疊加后的調(diào)制信號頻譜波形 39。) %調(diào)制 yv=y6.*yc。 %乘以載波進(jìn)行解調(diào) Ws=yv.^2。 p1=fcfm。 [k,Wn,beta,ftype]=kaiserord([p1 fc],[1 0],[ ],fs)。 %Fir 數(shù)字低通濾波 window=kaiser(k+1,beta)。 %使用 kaiser 窗函數(shù) b=fir1(k,Wn,ftype,window,39。noscale39。)。 %使用標(biāo)準(zhǔn)頻率響應(yīng)的加窗設(shè)計(jì)函數(shù) yt=filter(b,1,yv)。 yssdb=yt.*22。 figure(6) subplot(2,1,1),plot(
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1