clc;
clear?all;
Wp=2*pi*10;Ws=2*pi*20;Rp=0.5;Rs=40;fs=1000;W=2*pi*fs;???%構(gòu)造低通濾波器
Wn=0.06;
[ba]=butter（4Wn‘low‘）;
[hf]=freqz（ba1000fs）;
figure（2）;plot（fabs（h））;?xlabel（‘f/Hz‘）;ylabel（‘|H（jf）|‘）;
axis（[0?100?0?1.2]）;grid?on;
title（‘低通濾波器幅頻響應(yīng)‘）;

y=wgn（110240）;???%產(chǎn)生高斯白噪聲
t=0:1023;
gss=filter（bay）;?%濾波產(chǎn)生高斯色噪聲
figure（3）;
subplot（241）;plot（tgss）;
title（‘高斯色噪聲‘）;axis（[0?1024?-5?5]）;

gss1=mean（gss）;%均值
subplot（242）;plot（tgss1）;
title（‘高斯色噪聲均值‘）;axis（[0?1024?-1?1]）;

gss2=var（gss）;%方差
subplot（243）;plot（tgss2）;
title（‘高斯色噪聲方差‘）;axis（[0?1024?-0.5?1.5]）;

gss3=gss2+gss1.*gss1;%均方值
subplot（244）;plot（tgss3）;
title（‘高斯色噪聲均方值‘）;axis（[0?1024?-0.5?1.5]）;

[gss4lag]=xcorr（gss‘unbiased‘）??%自相關(guān)函數(shù)
subplot（245）;plot（laggss4）;
title（‘高斯色噪聲自相關(guān)函數(shù)‘）;axis（

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

?????文件????????721??2012-10-15?11:45??gauss.m

?????文件????????932??2012-10-15?11:48??Whitenoise.m

?????文件???????1274??2012-10-15?11:36??Colornoise.m

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

?????????????????2927????????????????????3