《数值分析上机(四)》由会员分享，可在线阅读，更多相关《数值分析上机(四)（4页珍藏版）》请在金锄头文库上搜索。
1、Author by prs,版权所有Ps :题目均来自数值分析第五版作者：李庆扬，王能超，易大义编出 版社：清华大学出版社误差分析问题：求下列方程的实根(1) x2 3x + 2 ex = 0(2) x3 +2x2+10x - 20 = 0要求：(1)设计一种不动点迭代法，要使迭代序列收敛，然后再用斯特芬森加速迭代，计算到I xk-xj | 10为止。(2)用牛顿迭代，同样计算到| xk-xy |1.0e-8x(k+1)=x(k)-fun1(x(k)/dfun1(x(k);des=abs(x(k+1)-x(k);k=k+1; endy=x(k); k=k; endfunction y,k= s
2、itefensen(x0,f)%SITEFENSEN Summary of this function goes here % Detailed explanation goes here%x0为初值，n为迭代次数，f为迭代函数x(1)=x0;des=1;k=1;while des1.0e-8 y(k)=f(x(k); z(k)=f(y(k);x(k+1)=x(k)-(y(k)-x(k)A2/(z(k)-2*y(k)+x(k);des=abs(x(k+1)-x(k);k=k+1;endy=x(k);k=k;end%fun的导数function y= dfun(x)%DFUN Summary o
3、f this function goes here% Detailed explanation goes herey=3*xA2+ 4*x+10;end%fun1的导数function y = dfun1( x )%DFUN1 Summary of this function goes here% Detailed explanation goes here y=2*x-exp(x)-3;endclearclc%不动点迭代法n=100;x0=0.5;% 初值k=0;des=1;while des1.0e-8x=(x0A2+2-exp(x0)/3;%2*log(x)+log(3)%x=(-x0A
4、3-10*x0+20)/2*(x0+eps);des=abs(x-x0);k=k+1;x0=x;end disp(不动点迭代解-) fprintf(%fn,x) disp(迭代次数-) fprintf(%dn,k) disp(误差-)fprintf(%fn,abs(0-fun1(x) %斯蒂芬森加速 f=(xA2+2-exp(x)/3;f=inline(f);yy,kk=sitefensen(0.5,f); disp(斯蒂芬森加速解-) fprintf(%fn,yy) disp(迭代次数-) fprintf(%dn,kk) disp(误差-) fprintf(%fn,abs(0-fun1(yy
5、) %牛顿迭代法 yy1,kk1=niudun(0.5);disp(牛顿迭代法解-) fprintf(%fn,yy1) disp(迭代次数-) fprintf(%dn,kk1) disp(误差-)fprintf(%fn,abs(0-fun1(yy1)第一个方程的运行结果如下： 不动点迭代解-0.257530迭代次数-14误差-0.000000斯蒂芬森加速解-0.257530迭代次数-5误差-0.000000牛顿迭代法解-0.257530迭代次数-5误差-0.000000结论：由上述三个表格可以看由去在迭代初值相同的情况 下，斯蒂芬森加速和牛顿加速迭代次数都明显少于不动点迭 代。但在迭代初值和真
6、实值相差较大的时候，牛顿法的迭代 次数和不动点差不多， 而斯蒂芬森加速法依然保持很好的收 敛快速性。止匕外，在初值取得再大时，不动点迭代呈发散状, 故在方法的选取中，我们应该选择斯蒂芬森加速迭代法。理制度1范围本标准规定了龙腾公司合同管理工作的管理机构、职责、合同的授权委托、洽谈、承办、会签、订F阅、履行和变更、终止及争议处理和合同管理的处罚、奖励;本标准适用于龙腾公司项目建设期间的各类合同管理工作，厂内各类合同的管理，厂内所属各具法人资格的部门，参照本标准执行。2规范性引用中华人民共和国合同法龙腾公司合同管理办法3定义、符号、缩F略语4职责4.1 总经理：龙腾公司经营管理的法定代表人。负责对
7、厂内各类合同管理工作实行统一领导。以法人代表名义或授权委托他人签订各类合法合同，并对电厂负责。4.2 工程部：是发电厂建设施工安装等工程合同签订管理部门；负责签订管理基建、安装、人工技术的工程合同。4.3 经营部：是合同签订管理部门，负责管理设备、材料、物资的订购合同。4.5 合同管理部门履行以下职责:4.5.1 建立健全合同管理办法并逐步完善规范;4.5.2 参与合同的洽谈、起草、审查、签约、变更、解除以及合同的签证、公证、调解、诉讼等活动，全程跟踪和检查合同的履行质量;4.5.3 审查、登记合同对方单位代表资格及单位资质，包括营业执照、经营范围、技术装备、信誉、越区域经营许可等证件及履约能
8、力（必要时要求对方提供担保）,检查合同的履行情况;4.5.4 保管法人代表授权委托书、合同专用章，并按编号归口使用;4.5.5 建立合同管理台帐，对合同文本资料进行编号统计管理;4.5.6 组织对法规、制度的学习和贯彻执行，定期向有关领导和部门报告工作;4.5.7 在总经理领导下，做好合同管理的其他工作,4.6 工程技术部：专职合同管理员及材料、燃料供应部兼职合同管理员履行以下职责:4.6.1 在主任领导下，做好本部门负责的各项合同的管理工作，负责保管-法人授权委托书”4.6.2 签订合同时，检查对方的有关证件，对合同文本内容依照法规进行检查，检查合同标的数量、金额、日期、地点、质量要求、安全责任、违约责任是否明确，并提出补充及修改意见。重大问题应及时向有关领导报告，提出解决方案;4.6.3 对专业对口的合同统一编号、登记、建立台帐，分类整理归档。对合同承办部门提供相关法规咨询和日常协作服务工作;4.6.4 工程技术部专职合同管理员负责收集整理各类合同，建立合同统计台帐，并负责
