 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 575 | 回復(fù): 2 | |||
旖旎落下金蟲 (小有名氣)
|
[求助]
求解非線性六元方程組，時(shí)間緊，自己來不及學(xué)了，麻煩大家?guī)兔Γ瑔栴}簡單，懸賞多。 已有1人參與
|
|
用MATLAB 非線性求解，不會(huì) >> syms qn dn dr pn pr cn cr sn sr hn hr t e a o A T u1 k Fn Fr zn zr >>o=0.8; >>a=0.1; >>cn=0.4; >>cr=0.2; >>k=2; >>T=1; >>hn=0.1; >>sn=0.1; >>hr=0.05; >>sr=0.05; >>A=0.1; >> dn=1-(pn-pr)/(1-o)>> dr=(o*pn-pr)/(o-o^2) >> Fn=20*(qn-dn) >> Fr=20*(qn*e*(a*t+1)-dr) >> zn=qn-dn >> zr=qn*e*(a*t+1)-dr >> eq1='-(pn+hn-cn-sn)*Fn+pn+hn-cn-e*(a*t+1)*(pr+hr+A*(T-t)-cn-sr)*Fr+e*(a*t+1)*(pr+hr+A*(T-t)-cn)+u1*(1-e*(a*t+1))=0' >> eq2='qn+(hn-hr)/(1-o)+(1+qn-(2*pn-pr-sn-cn+hn)/(1-o))*Fn+Fr*(pr+hr+A*(T-t)-cn-sr)/(1-o)+10*zn^2=0' >> eq3='e*qn*(a*t+1)+(hr-o*hn)/(o-o^2)+Fn*(pn+hn-cn-sn)/(1-o)+Fr*(o*pn-2*pr+cn+sr-A*(T-t)-hr-o*(1-o)*e*qn*(a*t+1))/(o-o^2)+10*0.05^2-10*zr^2=0' >> eq4='a*e*qn*(hr+pr+A*(T-2*t-1/a)-cn)+a*e*qn*(sr+cn-hr-A*(T-2*t-1/a)-pr)*Fr-A*(20*dr*zr+10*zr^2)-u1*qn*e*a=0' >> eq5='qn*(a*t+1)*(pr+hr+A*(T-t)-cn)-qn*(a*t+1)*(pr+hr+A*(T-t)-cn-sr)*Fr-2*k*e-u1*a*t*qn=0' >> eq6='qn*(1-e*(a*t+1))=0' 最終求解qn pn pr t e u1 是用fsolve算嗎，初始值不知道怎么賦，大致推遲大概是0.5 0.5 0.5 0.5 0.5 謝謝 |
回復(fù)此樓» 猜你喜歡
申博26年
已經(jīng)有4人回復(fù)
-
307求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿C9材料與化工專業(yè)總分300求調(diào)劑
已經(jīng)有5人回復(fù)
-
274求調(diào)劑
已經(jīng)有3人回復(fù)
-
299求調(diào)劑
已經(jīng)有7人回復(fù)
-
0703化學(xué)調(diào)劑，求導(dǎo)師收
已經(jīng)有4人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有3人回復(fù)
-
材料學(xué)碩，求調(diào)劑
已經(jīng)有6人回復(fù)
-
307求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩331求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助：非線性方程組的求解（郁悶中）
已經(jīng)有26人回復(fù)
- 數(shù)學(xué)小白求助，解超定非線性方程組，求大神幫忙
已經(jīng)有4人回復(fù)
- Matlab用牛頓法求解非線性方程組問題
已經(jīng)有8人回復(fù)
- 用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？
已經(jīng)有11人回復(fù)
- matlab如何求解一個(gè)非線性微分方程組
已經(jīng)有8人回復(fù)
- MATLAB求解非線性方程組
已經(jīng)有5人回復(fù)
- 關(guān)于matlab求解非線性指數(shù)方程組出現(xiàn)問題
已經(jīng)有3人回復(fù)
- matlab求解非線性方程組，錯(cuò)誤提示怎么解決
已經(jīng)有5人回復(fù)
- matlab求解非線性方程組，求助！
已經(jīng)有6人回復(fù)
- 求助！matlab用fsolve函數(shù)求解非線性方程組的問題！
已經(jīng)有19人回復(fù)
- 如何使用matlab求解非線性方程組的所有整數(shù)解？
已經(jīng)有9人回復(fù)
- matlab求解非線性方程組
已經(jīng)有16人回復(fù)
- 求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復(fù)
- 【求助】用mathematica 5.0求解一個(gè)非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復(fù)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
旖旎落下: 金幣+50, ★★★很有幫助, 好。 2014-01-09 10:23:14
感謝參與，應(yīng)助指數(shù) +1
旖旎落下: 金幣+50, ★★★很有幫助, 好。 2014-01-09 10:23:14
|
以下是如何求數(shù)值解。注意，為了簡化，讓優(yōu)化工具箱自行計(jì)算方程組的數(shù)值一階微分（雅可比矩陣）和二階微分（Hessian矩陣）。如果結(jié)果不理想，可以先調(diào)節(jié) options 中的參數(shù)。如果還不理想，可以利用符號(hào)運(yùn)算工具箱求出雅可比矩陣和 Hessian 矩陣。 1. 建立一個(gè)函數(shù)文件 func1.m function y=func1(x) %未知量 qn=x(1); pn=x(2); pr=x(3); t=x(4); e=x(5); u1=x(6); %常量 o=0.8; a=0.1; cn=0.4; cr=0.2; k=2; T=1; hn=0.1; sn=0.1; hr=0.05; sr=0.05; A=0.1; %簡化表達(dá)式 dn=1-(pn-pr)/(1-o); dr=(o*pn-pr)/(o-o^2); Fn=20*(qn-dn); Fr=20*(qn*e*(a*t+1)-dr); zn=qn-dn; zr=qn*e*(a*t+1)-dr; %非線性方程組 y = [-(pn+hn-cn-sn)*Fn+pn+hn-cn-e*(a*t+1)*(pr+hr+A*(T-t)-cn-sr)*Fr+e*(a*t+1)*(pr+hr+A*(T-t)-cn)+u1*(1-e*(a*t+1)); qn+(hn-hr)/(1-o)+(1+qn-(2*pn-pr-sn-cn+hn)/(1-o))*Fn+Fr*(pr+hr+A*(T-t)-cn-sr)/(1-o)+10*zn^2; e*qn*(a*t+1)+(hr-o*hn)/(o-o^2)+Fn*(pn+hn-cn-sn)/(1-o)+Fr*(o*pn-2*pr+cn+sr-A*(T-t)-hr-o*(1-o)*e*qn*(a*t+1))/(o-o^2)+10*0.05^2-10*zr^2; a*e*qn*(hr+pr+A*(T-2*t-1/a)-cn)+a*e*qn*(sr+cn-hr-A*(T-2*t-1/a)-pr)*Fr-A*(20*dr*zr+10*zr^2)-u1*qn*e*a; qn*(a*t+1)*(pr+hr+A*(T-t)-cn)-qn*(a*t+1)*(pr+hr+A*(T-t)-cn-sr)*Fr-2*k*e-u1*a*t*qn; qn*(1-e*(a*t+1))]; 2. 建立一個(gè)主程序 solve_6_unknowns.m %清屏，清工作區(qū) clc clear all %設(shè)置優(yōu)化算法參數(shù) maxiter = 20000; maxfuneval = length(x0)*maxiter; options = optimset('Display‘， ’off',... 'GradObj', 'off',... 'Hessian', 'off',... 'TolX', 1e-6,... 'TolFun', 1e-6,... 'MaxIter', maxiter,... 'MaxFunEvals', maxfuneval); %設(shè)置初始值 x0 = [0.5, 0.5, 0.5, 0.5, 0.5, 0.5]; %調(diào)用優(yōu)化函數(shù) [x, fval, exitflag, output] = fsolve(@func1, x0, options); %取結(jié)果 qn=x(1); pn=x(2); pr=x(3); t=x(4); e=x(5); u1=x(6); %顯示結(jié)果 disp(['Iterations: ',num2str(output.iterations)]) disp(['Func Evals: ', num2str(output.funcCount)]) disp(['Algorithm: ',output.algorithm]) disp(['exit flag = ',num2str(exitflag)]) disp(['error = ( ',num2str(fval','%-15.6e'),' )']) disp(['qn = ( ',num2str(qn,'%-15.6f'),' )']) disp(['pn = ( ',num2str(pn,'%-15.6f'),' )']) disp(['pr = ( ',num2str(pr,'%-15.6f'),' )']) disp(['t = ( ',num2str(t,'%-15.6f'),' )']) disp(['e = ( ',num2str(e,'%-15.6f'),' )']) disp(['u1 = ( ',num2str(u1,'%-15.6f'),' )']) |
2樓2014-01-08 22:57:57
旖旎落下
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 778.5
- 帖子: 53
- 在線: 34.4小時(shí)
- 蟲號(hào): 2510792
- 注冊(cè): 2013-06-17
- 性別: MM
- 專業(yè): 物流與供應(yīng)鏈管理
3樓2014-01-09 10:23:46
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 旖旎落下 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +5 | WW.' 2026-03-23 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +6 | 加一一九 2026-03-24 | 6/300 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +6 | 小羊睡著了? 2026-03-23 | 6/300 |
|
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +12 | 吃吃吃才有意義 2026-03-19 | 13/650 |
|
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
|
[考研] 07化學(xué)280分求調(diào)劑 +4 | 722865 2026-03-23 | 4/200 |
|
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[考研] 工科0856求調(diào)劑 +5 | 沐析汀汀 2026-03-21 | 5/250 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研] 招08考數(shù)學(xué) +6 | laoshidan 2026-03-20 | 14/700 |
|
|
|
[考研] 求助 +5 | 夢(mèng)里的無言 2026-03-21 | 6/300 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +9 | 簡單點(diǎn)0 2026-03-17 | 9/450 |
|
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|