- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第七章 多项式
第一节 多项式的创建
1、多项式由其系数向量唯一确定:
p=[1 2 3 4 5] 主要内容:
主要内容:
y=poly2sym(p) 第一节 多项式的创建
第一节 多项式的创建
disp(y) 第二节 多项式的运算
第二节 多项式的运算
2、求矩阵A的特征多项式:
a=[1 2 3;3 2 4;5 7 1] %方正
poly(a) %求a的特征多项式
poly(a(1,:))
%x^3-6x^2+11x-6(1,2,3为方程的根)
钟培华2006 -7
3、计算以p为系数的矩阵在X处的值:
Y=polyval(p,X)%当X是矩阵或者向量,将对每个
元素进行运算。
p=[1 0 0 1];a=1,b=[1 2];c=[1
2;1 2;1 3];
pola=polyval(p,a)%x^3+1在a的值
polb=polyval(p,b)%在b的各元素
的值
polb=polyval(p,c)
4、Y=polyvalm(p,x)%求以p为系数的多项矩阵:
Y=p(1)x^n+…+p(n-1)x+p(0)E,x为多项式。
钟培华2006 -7
5、求多项式x^4+3x^2+12x-7的根:
p=[1 0 3 12 -7]
roots(p) %求根
6、多项式的加减法对向量加减即可。
7、以a,b为系数的多项式的积:
a=[1 2 3 4];b=[4 5 6 7];
C=conv(a,b)
cx=poly2sym(C)
钟培华2006 -7
7.反褶积(解卷)和多项式除法运算
函数 deconv
格式 [q,r] = deconv(v,u) %多项式v除以多项
式u,返回商多项式q和余多项式r。
注意:v、u、q、r都是按降幂排列的多项式系数
向量。
3 2 2
(x +2x +3x +4)(10x +20x +30) 的卷积为
u = [1 2 3 4]
v = [10 20 30]
c = conv(u,v)
%c=[10 40 100 160 170 120]
钟培华2006 -7
则反褶积为
[q,r] = deconv(c,u)%q =[10 20 30]r
=[0 0 0 0 0]
8、求导与求积:
Polyder(p)%求以p为系数的矩阵的倒数(的系数)
Polyder(a,b)%求多项式的积的导数。
p=[3 1 8 8]; a=[9 4 9 4];
b=[8 4 6 7]
c=polyder(a)
w=polyder(a,b)
polyint(a)
钟培华2006 -7
9、[q,d]=Polyder(a,b)%求多项式的商的导数,
并以q/d的形式返回。
Polyint(p,k)%返回多项式p的积分,并以k为常
数项。
polyder([1 0 1],[1 0])%(x^2+1)/x
的倒数
%q=[1 0 -1],d=[1 0 0]
polyint([1 1 2 4 5],3)
ans =
0.2000 0.2500 0.6667 2.0000 5.0000
3.0000
钟培华2006 -7
文档评论(0)