 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 556 | 回復(fù): 2 | |||
wlx0914金蟲 (小有名氣)
|
[交流]
【求助】fminsearchi求解參數(shù)問題
|
|
已知數(shù)據(jù)T，xx, xb, 想要把xx表示成以下方程： lnxx=a+b/T+c*lnT+d*xb+e*xb^2+f*xb^3; 根據(jù)數(shù)據(jù)求方程參數(shù)a到e, 目標(biāo)函數(shù)為sum((T-TCAL).^2)最小，TCAL為計(jì)算值 自己編的程序有錯(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 292.15 0.016489969 0.399943158 297.09 0.017925164 0.399943158 300.79 0.019553686 0.399943158 305.09 0.0215646 0.399943158 309.99 0.024057562 0.399943158 314.39 0.026586429 0.399943158 318.29 0.029056351 0.399943158 322.49 0.03187043 0.399943157 326.69 0.035019471 0.399943158 330.92 0.038027106 0.399943158 335.29 0.041484872 0.399943158 340.5 0.045715199 0.399943158 344.8 0.04958667 0.399943157 348.09 0.052760705 0.399943158 ]; T1= ExpData(:,1); xx1= ExpData(:,2); xb1= ExpData(:,3); T=T1';xx=xx1';xb=xb1'; %用fminsearch()—求解非線性Nelder-Mead單純形法（非線性數(shù)據(jù)擬合）問題 A0= [-76.4362523947615 1917.34072551572 11.5199849169552 0.0858354183915545 2.96421571145088 -1.90592987130477 ]; options=optimset('MaxFunEvals',1e+008,'MaxIter',1e+008) A=fminsearch(@obj,A0,options,T,xx,xb) function f =Eqs(T,A,xx,xb) % 模型方程 f=xx-exp(A(1)+A(2)./T+A(3).*log(T)+A(4).*xb+A(5).*xb.^2+A(6).*xb.^3); function f=obj(A,T,xx,xb) %目標(biāo)函數(shù)方程 n=length(T); for i=1:n x=fzero(@Eqs,T(i),[],xx(i),xb(i),A(1),A(2),A(3),A(4),A(5),A(6)); Tcal(i)=x; end f=sum((T-TCAL).^2) |
回復(fù)此樓» 猜你喜歡
材料調(diào)劑
已經(jīng)有7人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有10人回復(fù)
-
289求調(diào)劑
已經(jīng)有15人回復(fù)
-
291 求調(diào)劑
已經(jīng)有7人回復(fù)
-
274求調(diào)劑
已經(jīng)有14人回復(fù)
-
北京林業(yè)大學(xué)碩導(dǎo)招生廣告
已經(jīng)有3人回復(fù)
-
309求調(diào)劑
已經(jīng)有5人回復(fù)
-
292求調(diào)劑
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿 西北大學(xué) 總分282 英語一62 求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 基金委網(wǎng)上報(bào)到問題
已經(jīng)有8人回復(fù)
- 求助幾個(gè)儀器分析中基本概念的問題
已經(jīng)有7人回復(fù)
- 手征參數(shù)與碳納米管螺旋角和直徑的關(guān)系
已經(jīng)有10人回復(fù)
- 偶聯(lián)劑570的用法問題
已經(jīng)有14人回復(fù)
- 急！ 關(guān)于Ms從Windows提交到Linux計(jì)算的問題？
已經(jīng)有3人回復(fù)
- 做消費(fèi)行為研究，兩個(gè)問題得分相乘得出來的分?jǐn)?shù)是不是連續(xù)變量？
已經(jīng)有4人回復(fù)
- R語言圖書-求解差分方程和并行計(jì)算！
已經(jīng)有43人回復(fù)
- 求助一個(gè)糾結(jié)的問題？？？？
已經(jīng)有4人回復(fù)
- 在職考博，有個(gè)問題咨詢，謝謝大家了
已經(jīng)有15人回復(fù)
- LC-MS 求解
已經(jīng)有23人回復(fù)
- 二苯醚回收問題 求解
已經(jīng)有6人回復(fù)
- 求解實(shí)驗(yàn)室發(fā)酵酸奶分層了怎么回事？
已經(jīng)有10人回復(fù)
- 物理化學(xué)問題，求解！
已經(jīng)有7人回復(fù)
- 質(zhì)量監(jiān)督系統(tǒng)考試 色弱可以報(bào)嗎求解
已經(jīng)有4人回復(fù)
- 二元冪函數(shù)參數(shù)求解
已經(jīng)有10人回復(fù)
- 【求助】fminsearchi求解參數(shù)
已經(jīng)有3人回復(fù)
- λh方程的參數(shù)如何求解（一天之內(nèi)回答加40個(gè)金幣）
已經(jīng)有2人回復(fù)
lidaxue
木蟲 (正式寫手)
之乎者也
- 應(yīng)助: 10 (幼兒園)
- 金幣: 1952.1
- 散金: 612
- 紅花: 2
- 帖子: 527
- 在線: 856.1小時(shí)
- 蟲號(hào): 259962
- 注冊(cè): 2006-06-17
- 性別: GG
- 專業(yè): 幾何學(xué)
|
這是運(yùn)行樓主的程序產(chǎn)生的錯(cuò)誤信息 ??? Error using ==> fzero at 317 FZERO cannot continue because user supplied function_handle ==> Eqs failed with the error below. Too many input arguments. Error in ==> obj at 4 x=fzero(@Eqs,T(i),[],xx(i),xb(i),A(1),A(2),A(3),A(4),A(5),A(6)); Error in ==> fminsearch at 205 fv(:,1) = funfcn(x,varargin{:}); Error in ==> myfun at 39 A=fminsearch(@obj,A0,options,T,xx,xb) |
2樓2010-11-10 17:17:33
lidaxue
木蟲 (正式寫手)
之乎者也
- 應(yīng)助: 10 (幼兒園)
- 金幣: 1952.1
- 散金: 612
- 紅花: 2
- 帖子: 527
- 在線: 856.1小時(shí)
- 蟲號(hào): 259962
- 注冊(cè): 2006-06-17
- 性別: GG
- 專業(yè): 幾何學(xué)
3樓2010-11-10 17:31:53
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wlx0914 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +6 | 呼呼？~+123456 2026-03-24 | 6/300 |
|
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 11/550 |
|
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +4 | 13831862839 2026-03-24 | 5/250 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +6 | Stella_Yao 2026-03-20 | 6/300 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
|
[考研]
 材料考研調(diào)劑生
+3
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +5 | leaflight 2026-03-24 | 5/250 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-24 | 5/250 |
|
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|