数学建模MATLAB简介第二部分微积分常用算法.docVIP

第二部分 微积分常用算法 一、极限、导数和极值 limit 求极限（符号运算） diff 求导数（符号运算） taylor 函数展开成幂级数（符号运算） symsum 级数求和（符号运算） fmin 求一元函数的极小值 fmins 求多元函数的极小值 例1、计算 syms n x y; limit 1+x/n ^n,n,inf symsum -1 ^n*x^n/n,n,1,inf diff sin x*y ,x,2 例2、求sinx在x 0点的幂级数展开式，到6次幂为止。 syms x; taylor sin x ,6,0 例3、求函数在[-6，18]中的极小值。 fun ‘x^3-5*x^2+3*x+2’; x fmin fun,-6,18 f eval fun eval 函数的功能就是将括号内的字符串视为语句并运行 比如 eval y1 sin 2 和语句y1 sin 2 等价 多在循环中使用，可以对多个名字有规则的变量或文件进行操作，比如 for x 1:5 eval [y,num2str x , ,num2str x^2 ,;] end在原点附近的极大值。 fun ’x 1 ^4+x 2 ^4-4*x 1 *x 2 -5’; x fmins fun,[0,0] f -eval fun 二、积分 int 符号积分 trapz 数值积分 例5、计算 syms x; t1 int exp -x *sin x t2 int exp -x *sin x ,0,1 t2 vpa t2 t3 int exp -x *sin x ,1,inf t3 vpa t3,5 在Matlab中，往往会出现数据精度不够，精确计算失败，或者大树吃掉小数的现象，比如3.1415926会用3.1416来计算，1在大数面前会当零处理，这样就需要自定义数据类型或者数据精度vpa（，）函数会很好的解决， 比如 3.1415926/2 1.5708 vpa 3.1415926/2,10 1.座机电话号码0 vpa 3.1415926/2,10 1.座机电话号码0 vpa 3.1415926/2 1.座机电话号机电话号码1285482 例6、计算 syms x y; iy int sqrt 1-x^2 ,y,-sqrt 1-x^2 ,sqrt 1-x^2 ; int iy,x,-1,1 例7、计算 clear; x -1:0.1:1; y exp -x.^2 ; trapz x,y 三、求解微分方程（组） dsolve 符号解析解 例8、求的通解。 s dsolve ‘Dy a*y+b’ %自变量缺省值为t 例9、求解 dsolve ‘Dy y-2*t/y’,’y 0 1’ 例10、求解方程 s dsolve ‘D2y cos 2*x -y’,’y 0 1’,’Dy 0 0’,’x’ ; simplify s 例11、求解方程组 s dsolve ‘Df f+g’,’Dg -f+g’,’f 0 1’,’g 0 2’ ; s.f, s.g