%白噪聲及有色噪聲序列的產生
clear?all;?close?all;

L=500;?%仿真長度
d=[1?-1.5?0.7?0.1];?c=[1?0.5?0.2];?%D、C多項式的系數（可用roots命令求其根）
nd=length（d）-1;?nc=length（c）-1;?%nd、nc為D、C的階次
xik=zeros（nc1）;?%白噪聲初值，相當于ξ（k-1）...ξ（k-nc）
ek=zeros（nd1）;?%有色噪聲初值
xi=randn（L1）;?%randn產生均值為0，方差為1的高斯隨機序列（白噪聲序列）

for?k=1:L
????e（k）=-d（2:nd+1）*ek+c*[xi（k）;xi