24小时热门版块排行榜

返回列表

wlx0914

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1187.3
散金: 76
红花: 4
帖子: 180
在线: 38.6小时
虫号: 903494
注册: 2009-11-15
性别: GG
专业: 化工热力学和基础数据

[交流] 【求助】fminsearchi求解参数问题

已知数据T，xx, xb, 想要把xx表示成以下方程：
lnxx=a+b/T+c*lnT+d*xb+e*xb^2+f*xb^3;
根据数据求方程参数a到e,
目标函数为sum((T-TCAL).^2)最小，TCAL为计算值
自己编的程序有错，请高手指教

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单纯形法（非线性数据拟合）问题
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)    %目标函数方程
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)

回复此楼

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有8人回复
最失望的一年已经有10人回复
求推荐英文EI期刊已经有5人回复
存款400万可以在学校里躺平吗已经有27人回复
请教限项目规定已经有4人回复
国自然申请面上模板最新2026版出了吗？已经有20人回复
26申博已经有3人回复
基金委咋了？2026年的指南还没有出来？已经有10人回复
基金申报已经有6人回复
疑惑？已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2010-11-08 20:30:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lidaxue

木虫 (正式写手)

之乎者也

应助: 10 (幼儿园)
金币: 1952.1
散金: 612
红花: 2
帖子: 527
在线: 856.1小时
虫号: 259962
注册: 2006-06-17
性别: GG
专业: 几何学

这是运行楼主的程序产生的错误信息
??? 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)