请教马尔可夫的matlab编程
设一阶马尔可夫过程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 编辑 ] 我修改的是这样，你看看！
a=0.05;c=0.2;T=1;
n(1)=1;
w=c*randn(1,1000);
for k=1:1000
n(k+1)=exp(-a*T)*n(k)+w(k);
end
plot(1:k,n(1:1000))
页:
[1]