| 查看: 632 | 回复: 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) |
回复此楼» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 基金委网上报到问题
已经有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
20