MATLAB環境下的正弦信號及高斯白噪聲仿真程序說明_第1頁
MATLAB環境下的正弦信號及高斯白噪聲仿真程序說明_第2頁
MATLAB環境下的正弦信號及高斯白噪聲仿真程序說明_第3頁
全文預覽已結束

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

1、姓名：朱奇峰 專業：電子與通信工程 方向：數字廣播電視技術 學號：103320430109033MATLAB環境下的正弦信號及高斯白噪聲仿真程序說明一、信號的產生及時域觀察1、設定正選信號的頻率為10HZ，抽樣頻率為100HZ；2、設定N(0,0.25)高斯白噪聲，及噪聲功率為0.25W；3、最后將噪聲疊加到正弦信號上，觀察其三者時域波形。二、信號頻譜及白噪聲功率譜的求解與觀察1、對原正弦信號直接進行FFT，得出其頻譜；2、求白噪聲的自相關函數，隨機序列自相關函數的無偏估計公式為： 對所求自相關函數進行FFT變換，求的白噪聲的功率譜函數。三、仿真結果：附源程序代碼：fs=100;fc=10;x

2、=(0:1/fs:2);n=201;y1=sin(2*pi*fc*x); %原正弦信號，頻率為10a=0;b=0.5; %均值為a，方差為b2subplot(3,2,1);plot(x,y1,'r');title('y=sin(20pi*x)');ylabel('y');xlabel('x/20pi');grid;y2=a+b*randn(1,n); %高斯白噪聲subplot(3,2,2);plot(x,y2,'r');title('N(0，0.25)的高斯白噪聲');ylabel('y&

3、#39;);xlabel('x/20pi');grid;y=y1+y2; %加入噪聲之后的信號subplot(3,2,3);plot(x,y,'r');title('疊加了高斯白噪聲的sinx');ylabel('y');xlabel('x/20pi');grid;FY=fft(y); %傅里葉變換得出頻譜函數FY1=fftshift(FY); %頻譜校正f=(0:200)*fs/n-fs/2;subplot(3,2,4);plot(f,abs(FY1),'r');title('函數頻譜圖&

4、#39;);ylabel('F(jw)');xlabel('w');grid;%求高斯白噪聲的自相關函數m=50;i=-0.49:1/fs:0.49;for j=1:m R(j)=sum(y2(1:n-j-1).*y2(j:199),2)/(n-j);%無偏自相關函數的估計 Rx(49+j)=R(j); Rx(51-j)=R(j);endsubplot(3,2,5);plot(i,Rx,'r');%plot(leg*1/fs,Rx,'r');title('白噪聲自相關函數圖');ylabel('Rx');xlabel('x');grid;Fy2=fft(Rx); %傅里葉變換得出白噪聲功率譜函數Fy21=fftshift(Fy2); %功率譜校正f=(0:98)*fs/99-fs/2;subplot(3,2,6);plot(f,abs(Fy21),'r');axis(-50 50 -0.5 1);titl

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間，僅對用戶上傳內容的表現方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容，請與我們聯系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論