matlab数值分析代码,数值分析matlab代码

《数值分析matlab代码》由会员分享，可在线阅读，更多相关《数值分析matlab代码(24页珍藏版)》请在人人文库网上搜索。

1、1、%用牛顿法求f(x)=x-sin x 的零点，e=10(-6)disp(牛顿法);i=1;n0=180;p0=pi/3;tol=10(-6);for i=1:n0p=p0-(p0-sin(p0)/(1-cos(p0);if abs(p-p0)0a=p;else b=p;endendif i=n0&(fp=0|(abs(b-a)/2)0disp(|p-p1|=)disp(abs(pp-p1)disp(用不动点迭代法求得方程的正根为)disp(pp);elsedisp(用不动点迭代法无法求出方程的正根);endbreak;elsep1=pp;endendif i=n0disp(n0)disp(。

2、次不动点迭代后无法求出方程的正根)end2、%用牛顿法求方程 x-ex+4=0的正根与负根，误差限是10-6 disp(牛顿法)n0=80;p0=1;for i=1:n0p=p0-(p0-exp(p0)+4)/(1-exp(p0);if abs(p-p0)bp=j; endendif a(p,i)=0disp(该方程组没有唯一解); endif i=pb=a(i,:);a(i,:)=a(p,:);a(p,:)=b; endfor j=(i+1):nm=a(j,i)/a(i,i);a(j,:)=a(j,:)-m*a(i,:);endendif a(n,n)=0disp(该方程组没有唯一解); e。

3、ndx(n)=a(n,n+1)/a(n,n);for i=(n-1):(-1):1sum=0;for j=(i+1):nsum=sum+a(i,j)*x(j); endx(i)=(a(i,n+1)-sum)/a(i,i);enddisp(方程组的解 X=);disp(x);3、disp(按比例选主元高斯消元法);a=1.19,2.11,-100,1,1.12;14.2,-0.122,12.2,-1,3.44;0,100,-99.9,1,2.15;15.3,0.11,-13.1,-1,4.16;for i=1:ns(i)=abs(a(i,1);for j=1:nif abs(a(i,j)s(i)。

4、s(i)=abs(a(i,j);endendif s(i)=0disp(该方程组没有唯一解);break;endendfor i=1:(n-1)b=abs(a(i,i)/s(i);for j=i:nif (abs(a(j,i)/s(j)bp=j; endendif a(p,i)=0disp(该方程组没有唯一解); endif i=pb=a(i,:);a(i,:)=a(p,:);a(p,:)=b; endfor j=(i+1):nm=a(j,i)/a(i,i);a(j,:)=a(j,:)-m*a(i,:);endendif a(n,n)=0disp(该方程组没有唯一解); endx(n)=a(n。

5、,n+1)/a(n,n);for i=(n-1):(-1):1sum=0;for j=(i+1):nsum=sum+a(i,j)*x(j); endx(i)=(a(i,n+1)-sum)/a(i,i);enddisp(方程组的解 X=);disp(x);题目：采用Jacobi方法和G-S方法求解线性方程组Ax=b，满足按范数误差在之内，A的元素为，并且b的元素为。%Jacobi 方法求方程组的解disp(Jacobi 方法)for i=1:80j=i;a(i,j)=2*i;endfor i=1:78j=i+2;a(i,j)=0.5*i;endfor i=3:80j=i-2;a(i,j)=0.5。

6、*i; endfor i=1:76j=i+4;a(i,j)=0.25*i;endfor i=5:80j=i-4;a(i,j)=0.25*i;endfor i=1:80b(i)=pi; endN=180;n=80;k=1;tol=10(-5);for i=1:80x0(i)=1; endwhile knormXnormX=abs(x(l);endend%求X0的范数normX0=0;for l=1:nif abs(x0(l)normX0normX0=abs(x0(l);endendif(abs(normX-normX0)Ndisp(k);disp(次迭代后无法求出方程组的解); end%Gauss-Seidel 方法求方程组的解disp(Gauss-Seidel 方法)N=80;n=80;k=1;tol=10(-5);for i=1:80x0(i)=1; endk=1;while knormXnormX=abs(x(l);endend%求X0的范数normX0=0;for l=1:nif abs(x0(l)normX0normX0=abs(x0(l);endendif(abs(normX-normX0)Ndisp(k);disp(次迭代后无法求出方程组的解); end。