| 查看: 1391 | 回复: 5 | |||
[交流]
【求助】最小二乘问题（急） 已有3人参与
|
|
我是个MATLAB新手，在应用编辑器编辑函数： function F = myfun(x) a_m=534.87,b_m=-301.56,a_n=534.86,b_n=-301.58,c_m=-0.15,d_m=-0.27,c_n=-0.15,d_n=-0.27,l=300; A=0.5*(exp(x(1)*l)+exp(-x(1)*l))*cos(x(2)*l); B=0.5*(exp(x(1)*l)-exp(-x(1)*l))*sin(x(2)*l); C=0.5*(exp(x(1)*l)-exp(-x(1)*l))*cos(x(2)*l); D=0.5*(exp(x(1)*l)+exp(-x(1)*l))*sin(x(2)*l); E=0.5*(exp(0.5*x(1)*l)+exp(-0.5*x(1)*l))*cos(0.5*x(2)*l); I=0.5*(exp(0.5*x(1)*l)-exp(-0.5*x(1)*l))*sin(0.5*x(2)*l); G=0.5*(exp(0.5*x(1)*l)-exp(-0.5*x(1)*l))*cos(0.5*x(2)*l); H=0.5*(exp(0.5*x(1)*l)+exp(-0.5*x(1)*l))*sin(0.5*x(2)*l); J=a_m*A-b_m*B-x(3)*c_m*C+x(4)*d_m*C+x(4)*c_m*D+x(3)*d_m*D; K=a_n*A-b_n*B-x(3)*c_n*C+x(4)*d_n*C+x(4)*c_n*D+x(3)*d_n*D; L=A*b_m+B*a_m-x(4)*c_m*C-x(3)*d_m*C-x(3)*c_m*D+x(4)*d_m*D; M=A*b_n+B*a_n-x(4)*c_n*C-x(3)*d_n*C-x(3)*c_n*D+x(4)*d_n*D; N=a_m*E-b_m*I-x(3)*c_m*G+x(4)*d_m*G+x(4)*c_m*H+x(3)*d_m*H; O=E*b_m+I*a_m-x(4)*c_m*G-x(3)*d_m*G-x(3)*c_m*H+x(4)*d_m*H; P=a_n*E-b_n*I-x(3)*c_n*G+x(4)*d_n*G+x(4)*c_n*H+x(3)*d_n*H; Q=E*b_n+I*a_n-x(4)*c_n*G-x(3)*d_n*G-x(3)*c_n*H+x(4)*d_n*H; F=[J*K-L*M-a_n*a_m+b_m*b_n; J*M+K*L-a_m*b_n-a_n*b_m; P*J-Q*L-N*a_n+O*b_n; P*L+Q*J-N*b_n-O*a_n] 然后在命令窗输入： x0=[1 2 3 4]; [x]=lsqnonlin(@myfun,x0) 结果出现这种情况： ??? Error using ==> eig NaN or Inf prevents convergence. Error in ==> C:\MATLAB6p5p1\toolbox\optim\private\trust.m On line 28 ==> [V,D] = eig(H); Error in ==> C:\MATLAB6p5p1\toolbox\optim\private\trdog.m On line 108 ==> [st,qpval,po,fcnt,lambda] = trust(rhs,MM,delta); Error in ==> C:\MATLAB6p5p1\toolbox\optim\private\snls.m On line 375 ==> [sx,snod,qp,posdef,pcgit,Z] = trdog(x,g,A,D,delta,dv,... Error in ==> C:\MATLAB6p5p1\toolbox\optim\private\lsqncommon.m On line 207 ==> [x,FVAL,LAMBDA,JACOB,EXITFLAG,OUTPUT,msg]=... Error in ==> C:\MATLAB6p5p1\toolbox\optim\lsqnonlin.m On line 122 ==> [x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ... 请问各位大虾，应该怎么解决。 以前用都可以，只是方程变了 [ Last edited by ajicsu on 2010-5-19 at 12:12 ] |
回复此楼» 猜你喜欢
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有11人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
2楼2010-05-19 12:32:10
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
3楼2010-05-19 12:35:51
4楼2010-05-19 14:17:55
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
5楼2010-05-19 20:19:51
6楼2010-05-21 18:36:55
30