5-符号运算.ppt

一元函数积分 鲁 坚 jianlu@szu.edu.cn 深圳大学 数学与计算科学学院 syms x n int(x^n) Ans= x^(n+1)/(n+1) int(1/(1+x^2)) Ans= atan(x) syms x int(x^7,0,1); Ans= 1/8 int(1/x,1,2) Ans= log(2) syms x a f=‘a*x*sin(x)’; F1=int(f) F2=int(f,a) F3=int(f,0,2*pi) 结果： F1=a*(sin(x)-x*cos(x)) F2=1/2*a^2*x*sin(x) F3=-2*pi*a syms x f=exp(-x^2) F1=int(f,x,-inf,inf) 结果： F1=pi^(1/2) 微分方程 符号微分方程求解 在MATLAB中，用大写字母D表示导数。例如，Dy表示y，D2y表示y，Dy(0)=5表示y(0)=5。D3y+D2y+Dy-x+5=0表示微分方程y+y+y-x+5=0。符号常微分方程求解可以通过函数dsolve来实现，其调用格式为： dsolve(e,c,v) 该函数求解常微分方程e在初值条件c下的特解。参数v描述方程中的自变量，省略时按缺省原则处理，若没有给出初值条件c，则求方程的通解。 dsolve在求常微分方程组时的调用格式为： dsolve(e1,e2,…,en,c1,…,cn,v1,…,vn) 该函数求解常微分方程组e1,…,en在初值条件c1,…,cn下的特解，若不给出初值条件，则求方程组的通解，v1,…,vn给出求解变量。 求 的通解 y=dsolve(‘Dy=a*y’) y= exp(a*t)*C1 求 的特解 y=dsolve(‘Dy=a*y’,’y(0)=1’) y= exp(a*t) 求 的二阶通解 y=dsolve(‘D2y=-a^2*y’) y= C1*cos(a*t)+C2*sin(a*t) 求 的二阶特解 y=dsolve(‘D2y=-a^2*y’,’y(0)=1,Dy(pi/a)=0’) y= cos(a*t) 求 的通解 y=dsolve(‘Dy-(x^2+y^2)/x^2/2’,’x’) y= x*(-log(x)+2+C1)/(-log(x)+C1) 求 的特解，y(2)=1 y=dsolve(‘Dy-x^2/(1+y^2)’,’y(2)=1’,’x’) 求 的通解 [x,y]=dsolve(‘Dx=4*x-2*y’,’Dy=2*x-y’,’t’) x = -1/3*C1+4/3*C1*exp(3*t)-2/3*C2*exp(3*t)+2/3*C2 y = 2/3*C1*exp(3*t)-2/3*C1+4/3*C2-1/3*C2*exp(3*t) 多元函数微分 多元函数定义 二元符号表达式函数定义法 例 定义二元函数 syms x y z=5*x^2+4*y^3+sin(x*y) M文件定义函数 function z=f1(x,y) z=5*x.^2+4*y.^3+sin(x.*y) 绘出上例函数图形。 x=linspace(-10,10,40); y=3*x; [x,y]=meshgrid(x,y); z=f1(x,y); surf(x,y,z) shading interp 结果 多元函数偏导数和高阶偏导数 五种格式 diff(z,x) %函数对x求偏导数 diff(z,y) %函数对y求篇导数 diff(z,x,2) %函数对x求2阶偏导数 diff(z,y,2) %函数对y求2阶偏导数 diff(diff(z,x),y) %函数对x再对y求二阶偏导数 syms x y z Z=x^2*y^3+sin(x*y) F1=diff(z,x) F2=diff(z,y) F3=diff(z,x,2) F4=diff(z,y,2) F5=diff(diff(z,x),y) x=1;y=pi;F6=eval(f1) 多元函数全微分 公式 格式 Dz=diff(z,x)*’dx’+diff(z,y)*’dy’ syms x y z=(x^2+y^2)*sin(x*y); dz=diff(z,x)*’dx’+diff(z,y)*’dy’ 结果： dz=(2*x*sin(x*y)+(x^2+y^2)*cos(x*y)*y)*dx+(2*y*sin(x*y)+(x^2+y^2)*cos(x*y)*x)*dy 重