实验4常微分方程积分与微分的运算-答案.docVIP

实验4 常微分方程、积分与微分的运算-答案 1、用solve函数求下列非线性方程组的解 [x,y] solve cos x +y*exp x -2 0,sin y +x*exp y -2 0 x .座机电话号码 y .座机电话号码0座机电话号码0座机电话号码668096 2、对于二阶微分方程 （1）利用ode45方法，求当，在时的数值图解。 （2）利用dsolve函数求当，时的特解，画出时的曲线，并与（1）中的数值图解作比较。 function dx ff t,x dx [x 2 ;-2*x 2 -x 1 +2*sin t ]; 建立调用函数 x0 [0 1]; [t,x] ode45 ff,[0,30],x0 plot t,x :,1 （2） 求的解： y dsolve D2y+2*Dy+y 2*sin t ,y 0 0,Dy 0 1,t y exp -t +2*exp -t *t-cos t 作曲线： t 0:0.1:30; y exp -t +2*exp -t .*t-cos t ; plot t,y 3、分别用Simpson法、 Newton-Cotes法、梯形法trapz以及符号积分函数int计算定积分。 function f ff x f sin x ; Simpson法： 在主窗口调用： [S,n] quad ff,0,pi S 2.0000 n 33 Newton-Cotes法： [S,n] quad8 ff,0,pi S 2.0000 n 18 trapz方法： x 0:pi/100:pi; y sin x ; z trapz x,y z 1.9998 符号积分int法： syms x int sin x ,0,pi ans 2 4、用符号积分int法求下列积分： （1） （2） （3） （4） syms x int 3-x^2 ^3 ans 27*x-1/7*x^7+9/5*x^5-9*x^3 （2） syms x t int 5*x*t/ 1+x^2 ,t ans 5/2*x*t^2/ 1+x^2 （3） syms x int abs 1-x ,1,2 ans 1/2 （4） syms x t int 4*x/t,t,2,sin x ans 4*log sin x *x-4*log 2 *x 5、求下列函数的导数。 （1），求 （2），求 ， syms t a d diff a*sin t /diff a*cos t d -cos t /sin t diff d /diff a*cos t ans - 1+cos t ^2/sin t ^2 /a/sin t （2） syms t f diff 1-sqrt t / 1+sqrt t f -1/2/t^ 1/2 / 1+t^ 1/2 -1/2* 1-t^ 1/2 / 1+t^ 1/2 ^2/t^ 1/2 t 4; f -1/2/t^ 1/2 / 1+t^ 1/2 -1/2* 1-t^ 1/2 / 1+t^ 1/2 ^2/t^ 1/2 f -0.0556