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

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????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