大连理矩阵上机作业.docxVIP

第一题Lagrange插值函数function y=lagrange(x0,y0,x);n=length(x0);m=length(x);for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s;endx0=[1:10];y0=[67.052,68.008,69.803,72.024,73.400,72.063,74.669,74.487,74.065,76.777];lagrange(x0,y0,17)ans=-1.9516e+12x=[1:0.1:10];x=x;plot(x0,y0,r);hold onplot(x,y,k);legend(原函数,拟合函数)拟合图像如下拟合函数出现了龙格现象，运用多项式进行插值拟合时，效果并不好，高次多项式会因为误差的不断累积，导致龙格现象的发生。第二题function fun =nihe(n)m=[67.052*10^6,68.008*10^6,69.803*10^6,72.024*10^6,73.400*10^6,72.063*10^6,74.669*10^6,74.487*10^6,74.065*10^6,76.777*10^6];w=[1,2,3,4,5,6,7,8,9,10];d1=0;d2=0;d3=0;y1=polyfit(m,w,1);y2=polyfit(m,w,2);y3=polyfit(m,w,3);y2=poly2sym(s2);y3=poly2sym(s3);y4=poly2sym(s4);f1=subs(y1,17);f2=subs(y2,17); f3=subs(y3,17);for p=1:10; d1=d1+(subs(y1,w(p))-m(p))^2; d2=d2+(subs(y2,w(p))-m(p))^2; d3=d3+(subs(y3,w(p))-m(p))^2;endd1=sqrt(d1);d2=sqrt(d2); d3=sqrt(d3);fun=[f1 f2 f3;d2 d3 d4];return;结果三次函数的均方误差最小，拟合的最好。函数function f=fun(x)syms aa=x;f=a*a*a+a*a+a-3;梯度函数function df=dfun(x)df=3*x*x+2*x+1;Newton法function result=didainewton(x0)k=0;xk=x0;xi=1;e0=abs(x0-xi);ek=e0;m=zeros(7,1);n=zeros(7,1);p=zeros(7,1);result=zeros(7,3);while k7 ak=feval(fun,xk); bk=feval(dfun,xk); xk=xk-ak/bk; e0=ek; k=k+1; m(k)=xk; ek=abs(m(k)-xi); jingdu=ek/(e0*e0); n(k)=ek; p(k)=jingdu;endresult=[m,n,p];return;计算结果iXiEi=|xi-x*|Ei/ Ei-120-0.71.7_12.62061.62060.560721.70840.70840.269831.20640.20640.411241.02420.02420.567351.00040.00040.653561.00000.00000.666571.00000.00000.6609Newton迭代法的收敛速度快，至少是平方收敛的，GAUSS消去法function x=DelGauss(N)%Gauss??è￥·¨syms M;M=N;a=zeros(M);b=ones(M,1);for i=1:Mfor j=1:M a(i,j)=1/(i+j-1);endend[n,m]=size(a);nb=length(b); det=1; x=zeros(n,1);for k=1:n-1 for i=k+1:n if a(k,k)==0returnend m=a(i,k)/a(k,k);for j=k+1:n a(i,j)=a(i,j)-m*a(k,j);end b(i)=b(i)-m*b(k);enddet=det*a(k,k); end det=det*a(n,n);for k=n:-1:1for j=k+