 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1193 | 回復(fù): 5 | |||
[交流]
【求助】最小二乘問題（急） 已有3人參與
|
|
我是個(gè)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] = ... 請(qǐng)問各位大蝦，應(yīng)該怎么解決。 以前用都可以，只是方程變了 [ Last edited by ajicsu on 2010-5-19 at 12:12 ] |
回復(fù)此樓» 猜你喜歡
求助
已經(jīng)有4人回復(fù)
-
0805材料320求調(diào)劑
已經(jīng)有3人回復(fù)
-
297求調(diào)劑
已經(jīng)有11人回復(fù)
-
268求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有10人回復(fù)
-
278求調(diào)劑
已經(jīng)有7人回復(fù)
-
26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑
已經(jīng)有5人回復(fù)
-
280求調(diào)劑
已經(jīng)有9人回復(fù)
-
279求調(diào)劑
已經(jīng)有5人回復(fù)
-
286求調(diào)劑
已經(jīng)有5人回復(fù)
2樓2010-05-19 12:32:10
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲號(hào): 496517
- 注冊(cè): 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
robert2020(金幣+1):鼓勵(lì)交流討論�？紤]可能是因?yàn)閿?shù)據(jù)太大（像運(yùn)行程序后得到的 F 值）導(dǎo)致。樓主最后想要的數(shù)據(jù)沒有出來(lái)。 2010-05-19 12:39:28
kuhailangyu(金幣+1):謝謝參與 2010-05-19 16:50:07
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
robert2020(金幣+1):鼓勵(lì)交流討論�？紤]可能是因?yàn)閿?shù)據(jù)太大（像運(yùn)行程序后得到的 F 值）導(dǎo)致。樓主最后想要的數(shù)據(jù)沒有出來(lái)。 2010-05-19 12:39:28
kuhailangyu(金幣+1):謝謝參與 2010-05-19 16:50:07
 這是啥數(shù)? F = 1.0e+265 * 1.568034000166487 -3.175955144234988 0.000000000000000 0.000000000000000 |
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小時(shí)
- 蟲號(hào): 496517
- 注冊(cè): 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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +8 | 楓橋ZL 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 299求調(diào)劑 +3 | 某某某某位 2026-03-21 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +4 | shxchem 2026-03-20 | 6/300 |
|
|
|
[考研] 336求調(diào)劑 +3 | rmc8866 2026-03-21 | 3/150 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開心 2026-03-16 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|