 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 574 | 回復(fù): 2 | ||
旖旎落下金蟲 (小有名氣)
|
[求助]
求解非線性六元方程組，時間緊，自己來不及學(xué)了，麻煩大家?guī)兔�，問題簡單，懸賞多。 已有1人參與
|
|
用MATLAB 非線性求解，不會 >> 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ù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(fù)
-
生物學(xué)學(xué)碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
-
求材料，環(huán)境專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
335求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿吉大化學(xué)322求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助：非線性方程組的求解（郁悶中）
已經(jīng)有26人回復(fù)
- 數(shù)學(xué)小白求助，解超定非線性方程組，求大神幫忙
已經(jīng)有4人回復(fù)
- Matlab用牛頓法求解非線性方程組問題
已經(jīng)有8人回復(fù)
- 用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？
已經(jīng)有11人回復(fù)
- matlab如何求解一個非線性微分方程組
已經(jīng)有8人回復(fù)
- MATLAB求解非線性方程組
已經(jīng)有5人回復(fù)
- 關(guān)于matlab求解非線性指數(shù)方程組出現(xiàn)問題
已經(jīng)有3人回復(fù)
- matlab求解非線性方程組，錯誤提示怎么解決
已經(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ù)
- 求高人指點用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復(fù)
- 【求助】用mathematica 5.0求解一個非線性方程組失敗，特發(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)化工具箱自行計算方程組的數(shù)值一階微分（雅可比矩陣）和二階微分（Hessian矩陣）。如果結(jié)果不理想，可以先調(diào)節(jié) options 中的參數(shù)。如果還不理想，可以利用符號運算工具箱求出雅可比矩陣和 Hessian 矩陣。 1. 建立一個函數(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. 建立一個主程序 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小時
- 蟲號: 2510792
- 注冊: 2013-06-17
- 性別: MM
- 專業(yè): 物流與供應(yīng)鏈管理
3樓2014-01-09 10:23:46
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 旖旎落下 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +11 | 手機(jī)用戶 2026-03-17 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +7 | vv迷 2026-03-22 | 7/350 |
|
|
|
[考研] 299求調(diào)劑 +5 | 某某某某位 2026-03-21 | 5/250 |
|
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +5 | 葵梓衛(wèi)隊 2026-03-18 | 7/350 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項目，四級已過 +6 | 幸運哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
|
[考研] 生物學(xué)調(diào)劑招人�。�！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|