资源描述：
《华工数值分析课程设计.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数值分析课程设计第一题：1.设计思路：我打算用选主元法，先算出每一列，然后把买一列加起来就是结果了。2.程序清单：functionx=mat(a,b,flag)ifnargin<3,flag=0;endn=length(b);a=[a,b];fori=1:(n-1)[ar,r]=max(abs(a(i:n,i)));r=r+i-1;ifr>it=a(i,:);a(i,:)=a(r,:);a(r,:)=t;enda((i+1):n,(i+1):(n+1))=a((i+1):n,(i+1):(n+1))-
2、a((i+1):n,i)/a(i,i)*a(i,(i+1):(n+1));a((i+1):n,i)=zeros(n-i,1);ifflag==0,a,endendx=zeros(n,1);x(n)=a(n,n+1)/a(n,n);fori=n-1:-1:1x(i)=(a(i,n+1)-a(i,(i+1):n)*x((i+1):n))/a(i,i);end3.流程图：我的例子是九阶对角矩阵对角元素为1,2,3,4,5,6,7,8,9,。流程如图所示通过对比知道，所编的程序算的结果和inv的结果一致。我的
3、方法时间约为0.003*9=0.027秒，而inv方法时间约为0.0031秒，明显inv方法快。我第二题：1.设计思路：我通过limit与求阶方法编程。2.程序清单：functionn=qiujie(x0)k=0;symsx;x1=0.99*x-x^2;h=(abs(x0-x1))/(abs(x0-x))^k;whilelimit(h,x,x0)==0k=k+1;h=(abs(x0-x1))/(abs(x0-x))^k;endn=kfunction[y,m]=bdd(x,eps)ifnargin==1
4、eps=1e-100;elseifnargin<1errorreturnendx1=0.99*x-x^2;n=1;while(abs(x1-x)>=1e-100)&&(n<1e100)x=x1;x1=0.99*x-x^2;n=n+1;endy=x1m=n1.流程图：由所得知，所求阶数为一阶。第三题：1.设计思路：按照题意用三次样条。从题意分析可知x=4.9与22.9的一阶导数均为0所以带入七组数据即可。2.程序清单：functions=sp(x,y,dx0,dxn)n=length(x)-1;h=di
5、ff(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));fork=2:n-1temp=a(k-1)/b(k-1);b(k)=b(k)-temp*c(k-1);u(k)=u(k)-temp*u(k-1);endm(n)=u(n-1)/b(n-
6、1);fork=n-2:-1:1m(k+1)=(u(k)-c(k)*m(k+2))/b(k);endm(1)=3*(d(1)-dx0)/h(1)-m(2)/2;m(n+1)=3*(dxn-d(n))/h(n)-m(n)/2;fork=0:n-1s(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);end1.流程图：可知所求函数为T=
7、-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.1，。。。，2的值进行对比。2.程序清单：function[x,
8、y]=eu(fname,xspan,y0,h)x=xspan(1):h:xspan(2);y(1)=y0;forn=1:length(x)-1y(n+1)=y(n)+h*feval(fname,x(n),y(n));endx=x';y=y';function[x,y]=eu2(fname,xspan,y0,h)x=xspan(1):h:xspan(2);y(1)=y0;forn=1:length(x)-1k1=feval(fname,x(n),y(n))
