Xiaoxing He分享 http://blog.sciencenet.cn/u/xiaoxinghe GNSS Time series analysis, Wireless positioning, 3D laser scan, etc

博文

多项式操作

已有 2540 次阅读 2012-4-8 15:58 |个人分类:Matlab|系统分类:科研笔记|关键词:学者| div, 系数, 多项式

//**************   根    **************
1 求线性方程的根:输入多项式系数:P=[*****],r=roots(P),规定多项式是行向量,根是列向量。。
2 根据根构造多项式:pp=poly(r),real(pp);因为MATLAB无隙地处理复数,当用根重组多项式时,如果一些根有虚部.则poly的结果有一些小的虚部,这是很普通的.只要使用函数real抽取实部.
***************   乘法   ************
函数conv支持多项式乘法(执行两个数组的卷积)
系数a=[***],b=[****],c=conv(a,b),d=a+b 

M文件代码:
————————————————
调用;test(c,d),c,d要先赋值多项式的系数。

%加法和减法(减法时,调用参数改为-d即可)

function p=mmpadd(a,b)
%  mmpadd Polynomial addition
% mmpadd(A,B) adds the polynomial A and B 
if nargin<2
    error('Not enough input arguments')
end
a=a(:).';    %make sure inputs are polynomial row vector 
%使输入的的为行向量
b=b(:).';
na=length(a);
nb=length(b);
p=[zeros(1,nb-na) a]+[zeros(1,na-nb) b]; 
%增加必要的系数00000
————————————————————————
除法:
     [q , r]=deconv(c , b)  %结果是b被c除,给出商多项式q和余数r,在现在情况下r是零因为b和q的乘积恰好是c

求导:polyder,h=polyder(g)
--------------**** 估值:****---------------------------------
函数polyval来完成
x=linspace(-1, 3) ;      %   choose 100 data points between -1and 3. 
p=[1  4  -7  -10] ;      %   uses polynomial p(x) = x^3+4x^2-7x-10 
v=polyval(p , x) ;  %计算计算x值上的p(x),把结果存在v里。然后用函数plot绘出结果
plot(x , v),title(' x^3+4x^2-7x-10 '),  xlabel(' x ') %绘制出结果
//*******有理多项式:
有理多项式由它们的分子多项式和分母多项式表示,对有理多项式进行运算的两个函数,residue和polyder.函数residue 执行部分分式展开



https://m.sciencenet.cn/blog-645611-556774.html

上一篇:数据分析函数
下一篇:多项式操作特性

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-13 06:05

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部