




已阅读5页,还剩8页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
数学实验报告实验序号: 03 日期: 2013 年 12 月 5 日班级应数二班姓名黄怡丹学号1101114037实验名称求代数方程的近似根问题背景描述:求代数方程的根是最常见的数学问题之一,当是一次多项式时,称为线性方程,否则称之为非线性方程当是非线性方程时,由于的多样性,尚无一般的解析解法可使用,但如果对任意的精度要求,能求出方程的近似根,则可以认为求根的计算问题已经解决,至少能满足实际要求本实验介绍一些求方程实根的近似值的有效方法,要求在使用这些方法前先确定求根区间a,b,或给出某根的近似值实验目的:1. 了解对分法、迭代法、牛顿切线法求方程近似根的基本过程. 2.求代数方程(组)的.。实验原理与数学模型:1.对分法基本思想:将有根区间进行对分,判断出解在某个分段内,然后再对该段对分,依次类推,直到满足给定的精度为止适用范围求有根区间内的 单根 或 奇重实根数学原理:介值定理 设在上连续,即 ,或,则根据连续函数的介值定理,在内至少存在一点 ,使具体步骤:(1) 令,计算;(2) 若,则是的根,停止计算,输出结果若 ,则令,若,则令,;,有、以及相应的(3) 若 (为预先给定的精度要求),退出计算,输出结果;反之,返回(1),重复(1),(2),(3)以上方法可得到每次缩小一半的区间序列,在中含有方程的根当区间长很小时,取其中点为根的近似值,显然有以上公式可用于估计对分次数2. 迭代法基本思想:构造的一个等价方程从某个近似根出发,计算,可得一个迭代序列,这种方法称为迭代法若 收敛,即,假设连续,有即可知,的极限是的根,也就是的根当然,若发散,迭代法就失败迭代过程收敛的常用判别标准:当根区间较小,且对某一,明显小于1时,则迭代收敛 迭代法的加速:(1)松弛法:若与同是的近似值,则是两个近似值的加权平均,其中称为权重,现通过确定看能否得到加速迭代方程是:其中,令,试确定:当时,有,即当,时,可望获得较好的加速效果,于是有松弛法:,(2) Altken方法:,是它的根,是其近似根设,因为,用差商近似代替,有 ,解出,得由此得出公式 ;,这就是Altken 公式。3. 牛顿(Newton)法(牛顿切线法)基本思想:是非线性方程,一般较难解决,多采用线性化方法记:是一次多项式,用作为的近似方程的解为 记为,一般地,记 即为牛顿法公式。实验所用的软件及版本:MatlabR2012b主要内容(要点):4.分别用对分法、普通迭代法、松弛迭代法、Altken 迭代法、牛顿切法线等5种方法,求方程 的正的近似根,(建议取 时间许可的情况下,可进一步考虑的情况)实验过程记录(含:基本步骤、主要程序清单及异常情况记录等):第三题1.对分法syms x fx;a=0;b=2;fx=0.5*x-sin(x);x=(a+b)/2;i=0;ffx=subs(fx,x,x);if ffx=0; disp(the root is:,num2str(x)%第一次the root is:少加单引号第二次改的时候只把the root is了单引号,把 :放在了单引号外面else disp(i ai bi f(xi) while abs(ffx)0.0001&ab; disp(num2str(i), ,num2str(a), ,num2str(b), ,num2str(ffx) %把ffx误写为fx fa=subs(fx,x,a);ffx=subs(fx,x,x); if fa*ffx ti1Error: File: ti1.m Line: 7 Column: 23Expression or statement is incorrect-possibly unbalanced (, , or . ti1Error: File: ti1.m Line: 7 Column: 25Expression or statement is incorrect-possibly unbalanced (, , or . ti1Error: File: ti1.m Line: 21 Column: 9The input character is not valid in MATLAB statements or expressions. ti1Error: File: ti1.m Line: 6 Column: 1At least one END is missing: the statement may begin here. 正确结果 ti1i ai bi f(xi)0 0 2 -0.341471 1 2 -0.341472 1.5 2 -0.247493 1.75 2 -0.108994 1.875 2 -0.0165865 1.875 1.9375 0.0352366 1.875 1.9063 0.00886397 1.8906 1.9063 -0.00397688 1.8906 1.8984 0.00241479 1.8945 1.8984 -0.0007882910 1.8945 1.8965 0.0008113911 1.8945 1.8955 1.1094e-05所求的解是:x=1.895020,迭代步数是:11/n2. 普通迭代法syms x fx gx;gx=sin(x)/0.5;fx=0.5*x-sin(x);disp(i x f(x)x=1.5;i=0;ffx=subs(fx,x,x);while abs(ffx)0.0001; disp(num2str(i), ,num2str(i), ,num2str(ffx); x=subs(gx,x,x);ffx=subs(fx,x,x);i=i+1;enddisp(num2str(i), ,num2str(x), ,num2str(ffx)fprintf(所求的解是:x=%f,迭代步数是:%d/n,x,i)异常记录 ti2Error: File: ti2.m Line: 2 Column: 21Unexpected MATLAB expression. ti2Error: File: ti2.m Line: 8 Column: 6The expression to the left of the equals sign is not a valid target for an assignment. ti2Undefined function or variable sinx.Error in ti2 (line 2)gx=sin(x)/0.5;fx=0.5*x-sinx; 异常结果 ti2i x f(x)0 0 -0.341210 0 -0.0864850 0 0.0507390 0 -0.0332380 0 0.0206770 0 -0.0133570 0 0.00844330 0 -0.0054160 0 0.00344310 0 -0.00220170 0 0.00140280 0 -0.000895840 0 0.000571250 0 -0.000364620 0 0.000232590 0 -0.000148420 1.8956 9.4692e-05所求的解是:x=1.895610,迭代步数是:0/n正确运行结果 ti2i x f(x)0 0 -0.247491 1 0.0861242 2 -0.0570583 3 0.0346854 4 -0.0225655 5 0.0141636 6 -0.00911597 7 0.00578028 8 -0.00370179 9 0.00235610 10 -0.001505511 11 0.0009596512 12 -0.0006126813 13 0.0003907614 14 -0.0002493915 15 0.0001590916 16 -0.0001015217 1.8956 6.477e-05所求的解是:x=1.895573,迭代步数是:17/n 3.松弛迭代法syms fx gx x dgx;gx=sin(x)*2;fx=0.5*x-sin(x);dgx=diff(gx,x);x=1.5;k=0;ggx=subs(gx,x,x);ffx=subs(fx,x,x);dgxx=subs(dgx,x,x);% 在写ffx=subs(fx,x,x)时少加了一个单引号disp(i x c)while abs(ffx)0.0001; c=1/(1-dgxx); %把dgxx写成dgx disp(num2str(i), ,num2str(x), ,num2str(c)%漏写单引号 x=(1-c)*x+c*ggx;i=i+1; ggx=subs(gx,x,x);ffx=subs(fx,x,x);dgxx=subs(dgx,x,x);end disp(num2str(i), ,num2str(x), ,num2str(c)fprintf(所求的解是:x=%f,迭代步数是:%d/n,x,i)异常记录 ti4Error: File: ti4.m Line: 4 Column: 32Unexpected MATLAB expression. ti4Error: File: ti4.m Line: 8 Column: 39Unexpected MATLAB expression. ti4Error: File: ti4.m Line: 8 Column: 40Unexpected MATLAB expression. ti4i x cUndefined function max for input arguments of type sym.Error in num2str (line 84) xmax = double(max(abs(x(:);运行结果 ti4i x c17 1.5 1.164818 2.0766 0.5078919 1.9105 0.6000920 1.8956 0.6103921 1.8955 0.61039所求的解是:x=1.895494,迭代步数是:21 4.altken法syms fx gx x ;gx=sin(x)*2;fx=0.5*x-sin(x);disp(i x x1 x2)x=1.5;i=0;ffx=subs(fx,x,x);while abs(ffx)0.0001; u=subs(gx,x,x);v=subs(gx,x,u); disp(num2str(i), ,num2str(x), ,num2str(u), ,num2str(v) x=v-(v-u)2/(v-2*u+x);i=i+1;ffx=subs(fx,x,x);%漏写了2end disp(num2str(i), ,num2str(x), ,num2str(v),num2str(v) % num2str拼错了 fprintf(所求的解是:x=%f,迭代步数是:%dn,x,i)异常记录 ti5Error: File: ti5.m Line: 9 Column: 13Unexpected MATLAB expression. ti5i x x1 x2Undefined function num2sr for input arguments of type double.Error in ti5 (line 8) disp(num2sr(i), ,num2str(x), ,num2str(u), ,num2str(v)运行结果 ti5i x x1 x20 1.5 1.995 1.82271 1.8672 1.9128 1.88422 1.8952 1.8957 1.89543 1.8955 1.8957 1.8954所求的解是:x=1.895494,迭代步数是:35.牛顿法syms x fx gx;fx=0.5*x-sin(x);gx=diff(fx,x);x1=0.5;x2=1.5;x3=3;i=0;disp(i x1 x2 x3)fx1=subs(fx,x,x1);fx2=subs(fx,x,x2);fx3=subs(fx,x,x3);gx1=subs(gx,x,x1);gx2=subs(gx,x,x2);gx3=subs(gx,x,x3);while abs(fx1)0.0001|abs(fx2)0.0001|abs(fx3)0.0001; disp(num2str(i), ,num2str(x1), ,num2str(x2), ,num2str(x3) x1=x1-fx1/gx1;x2=x2-fx2/gx2;x3=x3-fx3/gx3;i=i+1; fx1=subs(fx,x,x1);fx2=subs(fx,x,x2);fx3=subs(fx,x,x3); gx1=subs(gx,x,x1);gx2=subs(gx,x,x2);gx3=subs(gx,x,x3);end disp(num2str(i), ,num2str(x1), ,num2str(x2), ,num2str(x3) fprintf(所求的解是:x1=%f,x2=%f,x3=%f,迭代步数:%dn,x1,x2,x3,i)运行结果 ti6i x1 x2 x30 0.5 1.5 31 -0.10762 2.0766 2.0882 0.00083966 1.9105 1.91223 -3.9465e-10 1.8956 1.89574 0 1.8955 1.8955所求的解是:x1=0.000000,x2=1.895494,x3=1.895494,迭代步数:4 实验报告结果与实验总结实验报告结果:第三题1.对分法运行结果 ti1i ai bi f(xi)0 0 2 -0.341471 1 2 -0.341472 1.5 2 -0.247493 1.75 2 -0.108994 1.875 2 -0.0165865 1.875 1.9375 0.0352366 1.875 1.9063 0.00886397 1.8906 1.9063 -0.00397688 1.8906 1.8984 0.00241479 1.8945 1.8984 -0.0007882910 1.8945 1.8965 0.0008113911 1.8945 1.8955 1.1094e-05所求的解是:x=1.895020,迭代步数是:11/n3. 普通迭代法运行结果 ti2i x f(x)0 0 -0.247491 1 0.0861242 2 -0.0570583 3 0.0346854 4 -0.0225655 5 0.0141636 6 -0.00911597 7 0.00578028 8 -0.00370179 9 0.00235610 10 -0.001505511 11 0.0009596512 12 -0.0006126813 13 0.0003907614 14 -0.0002493915 15 0.00015909
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 钢筋工设备调试考核试卷及答案
- 沙坪坝防水防腐施工方案
- 油气管道腐蚀防护效果分析报告
- 陶瓷成型施釉工标准化作业考核试卷及答案
- 糖果成型工设备调试考核试卷及答案
- 公益讲座活动策划及执行方案
- 同辈心理咨询方案
- 店长人力资源配置分析报告
- 竹蔗网络营销策划方案
- 会计考试题库及答案6
- 新疆乌鲁木齐市2025年通信工程师考试高、中级(通信专业实务传输与接入)题库及答案
- 关于医院“十五五”发展规划(2026-2030)
- 中外航海文化知到课后答案智慧树章节测试答案2025年春中国人民解放军海军大连舰艇学院
- 部编高教版2023·职业模块 中职语文 2.《宁夏闽宁镇:昔日干沙滩今日金沙滩》 课件
- 数学人教A版(2019)必修第一册1.3集合的基本运算(共17张ppt)
- 塑料电镀原理与工艺介绍(12013-03-08)
- 北京理工大学应用光学课件(大全)李林
- 新概念英语第一册第四课(课堂PPT)
- 井下作业公司试气作业培训
- 化学元素周期表word(可打印)
- 室外监控系统防雷设计方案范本
评论
0/150
提交评论