解非线性方程实验.docVIP

实验?2??解非线性方程实验 ? ? ?? ??专业班级?数学122班????? ????学号?201212010215?????? ??? ?? ??姓名??拓林芳???????????? ????报告日期??2014年4月4日? . ? 实验类型：●验证性实验??○综合性实验??○设计性实验 实验目的：进一步熟练掌握解非线性方程的二分法算法、Steffensen算法、牛顿迭代法，提高编程能力和解算非线性方程问题的实践技能。 实验内容：用二分法算法(取[a,b]=[1,2])、Steffensen算法(取)、牛顿迭代法(取)解算非线性方程的根 实验原理???二分法算法 ????????? Steffensen算法 ??????????牛顿迭代法 ? 实验步骤 1?要求上机实验前先编写出程序代码 2?编辑录入程序 3?调试程序并记录调试过程中出现的问题及修改程序的过程 4?经反复调试后，运行程序并验证程序运行是否正确。 5?记录运行时的输入和输出。 实验思考?试验改变?Steffensen算法的迭代初值进行计算，试验改变?牛顿迭代算法的迭代初值进行计算，分别对比计算结果并作出分析总结 数学建模与应用：请建立数学模型并应用你的程序解决一个物理计算问题--木质球体浸入水中深度问题。假若球体是由密度为ρ=0.638且半径为r=10cm的长叶松木质球体构成的，问球体浸入水中深度是多少？又假若球体是由密度为ρ=0.638且三个长半轴分别为a=10cm,b=15cm,c=20cm的长叶松木质椭圆球体构成的，问球体浸入水中深度是多少？（建模可参考课本32页） 实验报告书写要求：根据实验情况和结果撰写并递交实验报告。实验报告应当有算法原理简介，算法流程图，程序代码，运行调试记录，计算结果阐释；对于数学建模的应用问题要有分析建模的过程、求解计算的过程、实际问题解答或阐释的内容。 实验总结(学会了......;?掌握了......;?训练了......;?发现了......;?今后学习中......有待提高。） 参考程序 1.?????二分法 1、实验原理： ??? ??二分法 2、实验步骤 ? 1建立函数 function?y=f3(x) y=x.^3-x-1 2建立函数 function[c,err,yc]=bisect(f,a,b,delta) ya=feval(f,a); yb=feval(f,b); if?ya*yb0 ? ??return end max1=1+round((log(b-a)-log(delta))/log(2)); for?k=1:max1 ? ? c=(a+b)/2; ? ? yc=feval(f,c) ? ??if?yc==0 ????? ? a=c; ????? ? b=c; ? ??elseif?yb*yc0 ????????? ? b=c; ????????? ? yb=yc; ????? ??else ????????? ? a=c; ????????? ? ya=yc; ????? ??end ????? ??if?b-adelta ????????? ??break ????? ??end ? ??end ? ? c=(a+b)/2; ? ? err=abs(b-a); ? ? yc=feval(f,c); 3调用函数 [c,err,yc]=bisect(f3,1,2,0.0001) y = ? ? -1 y = ?? ? 5 y = ? ? 0.8750 yc = ? ? 0.8750 y = ?? -0.2969 yc = ?? -0.2969 y = ? ? 0.2246 yc = ? ? 0.2246 y = ?? -0.0515 yc = ?? -0.0515 y = ? ? 0.0826 yc = ? ? 0.0826 y = ? ? 0.0146 yc = ? ? 0.0146 y = ?? -0.0187 yc = ?? -0.0187 y = ?? -0.0021 yc = ?? -0.0021 y = ? ? 0.0062 yc = ? ? 0.0062 y = ? ? 0.0020 yc = ? ? 0.0020 y = ?-4.6595e-005 yc = ?-4.6595e-005 y = ? 9.9479e-004 yc = ? 9.9479e-004 y = ? 4.7404e-004 yc = ? 4.7404e-004 y = ? 2.1371e-004 yc = ? 2.1371e-004 y = ? 8.3552e-005 c = ?? 1.3247 err = ? 6.1035e-005 yc = ? 8.3552e-005 2、牛顿迭代 ??1、实验原理