数值计算方法上机答案.pdf

《数值计算方法上机答案.pdf》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、本科实验报告课程名称：实验项目：实验地点：机房专业班级：采矿1206学号：2012002896学生姓名：康明月指导教师：2014年7月3日一、求非线性方程的根。1、求方程f(x)xcosx0在x1.5附近的是根，要求精度满足03xx10.(牛顿切线法)k1kf=inline('x-cos(x)');%f(x)df=inline('1+sin(x)');%f'(x)n=1;x0=input('x0=');del=input('del=');N=input('N=');fprintf('kx(k)');fprintf('%2d%f',

2、0,x0);F0=f(x0);dF0=df(x0);whilen

3、abs(F1)

4、nx0=1.5del=1e-4N=100kx(k)01.50000010.78447220.739519结果：0.739085322、求方程f(x)xx0.80在x1附近的是根，求出具有思维有效数字的根0近似值..(简单迭代法)clearclcphi=inline('(0.8+x^2)^(1/3)');%迭代函数x0=input('x0=');del=input('del=');N=input('N=');n=1;fprintf('%2d%f',0,x0);whilen

5、近似解=%f',x);returnendfprintf('%2d%f',n,x);n=n+1;x0=x;endfprintf('%fd次迭代后未达到精度要求.',N)x0=1del=1e-4N=10001.00000011.21644021.31611631.36300441.38518051.39568861.40067171.40303481.40415591.404687101.404939111.405059近似解=1.405116100.000000d次迭代后未达到精度要求.二、求解线性方程组（直接法或迭代法）22

6、13x1812113x111、2x8138x12310443x482(列主元素消元法)a=input('a=')%[2,2,1,-3,8;-2,1,-1,-3,1;8,-1,3,8,-1;10,4,4,3,8];[p,n]=size(a);forw=1:p[x,y]=find(a(w:p,w)==max(max(a(w:p,w))));q=a(w,:);a(w,:)=a(x,:);a(x,:)=q;endforj=1:(p-1)fori=(j

7、+1):pa(i,:)=a(j,j)/a(i,j).*a(i,:)-a(j,:);endendm=p;whilem>1s(m)=a(m,n);j=p;while((j>2)&(j>=m+1)&(j

8、071.0060134.21112.16.31.700.5130.511.5A4.5013.113.90761.70512.178.9990.1010.8120.0170.914.9180.11234.5521.8030.23122.3222543b精确解x,E1E6(误差限)。240.236429.3045117.8186function[x,det,flag]=Gauss(A,b)[n,m]

9、=size(A);nb=length(b);ifn~=merror('thero