中央财经大学MATLAB第4讲课余作业.docVIP

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

《第4讲 符号计算》课余作业 1、创建符号表达式f(x)=sin(x)+cos(x)-tan(x) syms x f= sin(x)+cos(x)-tan(x) f = cos(x) + sin(x) - tan(x) 2、计算习题1中的符号表达式在x=0, x=1, x=2*pi处的值 syms x f=sym(sin(x)+cos(x)-tan(x)); subs(f,x,[0 1 2*pi]) ans = 1.0000 -0.1756 1.0000 3、计算符号表达式f(x)=x+cos(x)-sin(x)在x=pi-1处的值,并将结果设置为以下5中精度,即分别为小数点后1位、2位、10位、20位和50位有效数字。 syms x f=x+cos(x)-sin(x); a=subs(f,x,pi-1); a1=vpa(a,1) % 在Matlab2009a版本中,此程序无效;在Matlab6.5中可用 fprintf(a1=\n%3.1f\n,a) % 输出一位有效数字 a2=vpa(a,2) a10=vpa(a,10) a20=vpa(a,20) a50=vpa(a,50) a1 = 0.75981936291375673508952104384662 a1= 0.8 a2 = 0.76 a10 = 0.7598193629 a20 = 0.75981936291375673509 a50 = 0.75981936291375673508952104384661652147769927978516 4、设x为符号变量,f(x)=x^4+2*x^2+1, g(x)=x^3+6*x^2+3*x+5,试进行如下运算: (1)f(x)+g(x) (2)f(x)*g(x) (3)对f(x)进行因式分解 (4)求g(x)的反函数 syms x f=x^4+2*x^2+1; g=x^3+6*x^2+3*x+5; f+g f*g horner(f) finverse(g) ans = x^4 + x^3 + 8*x^2 + 3*x + 6 ans = (x^4 + 2*x^2 + 1)*(x^3 + 6*x^2 + 3*x + 5) ans = x^2*(x^2 + 2) + 1 Warning: finverse(x^3 + 6*x^2 + 3*x + 5) is not unique. In sym.finverse at 46 ans = 3/(x/2 + ((x/2 - 15/2)^2 - 27)^(1/2) - 15/2)^(1/3) + (x/2 + ((x/2 - 15/2)^2 - 27)^(1/2) - 15/2)^(1/3) - 2 5、设x为符号变量,令y=x^4+2*x^2+1,z=y^2,w=sin(z),试求w关于x的符号表达式,并求当x=5时,w的小数点之后具有1位、2位、10位、20位和50位有效数字的数值解。 syms x y=x^4+2*x^2+1; z=y^2; w=sin(z) w0=subs(w,x,5) fprintf(w1=\n%3.1f\n,w0) vpa(w0,2) vpa(w0,10) vpa(w0,20) vpa(w0,50) w = sin((x^4 + 2*x^2 + 1)^2) w0 = -0.0673 w1= -0.1 ans = -0.067 ans = -0.06734017262 ans = -0.067340172619870686255 ans = -0.067340172619870686254728298081317916512489318847656 6、使用sym函数生成如下符号矩阵, a=sym([1/x,1/(x+1);1/(x+2) 1/(x+3)]), b=sym([x, 1; x+2,0]),试对着两个符号矩阵进行如下操作: (1)a-b (2)a*b (3)求a的倒置 (4)求a的行列式 (5)求b的逆 (6)求a的秩 (7)求a^3 a=sym([1/x,1/(x+1);1/(x+2) 1/(x+3)]); b=sym([x, 1; x+2,0]); a-b a*b a % 求a的倒置 det(a) % 求a的行列式 inv(b) % 求b的逆 rank(a) % 求a的秩 a^3 ans = [ 1/x - x, 1/(x + 1) - 1] [ 1/(x + 2) - x - 2, 1/(x + 3)] ans = [ (x + 2)/(x + 1) + 1, 1

文档评论(0)

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

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

1亿VIP精品文档

相关文档