24小时热门版块排行榜

返回列表

wlx0914

金虫 (小有名气)

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

[交流] 【求助】关于fminsearch和fsolve参数估计

传热实验的努塞尔数、雷诺数及普朗特数之间的拟合问题Nu=c1*Re^c2*Pr^c3.
已知nu=[1.127,2.416,2.205,2.312,1.484,6.038,7.325];
re=[100,200,300,500,100,700,800];
pr=[2,4,1,0.3,5,3,4];
用fminsearch求参数使计算值与实验Nu值差的平方和最小。下面是我的程序
function f=nurepr(A,nu,re,pr) %子程序
nu1=A(1)*re.^A(2).*pr.^A(3); %努塞尔数表达式
f=sum((nu-nu1).^2)；%目标函数
命令窗口
nu=[1.127,2.416,2.205,2.312,1.484,6.038,7.325];
re=[100,200,300,500,100,700,800];
pr=[2,4,1,0.3,5,3,4];
A0=[1,1,1]; %初值
A=fminsearch('nurepr',A0,[],nu,re,pr); %寻找f的最小值，并返回参数A
这样能够算出参数，现在想用fsolve函数解出nu1,然后再代入目标函数的方法求解，请高手指教

回复此楼