要估计y=a(5)+a(1).*x1+a(2).*x2.^2+a(3).*((x4.^a(4))*x3);
其中X4是一个n*n的矩阵
y=data(:,25);
n=length(y);
x1=data(:,2);
x2=data(:,4);
x3=data(:,3);
x4=data(1:20,5:24);
x=[x1 x2 x3 x4];
nlsfun=inline('a(5)+a(1).*x(:,1)+a(2).*x(:,2).^2+a(3).*((x(1:20,4:23).^a(4))*x(:,3)','a','x')
result=nlinfit(x,y,nlsfun,[0 0 0 1 1]);
result

程序跑不出来，请教高手正确程序。