Matlab 符号计算.pptVIP

  1. 1、本文档共61页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Matlab 符号计算.ppt

syms a b c d r = [a,b,c,d] A = [r.^0;r;r.^2;r.^4] B = det(A) C = factor(B) 2、expand expand(S) 把 S 展开为多项式形式,也可展开复合的三角函数、指数函数 syms x expand((x-2)*(x-4)) ans = x^2-6*x+8 syms x y; expand(cos(x+y)) ans = cos(x)*cos(y)-sin(x)*sin(y) syms a b expand(exp((a + b)^2)) ans = exp(a^2)*exp(a*b)^2*exp(b^2) 3、collect collect(S) 按 默认变量 对 S 合并同类项 collect(S,v) 按指定变量 v 对 S 合并同类项 syms x y; f = (x+y)*(x^2+y^2+1) f1 = collect(f) f2 = collect(f,y) f1 = x^3+y*x^2+(y^2+1)*x+y*(y^2+1) f2 = y^3+x*y^2+(x^2+1)*y+x*(x^2+1) 4、simplify simplify(S): 用内置的化简规则对 S 进行化简 syms x; f = sin(x)^2 + cos(x)^2 ; f1 = simplify(f) syms a b c; g = exp(c*log(sqrt(a+b))); g1 = simplify(g) f1 = 1 g1 = (a+b)^(1/2*c) 5、simple r = simple(S) 对 S 尝试多种不同的算法(包括simplify)进行简化,但不显示相应的结果,只是返回其中最简短的形式并赋值给变量 r [r, How] = simple(S) 对 S 尝试多种不同的算法(包括simplify)进行简化,但不显示相应的结果,只是返回其中最简短的形式并赋值给变量 r,同时返回字符串How ,记录得到 r 所使用的化简方法。 simple(S) 对 S 尝试多种不同的算法(包括simplify)进行简化,并显示相应的结果,最后返回其中最简短的形式并赋值给默认变量ans syms x; f = (1/x^3+6/x^2+12/x+8)^(1/3); f1 = simple(f) f2 = simple(f1) f3 = simple(f2) 多次使用 simple ,直到结果与前一次相同,以达到最简单的表示形式。 f1 = (2*x+1)/x f2 = 2+1/x f3 = 2+1/x 6、horner horner(P) 将多项式 P 改写为horner形式(这种形式可以快速计算多项式的值) syms x; f = x^4+2*x^3+4*x^2+x+1; g = horner(f) g = 1+(1+(4+(x+2)*x)*x)*x 7、numden [N,D] = numden(f) 返回值: N 为 f 的分子,D 为 f 的分母 syms x y; f = x/y + y/x; [N,D] = numden(f) N = x^2+y^2 D = y*x 四、求解符号方程 1、线性方程组的符号解 2、符号方程(组)求解 3、常微分方程(组)的求解 1、线性方程组的符号解 clc; clear all; A = sym([1,1,1;3,-1,6;0,1,3]) b = sym([a;7;4]) s = A\b s = -14/15+3/5*a -3/5+3/5*a 23/15-1/5*a 2、符号方程(组)求解 s = solve(eq,var) 求方程指定变量var的解 s = solve(eq) 求方程默认自变量的解 eq 是用字符串表示的方程,或符号表达式。 (1)用 字符串 创建方程 eq1 = a*x^2+b*x+c=d sx = solve(eq1) sa = solve(eq1,a) eq1 = a*x^2+b*x+c-d sx = solve(eq1) sa = solve(eq1,a) 若方程 eq 中不含等号,则表示解方程 eq = 0 sx = -1/2*(b-(b^2-4*a*c+4*a*d)^(1/2))/a -1/2*(b+(b^2-4*a*c+4*a*d)^(1/2))/a sx是2行1列的符号数组 sa = -(b*x+c-d)/x^2 以a为方程的变量求出的解 syms a b c d x eq1 = a*x^2+b*x+c-d; sx = solve(eq1) sa = solve(eq1,

文档评论(0)

fa159yd + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档