%白噪声及有色噪声序列的产生clearall;closeall;L=500;%仿真长度d=[1-1.50.70.1];c=[10.50.2];%D、C多项式的系数(可用roots命令求其根)nd=length(d)-1;nc=length(c)-1;%nd、nc为D、C的阶次xik=zeros(nc,1);%白噪声初值，相当于ξ(k-1)...ξ(k-nc)ek=zeros(nd,1);%有色噪声初值xi=randn(L,1);%randn产生均值为0，方差为1的高斯随机序列（白噪声序列）fork=1:Le(k)=-d(2:nd+1)*ek+c*[xi(k);xik];%产生有色噪声