clc;
clear?all;
Wp=2*pi*10;Ws=2*pi*20;Rp=0.5;Rs=40;fs=1000;W=2*pi*fs;???%構造低通濾波器
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=wgn（110240）;???%產生高斯白噪聲
t=0:1023;
gss=filter（bay）;?%濾波產生高斯色噪聲
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‘）??%自相關函數
subplot（245）;plot（laggss4）;
title（‘高斯色噪聲自相關函數‘）;axis（

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

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