單邊帶信號調制與解調matlab-預覽頁
【正文】 表示低通濾波器的截止頻率，混頻信號y中高于此頻率的分量將被濾除，得到解調信號；而“DEM=H.*Y。另外要注意對采樣頻率fs的選取。其他一些參數(shù)設置為：信號保持時間t0=，濾波器截止頻率fl=700Hz，頻率分辨率df=。因為輸入計算機中的連續(xù)信號被轉換成了離散的序列，而設其序列共有N個點，則兩根譜線間的距離＝fs/N，它也是頻率分辨率的一種度量。我在MATLAB中實現(xiàn)SSB解調的仿真就是基于以上的程序和參數(shù)設定，運行調用了plot函數(shù)的M文件，可以分別得出解調前后的頻域與時域對比圖，為了便于參照，將LSSB信號的時域與頻域圖也一同進行比較。三者的頻域圖比較如圖c所示：圖c調制信號、LSSB信號、解調信號的頻域對比圖由頻域圖可以明顯地看出，下邊帶信號將基帶信號的頻譜向兩側進行了搬移，然后去除了其中上邊帶的那一半。為了探究低通濾波器的設置對解調的影響，我將濾波器的截止頻率fl分別改為1000Hz和50Hz進行比較，根據(jù)結果來判斷原因。接下來將截止頻率fl設為50Hz，觀察在fl較低的情況下的解調輸出結果。由此可以估計出一種可能性：基帶信號的一部分丟失了，解調信號只是原基帶信號的一部分。這就是問題所在。因此為了保全基帶信號的完整性，fl不宜過低。所以開始所設的fl=700Hz是一個較為合適的值，解調結果也很接近基帶信號。在本程序中，實現(xiàn)調制過程使用語句“u=ssbmod(m,fc,fs)。解調后的信號基本與基帶信號相同，只是邊緣處有少許失真。這種方法不用自己設置濾波器，程序簡單，在實現(xiàn)解調時較為簡便。而且大部分噪聲為隨機的高斯白噪聲，所以在設計時引入噪聲，才能夠真正模擬實際中信號傳輸所遇到的問題，進而思考怎樣才能在接受端更好地恢復基帶信號。measured39。的功能是在疊加噪聲前先測量信號x的功率，再根據(jù)所測結果來產(chǎn)生所設信噪比的噪聲。measured39。)。 [S2,s2,dfl]=fft_seq(s2,ts,df)。再進一步觀察圖1中的時域變化：當疊加大噪聲時，解調波形出現(xiàn)明顯的紊亂，產(chǎn)生了很多不規(guī)則的毛刺雜波，與基帶信號相去甚遠；而疊加小噪聲時，波形圖只是局部稍有不平滑現(xiàn)象，總體還是與未加噪聲時的解調信號相同，出入不大。所以當信噪比低于一定大小時，會給解調信號帶來嚴重的失真，導致接收端無法正確地接收有用信號。這樣不利于我觀察頻譜特性，為了能讓頻譜圖居中顯示，我對程序設置進行分析，發(fā)現(xiàn)關鍵在于頻率矢量的設置上。開始我將截止頻率設置為fl=50Hz，但仿真結果顯示解調信號的幅度明顯小于基帶信號。對于其合適的參數(shù)我進行了多次測試，力求在可以濾去2高頻分量的情況下使解調信號損失最小。用數(shù)學表達式可以表示為 哪些語句是生成噪聲的？我用函數(shù)“awgn(x,snr,39。measured39。s1 = awgn(lssb,snr1,39。measured39。這些具體的問題在書上只有少量的理論介紹，要在實際運用中結合這些理論是一件不太容易的事。通過這次的實踐，我明白了要將理論與實際相結合的道理，盡管這個過程會有一些辛苦，但通過努力實現(xiàn)后，就能大大深化我對知識的理解程度，增長實踐經(jīng)驗。 %信號保持時間ts=1/2500。 %濾波器的截止頻率df=。 %時間矢量m=cos(2*pi*50.*t)。 %下邊帶信號s1 = awgn(lssb,snr1,39。measured39。 %加大噪聲后混頻y2=s2.*c。 %傅里葉變換求下邊帶信號頻譜LSSB=LSSB/fs。 %求下邊帶信號加小噪聲的頻譜S2=S2/fs。 %傅里葉變換得出混頻頻譜Y1=Y1/fs。 %設計濾波器f=[0:dfl:dfl*(length(m)1)]fs/2。DEM=H.*Y。 %濾波器輸出(即解調后)信號dem1=real(ifft(DEM1))*fs。plot(t,m(1:length(t)))。subplot(3,1,2)。)。解調信號39。plot(f,abs(fftshift(M)))。subplot(3,1,2)。)。解調信號頻譜39。plot(t,lssb(1:length(t)))。subplot(3,2,2)。)。加大噪聲后下邊帶信號39。 %顯示加大噪聲后下邊帶信號頻譜title(39。plot(t,s2(1:length(t)))。subplot(3,2,6)。)。 %顯示解調信號title(39。plot(t,dem1(1:length(t)))。subplot(3,1,3)。)。 %顯示解調信號頻譜title(39。plot(f,abs(fftshift(DEM1)))。subplot(3,1,3)。)。endn2=length(m)。 %對m補零df=fs/n。 %此為載波頻率fs=1/ts。 %原始信號u=ssbmod(m,fc,fs)。 %縮放[U,u,dfl]=fft_seq(u,ts,df)。 %縮放f=[0:df:df*(length(m)1)]fs/2。 %顯示調制信號title(39。plot(t,u(1:length(t)))。subplot(3,1,3)。)。 %顯示調制信號頻譜title(39。plot(f,abs(fftshift(U)))。subplot(3,1,3)。
公司管理相關推薦
鄂ICP備17016276號-1