 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 931 | 回復(fù): 0 | ||
jjliu2012新蟲(chóng) (小有名氣)
|
[求助]
【求助 一個(gè)關(guān)于最優(yōu)控制的matlab問(wèn)題（新手，請(qǐng)大家多多幫助）
|
|
Sample Text 最速下降法求最優(yōu)控制u的matlab算法，各位大俠能否把這個(gè)主循環(huán)補(bǔ)充一下，不知道怎么回事，運(yùn)行不了 。是附件里文獻(xiàn)的例子2，在第7頁(yè)，多謝 %%%%%%%%%%%%%%%%%%%%%% for i = 1:max_iteration % 1) start with assumed control u and move forward [Tx,X] = ode45(@(t,x) stateEq(t,x,u,Tu), [t0 tf], ... initx, options); % 2) Move backward to get the trajectory of costates x1 = X(:,1); x2 = X(:,2); [Tp,P] = ode45(@(t,p) costateEq(t,p,u,Tu,x1,x2,Tx), ... [tf t0], initp, options); p1 = P(:,1); % Important: costate is stored in reverse order. The dimension of % costates may also different from dimension of states % Use interploate to make sure x and p is aligned along the time axis p1 = interp1(Tp,p1,Tx); % Calculate deltaH with x1(t), x2(t), p1(t), p2(t) dH = pH(x1,p1,Tx,u,Tu); H_Norm = dH'*dH; % Calculate the cost function J(i,1) = tf*(((x1')*x1 + (x2')*x2)/length(Tx) + ... 0.1*(u*u')/length(Tu)); % if dH/du < epslon, exit if H_Norm < eps % Display final cost J(i,1) break; else % adjust control for next iteration u_old = u; u = AdjControl(dH,Tx,u_old,Tu,step); end; end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % State equations function dx = stateEq(t,x,u,Tu) dx = zeros(2,1); u = interp1(Tu,u,t); % Interploate the control at time t dx(1) = -2*(x(1) + 0.25) + (x(2) + 0.5)*exp(25*x(1)/(x(1)+2)) ... - (x(1) + 0.25).*u; dx(2) = 0.5 - x(2) -(x(2) + 0.5)*exp(25*x(1)/(x(1)+2)); %%%%%%%%%%%%%%%%%%%%%%% % Costate equations function dp = costateEq(t,p,u,Tu,x1,x2,xt) dp = zeros(2,1); x1 = interp1(xt,x1,t); % Interploate the state varialbes x2 = interp1(xt,x2,t); u = interp1(Tu,u,t); % Interploate the control dp(1) = p(1).*(u + exp((25*x1)/(x1 + 2)).*((25*x1)/(x1 + 2)^2 - ... 25/(x1 + 2))*(x2 + 1/2) + 2) - ... 2*x1 - p(2).*exp((25*x1)/(x1 + 2))*((25*x1)/(x1 + 2)^2 - ... 25/(x1 + 2))*(x2 + 1/2); dp(2) = p(2).*(exp((25*x1)/(x1 + 2)) + 1) - ... p(1).*exp((25*x1)/(x1 + 2)) - 2*x2; %%%%%%%%%%%%%%%%%%%%%%%%%%% % Partial derivative of H with respect to u function dH = pH(x1,p1,tx,u,Tu) % interploate the control u = interp1(Tu,u,tx); R = 0.1; dH = 2*R*u - p1.*(x1 + 0.25); %%%%%%%%%%%%%%%%%%%%%%%%% % Adjust the control function u_new = AdjControl(pH,tx,u,tu,step) % interploate dH/du pH = interp1(tx,pH,tu); u_new = u - step*pH; %%%%%%%%%%%%%%%% |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : SolvingoptimalcontrolproblemswithMATLAB.pdf
2013-09-01 16:16:28, 1.95 M
» 猜你喜歡
323求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有4人回復(fù)
-
352求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有8人回復(fù)
-
277材料科學(xué)與工程080500求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有18人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
280分求調(diào)劑 一志愿085802
已經(jīng)有7人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有3人回復(fù)
-
263求調(diào)劑
已經(jīng)有4人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
關(guān)于Matlab求孔隙率(找最佳閾值)
已經(jīng)有6人回復(fù)
- MATLAB將三維圖和其中一個(gè)面的等高線畫(huà)出來(lái)的問(wèn)題，謝謝了。
已經(jīng)有7人回復(fù)
- 求助一個(gè)matlab關(guān)于分類討論的問(wèn)題
已經(jīng)有9人回復(fù)
- 請(qǐng)教大家一個(gè)反卷積的問(wèn)題~~幫忙推導(dǎo)或者matlab編程計(jì)算~~
已經(jīng)有5人回復(fù)
- 一篇關(guān)于最優(yōu)控制的matlab tutorial
已經(jīng)有215人回復(fù)
- 菜鳥(niǎo)求助~matlab怎么樣讓一個(gè)1xn的向量中的每個(gè)元素是個(gè)矩陣？
已經(jīng)有11人回復(fù)
- 求助：能不能幫忙寫(xiě)一個(gè)用matlab求一個(gè)函數(shù)的傅里葉變換的程序
已經(jīng)有16人回復(fù)
- 關(guān)于matlab 2010b 的一個(gè)小問(wèn)題
已經(jīng)有7人回復(fù)
- 【求助】matlab求解一個(gè)簡(jiǎn)單的積分（急！�。。�
已經(jīng)有7人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題
已經(jīng)有21人回復(fù)
- 【求助】一個(gè)matlab窗口占25%CPU，如何提高？
已經(jīng)有11人回復(fù)
- 【求助】請(qǐng)教一個(gè)matlab程序，謝謝幫忙！！
已經(jīng)有15人回復(fù)
- 【求助】如何在MATLAB中作出一個(gè)隨機(jī)平面，急急急�。�！
已經(jīng)有3人回復(fù)
- 【求助】matlab中如何將程序的結(jié)果輸出到一個(gè)文件中
已經(jīng)有11人回復(fù)
- 【求助】MATLAB多個(gè)mat中 數(shù)據(jù)畫(huà)圖問(wèn)題
已經(jīng)有7人回復(fù)
- 【求助】matlab 二次規(guī)劃的優(yōu)化的問(wèn)題
已經(jīng)有4人回復(fù)
- 【求助】問(wèn)一個(gè)關(guān)于matlab中l(wèi)egend的問(wèn)題
已經(jīng)有11人回復(fù)
- 【求助】matlab中保存計(jì)算數(shù)據(jù)的問(wèn)題
已經(jīng)有8人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 jjliu2012 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
|
[考研] 招08考數(shù)學(xué) +4 | laoshidan 2026-03-20 | 8/400 |
|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|