| 查看: 1462 | 回復(fù): 12 | |||
chain416新蟲 (初入文壇)
|
[求助]
如何使用fsolve求解經(jīng)過微分計算生成的函數(shù)的解?不是微分方程 已有1人參與
|
||
|
標(biāo)題可能不清楚,意思就是fsolve求解一個方程我已經(jīng)會了,但是方程必須是固定已知的,而不能是通過微分計算計算出來的。 比方說: x=fsolve('myfun',[0.5],optimset('Display','off')) function q=myfun(p) alpha1 = p; q =25*exp(5/(2*alpha1)) - 100/(9*exp(5/(2*(alpha1 - 1)))) - 250/(9*exp(5/(2*(alpha1 - 1)))*(alpha1 - 1)) - (125*exp(5/(2*alpha1)))/(2*alpha1) - 125/9; 是可以計算出結(jié)果的,而 x=fsolve('myfun',[0.5],optimset('Display','off')) function q=myfun(p) RT=2.5; h1=0.2; h2=0.3; r1=abs(h1)^2; r2=abs(h2)^2; [P1, P3, beta1]=solve ('alpha1*log(1+P1*r1/alpha1) = RT','beta1*log(1+P3*r2/beta1) = RT','alpha1 + beta1 = 1','P1','P3','beta1'); P1 = subs(P1); P3 = subs(P3); y = P1 + P3; q = diff (y); alpha1 = p; 就會報錯,因?yàn)閝是經(jīng)過微分計算所得,alpha1總是無法將數(shù)值解帶入。 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +31 |

新蟲 (初入文壇)
新蟲 (初入文壇)
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +31 |

新蟲 (初入文壇)
送紅花一朵 |
'alpha1*log(1+P1*r1/alpha1) = RT','beta1*log(1+P3*r2/beta1) = RT','alpha1 + beta1 = 1' 1 這三個等式中有4個未知量,所以假設(shè)alpha1已知,用alpha1來表示另外3個未知量(此時用solve來求解一切正常):'P1','P3','beta1'。 2 然后目標(biāo)函數(shù)P+P3對alpha1進(jìn)行求導(dǎo),得出P1+P3對alpha1的一階導(dǎo)數(shù)。 3 最后再求導(dǎo)數(shù)等于0時,alpha1的數(shù)值解(此時用fsolve無法完成求解)。 請幫忙看下,謝謝。! |
新蟲 (初入文壇)
新蟲 (初入文壇)
新蟲 (初入文壇)
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研化學(xué)308分求調(diào)劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 350求調(diào)劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
[考研]
|
pk3725069 2026-03-19 | 14/700 |
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 3/150 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考研] 一志愿華中科技大學(xué)071000,求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[論文投稿]
申請回稿延期一個月,編輯同意了。但系統(tǒng)上的時間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|