2014重庆大学数学实验报告.docVIP

重 庆 大 学 学 生 实 验 报 告 实验课程名称 开课实验室 学 院 年级 专业班 学 生 姓 名 学 号 开 课 时 间 至 学年第 学期 总 成 绩 教师签名 数 学 与 统 计 学 院 制 开课学院、实验室： DS1422 实验时间 ：2015年 3 月 20 日 课程 名称 数学实验 实验项目 名 称 方程模型及其求解算法 实验项目类型 验证 演示 综合 设计 其他 指导 教师 肖剑 成 绩 实验目的 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。 实验内容 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。 基础实验 一、问题重述 1．用图形放大法求解方程 x sin(x) = 1. 并观察该方程有多少个根。 2．将方程x5 +5x3- 2x + 1 = 0 改写成各种等价的形式进行迭代，观察迭代是否收敛，并给出解释。 3．求解下列方程组 直接使用MATLAB命令：solve()和fsolve()对方程组求解。 4．编写用二分法求方程根的函数M文件。 5. 设非线性方程组为 其中已知，随机产生数据后，用fsolve解这个方程组。 6. 使用fsolve计算方程组的解时，为验证初值是否对解有影响，采用随机产生的100组随机数作为初始值，依次进行求解。 二、实验过程（一般应包括实验原理或问题分析，变量说明、程序、调试情况记录、图表等， 实验结果及分析） 1.解答： 程序： x=-9:0.01:9; y=x.*sin(x)-1; plot(x,y),hold on; line([-9,9],[0,0]) 运行结果： 再放大区间：程序： x=-50:0.01:50; y=x.*sin(x)-1; plot(x,y),hold on; line([-50,50],[0,0]) 运行结果： 综上：该方程有无数多个根。 2.解答： 程序： x=2;y=2;z=2; for k=1:30; x=x.^5+5*x.^3+1; y=-x.^3/5+2/(5*x)-1/(5*x.^2); z=-5/x+2/(x.^3)-1/(x.^4) x,y,z; end 运行结果： x = 73 y = -7.7803e+04 z = -2.4096e-09 x = 2.0750e+09 y = -1.7869e+27 z = -1.2998e-46 x = 3.8469e+46 y = -1.1386e+139 z = -5.9351e-233 x = 8.4244e+232 y = -Inf z = 0 迭代失败，使用快速迭代法； 程序： x=2;y=2;z=2; for k=1:30 x=(-4*x^5-10*x^3+1)/(2-5*x^4-15*x^2); y=(2*y^6+4*y^2-3*y)/(5*y^3+3*y^5+2*y-2); z=(8*z^2-2*z)/(z^5+5*z^3+6*z-1); x,y,z; end 运行结果： x = 1.5000 y = 1 x = 1.1062 y = 0.3750 x = 0.8037 y = 0.5777 x = 0.5660 y = -1.0368 x = 0.3151 y = -0.7473 x =