Matlab_ch8_符号运算.ppt

8. Matlab的符号运算 调用格式： s = solve( f1 , f2 , ... , fN , v1 , v2 , ... , vN) 求解由 f1 , f2 , ... , fN 确定的方程组关于 v1 , v2 , ... , vN 的解 s 是一个架构数组，如果要显示(引用)每个解，使用 s.v1，s.v2 若没有提供v1 , v2 , ... , vN ，Matlab会自动确定默认变量。 g1=sym(x+y+z=1) , g2=sym(x-y+z=2), g3=sym(2*x-y-z=1) [x y z]=solve(g1,g2,g3) z = 5/6, y = -1/2, x = 2/3 y = dsolve(eq1,eq2, ... ,cond1,cond2, ... ,v) y ：输出； eq1、eq2、...：微分方程； cond1、cond2、...：初值条件； v ：自变量。 只有很少一部分微分方程（组）能求出解析解。大部分微分方程（组）只能利用数值方法求数值解。 如果省略初值条件，则表示求通解； 如果省略自变量，则默认自变量为 t dsolve(Dy=2*x,x); dy/dx = 2x dsolve(Dy=2*x); dy/dt = 2x 若找不到解析解，则返回其积分形式。 微分方程中用 D 表示对 自变量 的导数，如： Dy y；D2y y；D3y y 或 或 或 y的一阶导数—— Dy y的二阶导数—— D2y y的 n 阶导数—— Dny [y1,y2…]=dsolve(x1,x2,…xn) —— 返回微分方程的解 eq1 = Dy+2*x*y=x*exp(-x^2) y = dsolve(eq1,x) eq1 = x*Dy+y-exp(x)=0 cond1 = y(1)=2*exp(1) y = dsolve(eq1,cond1,x) y=dsolve(D2y+2*Dy+2*y=0,y(0)=1,Dy(0)=0) ans = exp(-x)*cos(x)+exp(-x)*sin(x) ezplot(y) —— 方程解y(t)的时间曲线图 求该方程的解 roots(p) solve(p) p=[1 -4 9 -10] r=roots(p) 或 s1=sym(x^3-4*x^2+9*x-10); solve(s1) 例: 求方程 的所有根 [x,y]= solve(x+y=1,x-11*y=5,x,y) 求极限 syms x; y1=(1+4*x)^(1/x); y2=(exp(x)-1) /x; limit(y1,x,0) limit(y2,x,0) syms x; f=3*x^3+5*x+1; diff(f,1) syms x y; z=x^2*sin(2*y); B=diff(z,x) syms x ; y=(x^2+sin(x))/(1+x^2); int(y,x,-1,1) syms x y ; int(int(x^2+y^2+1,y,x,x+1),x,0,1) dsolve(Dy=1+y^2, x) dsolve(Dy=1+y^2, y(0)=1 ,x) dsolve(x^2*D2y+x*Dy+(x^2-1/2)*y=0,y(pi/2)=2,Dy(pi/2)=-2/pi,x) * syms x a t h; limit(sin(x)/x) = 1 limit(1/x,x,0,right) = inf limit(1/x,x,0,left) = -inf limit((sin(x+h)-sin(x))/h,h,0) = cos(x) v = [(1 + a/x)^x, exp(-x)]; limit(v,x,inf,left) = [exp(a), 0] MATLAB 程序设计 符号运算的基本知识 符号运算 在进行符号运算时，必须先定义基本的符号量，可以是符号常量、符号变量。 含有符号对象的表达式称为符号表达式。 符号矩阵/数组：元素为符号表达式的矩阵/数组。 Matlab 的符号数学工具箱可以完成几乎所有的符号运算功能。主要