实验2方程模型及其求解算法.docVIP

实验2 方程模型及其求解算法 一、实验目的及意义 [1] 复习求解方程及方程组的基本原理和方法； [2] 掌握迭代算法； [3] 熟悉MATLAB软件编程环境；掌握MATLAB编程语句(特别是循环、条件、控制等语句)； [4] 通过范例展现求解实际问题的初步建模过程； 通过该实验的学习，复习和归纳方程求解或方程组求解的各种数值解法（简单迭代法、二分法、牛顿法、割线法等），初步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。 二、实验内容 1．方程求解和方程组的各种数值解法练习 2．直接使用MATLAB命令对方程和方程组进行求解练习 3．针对实际问题，试建立数学模型，并求解。 三、实验步骤 1．开启软件平台——MATLAB，开启MATLAB编辑窗口； 2．根据各种数值解法步骤编写M文件 3．保存文件并运行； 4．观察运行结果(数值或图形)； 5．根据观察到的结果写出实验报告，并浅谈学习心得体会。 四、实验要求与任务 基础实验 用图形放大法求解方程 x sin(x) = 1. 并观察该方程有多少个根。 用matlab实现： syms x y x=-3:0.1:3 y=x.*sin(x)-1 plot(x,y) grid ezplot(x*sin(x)-1,[-100,100]);grid 如图： 放大图像求得两个个解为-1.11和1.11 所以有无数个解。 2．将方程x5 +5x3- 2x + 1 = 0 改写成各种等价的形式进行迭代，观察迭代是否收敛，并给出解释。 首先先观察其本身图像： x1=-2:0.01:2; x2=-0.8:0.01:0.7; y1=x1.^5 +5*x1.^3-2*x1+1; y2=x2.^5 +5*x2.^3-2*x2+1; subplot(2,2,1),plot(x1,y1), title(原图) ,grid on, subplot(2,2,2),plot(x4,y4), title(放大) ,grid on 所以解在-0.76—-0.77之间 构造迭代函数： 用matlab实现如下： x(1)=-0.77;y(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); end k=1:20 subplot(2,2,1), plot(k,x), title(f(1)); subplot(2,2,2), plot(k,y), title(f(2)); 如图： 经观察发现不收敛。 设定迭代函数为： 初始值设置为0.75 观察发现收敛。 3．求解下列方程组 直接使用MATLAB命令：solve()和fsolve()对方程组求解。 (1) syms x y [x,y]=solve(2*x-y-exp(-x)=0,2*y -x -exp(-y)=0,x,y) x = 0.5671 y = 0.5671 (2) syms x1 x2 x3 syms x1 x2 x3 [x1,x2,x3]=solve(x1^2-5*x2^2+7*x3^2+12,3*x1*x2+x1*x3-11*x1=0,2*x2*x3+40*x1=0,x1,x2,x3); x1=double(x1), x2=double(x2), x3=double(x3) x1 =1.0e+002 * 0 0 0 0 0.0100 -0.0031 -3.8701 - 0.3270i -3.8701 + 0.3270i x2 = -1.5492 1.5492 0 0 5.0000 2.9579 -0.3123 +50.8065i -0.3123 -50.8065i x3 =1.0e+002 * 0 0 0 + 0.0131i 0 - 0.0131i -0.0400 0.0213 0.1194 - 1.5242i 0.1194 + 1.5242i 4．编写用二分