用MATLAB求解高等数学问题2.docVIP

实验三：用MATLAB求解高等数学问题-2 姓名： 学号： 专业： 成绩： 一、 实验名称：用MATLAB求解高等数学基本计算问题-2 二、 实验目的： 学会利用MATLAB求解有关非线性方程（组）符号解问题； 学会利用MATLAB求解有关非线性方程（组）数值解问题； 学会利用MATLAB求解梯度向量并画图； 学会利用MATLAB求解数值积分问题； 学会利用MATLAB解决微分方程（组）符号求解问题； 学会利用MATLAB解决微分方程（组）数值求解问题。 三、 实验准备： 复习高等数学中相关主题的理论知识和方法。 四、 实验内容 求下列方程的精确解（符号解）； （1） x3-5x2 +3x + 6 = 0 ； （2） x4 + x2 - a = 0. 求方程x4 = 2在限值条件（一2x2）下的根； 画出下面两个椭圆的图形，并求出它们所有的交点坐标； (x-2)2 + (j + 2x-3)2=5 18(x-3)2 +j2 = 36 作图表示函数z =亦一宀（-1 xl,O j2）沿兀轴方向的梯度； 求函数/（u,z） =兀+ 2y + 3z在平面兀一y + z = 1与柱面x2 + j2 = 1的交线上的最大值; 用梯形求积公式和辛普森（抛物线）公式计算积分 并估计误差。 7.若用复化梯形求积公式求e~xdx的近似值，问要将积分区间［0,1］分成多少等份才能保证计算结果 有四位有效数字？若用复化抛物线求积公式呢？ 提示：（1）复化梯形公式，把区间［a,〃］等分为〃等分，节点 ya + kh, TOC \o 1-5 \h \z \o Current Document \h b _ a h k = 0,l,???,〃，h =——，则积分 / - -［/（a） + f（b）4- 2V （a + kh）］； n 2 （2）复化抛物线公式，把区间［a,b］等分为2加等分，力=2二纟， h m m-I 则积分 1 = 一 f(a) + f(b) + + (2k 一 1)创+ 2》(a + 2kh) \o Current Document \h 3 r=i o,y求解下列微分方程的解： o,y ⑴ /-4y + 13j = O,j|x=o = 0,y\x=Q = 6； (2) y^-yf-2y = (l-2x)ex. 解下列微分方程： x = 2x + 3j,/ = 2x + x(0) = -2,j(0) = 2.0vfvlO,作相平面图; 提示：相平面图就是画在一个平面内，即横轴为工，纵轴为y. /-0.01(y)2 + 2j = sinr, j(0) = 0,/(0) = 1,0 Z 5 ,作丿的图。 d2x2(x + //) /z(x-2) d2x 2(x + //) /z(x-2) TOC \o 1-5 \h \z \o Current Document \h = 2^x-. , dt rX r2 \o Current Document \h =_2竺+厂绰-够dt r23 其中，“ = l/82?45,兄=1 一“,厶=J(x~F“)2 + y2，E = ,试在初值 x(0) = 1.2,x(0) = 0,j(0) = 0,y(0) = —10/100下求解,并绘制飞船轨迹图。 五、程序及实验结果 1. (1)程序： syms x; f=symC x 飞-5*x2+x+6二O); y=sol ve (f, x) 输出结果：y二 22/(9*(43/54 + (50371/2)*1)/6)71/3)) + (43/54 + (50371/2)*1)/6)71/3) + 5/3 5/3 - (43/54 + (50371/2)*1)/6)71/3)/2 - 11/(9*(43/54 + (503 (l/2)*i)/6厂(1/3))- (3, (1/2)*(22/(9*(43/54 + (50371/2)*1)/6) (1/3)) - ((50371/2)*1)/6 + 43/54厂(1/3))*i)/2 5/3 - (43/54 + (50371/2)*i)/6) 71/3)/2 - 11/(9*(43/54 + (503, (l/2)*i)/6厂(1/3)) + (3八(1/2)*(22/(9*(43/54 + (503, (l/2)*i)/6厂(1/3)) - ((50371/2)*1)/6 + 43/54厂(1/3) )*i)/2 (2)程序： 输出结果： 程序： syms x a f=sym(，x4+x“2-8=0) y二solve (f, x) 输出结果：f 二x4 + x^2 - a = 0 y 二(—(4*a + 0 71/2)/2 一 1/2厂(1/2) ((4*a