原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
matlab语言把多项式表达成一个行向量，该向量中的元素是按多项式降幂排列的。f(x)=anxn+an-1xn-1+……+loa0可用行向量p=[anan-1……a1+a0]表示poly——产生特征多项式系数向量特征多项式一定是n+1维的特征多项式第一个元素一定是1四、多项式运算第29页，共56页，星期日，2025年，2月5日例:a=[123;456;780];p=poly(a)p=1.00-6.00-72.00-27.00p是多项式p(x)=x3-6x2-72x-27的matlab描述方法，我们可用：p1=poly2str(p,‘x’)—函数文件，显示数学多项式的形式p1=x^3-6x^2-72x-27第30页，共56页，星期日，2025年，2月5日2.roots——求多项式的根a=[123;456;780];p=poly(a)p=1.00-6.00-72.00-27.00r=roots(p)r=12.12-5.73——显然r是矩阵a的特征值-0.39第31页，共56页，星期日，2025年，2月5日当然我们可用poly令其返回多项式形式p2=poly(r)p2=1.00-6.00-72.00-27.00matlab规定多项式系数向量用行向量表示，一组根用列向量表示。第32页，共56页，星期日，2025年，2月5日3.conv,convs多项式乘运算例:a(x)=x2+2x+3;b(x)=4x2+5x+6;c=(x2+2x+3)(4x2+5x+6)a=[123];b=[456];c=conv(a,b)c=conv([123],[456])c=4.0013.0028.0027.0018.00p=poly2str(c,x)p=4x^4+13x^3+28x^2+27x+18第33页，共56页，星期日，2025年，2月5日4.deconv多项式除运算a=[123];c=[4.0013.0028.0027.0018.00]d=deconv(c,a)d=4.005.006.00[d,r]=deconv(c,a)余数c除a后的整数第34页，共56页，星期日，2025年，2月5日5.多项式微分matlab提供了polyder函数多项式的微分。命令格式：polyder(p):求p的微分polyder(a,b):求多项式a,b乘积的微分[p,q]=polyder(a,b):求多项式a,b商的微分例：a=[12345];poly2str(a,x)ans=x^4+2x^3+3x^2+4x+5b=polyder(a)b=4664poly2str(b,x)ans=4x^3+6x^2+6x+4第35页，共56页，星期日，2025年，2月5日五、代数方程组求解matlab中有两种除运算左除和右除。对于方程ax=b，a为an×m矩阵，有三种情况：?当n=m时，此方程成为“恰定”方程?当nm时，此方程成为“超定”方程?当nm时，此方程成为“欠定”方程matlab定义的除运算可以很方便地解上述三种方程第36页，共56页，星期日，2025年，2月5日1.恰定方程组的解方程ax=b(a为非奇异)x=a-1b矩阵逆两种解:x=inv(a)?b—采用求逆运算解方程x=a\b—采用左除运算解方程第37页，共56页，星期日，2025年，2月5日方程ax=ba=[12;23];b=[8;13];?x=inv(a)*b?x=a\bx=x=2.002.003.003.00=ax=
文档评论（0）