数学实验1方程求解.docVIP

湖 南 商 学 院 学 生 实 验 报 告 实验名称 方程求解 课程名称 常用算法的matlab设计 姓 名 符洋 班 级 信科1101 时 间 2013.9.27 成 绩 教师签名 课程 名称 常用算法的Matlab设计 实验项目 名 称 方程求解 实验项目类型 验证 演示 综合 设计 其他 时间 成 绩 √ √ 实验目的 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句 特别是循环、条件、控制等语句 ； [4] 通过范例展现求解实际问题的初步建模过程； 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。x 2*pi:2*pi/100:8*pi; y1 sin x ; y2 1./ x+eps ; plot x,y1,x,y2 结果： 图1 图像分析： 由图1可以看出y1和y2有无数组解，且呈现周期分布，每2*pi的长度内会出现两个交点，现在将区间缩小至[2*pi’4*pi];得到的图像如下 图2 再将区间缩小至[9,9.5],得到的图形如下 图3 由图1、图2、图3可以看出x sin x 1有无数个根，且呈周期出现，根的值为X 9.31+k*pi; （K为任意整数） 2．将方程x5 +5x3- 2x + 1 0 改写成各种等价的形式进行迭代，观察迭代是否收敛，并给出解释。x 1;y 1;z 1; for k 1:10 x x^5+5*x^3+1 ^2; y 2*y-1-5*y^3 ^ 1/5 ; z 2*z-1-z^5 /5 ^ 1/3 ; h [k,x,y,z] end 实验结果： h 1.0000 49.0000 1.0675 + 0.7756i 0 h 1.0e+016 * 0.0000 8.0125 0.0000 - 0.0000i 0.0000 + 0.0000i h 1.0e+169 * 0.0000 1.0906 0.0000 + 0.0000i 0.0000 + 0.0000i h 4.0000 Inf 1.8323 - 0.6099i 0.4701 + 0.2548i h 5.0000 Inf 1.8377 + 0.8353i 0.3935 + 0.2364i h 6.0000 Inf 1.9443 - 0.7354i 0.3686 + 0.2837i h 7.0000 Inf 1.9458 + 0.8285i 0.3947 + 0.3044i h 8.0000 Inf 1.9884 - 0.7881i 0.4103 + 0.2922i h 9.0000 Inf 1.9889 + 0.8249i 0.4057 + 0.2811i h 10.0000 Inf 2.0055 - 0.8091i 0.3983 + 0.2821i 结果分析： 由结果可以看出，采用x x^5+5*x^3+1 ^2的结果是发散的 而y 2*y-1-5*y^3 ^ 1/5 ;和 z 2*z-1-z^5 /5 ^ 1/3 ;两种方式迭代结果收敛，且第三种比第二种迭代收敛速度快。 3．求解下列方程组 function f died x f x 1 +x 2 -exp x 1 -exp x 2 ; 在命令窗口输入以下指令： options optimset Jacobian,off,Display,iter ; x fsolve died,[0.5,0.5],options 运行结果 Iteration Func-count Residual Step-size derivative 0 3 5.27824 1 9 4.00118 0.894 -0.0745 Conditioning of Gradient Poor - Switching To LM method 2 16 4 1.27 -1.61e-009 2.55741 3 17 4 1 -2.78e-015 2.55795 x 1.0e-007 * 0.1729 0.1729 结果分析 由运行结果可知x1 1.729*10^ -8 ; x2 1.729*10^ -8 2 实验代码： [x1,x2,x3] solve x1^2-5*x2^2+7*x3^2+12 0,3*x1*x2+x1*x3-11*x1 0,2*x2*x3+40*x1 0,x1,x2,x3 实验结果： double x1 ans 1.0e+002 * 0 0 0 0 0.0100 -0.0031