5 符号运算解析.pptVIP

f3=[x,y^2;x*sin(y),exp(-x)*log(y)] dfdx=diff(f3) dfdx = [ 1, 0] [ sin(y), -exp(-x)*log(y)] dfdx2=diff(f3,x,2) dfdx2 = [ 0, 0] [ 0, exp(-x)*log(y)] dfdxdy=diff(diff(f3,x),y) dfdxdy = [ 0, 0] [ cos(y), -exp(-x)/y] * 【例5-20】 使用diff指令求函数的导数。 【例5-21】 求多元向量函数的Jacobian矩阵。 syms x y; f=[x^2+y^2;sin(x)*log(y);x*exp(-y)]; jacobian(f,[x y]) ans = [ 2*x, 2*y] [ cos(x)*log(y), sin(x)/y] [ exp(-y), -x*exp(-y)] * 3. 符号积分运算 符号积分指令的调用格式为 int(f,?v)和int(f,?v,?a,?b) int(f,?v)用来计算函数f对变量v的不定积分， int(f,?v,?a,?b)用来计算f对变量v的定积分，其中a、b分别为积分的下、上限。 注：和符号微分一样，如果f是矩阵，则积分运算将对矩阵中的元素逐个进行。 * 【例5-22】 使用int指令求函数的积分。 syms x u t; int((x+sin(x))/(1+cos(x))) %求不定积分 ans = x*tan(1/2*x) int(cos(x)*exp(x),0,1) %求定积分 ans = 1/2*cos(1)*exp(1)+1/2*exp(1)*sin(1)-1/2 int(u*sin(x),x,1,sin(t)) ans = -u*cos(sin(t))+u*cos(1) * 【例5-22】 使用int指令求函数的积分。 f=[u*x^2,t*sin(x);2*x,x*exp(x)]; int(f) ans = [ 1/3*u*x^3, -t*cos(x)] [ x^2, x*exp(x)-exp(x)] * 5.3 符号方程的求解 1. 符号代数方程组的求解 ① 线性方程组的解 矩阵计算是求解线性方程组的一种简便有效 的方法。在MATLAB中，不管数据对象是数值 还是符号，实现矩阵运算的指令形式几乎完全相 同。因此，对于线性方程组符号解的问题，可以 套用求数值解的指令编写方法进行求解。 * 【例5-28】 求线性方程组 的解。 为便于观察，先把方程组写成矩阵形式。 简记为AX=B的形式。 * A=sym([1 1 -1 -1;1 1/2 1/3 1/4;1/4 -1 1 0;1/2 0 -1 1]); B=sym([0;12;6;-6]); X=A\B %矩阵B除以矩阵A X = 8 2 6 4 即方程组的解为a=8，b=2，c=6，d=4。 * ② 非线性方程组的解 对于非线性方程组，通常使用solve指令来进 行求解，一般调用格式如下 S=solve(e1,?e2,?…,?en,?s1,?s2,?…sn) 在该指令中，‘e1’,‘e2’,…,‘en’是字符串表达的 方程组或字符串表达式，‘s1’,?‘s2’,?…‘sn’是字符 串表达的求解变量名，S是一个构架数组，如果 要显示求解结果，需要采用S.s1，S.s2,…， S.sn的援引方式。 * 【例5-29】 求方程组 关于x、y的解。 S=solve(a^2*x+b*y+c=0,x+y+a+b+c=0,x,y)； S.x % 显示x的解 ans = (c*b-c+a*b+b^2)/(-b+a^2) S.y % 显示y的解 ans = -(a^2*c-c+a^3+a^2*b)/(-b+a^2) * 2 . 符号微分方程(组)的求解 求解符号微分方程(组)最常用的调用格式如下 S=dsolve(e1,?e2,?...,?c1,?c2...,?v) 指令中 e表示微分方程， c表示初始条件或边界条件， v表示指定的独立变量。 其中微分方程是必不可少的输入内容，初始条件和独立变量是否省略可视需要而定。 注意：如果没有对独立变量进行定义，则默认小写字母t为独立变量。 * 微分方程和初始