//************** 根 **************
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
上一篇:
数据分析函数下一篇:
多项式操作特性