clear?all;
clc
load（‘chirp‘‘Fs‘‘y‘）;
p=audioplayer（yFs）;
play（p）;
subplot（211）;
plot（y）;
title（‘原始信號‘）
Py=fft（y）;
subplot（212）;
plot（abs（Py））;
title（‘原始信號頻譜‘）?????%原始信號時域圖、頻域圖及信號的播放
disp（‘按任意鍵繼續‘）;
pause（）;??????????????????%按任意鍵程序繼續執行
N=length（y）;
%?sigma=0.1;?			??%控制噪聲強度
%?z=zeros（N1）;
%?v=randn（N1）*sigma;?????%產生噪聲
%?z=y+v;
z=awgn（y30）;?????????????%將高斯白色噪聲添加到信號中.產生高斯白噪聲加噪信號z
??????????????????????????%標量30指定了每一個采樣點信號與噪聲的比率，單位為dB
q=audioplayer（zFs）;
play（q）;
figure（）;?????????????????%創建新圖像窗口
subplot（211）;
plot（z）;
title（‘觀測信號‘）
Pz=fftn（z）;
subplot（212）;
plot（abs（Pz））;
title（‘觀測信號的頻譜‘）?????%觀測信號時域圖、頻域圖及信號的播放
%加噪后均方誤差
err1=mean（（z-y）.^2）;???????????%對（z-y）列計算（a^2+b^2+...n^2）/n----
disp（[‘加噪后均方誤差為‘num2str（err1）]）;??????%把數值轉換成字符串
disp（‘按任

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????25276??2004-03-14?10:31??1\chirp.mat

?????文件???????1688??2018-12-04?19:02??1\winner.m

?????目錄??????????0??2018-12-28?15:34??1

-----------?---------??----------?-----??----

????????????????26964????????????????????3