第二章 matlab基础2.pptVIP

* 1、符号级数求和 MATLAB中利用symsum函数用于对通项为符号表达式s的级数求部分和或和。该函数的一般引用格式为： Symsum(s, n, a, b) 级数通项的符号表达式 从第a项开始求和 求和至第b项，如果求所有项的和，b可取inf 三、符号级数 求和变量单变量时可缺省 * 例：求级数 三、符号级数 syms k n symsum(1/k^2,1,10) ans = 1968329/1270080 ans1=double(ans) ans1=1.5498 symsum(1/k^2,1,inf) ans = 1/6*pi^2 * 泰勒级数 三、符号级数 MATLAB中提供了将函数展开为幂级数的函数taylor，其调用格式为： taylor(f, x, a, k) 展开函数的符号表达式 在x=a处展开，缺省时a=0 将f展开到k-1阶，缺省时k=6 自变量可缺省 * 例 将多项式p(x)=1+3x+5x2-2x3表示成x+1的多项式。 三、符号级数 x=sym(x); p=1+3*x+5*x^2-2*x^3; f=taylor(p,x,-1,4) f = -8-13*x+11*(1+x)^2-2*(1+x)^3 ans= 1-x+1/2*x^2-1/6*x^3+1/24*x^4-1/120*x^5 例 x=sym(x); taylor(exp(-x)) * MATLAB中展开泰勒级数的可视化工具——taylortool 三、符号级数 在命令窗口中输入taylortool，回车即可。 * 四、代数方程的符号求解 在MATLAB中，求解用符号表达式表示的代数方程可由函数solve实现，其调用格式为： solve(s1, s2, …, sn, v1, v2 ,…, vn) n个方程构成的方程组，方程都整理成等号右端为0的形式 n个未知数 注：当求解方程组时，solve命令将方程组的解以结构体的形式存储。 * 例 求解代数方程：ax2-bx-6=0 四、代数方程的符号求解 syms a b x solve(a*x^2-b*x-6) ans = [ 1/2/a*(b+(b^2+24*a)^(1/2))] [ 1/2/a*(b-(b^2+24*a)^(1/2))] * 四、代数方程的符号求解 注：当使用solve命令求解方程组时，解将存储在返回的数组或结构体中，如果用数组存储，所有的解必须按照字典顺序排列，否则会返回错误结果。 ， [x,y,z,w]=solve(x+y+z+w-10,2*x+y+z+w-1,x+2*y+z+w-12, ‘x+y+z+2*w-14’) 例如： MATLAB解得： x =4，y =-9，z =2，w =13， 应改为： [w,x,y,z]=solve(x+y+z+w-10,2*x+y+z+w-1,x+2*y+z+w-12, ‘x+y+z+2*w-14’) MATLAB解得： w =4，x =-9，y =2，z =13, 这才是正确解 显然错误 * 例 求方程组 四、代数方程的符号求解 syms u y v z w; S=solve(u*y^2+v*z+w=0,y+z+w=0,y,z); disp(S.y),disp(S.y),disp(S.z),disp(S.z) ， 关于y，z的解 S.y [ -1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))-w] [ -1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))-w] S.z [ 1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))] [ 1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))] * 1、替代 五、替代与转换 符号计算的结果经常需要进行变量代换、参数赋值等操作，MATLAB提供了用于替代的命令subs，具体格式为： subs(ES, old, new) 被执行替代操作的符号表达式 将ES中的old全部替换为new ES中被替换的符号变量old，如缺省按默认变量执行 注：如果一次执行多对替代，可分别将new和old存储在两个大小相同的数组中，按照对应位置进行替代。 * 例： 五、替代与转换 syms a x; f=a*sin(x)+5; f1=subs(f, ’sin(x)’,sym(y)) f1=a*y+5 f2=subs(f, {a,x},{2, sym(pi/3)})