| 查看: 604 | 回复: 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) |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
听说今天filecode变了
已经有33人回复
-
基金中了
已经有6人回复
-
filecode
已经有14人回复
-
固定端突然变了，今天
已经有7人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 基金委网上报到问题
已经有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
50