 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1191 | 回復(fù): 5 | |||
[交流]
【求助】最小二乘問題（急） 已有3人參與
|
|
我是個MATLAB新手，在應(yīng)用編輯器編輯函數(shù)： 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) 結(jié)果出現(xiàn)這種情況： ??? 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] = ... 請問各位大蝦，應(yīng)該怎么解決。 以前用都可以，只是方程變了 [ Last edited by ajicsu on 2010-5-19 at 12:12 ] |
回復(fù)此樓» 猜你喜歡
材料與化工（0856）304求 B區(qū) 調(diào)劑
已經(jīng)有3人回復(fù)
-
268求調(diào)劑
已經(jīng)有8人回復(fù)
-
265求調(diào)劑
已經(jīng)有14人回復(fù)
-
350 本科985求調(diào)劑，求老登收留
已經(jīng)有3人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑
已經(jīng)有9人回復(fù)
-
初試 317
已經(jīng)有3人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有15人回復(fù)
-
286求調(diào)劑
已經(jīng)有4人回復(fù)
-
279求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
2樓2010-05-19 12:32:10
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時
- 蟲號: 496517
- 注冊: 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
3樓2010-05-19 12:35:51
4樓2010-05-19 14:17:55
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時
- 蟲號: 496517
- 注冊: 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
5樓2010-05-19 20:19:51
6樓2010-05-21 18:36:55
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ajicsu 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +6 | 煙火先于春 2026-03-17 | 6/300 |
|
|
[考研]
 271材料工程求調(diào)劑
+8
|
.6lL 2026-03-18 | 8/400 |
|
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 5/250 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|