MATLAB源代码

埃特金插值

function f=Atken(x,y,x0) syms t;

if(length(x)==length(y)) n=length(x); else

disp; return; end

y1(1:n)=t; for(i=1:n-1) for(j=i+1:n)

y1(j)=y(j)*(t-x(i))/(x(j)-x(i))+y(i)*(t-x(j))/(x(i)-x(j)); end y=y1

simplify(y1); end

if(nargin==3)

f=subs(y1(n),'t',x0); else

simplify(y1(n)); f=collect(y1(n)); f=vpa(f,6); end

特征多项式法

function l=Chapoly(A)

%特征多项式法求矩阵特征值 %已知矩阵:A

%求得的矩阵特征值:l syms t; N=size(A); n=N(1,1);

y=det(A-t*eye (n,n)); l=solve(y);

l=vpa(1,5); % function f = Chebyshev(y,k,x0) syms t;

T(1:k+1) = t; T(1) = 1; T(2) = t;

Word文档免费下载Word文档免费下载:MATLAB源代码 (共29页,当前第1页)

MATLAB源代码相关文档

最新文档

返回顶部