| 查看: 610 | 回复: 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) |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 基金委网上报到问题
已经有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
