Matlab4-符号计算概述.ppt

ezsurf(x,y,z) ezsurfh(x,y,z,[smin,smax,tmin,tmax]) ezsurf(x,y,z,[a,b]) x = cos(s)*cos(t) y = cos(s)*sin(t) z = sin(s) ezsurf(x,y,z,[0,pi/2,0,3*pi/2]) （2）曲面由参数方程x=x(s,t),y=y(s,t),z=z(s,t)确定 * Matlab 演示 * Matlab 演示 5、一元函数的泰勒级数展开 taylor(f) taylor(f,n) taylor(f,n,v) taylor(f,n,v,a) 求f关于默认变量的5阶麦克劳林展开 求f关于默认变量的n-1阶麦克劳林展开 求f关于变量v的n-1阶麦克劳林展开 求f在v=a处的n-1阶泰勒展开式 求函数y=cosx在x=0点处的5阶taylor展开式 及在x=pi/3处的6阶taylor展开式 。 syms x y1=taylor(cos(x)) y2=taylor(cos(x),7,pi/3) 6、复合函数 7、反函数 求反函数 g = finverse(f) g = finverse(f,v) 四、求解符号方程 1、线性方程组的符号解 2、符号(代数)方程求解 3、常微分方程（组）的求解 1、线性方程组的符号解 clc; clear all; A = sym([1,1,1;3,-1,6;0,1,3]) b = sym([1;7;4]) x = A\b 2、符号(代数)方程求解 s = solve(f,v)：求方程关于指定自变量的解 s = solve(f)：求方程关于默认自变量的解 f 是用字符串表示的方程，或符号表达式。 eq1 = a*x^2+b*x+c=d sx = solve(eq1) sa = solve(eq1,a) eq2 = a*x^2+b*x+c-d sx = solve(eq2) sa = solve(eq2,a) A、直接用 字符串 创建方程 若方程 f 中不含等号，则表示解方程 f = 0 syms a b c d x eq1 = a*x^2+b*x+c-d; sx = solve(eq1) sa = solve(eq1,a) B、先定义符号变量，再创建符号表达式 符号表达式不能写成： eq1 = a*x^2+b*x+c=d 解方程 syms x; f = x^3-3*x+1; s = solve(f,x) s = solve(x^3-3*x+1,x) s = solve(x^3-3*x+1=0,x) s = solve( f1 , f2 , ... , fN , v1 , v2 , ... , vN) 求解由 f1 , f2 , ... , fN 确定的方程组关于 v1 , v2 , ... , vN 的解 solve 也可以用来解方程组 s 是一个架构数组，如果要显示(引用)每个解，使用 s.v1，s.v2 若没有提供v1 , v2 , ... , vN ，Matlab会自动确定默认变量。 例：解方程组 eq1 = x+2*y-z=27’ eq2 = x+z=3 eq3 = x^2+3*y^2=28’ s = solve(eq1,eq2,eq3) s.x s.y s.z 解方程组 clc； clear all; syms u v w y z eq1 = u*y^2 + v*z + w; eq2 = y + z + w; s = solve(eq1,eq2,y,z) 范德瓦尔斯气体的三个临界参量 clc; clear all; f = (P+a/V^2)*(V-b)=R*T D1 = diff(f,V); D2 = diff(D1,V); S = solve(f,D1,D2, P,T,V); Pc = S.P Vc = S.V Tc = S.T clc; clear all; syms a b R T V P = R*T/(V-b) - a/V^2 D1 = diff(P,V) D2 = diff(D1,V) S = solve(D1,D2,V,T) Vc = S.V Tc = S.T Pc = subs(P,[V,T],[S.V,S.T]) s = solve((x+2)^x=2) solve 在得不到解析解时，会给出数值解 3、常微分方程（组）的求解 y = dsolve(eq1,eq2, ... ,cond1,cond2, ... ,v) y 为输出， eq1、eq2、...为微分方程， cond1、cond2、...为初值条件， v 为自变量。 只有很少一部分微分方程（组）能求出解析解。大部分微分方程（组）只能利用数值方法求数值解。 或 或 或