计算方法上机题问题详解

《计算方法上机题问题详解》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、标准文档2.用下列方法求方程e^x+10x-2=0的近似根，要求误差不超过5*10的负4次方，并比较计算量（1）二分法（局部，大图不太看得清，故后面两小题都用局部截图）实用文案标准文档（2）迭代法（3）牛顿法顺序消元法#include#include实用文案标准文档#includeintmain(){intN=4,i,j,p,q,k;doublem;doublea[4][5];doublex1,x2,x3,x4;for(i=0;i

2、][j]);for(p=0;p

3、,x4);scanf("%lf",&a[i][j]);（这一步只是为了看到运行的结果）}运行结果实用文案标准文档列主元消元法function[x,det,flag]=Gauss(A,b)[n,m]=size(A);nb=length(b);flag='OK';det=1;x=zeros(n,1);fork=1:n-1max1=0;fori=k:nifabs(A(i,k))>max1max1=abs(A(i,k));r=i;endendifmax1<1e-10flag='failure';return;实用文案标准文档endifr>kforj=k:nz=A(k,j)

4、;A(k,j)=A(r,j);A(r,j)=z;endz=b(k);b(k)=b(r);b(r)=z;det=-det;endfori=k+1:nm=A(i,k)/A(k,k);forj=k+1:nA(i,j)=A(i,j)-m*A(k,j);endb(i)=b(i)-m*b(k);enddet=det*A(k,k);enddet=det*A(n,n)ifabs(A(n,n))<1e-10flag='failure';return;endx(n)=b(n)/A(n,n);实用文案标准文档fork=n-1:-1：1forj=k+1:nb(k)=b(k)-A(k,j)

5、*x(j);endx(k)=b(k)/A(k,k);end运行结果：雅可比迭代法实用文案标准文档functiony=jacobi(a,b,x0)D=diag(diag(a));U=-triu(a,1);L=-tril(a,-1);B=D(L+U);f=Db;y=B*x0+f;n=1;whilenorm(y-x0)>1e-4x0=y;y=B*x0+f;n=n+1;endyn实用文案标准文档高斯赛德尔迭代法functiony=seidel(a,b,x0)D=diag(diag(a));U=-triu(a,1);L=-tril(a,-1);G=(D-L)U;f=(

6、D-L)b;y=G*x0+f;n=1;whilenorm(y-x0)>10^(-4)x0=y;y=G*x0+f;n=n+1;endyn实用文案标准文档SOR迭代法functiony=sor(a,b,w,x0)D=diag(diag(a));U=-triu(a,1);L=-tril(a,-1);实用文案标准文档lw=(D-w*L)((1-w)*D+w*U);f=(D-w*L)b*w;y=lw*x0+f;n=1;whilenorm(y-x0)>10^(-4)x0=y;y=lw*x0+f;n=n+1;endyn实用文案标准文档实用文案标准文档1.分段线性插值：fu

7、nctiony=fdxx(x0,y0,x)p=length(y0);n=length(x0);m=length(x);实用文案标准文档fori=1:mz=x(i);forj=1:n-1ifz

8、结果0.3