characterfunction.m

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

function characterfunction
%程序C7_1  插值法
n=4;                     %%%%%%% 输入矩阵大小
A=[4,3,2,1;3,3,2,1;2,2,2,1;1,1,1,1];   %%%%%%% 输入矩阵
h=0.5;                     %%%%%%% 节点步长
for k=1:n+1
    t(k)=(k-1)*h;
    f_{1}(k)=det(A-t(k)*eye(n));    %%% 节点数值
end 
for k=2:n+1
    for i=1:(n+2-k)
        f_{k}(i)=f_{k-1}(i+1)-f_{k-1}(i);   %%%% 向前差分
    end
end
M=markov(n);
T=zeros(1,n+1);  %%%%%  T即输出多项式的系数
T(1)=f_{1}(1);
for i=2:n+1
    T(i)=f_{i}(1)*M(i-1,i-1);
      for k=i:n   %内循环
       T(i)=T(i)+f_{k+1}(1)*M(k,i-1);
      end
end
 for k=2:n+1
     T(k)=T(k)/(h^(k-1));     %%%%换回原来变量 s=(t-t(0))/h  set:t(0)=0
 end
 T
 eig(M)