 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1681 | 回復(fù): 8 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
matangtang木蟲(chóng) (正式寫手)
|
[求助]
matlab優(yōu)化報(bào)錯(cuò) 已有1人參與
|
||||
|
程序在附件中，其中dynamics_inverse_qiujie.m為優(yōu)化主程序，biaozhun_function.m為優(yōu)化目標(biāo)函數(shù)，yueshu.m為約束函數(shù)。 程序里面的初值a0已經(jīng)滿足我的需要了（是手調(diào)得到的），但仍然存在一個(gè)錯(cuò)誤。如圖片所示 yueshu.m中計(jì)算出來(lái)的g的個(gè)數(shù)本來(lái)是10個(gè)，最后一次計(jì)算變成了11個(gè)，如下圖所示。 還請(qǐng)高手解答 因?yàn)椴荒茉谠N上編輯，所以又重新發(fā)了這個(gè)求助帖，原帖在。如果能解決問(wèn)題，金幣可以都奉送哈 http://www.gaoyang168.com/bbs/viewthread.php?tid=8413221&page=2 1.png 2.png |
回復(fù)此樓» 猜你喜歡
287求調(diào)劑
已經(jīng)有7人回復(fù)
-
325求調(diào)劑
已經(jīng)有5人回復(fù)
-
343求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑推薦 材料 304
已經(jīng)有4人回復(fù)
-
316求調(diào)劑
已經(jīng)有4人回復(fù)
-
351求調(diào)劑
已經(jīng)有4人回復(fù)
-
349求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料學(xué)碩333求調(diào)劑
已經(jīng)有7人回復(fù)
-
342求調(diào)劑
已經(jīng)有3人回復(fù)
-
調(diào)劑求收留
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 使用matlab最優(yōu)化方法擬合獲得多個(gè)動(dòng)力學(xué)參數(shù)中的問(wèn)題
已經(jīng)有4人回復(fù)
- 菜鳥(niǎo)求助，老師給了一個(gè)分子讓優(yōu)化，可總是錯(cuò)誤，求大神指點(diǎn)
已經(jīng)有22人回復(fù)
- Matlab心得及學(xué)習(xí)方法（不斷更新）
已經(jīng)有26人回復(fù)
- matlab與comsol斷開(kāi)連接
已經(jīng)有10人回復(fù)
- matlab求函數(shù)積分出錯(cuò)
已經(jīng)有12人回復(fù)
- workbench計(jì)算報(bào)錯(cuò)
已經(jīng)有4人回復(fù)
- matlab一運(yùn)行就出現(xiàn)的問(wèn)題
已經(jīng)有7人回復(fù)
- 求教matlab線性優(yōu)化求最小值的問(wèn)題
已經(jīng)有7人回復(fù)
- 請(qǐng)教MATLAB中的LMI求解問(wèn)題，急�。�！
已經(jīng)有7人回復(fù)
- 求助：遺傳算法的matlab編程求解資源優(yōu)化問(wèn)題
已經(jīng)有7人回復(fù)
- 求助：資源優(yōu)化問(wèn)題的遺傳算法matlab編程
已經(jīng)有5人回復(fù)
- 求助matlab問(wèn)題
已經(jīng)有5人回復(fù)
- 求助有關(guān)Matlab有約束非線性最優(yōu)化問(wèn)題
已經(jīng)有9人回復(fù)
- Matlab兩個(gè)小問(wèn)題
已經(jīng)有10人回復(fù)
- 【求助】拜求高手解釋下matlab的報(bào)錯(cuò)-----“out of memory”
已經(jīng)有13人回復(fù)
- 【求助】matlab 二次規(guī)劃的優(yōu)化的問(wèn)題
已經(jīng)有4人回復(fù)
- 【分享】推薦一本好書《精通MATLAB最優(yōu)化計(jì)算》
已經(jīng)有880人回復(fù)
月只藍(lán)
主管區(qū)長(zhǎng) (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68127.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲(chóng)號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
matangtang: 金幣+50, ★★★★★最佳答案 2015-01-15 10:41:09
感謝參與，應(yīng)助指數(shù) +1
matangtang: 金幣+50, ★★★★★最佳答案 2015-01-15 10:41:09
|
您好！ 我大致看了一下，發(fā)現(xiàn)主要問(wèn)題在于 yueshu.m 中 變量g 的維數(shù)問(wèn)題。 按原來(lái)的程序： g=abs(double(z))-0.99; 那么 變量z 的維數(shù)就決定了 g 的維數(shù) 而 z 是用solve 函數(shù)解方程組得到的： [A B C Q H z]=solve('A=1-1.1/z','B=1+1.2/z','C=1+0.2/z',['Q=((',num2str(a(9)),')*(z-((',num2str(a(1)),')+(',num2str(a(2)),')*i))*(z-((',num2str(a(3)),')+(',num2str(a(4)),')*i)))/((z-((',num2str(a(5)),')+(',num2str(a(6)),')*i))*(z-((',num2str(a(7)),')+(',num2str(a(8)),')*i)))-B*C'],'H=z*Q-B*C-A*Q+z*B*C','z^2-(C*(Q+B*C))*z/H-((Q+B*C)*(B*C^2-z*Q*C+A*Q*C-z*B*C^2))/H^2','A,B,C,Q,H,z') 由于solve 函數(shù)會(huì)求解方程組所有的解，那么問(wèn)題就來(lái)了，如果原方程組是多解的，而且解的個(gè)數(shù)不唯一，就會(huì)引起z維數(shù)的變化 通過(guò)跟蹤 優(yōu)化過(guò)程中 z 維數(shù)的變化可見(jiàn)（以下以index_xy表示z的維數(shù)）： index_xy = 6 1 index_xy = 6 1 index_xy = 6 1 index_xy = 6 1 index_xy = 6 1 index_xy = 6 1 index_xy = 11 1 可見(jiàn)第7次，z 的維數(shù)為11。 為了進(jìn)一步驗(yàn)證，變量g和z 的維數(shù)就是調(diào)通程序的關(guān)鍵，嘗試將變量z 改成以下形式： g=abs(double(z(1:6)))-0.99; 即強(qiáng)制 g和z 維數(shù)=6 最終程序可以跑通的，結(jié)果見(jiàn)附圖1。當(dāng)然結(jié)果肯定不好，不過(guò)這樣確實(shí)說(shuō)明，調(diào)整 g和z 維數(shù)，確實(shí)可以跑通程序。 基本上找到了原因，那么怎么解決呢？ 由于solve 函數(shù)會(huì)求解方程組所有的解，解的個(gè)數(shù)會(huì)隨著方程的變化而不同，所以是否可以給 z 中的元素加個(gè)約束，只取出其中合適的6個(gè)，來(lái)代入優(yōu)化過(guò)程�；蛘�，不用 solve函數(shù)，用 fsolve函數(shù)代入不同的初值，求取不同的解。需要說(shuō)明的是，以上只是我初步的想法，我也沒(méi)有想到如何用代碼具體實(shí)現(xiàn)。 附圖1.png |
4樓2015-01-15 09:33:08
NicoleLas
木蟲(chóng) (正式寫手)
- 應(yīng)助: 66 (初中生)
- 金幣: 4288.7
- 紅花: 26
- 帖子: 514
- 在線: 177.2小時(shí)
- 蟲(chóng)號(hào): 2872321
- 注冊(cè): 2013-12-15
- 性別: GG
- 專業(yè): 計(jì)算機(jī)軟件
2樓2015-01-14 23:02:25
matangtang
木蟲(chóng) (正式寫手)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 3131.5
- 散金: 856
- 紅花: 6
- 帖子: 614
- 在線: 278.4小時(shí)
- 蟲(chóng)號(hào): 1871554
- 注冊(cè): 2012-06-27
- 性別: GG
- 專業(yè): 制造系統(tǒng)與自動(dòng)化
3樓2015-01-15 00:33:35
matangtang
木蟲(chóng) (正式寫手)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 3131.5
- 散金: 856
- 紅花: 6
- 帖子: 614
- 在線: 278.4小時(shí)
- 蟲(chóng)號(hào): 1871554
- 注冊(cè): 2012-06-27
- 性別: GG
- 專業(yè): 制造系統(tǒng)與自動(dòng)化
5樓2015-01-15 10:48:28
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +3 | Li李魚(yú) 2026-03-26 | 3/150 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +9 | A-哆啦Z夢(mèng) 2026-03-23 | 15/750 |
|
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap－ 2026-03-22 | 5/250 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +5 | leaflight 2026-03-24 | 5/250 |
|
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +4 | 幸運(yùn)的醬醬 2026-03-22 | 5/250 |
|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|