| 3 | 1/1 | 返回列表 |
| 查看: 555 | 回復(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ò),請高手指教 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) |
木蟲 (正式寫手)
之乎者也
|
這是運(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) |

木蟲 (正式寫手)
之乎者也

| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +17 | 陌の森林 2026-03-18 | 17/850 |
|
|
[考研] 289求調(diào)劑 +9 | 懷瑾握瑜l 2026-03-20 | 9/450 |
|
|
[考研] 318求調(diào)劑 +5 | plum李子 2026-03-21 | 8/400 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +9 | 苑豆豆 2026-03-20 | 9/450 |
|
|
[考研] 求調(diào)劑 +6 | 研研,接電話 2026-03-24 | 7/350 |
|
|
[考研] 材料292調(diào)劑 +8 | 橘頌思美人 2026-03-23 | 8/400 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +3 | 你好明天你好 2026-03-23 | 3/150 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[考研]
求調(diào)劑材料學(xué)碩080500,總分289分
5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 315分,誠求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|