东华大学matlab数学实验-第四次答案.doc

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
东华大学matlab数学实验-第四次答案

1 p=[1 1 1];x=roots(p) x = -0.5000 + 0.8660i -0.5000 - 0.8660i polyval(p,x) ans = 1.0e-15 * 0.3331 0.3331 clear; p=[3 0 -4 2 -1];x=roots(p) x = -1.4058 + 0.0000i 1.0000 + 0.0000i 0.2029 + 0.4427i 0.2029 - 0.4427i polyval(p,x) ans = 1.0e-13 * -0.2498 + 0.0000i -0.0266 + 0.0000i 0.0000 + 0.0072i 0.0000 - 0.0072i clear; p=[5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6 8 0 0 0 -5 0 0];x=roots(p) x = 0.0000 + 0.0000i 0.0000 + 0.0000i 0.9768 + 0.0000i 0.9388 + 0.2682i 0.9388 - 0.2682i 0.8554 + 0.5363i 0.8554 - 0.5363i 0.6615 + 0.8064i 0.6615 - 0.8064i 0.3516 + 0.9878i 0.3516 - 0.9878i -0.0345 + 1.0150i -0.0345 - 1.0150i -0.4609 + 0.9458i -0.4609 - 0.9458i -0.1150 + 0.8340i -0.1150 - 0.8340i -0.7821 + 0.7376i -0.7821 - 0.7376i -0.9859 + 0.4106i -0.9859 - 0.4106i -1.0416 + 0.0000i -0.7927 + 0.0000i polyval(p,x) ans = 1.0e-12 * 0.0000 + 0.0000i 0.0000 + 0.0000i 0.1305 + 0.0000i 0.0954 + 0.0245i 0.0954 - 0.0245i 0.0796 + 0.0419i 0.0796 - 0.0419i -0.0643 + 0.0680i -0.0643 - 0.0680i 0.1620 + 0.0550i 0.1620 - 0.0550i -0.0120 + 0.0557i -0.0120 - 0.0557i -0.0978 + 0.3836i -0.0978 - 0.3836i 0.0012 + 0.0121i 0.0012 - 0.0121i 0.0560 - 0.1047i 0.0560 + 0.1047i 0.2392 + 0.0861i 0.2392 - 0.0861i -0.0636 + 0.0000i 0.0011 + 0.0000i p1=[2 3]; p2=conv(p1,p1); p=conv(p2,p1) p = 8 36 54 27 p3=[8 36 54 23];x=roots(p3) x = -1.8969 + 0.6874i -1.8969 - 0.6874i -0.7063 + 0.0000i polyval(p3,x) ans = 1.0e-13 * 0.1066 + 0.0799i 0.1066 - 0.0799i 0.0355 + 0.0000i 2 clear; fun=@ (x)(x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-x/2); fplot(fun,[0,10]) fzero(fun,2) ans = 2.1155 3 fun=@ (x)[7/10*sin(x(1))+2/10*cos(x(2))-x(1),7/10*cos(x(1))-2/10*sin(x(2))-x(2)]; [x,f,h]=fsolve(fun,[0,0]) Equation solved. fsolve completed because the vector of function values is near zero as measured by the default value of the fun

文档评论(0)

shuwkb + 关注
实名认证
内容提供者

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

1亿VIP精品文档

相关文档