matlab中求多项式 1-1 2 1 3......

interp?

symssps=((s^2+1))^3*(s+5)^2*(s^4+4*s^2+7)ps1=expand(ps)结果：ps=(s^2+1)^3*(s+5)^2*(s^4+4*s^2+7)ps1=175+

倒是可以实现,但系数怎么存放好呢,是弄成像你写得这样是不行的[230,412,031]-------要是一行存放的话,只能借助cell了clearall;clc;symsxa=[2*x^2+3*x,4

symsabcdx;p=a*x^2+c*x+c+d*x;t=coeffs(p,x);t(2)%%%%输出结果=c+d即为所得.coeffs(p,x)的结果是按照变量的幂来排列的.如上t(1)为常系数c

x=1978:2002;y=[272.81,286.43,311.89,324.76,337.07,351.81,390.85,466.75,490.86,545.46,648.3,696.54,78

那你就都改成columnvectorsz=[-3;-2.5];p=[-5;-2;-15];

多项式的除法用deconv（)；多项式的乘法用conv()；例如上面的公式:a=[1011];b=[201];[q,r]=deconv(a,b)q=0.50000r=000.50001.0000r表示

clc;clearsymsa5b1b2b3b4b5zf=sym('-a5*b1*z^16+a5*b4*z^2*b3*b2*b1-a5*b3*z^12*b1-a5*b3*z^20*b2+a5*b2*z^

记该函数为f,用eval语句完整代码如下symsxx;f=-(390625*xx*(xx-1)*(xx+1)*(xx-1/5)*(xx+1/5)*(xx-2/5)*(xx+2/5)*(xx-3/5)*

poly2sym([-0.00140.0360-0.36681.8817-4.98835.9183-1.21980.0770])即可得到该多项式就像一个指数函数如果用多项式表示的话需要无穷多个多项式才

limit((x-3)./abs(x-3),x,3,'left')ans=-1limit((x-3)./abs(x-3),x,3,'right')ans=+1

matlab里面的矩阵就可以实现多项式的加减的,把多项式的系数赋给矩阵的每一个元素,次数不一样的话以高次的为准,没有高次的话就补0；然后用矩阵直接相加就可以啦.不知道有没有说清楚.

程序改成如下：A=[13631];s=poly2str(A,'x');x=unifrnd(0,5,150,1)%产生150个[05]上的均匀分布随机数y=polyval(A,x)%5阶多项式拟合p=p

套公式,代数!再问：呵呵呵，无语再答：那你要给个具体数据吧！

polyder的对应积分函数为polyint

1.根据实际系统理论模型的阶数拟合.2.当然一般情况下,1的条件并不知道,这时可以采用尝试的方法.3.还可以使用样条拟合（这个能基本保证所有点在曲线上）,就是spline

p=[13-5-6];a=roots(p)';A=blkdiag(a(1),a(2),a(3))先求出特征值,然后以这些特征值为对角线元素的矩阵就是所求

x=[];y=[];F=@(p,a)p(1)*a+p(2)+p(3)*sin(p(4)*a+p(5));p=lsqcurvefit(F,[11111],x,y)%p即为所拟合函数系数,分别为a,b,c

poly2sym(poly([x1,x2,x3,...]))x1,x2,x3...是根.

样条是分段函数,多项式的个数和分的段数有关的,这个段数应该是不定的,所以不太好弄.那么：你想要一个三次多项式的话,就用polyfit,得到多项式的系数如果你想三次样条插值,就用spline或者inte