MATLAB课件CSY..ppt

例6.9 计算3阶范得蒙德矩阵行列式的值 syms a b c; U=[a b c]; A=[[1 1 1];U;U.^2];% det(A) %... 例6.10 建立x、y的一般二元函数 syms x y; f=sym(‘f(x,y)’) 6.2.2 化简 1.因式分解 factor(s) 例6.11 对表达式f=x12-1进行因式分解 syms x; f=x^12-1; factor(f) 6.2.2 化简 2.展开 expand(s) 例6.12 展开(x-1)12 syms x; f=sym(‘(x-1)^12’); expand(f) 6.2.2 化简 3.同类项合并 collect(s,n) n? 例6.13对于表达式f=x(x(x-6)+12)t, 分别将自变量x和t的同类项合并 syms x t;f=x*(x*(x-6)+12)*t; collect(f) ％缺省按x合并同类项 collect(f,t) ％按t合并同类项 6.2.2 化简 4.化简 simplify(s) 例6.14对表达式(x-1)3/(x-1)进行化简 syms x; f=sym(‘(x-1)^3/(x-1); simplify(f) simple(s) 函数simple通过对表达式尝试不同的算法进行化简，以寻求符号表达式S的最简形式。 simple函数综合使用了下列化简方法： simplify radsimp combine collect factor convert 6.2.2 化简 5.分式通分 [n,d]=numden(s) 例6.15对表达式y=x/y+y/x,进行通分 syms x y; f=x/y+y/x; [n,d]=numden(f) §6.2 化简与替换 6.2.3 符号表达式的替换 [Y,SIGMA]=subexpr(S,SIGMA) R=subs(S,old,nem) 目的：通过符号替换使表达式的输出形式 简化，以得到一个简单的表达式。 例6.16求解并化简三次方程 x3+ax+1=0解 t=solve(x^3+a*x+1=0) [r,s]=subexpr(t,s) 例6.17分别用新变量替代表达式a+b和cosa+sinb中的变量 syms a b subs(a+b,a,4) subs(cos(a)+sin(b),{a,b},{sym(alpha),2}) §6.3 符号微积分 6. 3.1 符号积分 limit(F,x,a) limit(F, a) limit(F) limit(F,x,a,’right’) limit(F,x,a,’left’) §6.3 符号微积分 6. 3.2 符号微分 diff(S) diff(S, v) diff(S,n) 例6.20分别求表达式f=xx的导数和三次导数 syms x; f=x^x; diff(f) diff(f,3) 例6.21 求下列函数的导数 求 求 syms x y; f=sqrt(1+exp(x)); diff(f) %…… syms x y; f=x*cos(x); diff(f,x,2) %…… diff(f,x,3) %…… 例6.22 求曲线y=x3＋3x－2上与直线y=4x-1平行的切线的切点 x=sym(‘x’); y=x^3+3*x-2; %定义曲线 f=diff(y); %对曲线求导 g=f-4; solve(g) %求方程的根，即求曲线何处 的导数为4 §6.3 符号微积分 6. 3.3 符号积分 int(S) int(S,v) int(S,a,b) int(S,v,a,b) 例6.23 求下列不定积分 （1） （2） sym(‘x’); f=(3-x^2)^3; int(f) %求不定积分 sym(‘x’); f=sqrt(x^3+x^4); g=int(f) %求不定积分 simple(g) %化简结果 例6.24 求 的不定积分。 int(‘log(x)/(exp(x^2)’) Waring:Explicit integral could not be found. In E:\MATLAB\toolbox\symbolic\…… In E:\MATLAB\toolbox\symbolic\…… ans= int(log(x)/exp(x^2),x) 说明：当找不到原函数时，返回未经计算的函数！ §6.3 符号微积分 6. 3.4 积分变换 1 概念 通过积分运算，把一个函数f变成另外一个函数