 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 396 | 回復(fù): 3 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
wlx0914金蟲(chóng) (小有名氣)
|
[交流]
【求助】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ù)此樓» 猜你喜歡
面上5B能上會(huì)嗎？
已經(jīng)有4人回復(fù)
-
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）
已經(jīng)有9人回復(fù)
-
322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。
已經(jīng)有4人回復(fù)
-
085600 材料與化工 329分求調(diào)劑
已經(jīng)有14人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有5人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有21人回復(fù)
-
26申博自薦
已經(jīng)有5人回復(fù)
-
296求調(diào)劑
已經(jīng)有5人回復(fù)
-
總分293求調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】fminsearchi求解參數(shù)問(wèn)題
已經(jīng)有2人回復(fù)
change0618
鐵桿木蟲(chóng) (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲(chóng)號(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
change0618
鐵桿木蟲(chóng) (著名寫手)
方丈大師
- 應(yīng)助: 44 (小學(xué)生)
- 金幣: 17724.5
- 紅花: 17
- 帖子: 2413
- 在線: 546.7小時(shí)
- 蟲(chóng)號(hào): 496517
- 注冊(cè): 2008-01-19
- 專業(yè): 化學(xué)反應(yīng)工程
2樓2010-11-08 10:45:31
wlx0914
金蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1187.3
- 散金: 76
- 紅花: 4
- 帖子: 180
- 在線: 38.6小時(shí)
- 蟲(chóng)號(hào): 903494
- 注冊(cè): 2009-11-15
- 性別: GG
- 專業(yè): 化工熱力學(xué)和基礎(chǔ)數(shù)據(jù)
4樓2010-11-08 14:35:43
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。 +4 | XX小鄧 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考博] 26申博自薦 +5 | whh869393 2026-03-24 | 5/250 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +12 | 阿斯蒂芬2004 2026-03-28 | 12/600 |
|
|
|
[考研] 332求調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研]
 材料與化工考研調(diào)劑
+17
|
孅華 2026-03-22 | 17/850 |
|
|
|
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研]
 085404電子信息284分求調(diào)劑
+4
|
13659058978 2026-03-24 | 4/200 |
|
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|