MATLAB源代码.docVIP

埃特金插值 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; c(1:k+1) = 0.0; c(1)=int(subs(y,findsym(sym(y)),sym(t))*T(1)/sqrt(1-t^2),t,-1,1)/pi; c(2)=2*int(subs(y,findsym(sym(y)),sym(t))*T(2)/sqrt(1-t^2),t,-1,1)/pi; f = c(1)+c(2)*t; for i=3:k+1 T(i) = 2*t*T(i-1)-T(i-2); c(i) = 2*int(subs(y,findsym(sym(y)),sym(t))*T(i)/sqrt(1-t^2),t,-1,1)/2; f = f + c(i)*T(i); f = vpa(f,6); if(i==k+1) if(nargin == 3) f = subs(f,t,x0); else f = vpa(f,6); end end end 弦割法求解非线性方程 function [x k t]=ChordsecantToEquation(f,x0,x1,eps) %弦割法求解非线性方程 %[x k t]=ChordsecantToEquation(f,x0,x1,eps) %x:近似解 %k:迭代次数 %t:运算时间 %f:原函数，定义为内联函数 %x0,x1:初始值 %eps:误差限 % %应用举例： %f=inline(x^3+4*x^2-10); %x=ChordsecantToEquation(f,1,2,0.5e-6) %[x k]=ChordsecantToEquation(f,1,2,0.5e-6) %[x k t]=ChordsecantToEquation(f,1,2,0.5e-6) %函数的最后一个参数也可以不写，默认情况下，eps=0.5e-6 %[x k t]=ChordsecantToEquation(f,1,2) if nargin==3 eps=0.5e-6; end tic; k=0; while 1 x=x1-f(x1)*(x1-x0)./(f(x1)-f(x0)); k=k+1; if abs(x-x1) eps || k 30 break; end x0=x1; x1=x; end t=toc; if k = 30 disp(迭代次数太多。); x=0; t=0; end 复合公式 function [I,step] = CombineTraprl(f,a,b,eps) %a,b 积分上下限 %eps 精度 %I 积分结果 %step 积分的子区间数 if(nargin ==3) eps=1.0e-4; end n=1; h=(b-a)/2; I1=0; I2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h; while abs(I2-I1)eps n=n+1; h=(b-a)/n; I1=I2; I2=0; for i=0:n-1 x=a+h*i; x1=x+h; I2=I2+(h/2)*(subs(sym(f),findsym(sym(f)),x)+subs(sym(f),findsym(sym(f)),x1)); end end I=I2; step=n; function f = DCS(x,y