 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (6218)
- >考研 (2783)
- >導(dǎo)師招生 (2374)
- >文獻(xiàn)求助 (458)
- >碩博家園 (307)
- >休閑灌水 (250)
- >蟲(chóng)友互識(shí) (214)
- >考博 (179)
- >論文投稿 (106)
- >招聘信息布告欄 (102)
- >博后之家 (74)
- >基金申請(qǐng) (62)
- >綠色求助(高懸賞) (47)
- >攻關(guān)文獻(xiàn)(高獎(jiǎng)勵(lì)) (44)
- >找工作 (29)
- >公派出國(guó) (22)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 765 | 回復(fù): 3 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
niexianling兌換貴賓
|
[求助]
目標(biāo)函數(shù)出現(xiàn)錯(cuò)誤
|
|||
| 請(qǐng)各位給我看一下目標(biāo)函數(shù)lmsob是哪里的錯(cuò)誤以及怎樣修改，OBJVa.m的運(yùn)行總是提示lmsob目標(biāo)函數(shù)式編寫不正確，請(qǐng)具體詳述一下，急。 |
回復(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 : OBJVa.m
- 附件 2 : lmsob.m
2013-05-27 17:00:52, 1.44 K
2013-05-27 17:02:05, 2.16 K
» 猜你喜歡
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有65人回復(fù)
-
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有28人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有24人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 幫忙給看看這個(gè)matlab優(yōu)化函數(shù) 問(wèn)題
已經(jīng)有8人回復(fù)
- matlab程序問(wèn)題，總是出錯(cuò)，麻煩大家看看問(wèn)題在那里？
已經(jīng)有6人回復(fù)
- 一個(gè)區(qū)間搜索的c語(yǔ)言程序 出問(wèn)題了
已經(jīng)有8人回復(fù)
- 【求助】幫忙分析一下我的錯(cuò)誤【已完結(jié)】
已經(jīng)有6人回復(fù)
- 【求助】VS2005數(shù)值計(jì)算程序調(diào)試【急求�。。。〖咏饚拧�
已經(jīng)有9人回復(fù)
- 【求助】請(qǐng)教：隱式目標(biāo)函數(shù)和約束條件非線性規(guī)劃問(wèn)題
已經(jīng)有5人回復(fù)
月只藍(lán)
無(wú)蟲(chóng)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲(chóng)號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
3樓2013-05-27 21:14:47
月只藍(lán)
無(wú)蟲(chóng)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.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
niexianling: 金幣+100, ★★★很有幫助, 非常感謝您的指點(diǎn)，很專業(yè) 2013-05-28 21:11:15
感謝參與，應(yīng)助指數(shù) +1
niexianling: 金幣+100, ★★★很有幫助, 非常感謝您的指點(diǎn)，很專業(yè) 2013-05-28 21:11:15
|
其實(shí)，MATLAB提示你關(guān)于目標(biāo)函數(shù)的問(wèn)題，只是你程序存在的幾個(gè)問(wèn)題中的一個(gè)，解釋如下： 其一，目標(biāo)函數(shù)的問(wèn)題，你函數(shù)的命名有問(wèn)題，以致于兩個(gè)m文件邏輯關(guān)系混淆不清，其實(shí)m文件的命名有相當(dāng)?shù)淖杂�，后面的程序中的命名就是一個(gè)例子。 其二，定積分的求解定義自變量的問(wèn)題，int((1/(2*pi)^0.5)*exp((-1/2)*x_1^2),x_1,-inf,z_1)和int((1/(2*pi)^0.5)*exp((-1/2)*x_2^2),x_2,-inf,z_2)中的自變量x_1和x_2須用syms x_1 x_2來(lái)定義，不然會(huì)出現(xiàn)“??? Undefined function or variable 'x_1'”的錯(cuò)誤提示。 其三，定積分的積分積分區(qū)間問(wèn)題，t是一個(gè)51*1維的行向量，所以z_1和z_2也是51*1維的行向量，你想用int((1/(2*pi)^0.5)*exp((-1/2)*x_1^2),x_1,-inf,z_1)和int((1/(2*pi)^0.5)*exp((-1/2)*x_2^2),x_2,-inf,z_2)求兩次定積分，每次輸出51個(gè)結(jié)果，但是實(shí)際你help int 一下，就知道，int不能同時(shí)求多個(gè)積分區(qū)間的結(jié)果，至少用你這些的代碼是不可以的，你這樣的代碼系統(tǒng)會(huì)提示你： ??? Error using ==> maplemex Error, (in int) wrong number (or type) of arguments: invalid options or option values passed to definite integration. 這51個(gè)定積分單獨(dú)求出來(lái)，再組成一個(gè)向量就可以了，見(jiàn)后面的程序。 其四，在你的命名為lmsob的m文件，最后一行的end是完全沒(méi)必要的，好吧，這問(wèn)題不大，最大的問(wèn)題在于，倒數(shù)第二行的“OBJVa=c1*(A10/0.15*exp(0.15^2/2*0.15^2-(t-c2)/0.15)*int((1/(2*pi)^0.5)*exp((-1/2)*x_1^2),x_1,-inf,z_1))-c3*(A20/0.16*exp(0.17^2/2*0.16^2-(t-c4)/0.16)*int((1/(2*pi)^0.5)*exp((-1/2)*x_2^2),x_2,-inf,z_2));”，這個(gè)式子有點(diǎn)常，其實(shí)是由基本的三項(xiàng)： 一個(gè)列向量 乘以 一個(gè)行向量 乘以 一個(gè)行向量；先講一下，MATLAB中，一個(gè)行(列)向量 乘以一個(gè)行(列)向量會(huì)因維數(shù)不同，不能計(jì)算；一個(gè)行(列)向量 點(diǎn)乘(注意是點(diǎn)乘，“.*”)一個(gè)行(列)向量會(huì)得到一個(gè)新的行(列)向量，由于用到點(diǎn)乘，其實(shí)此時(shí)已經(jīng)是數(shù)組之間的運(yùn)算了，既然是數(shù)組之間的運(yùn)算，一個(gè)一個(gè)行(列)向量 點(diǎn)乘 一個(gè)列(行)向量，也是不允許的，不能計(jì)算；一個(gè)行向量乘以一個(gè)列向量會(huì)得到一個(gè)數(shù)，而反過(guò)來(lái)，即一個(gè)列向量乘以一個(gè)行向量會(huì)得到一個(gè)矩陣；現(xiàn)在問(wèn)題來(lái)了，你定義“OBJVa=……”中的第一項(xiàng)就是一個(gè)列向量，按照上面講的MATLAB矩陣數(shù)組運(yùn)算規(guī)律，這樣的計(jì)算結(jié)果只能一個(gè)矩陣或者一個(gè)列向量，而你最終的是想獲得20個(gè)數(shù)字的結(jié)果吧？如果最終是想獲得20個(gè)純數(shù)字的結(jié)果，參見(jiàn)后附的程序，如果想獲得是20個(gè)1*51維的列向量或者20個(gè)51*51維的矩陣的話，按照上面運(yùn)算規(guī)律修改程序即可。 %------------------------------第一個(gè)m文件---------------------------------------------------- function jisuanjieguo clear all;clc m=[0.481736644302983,2.94957203824238,1.09399774122023,3.75750613928427;1.06135672641442,2.91561309396085,1.33234054044775,3.89003886226546;1.21219908647450,3.15663543380302,0.881394167799156,3.73019860286579;1.38765080523568,3.23114464868989,1.27027161204492,4.03792623322128;0.359879059008655,3.05261473905062,2.01557021271726,3.96040054359488;0.502647855518203,3.29696111389266,1.39989825391603,3.96714112009155;0.312419499606609,3.10503492835515,2.18019902553465,3.71961900674725;0.105130356626851,2.99769305962854,2.41764853243688,4.00525331998188;0.544705219083041,3.04416613022435,2.24589323195766,3.90189256848580;0.128231451255275,2.94842190591994,1.13550109520063,3.74356865269532;0.763477472569916,2.98753737214792,0.545424442886775,3.96566006246573;1.06910543308776,2.86127649428987,0.126812771046420,3.78250244379277;1.49020841141549,3.14148153446344,1.75329214190687,3.86275707507808;0.0715845332951863,3.01569654054312,2.59471395951649,3.99104308227833;1.01382175380874,3.13320506401545,0.120408361347543,3.71376487137305;0.322224666332880,3.27335860572682,0.980276665474573,3.70803137591493;0.697928896073242,2.95730682116205,0.811642925398756,4.06080394821544;0.832650233602747,3.01431847984169,1.30039789571561,3.87632882721789;0.862088937844217,2.85931478435019,2.47468028467206,4.00851345874163;1.27468907746227,3.19224757408864,1.69120584393105,3.68661850606776;]; for i=1:20 OBJVa(i)=lmsob(m(i,1),m(i,2),m(i,3),m(i,4)) %計(jì)算目標(biāo)函數(shù)值 end OBJVa %-------------------------------------------------------------------------------------------------- %-------------------------------第二個(gè)m文件------------------------------------------------ function y=lmsob(c1,c2,c3,c4) %A10,A20都是一個(gè)51*1的列向量 syms x_1 x_2 A10=[1.65812050716062;1.37842188749411;1.32709929007728;1.75367065593514;2.51304929550515;2.79491385756236;1.07714921956842;2.66428084508187;2.13827457310178;1.94620054796526;2.35415605309299;1.31204317847664;1.51660859151163;2.20195161465751;2.10310222972196;1.76459748362657;2.06208279629935;2.07424058256854;1.07203440245290;2.65833257485091;1.48405130710365;1.57706857848243;2.87883422699138;1.90606734046389;2.07623645980896;2.27163482373316;1.77192220104947;2.11019745913312;1.74570261018346;1.22099926156193;2.07149834266615;1.09073995829449;2.35149412200773;1.98998423917103;1.71569344638625;1.14593753460427;2.04827142429856;1.75568474816854;2.54626241305846;2.37133833985307;1.79452517173189;1.94568581421309;1.66056195248652;1.51374245529544;2.94772481165852;1.47245169009091;2.92065621586633;2.11476472478847;2.33311868995606;1.75906126686106;1.57921864637486]; %第1個(gè)組分在51個(gè)波長(zhǎng)下的峰面積，51*1 A20=[4.55384241423182;4.89249075634115;4.85541503507305;5.19247597498081;5.45249297888724;4.28670916677997;5.93827455871261;4.07516990748543;5.78221317419875;4.36184853203795;5.62184092201618;5.22656041883435;4.55045860649714;4.61744333019453;5.08623878201647;5.77281803917650;4.61725879308014;4.98069032196413;5.94471324434871;5.51400301449397;4.39599644173578;5.98357804702743;5.20681939885120;4.41107033730235;4.77944226780546;5.75517971730936;5.74278513151668;5.61881126026832;5.54243601448795;4.39187816797730;5.74092819249945;5.09517257848863;5.38104415249068;4.50887012358101;5.34275530560960;5.79863621461785;5.01686763696078;5.97882849964250;4.85520498607941;4.48592494593789;5.69347349991529;4.81358361405434;4.99838477749102;5.33680256496134;5.28163362243688;5.61896565010412;5.24360419620764;5.91554962521572;4.98686106806965;5.96872459187650;4.30896308164311]; %第2個(gè)組分在51個(gè)波長(zhǎng)下的峰面積 t=0:0.2:10; %1*51的行向量 z_1=(t-c2)/0.15-(0.15/0.15); %1*51的行向量 z_2=(t-c4)/0.17-(0.17/0.16); %1*51的行向量 for i=1:51 int_1(i)=int((1/(2*pi)^0.5)*exp((-1/2)*x_1^2),x_1,-inf,z_1(i)); int_2(i)=int((1/(2*pi)^0.5)*exp((-1/2)*x_2^2),x_2,-inf,z_2(i)); end y=c1* A10'/0.15 *(( exp(0.15^2/2*0.15^2-(t-c2)/0.15) .* int_1) )'-... c3* A20'/0.16 *(( exp(0.17^2/2*0.16^2-(t-c4)/0.16) .* int_2) )'; %---------------------------------------------------------------------------------------------- |
2樓2013-05-27 21:01:47
月只藍(lán)
新蟲(chóng)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲(chóng)號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
4樓2013-05-27 22:48:41
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +7 | 太容易1018 2026-04-01 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 324求調(diào)劑 +10 | hanamiko 2026-03-26 | 12/600 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +5 | 哇呼哼呼哼 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +11 | 傾____蕭 2026-03-27 | 12/600 |
|
|
|
[論文投稿]
chinese chemical letters英文版投稿求助
 120+4
