白噪声及有色噪声序列的产生

%白噪声及有色噪声序列的产生
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