目录

%作者:John D'Errico%发行版本:2.0%发布日期:2006年8月8日%以下是polyfitn、polyvaln、poly2sympoly、Poly2Sym的示例用法

将一维模型拟合到cos(x)。我们只需要偶数订单条款。

x=-2:1:2;y=cos(x);p=polyfitn(x,y,'常数x^2 x^4 x^6')如果存在(“sympoly”) == 2%转化为一个sympoly。如果没有别的,这是显示模型的好方法。多项式2sympoly(p)终止如果存在(“sym”) == 2%向符号形式的转换。它也很好。多项式(p)终止%在某些点上评估回归模型polyvaln(p[0.51])
p=ModelTerms:[4x1 double]系数:[0.99996-0.49968 0.041242-0.0012079]参数var:[1.2876e-10 1.084e-09 4.6603e-10 1.3903e-11]参数std:[1.1347e-05 3.2925e-05 2.1588e-05 3.7286e-06]R2:1 RMSE:3.1468e-05变量名:{'x'}标量共聚合对象0.99996-0.49968^2+0.042*1240*0.9976^0.990

二维曲面模型,所有项都达到三阶。

%使用大量数据。n=1000;x=兰特(n,2);y=exp(和(x,2))+randn(n,1)/100;p=聚合体(x,y,3)如果存在(“sympoly”)==2多项式2sympoly(p)终止如果存在(“sym”)==2多项式(p)终止%在网格和绘图上进行评估:[xg,yg]=meshgrid(0.05:1);zg=polyvaln(p,[xg(:),yg(:)];冲浪(xg、yg、重塑(zg、尺寸(xg))保持在…上图3(x(:,1),x(:,2),y,“哦”)持有
p=ModelTerms:[10x2 double]系数:[1x10 double]参数var:[1x10 double]参数std:[1x10 double]R2:0.99992 RMSE:0.011198变量名:{}一个标量sypoly对象0.50402*X1^3+1.3927*X1^2*X2-0.017552*X1^2+1.3798*X1*X2^2+0.081473*X1+0.43959*X2^3+0.087656*X2^2+1.2254*X2+0.96196

二维线性模型,但没有常数项

uv=兰特(100,2);w=sin(和(uv,2));p=聚乙烯纤维(紫外线,w,“u,v”);如果存在(“sympoly”)==2多项式2sympoly(p)终止如果存在(“sym”)==2多项式(p)终止
标量共多边形对象0.76416*u+0.70472*v

具有各种指数的模型,而不是所有正整数。

%注:只有1个变量时,x&y可以是行向量或列向量。x=1:10;y=3+2./x+sqrt(x)+randn(尺寸(x))/100;p=聚管件(x,y,'常数x^-1 x^0.5');如果存在(“sympoly”)==2多项式2sympoly(p)终止如果存在(“sym”)==2多项式(p)终止席=1∶1:10;yi=polyvaln(p,xi);图(x,y,“罗”,席,一,“b-”)
标量共多边形对象2.9805+2.0448*x^-1+1.0041*x^0.5