| 3 | 1/1 | 返回列表 |
| 查看: 495 | 回復(fù): 2 | |||
amorphous鐵桿木蟲 (著名寫手)
小木蟲之超級奶爸
|
[求助]
matlab解非線性方程組,求助!
|
|
幫忙看一下下列程序,如何求得實(shí)數(shù)解. function q=zxxfun(n) logkp=[11.3,6.2,14.8,12.4,19.5,17.1,22.5,2.4,9.8]; p=1; p0=1; nn10=1; nn20=2; q(1)=logkp(1)-(log10(n(4))+log10(n(5))+2*log10(sum(n))-3*log10(n(2))-log10(n(1))-2*log10(p/p0)); q(2)=logkp(2)-(2*log10(n(4))+log10(n(6))+3*log10(sum(n))-4*log10(n(2))-2*log10(n(1))-3*log10(p/p0)); q(3)=logkp(3)-(2*log10(n(4))+log10(n(7))+4*log10(sum(n))-5*log10(n(2))-2*log10(n(1))-4*log10(p/p0)); q(4)=logkp(4)-(3*log10(n(4))+log10(n(8))+5*log10(sum(n))-6*log10(n(2))-3*log10(n(1))-5*log10(p/p0)); q(5)=logkp(5)-(3*log10(n(4))+log10(n(9))+6*log10(sum(n))-7*log10(n(2))-3*log10(n(1))-6*log10(p/p0)); q(6)=logkp(6)-(4*log10(n(4))+log10(n(10))+7*log10(sum(n))-8*log10(n(2))-4*log10(n(1))-7*log10(p/p0)); q(7)=logkp(7)-(4*log10(n(4))+log10(n(11))+8*log10(sum(n))-9*log10(n(2))-4*log10(n(1))-8*log10(p/p0)); q(8)=logkp(8)-(log10(n(2))+log10(n(3))-log10(n(1))-log10(n(4))); q(9)=logkp(9)-(log10(n(3))+log10(sum(n))-2*log10(n(1))-log10(p/p0)); q(10)=nn10-(n(1)+2*n(3)+n(4)); q(11)=2*nn20-(2*n(2)+2*n(4)+4*n(5)+4*n(6)+6*n(7)+6*n(8)+8*n(9)+8*n(10)+10*n(11)); clear options=optimset('Display','iter'); n=fsolve(@zxxfun,[0.5 0.5 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1,0.1],options); |
鐵桿木蟲 (職業(yè)作家)
|
不知道理解對否,用1stOpt計算: n4: 0.818320026230474 n5: 0.585572164827173 n1: 4.65521821183815E-6 n2: 0.0105340963919806 n3: 0.0908376592736609 n6: 1.85174177533351E-13 n7: 5.15897663173426E-7 n8: 8.17642694378187E-17 n9: 7.20355202419344E-12 n10: 1.14168605646827E-21 n11: 2.00691894813261E-18 |
榮譽(yù)版主 (著名寫手)
彩色掛圖
![]() |
專家經(jīng)驗(yàn): +2 |

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工學(xué)調(diào)劑 +7 | 用戶573181 2026-03-20 | 12/600 |
|
|---|---|---|---|---|
|
[考研]
|
13659058978 2026-03-24 | 3/150 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +6 | holy J 2026-03-21 | 6/300 |
|
|
[考研] 081700 調(diào)劑 267分 +8 | 迷人的哈哈 2026-03-23 | 8/400 |
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 材料專碩英一數(shù)二306 +8 | z1z2z3879 2026-03-18 | 8/400 |
|
|
[考研] 生物學(xué)一志愿985,分?jǐn)?shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-17 | 10/500 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
[考研] 考研調(diào)劑 +4 | 來好運(yùn)來來來 2026-03-21 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|