MATLAB - 04 - 符号计算.pptVIP

2、线性代数方程组的符号解 3、常微分方程（组）的求解 Matlab中，常微分方程： y’’+2y’=x 表示成： D2y + 2Dy = x D表示对变量求导数 Dy表示对y求1阶导数 Dny表示对y求n阶导数 dsolve * Matlab 演示 * Matlab 演示 * 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] 分式通分：numden [N,D]=numden(f): N 为通分后的分子，D 为通分后的分母 syms x y; f=x/y+y/x; [N,D]=numden(f) [n,d]=numden(sym(112/1024)) 嵌套形式的多项式：horner horner 多项式：嵌套形式的多项式 syms x; f=x^4+2*x^3+4*x^2+x+1; g=horner(f) 例： 4、subexpr命令 subexpr命令可以把复杂的表达式中所含的多个相同的子表达式用一个符号代替，使其表达更简洁，用法是： [RS,ssub]=subexpr(S,ssub) 其中： S是原表达式 ssub是符号变量 RS是结果表达是 注意：被代替的子表达式是matlab自动寻找的，并不肯定会得到简洁的结果，只有比较长的子表达式才会被替代。 5、pretty命令 pretty(s)：将s转变为书写格式，转换过程中，matlab不会对s进行化简 6、符号数值精度控制 double(x):把符号常数转化为16位相对精度的浮点数值对象。 digits(n)：设置数值计算以n位相对精度进行 xs＝vpa(x)：在digits指定的精度下，给出x的数值型符号结果xs xs=vpa(x,n)：在n位相对精度下，给出x的数值型符号结果xs 三、常用符号运算功能的实现 1、极限 2、微分 3、积分 4、级数求和 5、一元函数的泰勒级数展开 6、反函数 7、复合函数 1、极限 F：符号表达式 x：变量 默认变量：没有指定变量时，由findsym确定 limit(F,x,a) ：求x →a时F的极限 limit(F,a)：求默认变量→a时F的极限 limit(F)：求默认变量→0时F的极限 limit(F,x,a,’right’)：求x →a+0时F的极限（右极限） limit(F,x,a,’left’)：求x →a-0时F的极限（左极限） 例：求下列各函数式的极限 2、微分 diff(f,v,n) 求f对变量v的n阶导数 diff(f,n) 求f对默认变量的n阶导数 diff(f,v) 求f对变量v的导数 diff(f) 求f对默认变量的导数 例： diff命令还可以求插分 diff(X) 若X为一矢量则diff(X)的结果为插分，即： [X(2)-X(1), X(3)-X(2), ……, X(n)-X(n-1)] 若X为一矩阵则diff(X)的结果式对矩阵X各列求插分所得到的矩阵。 3、积分 求不定积分 int(f,v) inf(f) 求定积分与广义积分 int(f,v,a,b) int(f,a,b) a是积分下限，b是积分上限,当a或b取inf（或-inf）时，计算的就是广义积分。 例： 例：求下面的积分，给出50位精度的数值 4、级数求和 symsum(s,k)： 求k从0开始 到k-1为止s的前k项和 symsum(s) 功能同上，只不过变量为默认变量 symsum(s,v,a,b) 求自变量v从a变化到b时s的和 symsum(s,a,b) 功能同上，只不过变量为默认变量 例：求下列级数的和 5、一元函数的泰勒级数展开 taylor(f) 求f关于默认变量的5阶麦克劳林展开 taylor(f,n) 求f关于默认变量的n-1阶麦克劳林展开 taylor(f,n,v) 求f关于变量v的n-1阶麦克劳林展开 taylor(f,n,v,a) 求f在v=a处的n-1阶泰勒展开式 求泰勒展开式 多元函数的泰勒级数展开 mtaylor(f,v) mtaylor(f,v,n) mtaylor(f,v,n,w) f：待完全展开的代数表达式