| 查看: 1382 | 回复: 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 ] |
回复此楼» 猜你喜欢
能否退出参与的面上项目解除限项
已经有18人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有3人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
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
