已知X,Y,以及X,Y所对应的矩阵Z,各位大神怎么利用MATLAB编程(不用工具箱)来实现Z=f(X,Y)的拟合,需要拟合的函数形式为Z=p0+p1*X+p2*X^2+p3*Y+p4*Y^2+p5*X*Y,有会的大神请联系我,QQ514157119。
pp=nlinfit(xy,z,fun,[1 1 1 1 1 1]);里边的[1 1 1 1 1 1]是什么意思?程序里边也需要这么写么?
运行结果最后报错是怎么回事?
[1 1 1 1 1 1]是给p0~p5的6个初值,可以改变,一定要有的
出现错误主要是,你没有用z=Z(:);将z的数据变成数组
如果你的原始数据是X,Y,Z都是大小为mxn的矩阵记得先运行
xy=[X(:) Y(:)];
z=Z(:);
这时候xy是(mxn)x2的矩阵,z是(mxn)x1
如果过你原始的数据X是mx1,Y是nx1,Z是mxn,那么还需要将X,Y扩展为mxn的矩阵
我的原始数据X是20X1的,Y是12X1的,这个要怎么扩展成20X12的啊?能不能加我q下。
分给你了
等下把我写的程序发给你,你帮我改改,谢谢了
如果你的z是20x12
那么
[yy xx]=meshgrid(Y,X);
xy=[xx(:) yy(:)];
z=Z(:);
如果你的z是12x20
[xx yy]=meshgrid(X,Y);
xy=[xx(:) yy(:)];
z=Z(:);
这个时候xy应该是240x2,z应该是240x1
这个时候再用
fun=@(p,t) p(1)+p(2)*t(:,1)+p(3)*t(:,1).^2+p(4)*t(:,2)+p(5)*t(:,2).^2+p(6).*t(:,1).*t(:,2);
pp=nlinfit(xy,z,fun,[1 1 1 1 1 1]);
拟合