 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 577 | 回復: 2 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
旖旎落下金蟲 (小有名氣)
|
[求助]
求解非線性六元方程組，時間緊，自己來不及學了，麻煩大家?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 謝謝 |
回復此樓» 猜你喜歡
0805求調劑
已經(jīng)有5人回復
-
342求調劑
已經(jīng)有7人回復
-
調劑
已經(jīng)有6人回復
-
一志愿北京科技大學085601材料工程英一數(shù)二初試總分335求調劑
已經(jīng)有3人回復
-
273求調劑
已經(jīng)有4人回復
-
339求調劑
已經(jīng)有3人回復
-
267求調劑
已經(jīng)有7人回復
-
070300化學279求調劑
已經(jīng)有11人回復
-
318求調劑
已經(jīng)有3人回復
-
309分085801求調劑
已經(jīng)有4人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求助：非線性方程組的求解（郁悶中）
已經(jīng)有26人回復
- 數(shù)學小白求助，解超定非線性方程組，求大神幫忙
已經(jīng)有4人回復
- Matlab用牛頓法求解非線性方程組問題
已經(jīng)有8人回復
- 用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？
已經(jīng)有11人回復
- matlab如何求解一個非線性微分方程組
已經(jīng)有8人回復
- MATLAB求解非線性方程組
已經(jīng)有5人回復
- 關于matlab求解非線性指數(shù)方程組出現(xiàn)問題
已經(jīng)有3人回復
- matlab求解非線性方程組，錯誤提示怎么解決
已經(jīng)有5人回復
- matlab求解非線性方程組，求助！
已經(jīng)有6人回復
- 求助！matlab用fsolve函數(shù)求解非線性方程組的問題！
已經(jīng)有19人回復
- 如何使用matlab求解非線性方程組的所有整數(shù)解？
已經(jīng)有9人回復
- matlab求解非線性方程組
已經(jīng)有16人回復
- 求高人指點用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復
- 【求助】用mathematica 5.0求解一個非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復
旖旎落下
金蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 778.5
- 帖子: 53
- 在線: 34.4小時
- 蟲號: 2510792
- 注冊: 2013-06-17
- 性別: MM
- 專業(yè): 物流與供應鏈管理
3樓2014-01-09 10:23:46
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
旖旎落下: 金幣+50, ★★★很有幫助, 好。 2014-01-09 10:23:14
感謝參與，應助指數(shù) +1
旖旎落下: 金幣+50, ★★★很有幫助, 好。 2014-01-09 10:23:14
|
以下是如何求數(shù)值解。注意，為了簡化，讓優(yōu)化工具箱自行計算方程組的數(shù)值一階微分（雅可比矩陣）和二階微分（Hessian矩陣）。如果結果不理想，可以先調節(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; %簡化表達式 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 %設置優(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); %設置初始值 x0 = [0.5, 0.5, 0.5, 0.5, 0.5, 0.5]; %調用優(yōu)化函數(shù) [x, fval, exitflag, output] = fsolve(@func1, x0, options); %取結果 qn=x(1); pn=x(2); pr=x(3); t=x(4); e=x(5); u1=x(6); %顯示結果 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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311（085601）求調劑 +9 | liziyeyeye 2026-03-28 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 考研調劑 +9 | 小蠟新筆 2026-03-29 | 10/500 |
|
|
|
[考研] 材料科學與工程求調劑 +13 | 深V宿舍吧 2026-03-29 | 13/650 |
|
|
|
[碩博家園] 求調劑 有機化學考研356分 +11 | Nadiums 2026-03-25 | 12/600 |
|
|
|
[考研] 254材料與化工求調劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
|
[考研] 08開頭看過來！��！ +3 | wwwwffffff 2026-03-31 | 5/250 |
|
|
|
[考研] 一志愿211，335分，0856，求調劑院校和導師 +10 | 傾____蕭 2026-03-27 | 11/550 |
|
|
|
[考研] 生物考研337分求調劑 +4 | cgxin 2026-03-30 | 6/300 |
|
|
|
[考研] 一志愿西電085401數(shù)一英一299求調劑 六級521 +4 | 愛吃大鴨梨 2026-03-31 | 4/200 |
|
|
|
[考研] 297 地理學070500 復試求調劑 +3 | 小圓圈圈ooo 2026-03-30 | 3/150 |
|
|
|
[考研] 317求調劑 +8 | 十閑wx 2026-03-24 | 8/400 |
|
|
|
[考研] 310求調劑 +10 | 爭取九點睡 2026-03-30 | 10/500 |
|
|
|
[考研] 求調劑，一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
|
[考研] 求調劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
|
[考研] 327求調劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
|
[考研] 298求調劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[考研] 081200-11408-276學碩求調劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
|
[考研] 085602 化工專碩 338分 求調劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
|
[考研] 321求調劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿天津大學339材料與化工求調劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|