 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 395 | 回復(fù): 3 | |||
wlx0914金蟲 (小有名氣)
|
[交流]
【求助】fminsearchi求解參數(shù)
|
|
根據(jù)xx和T的數(shù)據(jù)，關(guān)聯(lián)方程參數(shù)a,b,c ，方程為 ln xx= a+b/T+c*ln T 下面是我的程序，提示有錯(cuò)，請(qǐng)高手指教 function y=myfun() clear all clc ExpData =[292.15 0.00848105 0.200113301 296.09 0.009329695 0.200113301 300.09 0.010580997 0.200113301 304.09 0.011993794 0.200113301 308.29 0.013486483 0.200113301 312.49 0.015286746 0.200113301 316.9 0.017209812 0.2001133 321.39 0.019224455 0.200113301 325.59 0.021446333 0.2001133 329.89 0.023818241 0.200113301 334.13 0.026117718 0.200113301 338.08 0.028576812 0.200113301 341.98 0.031254491 0.200113301 345.86 0.033917449 0.200113301 ]; T1= ExpData(:,1); xx1= ExpData(:,2); T=T1';xx=xx1'; %用fminsearch()—求解非線性Nelder-Mead單純形法（非線性數(shù)據(jù)擬合）問(wèn)題 A0= [1.22 1 1]; A=fminsearch(@Eqs,A0,[],T,xx) function f =Eqs(A,T,xx) % 模型方程 x1=exp(A(1)+A(2)/T+A(3)*log(T)); f=sum((xx-x1).^2) |
回復(fù)此樓» 猜你喜歡
環(huán)境工程 085701，267求調(diào)劑
已經(jīng)有5人回復(fù)
-
275求調(diào)劑
已經(jīng)有17人回復(fù)
-
334分 一志愿武理 材料求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑
已經(jīng)有15人回復(fù)
-
085600 286分 材料求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有10人回復(fù)
-
南京大學(xué)化學(xué)調(diào)劑
已經(jīng)有9人回復(fù)
-
085601材料工程找調(diào)劑
已經(jīng)有13人回復(fù)
-
考研調(diào)劑
已經(jīng)有6人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有13人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】fminsearchi求解參數(shù)問(wèn)題
已經(jīng)有2人回復(fù)
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲號(hào): 496517
- 注冊(cè): 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
2樓2010-11-08 10:45:31
change0618
鐵桿木蟲 (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲號(hào): 496517
- 注冊(cè): 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
★ ★ ★
wlx0914(金幣+15): 2010-11-08 14:36:07
robert2020(金幣+3):辛苦了！ 2010-11-08 23:00:36
wlx0914(金幣+15): 2010-11-08 14:36:07
robert2020(金幣+3):辛苦了！ 2010-11-08 23:00:36
|
function y=myfun() clear all clc ExpData =[292.15 0.00848105 0.200113301 296.09 0.009329695 0.200113301 300.09 0.010580997 0.200113301 304.09 0.011993794 0.200113301 308.29 0.013486483 0.200113301 312.49 0.015286746 0.200113301 316.9 0.017209812 0.2001133 321.39 0.019224455 0.200113301 325.59 0.021446333 0.2001133 329.89 0.023818241 0.200113301 334.13 0.026117718 0.200113301 338.08 0.028576812 0.200113301 341.98 0.031254491 0.200113301 345.86 0.033917449 0.200113301 ]; T1= ExpData(:,1); xx1= ExpData(:,2); T=T1';xx=xx1'; %用fminsearch()—求解非線性Nelder-Mead單純形法（非線性數(shù)據(jù)擬合）問(wèn)題 A0= [1.22 1 1]; A1=fminsearch(@Eqs,A0,[],T,xx); fprintf('\n\tA1 = \t%.4f\t%.4f\t%.4f\n',A1) options=optimset('TolFun',1e-10,'display','off'); A2= nlinfit(T,xx,@model,A0,options); TT=290:346; XX1=model(A1,TT); XX2=model(A2,TT); figure(1) plot(T,xx,'ro',TT,XX1,'b') grid on; figure(2) plot(T,xx,'ro',TT,XX2,'b') grid on; fprintf('\n\tA2 = \t%.4f\t%.4f\t%.4f\n',A2) function y=model(A,T) y=exp(A(1)+A(2)./T+A(3)*log(T)); function f =Eqs(A,T,xx) % 模型方程 x1=model(A,T); f=sum((xx-x1).^2) ; |
3樓2010-11-08 10:52:21
wlx0914
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1187.3
- 散金: 76
- 紅花: 4
- 帖子: 180
- 在線: 38.6小時(shí)
- 蟲號(hào): 903494
- 注冊(cè): 2009-11-15
- 性別: GG
- 專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)
4樓2010-11-08 14:35:43
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wlx0914 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 297求調(diào)劑 +11 | 田洪有 2026-03-26 | 11/550 |
|
|---|---|---|---|---|
|
[基金申請(qǐng)] 面上5B能上會(huì)嗎？ +3 | redcom 2026-03-29 | 3/150 |
|
|
|
[考研]
 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。
+10
|
nnnnnnn5 2026-03-25 | 10/500 |
|
|
|
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +8 | Koxui 2026-03-28 | 8/400 |
|
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 279 分 求調(diào)劑 +4 | 睡個(gè)好覺(jué)_16 2026-03-24 | 4/200 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪！？！ 2026-03-25 | 7/350 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|