请教马尔可夫的matlab编程

yr13236

设一阶马尔可夫过程n‘（t）＝-a*n（t)+w（t），a为一阶马尔可夫常数的倒数，w为高斯白噪声，强度为c。离散化后方程为n(k+1)＝exp(-aT)*n(k)＋w(k)，T为采样周期，那么是否一阶马尔可夫过程的程序可以编为：
n(1)=0;
T=1;
w=c*randn(N,1);
for k=1:N-1
    n(k+1)＝exp(-aT)*n(k)＋w(k)
end
plot(k,n);
请各位高手看看对不对！谢谢！

[ 本帖最后由 yr13236 于 2006-10-26 17:55 编辑 ]

lidavc

我修改的是这样，你看看！

1. a=0.05;c=0.2;T=1;
   
2. n(1)=1;
   
3. w=c*randn(1,1000);
   
4. for k=1:1000
   
5.     n(k+1)=exp(-a*T)*n(k)+w(k);
   
6. end
   
7. plot(1:k,n(1:1000))

复制代码