| 查看: 607 | 回复: 2 | |||
wlx0914金虫 (小有名气)
|
[交流]
【求助】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) |
回复此楼» 猜你喜欢
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有11人回复
-
重要来源：本周末出结果
已经有5人回复
-
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 基金委网上报到问题
已经有8人回复
- 求助几个仪器分析中基本概念的问题
已经有7人回复
- 手征参数与碳纳米管螺旋角和直径的关系
已经有10人回复
- 偶联剂570的用法问题
已经有14人回复
- 急！ 关于Ms从Windows提交到Linux计算的问题？
已经有3人回复
- 做消费行为研究，两个问题得分相乘得出来的分数是不是连续变量？
已经有4人回复
- R语言图书-求解差分方程和并行计算！
已经有43人回复
- 求助一个纠结的问题？？？？
已经有4人回复
- 在职考博，有个问题咨询，谢谢大家了
已经有15人回复
- LC-MS 求解
已经有23人回复
- 二苯醚回收问题 求解
已经有6人回复
- 求解实验室发酵酸奶分层了怎么回事？
已经有10人回复
- 物理化学问题，求解！
已经有7人回复
- 质量监督系统考试 色弱可以报吗求解
已经有4人回复
- 二元幂函数参数求解
已经有10人回复
- 【求助】fminsearchi求解参数
已经有3人回复
- λh方程的参数如何求解（一天之内回答加40个金币）
已经有2人回复
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) |
2楼2010-11-10 17:17:33
lidaxue
木虫 (正式写手)
之乎者也
- 应助: 10 (幼儿园)
- 金币: 1952.1
- 散金: 612
- 红花: 2
- 帖子: 527
- 在线: 856.1小时
- 虫号: 259962
- 注册: 2006-06-17
- 性别: GG
- 专业: 几何学
3楼2010-11-10 17:31:53