|
Yishengeryi 2026-03-30 | 5/250 |
|
|
|
[考研] 材料調(diào)劑 +10 | Eujd1 2026-03-31 | 11/550 |
|
|
|
[考研] 085410人工智能 初試316分 求調(diào)劑 +3 | 殘星拂曙 2026-03-31 | 3/150 |
|
|
|
[碩博家園] 博一被送出聯(lián)培感覺(jué)不適應(yīng)怎么辦 +3 | 全村的狗 2026-03-31 | 3/150 |
|
|
|
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +4 | 1段時(shí)間 2026-03-31 | 5/250 |
|
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +13 | 深V宿舍吧 2026-03-29 | 13/650 |
|
|
|
[考研] 070300化學(xué)354求調(diào)劑 +15 | 101次希望 2026-03-28 | 15/750 |
|
|
|
[考研] 求調(diào)劑 +8 | 11ggg 2026-03-30 | 8/400 |
|
|
|
[考研] 287求調(diào)劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +6 | 星沉uu 2026-03-26 | 7/350 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11 | 景隨風(fēng) 2026-03-29 | 16/800 |
|
|
|
[考研] 262求調(diào)劑 +7 | ZZ..000 2026-03-30 | 8/400 |
|
|
|
[考研] 328求調(diào)劑 +8 | 嗯滴的基本都 2026-03-27 | 8/400 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|