[求助]求大侠指导这个nlinfit拟合可行否

hrainbow

x=[t1 t2 t3 t4 t5  t6  th(4:57,2) const];   

aveW=sum(th(4:57,1))/nx;aveS=sum(th(4:57,2))/nx;     

b0=[ aveW^2 1 0.5 aveW  aveS 1];

myfunc=inline('beta(1)*(x(:,1)+x(:,2)*beta(2)+x(:,3)*beta(2)^2+x(:,4)*beta(3)+x(:,5)*beta(3)^2+x(:,6)*beta(3)^3)^2  +   beta(4)*(x(:,1)+x(:,2)*beta(2)+x(:,3)*beta(2)^2+x(:,4)*beta(3)+x(:,5)*beta(3)^2+x(:,6)*beta(3)^3)  +   beta(5)x(:,7)  +  beta(6)*x(:,8)','beta','x');

b=nlinfit(x,p,myfunc,b0);   

t1，t2等都是54*1的列向量，我试图用ninfit来拟合 结果报错如下

??? Error using ==> nlinfit
The inline model function generated the following error:
Error using ==> inlineeval
Error in inline expression ==> beta(1)*(x(:,1)+x(:,2)*beta(2)+x(:,3)*beta(2)^2+x(:,4)*beta(3)+x(:,5)*beta(3)^2+x(:,6)*beta(3)^3)^2+beta(4)*(x(:,1)+x(:,2)*beta(2)+x(:,3)*beta(2)^2+x(:,4)*beta(3)+x(:,5)*beta(3)^2+x(:,6)*beta(3)^3)+beta(5)x(:,7)+beta(6)*x(:,8)
??? Error: Unexpected MATLAB expression.
Error in ==> erci at 40
b=nlinfit(x,p,myfunc,b0);  

要拟合的方程有点长，求大虾指导我哪里弄错了？或者说nlinfit能否做这样的拟合？不行的话matlab还要其他什么命令能做到的？
再问一件事情呵呵，在这里常看到有什么说拟合用1stOpt拟合，这个软件貌似都是收费的吗？
式子有点长，问题也有点多，请各位多多指教啦，最近要用到刚这个头都大了，先谢过了

friendchj

beta(5)x(:,7)之间少了个“*”.另外用inline时不知道x(:,1)这种形式是合适的吗？

hrainbow

x的形式是可以的 前面做一次拟合时能行的
现在加了*号了 编译还有问题
??? Error using ==> nlinfit
The inline model function generated the following error:
Error using ==> inlineeval
Error in inline expression ==> beta(1)*(x(:,1)+x(:,2)*beta(2)+x(:,3)*beta(2)^2+x(:,4)*beta(3)+x(:,5)*beta(3)^2+x(:,6)*beta(3)^3)^2+beta(4)*(x(:,1)+x(:,2)*beta(2)+x(:,3)*beta(2)^2+x(:,4)*beta(3)+x(:,5)*beta(3)^2+x(:,6)*beta(3)^3)+beta(5)*x(:,7)+beta(6)*x(:,8)
??? Error using ==> mpower
Matrix must be square.
Error in ==> erci at 40
b=nlinfit(x,p,myfunc,b0);  
谢谢啦呵呵我自己都没看出来

hrainbow

再次谢谢啦，拟合出来了呵呵，大括号后面的^2应该是.^2，
原来这样也能拟合的，我还担心不行呢，看了我得细心点啊

hrainbow

再来问下matlab里边单行程序分行写是这样吗
myfunc=inline('beta(1)*(x(:,1)+x(:,2)*beta(2)+x(:,3)*beta(2)^2+x(:,4)*beta(3)+x(:,5)*beta(3)^2+x(:,6)*beta(3)^3)^2 ...
+   beta(4)*(x(:,1)+x(:,2)*beta(2)+x(:,3)*beta(2)^2+x(:,4)*beta(3)+x(:,5)*beta(3)^2+x(:,6)*beta(3)^3) ...
+   beta(5)x(:,7)  +  beta(6)*x(:,8)','beta','x');

我编译不行啊？

friendchj

myfunc=inline(['beta(1)*(x(:,1)+x(:,2)*beta(2)+x(:,3)*beta(2)^2+x(:,4)*beta(3)+x(:,5)*beta(3)^2+x(:,6)*beta(3)^3)^2', ...
'+   beta(4)*(x(:,1)+x(:,2)*beta(2)+x(:,3)*beta(2)^2+x(:,4)*beta(3)+x(:,5)*beta(3)^2+x(:,6)*beta(3)^3)', ...
'+   beta(5)x(:,7)  +  beta(6)*x(:,8)'],'beta','x');
这样试试