| 5 | 1/1 | 返回列表 |
| 查看: 1214 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[交流]
【求助】最小二乘問題(急) 已有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ā)表 | |
|---|---|---|---|---|
|
[考研] 296材料專碩求調(diào)劑 +11 | 202451007219 2026-04-02 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 302分初試 0854 +5 | 伶可樂 2026-04-02 | 5/250 |
|
|
[考研] 346求調(diào)劑 +5 | 鄭誠樂 2026-04-02 | 5/250 |
|
|
[基金申請]
請問共同通訊和共同一作的認(rèn)可度問題
10+4
|
psa1234 2026-04-01 | 8/400 |
|
|
[考研] 學(xué)碩機(jī)械工程303求調(diào)劑 +4 | 無名所以叫吳明 2026-03-30 | 4/200 |
|
|
[考研] 材料科學(xué)與工程調(diào)劑 +18 | 深V宿舍吧 2026-03-30 | 19/950 |
|
|
[考研] 324求調(diào)劑 +5 | 想上學(xué)求調(diào) 2026-04-01 | 6/300 |
|
|
[考研]
|
廈大化工 2026-04-01 | 5/250 |
|
|
[考研] 343求調(diào)劑085601 +4 | 要努力學(xué)習(xí)x 2026-03-29 | 5/250 |
|
|
[考研] 283求調(diào)劑 +9 | A child 2026-03-28 | 9/450 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-29 | 10/500 |
|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +6 | zzll03 2026-03-31 | 6/300 |
|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
[考研] 313求調(diào)劑 +6 | 賣個關(guān)子吧 2026-03-31 | 6/300 |
|
|
[有機(jī)交流]
10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué),化工,生物,生醫(yī)工等專業(yè)優(yōu)秀學(xué)子加入! +3 | llss0711 2026-03-28 | 6/300 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|