markov.m

来自「利用MATLAB进行程序C7_3 QR方法的算法」· M 代码 · 共 31 行

function M=markov(n)
%  算马尔可夫系数c_{n,i}
jc(1)=1;
for k=2:n
    jc(k)=jc(k-1)*k;  %算阶乘
end
a_{1}=[1,0];
b_{1}=a_{1};
for j=2:n
    a_{j}=[1,-j+1];
end
MM=zeros(n,n);
for k=2:n
    b_{k}=conv(b_{k-1},a_{k});
end
for ii=1:n
    for jj=1:ii
      MM(ii,jj)=b_{ii}(jj)/jc(ii);
    end
end
for i=1:n
    for j=1:i
      M(i,j)=MM(i,i+1-j);
    end
end