《重庆大学 数学实验报告 机械二零一六年机自 实验二》.docVIP

重 庆 大 学 学 生 实 验 报 告 实验课程名称 数学实验 开课实验室 DS1401 学 院 机械工程学院 年级 09级 专业班 机自6 学 生 姓 名 学 号 2009 开 课 时 间 2010 至 2011 学年第 2 学期 总 成 绩 教师签名 实验2 方程模型及其求解算法 一、实验目的及意义 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。 二、实验内容 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。 三、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 3．保存文件并运行； 4．观察运行结果(数值或图形)； 5．根据观察到的结果写出实验报告，并浅谈学习心得体会。 四、实验要求与任务 基础实验 1．用图形放大法求解方程 x sin(x) = 1. 并观察该方程有多少个根。 程序：x=linspace(-10,10,100); y=1./(sin(x)+eps); y1=zeros(size(x)); plot(x,y,x,y1) 结果： 结论：由图可知方程有无穷多个根 2．将方程x5 +5x3- 2x + 1 = 0 改写成各种等价的形式进行迭代，观察迭代是否收敛，并给出解释。 程序：（1）画图 x1=-6:0.01:6; x2=-3:0.01:3; x3=-1:0.01:1; x4=-0.8:0.01:-0.75; y1=x1.^5 +5*x1.^3-2*x1+1; y2=x2.^5 +5*x2.^3-2*x2+1; y3=x3.^5 +5*x3.^3-2*x3+1; y4=x4.^5 +5*x4.^3-2*x4+1; subplot(2,2,1),plot(x1,y1) ,title(图 (1)) ,grid on, subplot(2,2,2),plot(x2,y2) ,title(图 (2)),grid on, subplot(2,2,3),plot(x3,y3) ,title(图 (3)),grid on, subplot(2,2,4),plot(x4,y4) ,title(图 (4)) ,grid on, 有图可知x的初值在（-0.78，0.76）之间 （2）构造迭代函数 (3)迭代 设定初值 函数 用MATLAB编程： x(1)=-0.77;y(1)=-0.77;z(1)=-0.77; for i=2:20 x(i)=(x(i-1)^5+5*x(i-1)^3+1)/2; y(i)=(-y(i-1)^5+2*y(i-1)-1)/(5*y(i-1)^2); z(i)=(-5*z(i-1)^3+2*z(i-1)-1)/(z(i-1)^4); end n=1:20 subplot(221), plot(n,x), title(f(1)); subplot(222), plot(n,y), title(f(2)); subplot(223), plot(n,z), title(f(3)); 结果显示： 如图显示，显然不收敛。 （4）利用加速迭代收敛法变形后： 再次迭代： 设定初值 函数 用MATLAB编程： x(1)=-0.77;y(1)=-0.77;z(1)=-0.77; for i=2:30 x(i)=(-4*x(i-1)^5-10*x(i-1)^3+1)/(2-5*x(i-1)^4-15*x(i-1)^2); y(i)=(2*y(i-1)^6+4*y(i-1)^2-3*y(i-1))/(5*y(i-1)^3+3*y(i-1)^5+2*y(i-1)-2); z(i)=(8*z(i-1)^2-2*z(i-1))/(z(i-1)^5+5*z(i-1)^3+6*z(i-1)-1); end n=1:30 subplot(2