目录
将一维模型拟合到cos(x)。我们只需要偶数订单条款。
x=-2:1:2;y=cos(x);p=polyfitn(x,y,'常数x^2 x^4 x^6')如果存在(“sympoly”) == 2多项式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
具有各种指数的模型,而不是所有正整数。
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