华工数值分析课程设计.docxVIP

数值分析课程设计 第一题： 设计思路：我打算用选主元法，先算出每一列，然后把买一列加起来就是结果了。 2.程序清单：function x=mat(a,b,flag) if nargin3,flag=0;end n=length(b);a=[a,b]; for i=1:(n-1) [ar,r]=max(abs(a(i:n,i)));r=r+i-1; if ri t=a(i,:);a(i,:)=a(r,:);a(r,:)=t; end a((i+1):n,(i+1):(n+1))=a((i+1):n,(i+1):(n+1))-a((i+1):n,i)/a(i,i)*a(i,(i+1):(n+1)); a((i+1):n,i)=zeros(n-i,1); if flag==0,a,end end x=zeros(n,1);x(n)=a(n,n+1)/a(n,n); for i=n-1:-1:1 x(i)=(a(i,n+1)-a(i,(i+1):n)*x((i+1):n))/a(i,i); end 3.流程图：我的例子是九阶对角矩阵对角元素为1,2,3,4,5,6,7,8,9,。流程如图所示 通过对比知道，所编的程序算的结果和inv的结果一致。 我的方法时间约为0.003*9=0.027秒，而inv方法时间约为0.0031秒，明显inv方法快。 第二题： 设计思路：我通过limit与求阶方法编程。 程序清单： function n=qiujie(x0) k=0; syms x; x1=0.99*x-x^2; h=(abs(x0-x1))/(abs(x0-x))^k; while limit(h,x,x0)==0 k=k+1; h=(abs(x0-x1))/(abs(x0-x))^k; end n=k function [y,m]=bdd(x,eps) if nargin==1 eps=1e-100; elseif nargin1 error return end x1=0.99*x-x^2; n=1; while (abs(x1-x)=1e-100)(n1e100) x=x1; x1=0.99*x-x^2; n=n+1; end y=x1 m=n 流程图：由所得知，所求阶数为一阶。 第三题： 设计思路：按照题意用三次样条。从题意分析可知x=4.9与22.9的一阶导数均为0所以带入七组数据即可。 程序清单： function s=sp(x,y,dx0,dxn) n=length(x)-1;h=diff(x);d=diff(y)./h;a=h(2:n-1); b=2*(h(1:n-1)+h(2:n));c=h(2:n);u=6*diff(d); b(1)=b(1)-h(1)/2; u(1)=u(1)-3*(d(1)-dx0); b(n-1)=b(n-1)-h(n)/2; u(n-1)=u(n-1)-3*(dxn-d(n)); for k=2:n-1 temp=a(k-1)/b(k-1); b(k)=b(k)-temp*c(k-1); u(k)=u(k)-temp*u(k-1); end m(n)=u(n-1)/b(n-1); for k=n-2:-1:1 m(k+1)=(u(k)-c(k)*m(k+2))/b(k); end m(1)=3*(d(1)-dx0)/h(1)-m(2)/2; m(n+1)=3*(dxn-d(n))/h(n)-m(n)/2; for k=0:n-1 s(k+1,1)=(m(k+2)-m(k+1))/(6*h(k+1)); s(k+1,2)=m(k+1)/2; s(k+1,3)=d(k+1)-h(k+1)*(2*m(k+1)+m(k+2))/6; s(k+1,4)=y(k+1); end 流程图： 可知所求函数为 T=-0.008x^3-0.0911x^2+22.8, [4.9,9.1] T=0.0291x^3-0.1919x^2-1,1888x+20.6, [9.1,13.7] T=-0.0157x^3+0.2093x^2-1,1091x+13.9, [13.7,18.3] T=0.0039x^3-0.0071x^2-1,1793x+11.7, [13.7,18.3] 第五题： 设计思路：我设计了精确值与题意的要求，并根据结果进行比较，只对有相同的数值即对1,1.1，。。。，2的值进行对比。 程序