20115394王福临数学实验作业[解方程].docVIP

重 庆 大 学 学 生 实 验 报 告 实验课程名称 数学实验 开课实验室 学 院 计算机学院 年级 大二 专业 计科六班 学 生 姓 名 王福临 学 号 开 课 时 间 至 学年第 学期 总 成 绩教师签名 数 理 学 院 制 开课学院、实验室： DS1401 实验时间 ： 2013 年 3月 24 日 课程 名称数学实验实验项目 名 称 飞机如何定价——方程求解实验项目类型验证演示综合设计其他指导 教师陈道县成 绩实验目的 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。 基础实验 一、实验内容 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。 二、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 3．保存文件并运行； 4．观察运行结果(数值或图形)； 5．根据观察到的结果写出实验报告，并浅谈学习心得体会。 三、实验要求与任务 1．用图形放大法求解方程 x sin(x) = 1. 并观察该方程有多少个根。 2．将方程x5 +5x3- 2x + 1 = 0 改写成各种等价的形式进行迭代，观察迭代是否收敛，并给出解释。 3．用MATLAB命令求下列方程的根 1）e-3xsin(4x+2)+4 e-0.5xcos(2x) =0.5 2） 应用实验 1. 炮弹发射角的问题 炮弹发射视为斜抛运动，已知初始速度为200 m/s，问要击中水平距离360m、垂直距离160m 的目标，当忽略空气阻力时，发射角应多大？此时炮弹的运行轨迹如何？试进行动态模拟。 进一步思考：如果要考虑水平方向的阻力，且设阻力与（水平方向）速度成正比，系数为 0.1（1/s），结果又如何？此时炮弹的运行轨迹如何？试进行动态模拟。 实验过程及结果 用图形放大法求解方程 x sin(x) = 1. 并观察该方程有多少个根。 程序： x=-100:0.01:100; y=x.*sin(x); plot(x,y),hold on, line([-100,100],[0,0]) 运行程序： 逐次缩小区间，观察一个根在-4到-2 取x=-10到10 取x=-10到0 取x=-4到0 取x=-4到-2 结果：由图可知，该方程有无数组解 将方程x5 +5x3- 2x + 1 = 0 改写成各种等价的形式进行迭代，观察迭代是否收敛，并给出解释。 解：①迭代函数为，算法设计为： x1=0; x2=(x1^5+5*x1^3+1)/2; while abs(x1-x2)10^(-5) x1=x2; x2=(x1^5+5*x1^3+1)/2; end x1 输出结果为： x1 = Inf 因此x=j(x)迭代不收敛，则不直接使用j(x)迭代，用加速迭代函数 ， 算法设计为： x1=0; x2=(-4*x1^5-10*x1^3+1)/(-5*x1^4-15*x1^2+2); while abs(x1-x2)10^(-5) x1=x2; x2=(-4*x1^5-10*x1^3+1)/(-5*x1^4-15*x1^2+2); end x1 输出结果为： x1 = -0.7685 ②迭代函数为，算法设计为： x1=1; x2=((2*x1-x1^5-1)/5)^(1/3); while abs(x1-x2)10^(-5) x1=x2;