| 查看: 520 | 回复: 1 | |||
[交流]
【求助】matlab 拟合问题
|
|
自变量t是时间,一维数组; 因变量c是浓度,包含2种,一个是料液相浓度,另一个是接受相浓度,二维数组。 c1=g(K1,K2,t), c2=G(K1,K2,t); g和G是不同的函数。 现在有实验数据t和c=[c1 c2]; 请问如何在matlab中编程对参数K1,K2进行非线性拟合? 注意,K1,K2在c1,c2中都出现了,所以2个函数需要同时拟合。我在用nlinfit函数时发现好像只能对一个因变量进行拟合。 敬请赐教!谢谢 |
» 猜你喜欢
Sci. Bull. 悲剧经验
已经有5人回复
国自然评分
已经有4人回复
今年审到国自然15份,谈谈感受
已经有19人回复
提交了我也来说说感想
已经有4人回复
面上本子正文33页,违规吗?会被低分嘛?
已经有11人回复
评审有感
已经有22人回复
反应很差,大量原料没有反应
已经有6人回复
如果工作了想读博,可以边工作边读全日制嘛?
已经有5人回复
如何实现卤原子转化
已经有8人回复
河北省自然科学基金
已经有3人回复
» 本主题相关价值贴推荐,对您同样有帮助:
请教matlab反应动力学参数估计遇到的问题,谢谢
已经有15人回复
matlab拟合方程参数时初值的选择
已经有15人回复
matlab拟合求值
已经有17人回复
在matlab中拟合曲面的表达式
已经有29人回复
matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
matlab拟合拟均相动力学参数
已经有6人回复
需要用matlab拟合数据的外行人急需帮助
已经有10人回复
matlab非线性参数拟合问题
已经有7人回复
【求助】matlab曲面拟合表达式
已经有10人回复
【求助】使用Matlab拟合反应动力学方程问题
已经有7人回复
» 抢金币啦!回帖就可以得到:
祈福,面上申请连续三年了,希望今年面上能中!!!
+1/953
12年,离婚了。
+1/298
坐标深圳,愿遇良缘
+1/224
山东征女友,坐标济南
+1/162
海南医科大学肖华教授(国家级青年人才)招收博士生1名---报名从速
+4/92
ASTM F1980-21
+1/90
2026江西师范大学人文地理与低碳发展拟招1名“申请审核制”学术博士研究生(支持出国)
+1/84
钯炭、铂炭、钌炭等加氢催化剂
+1/80
北京科技大学能源与环境工程学院王庆功教授团队,招收全日制博士生1-2名,5月21日前
+1/79
电力全国重点实验室双一流A类长江学者团队补招26年全日制学术博士1名(5月19日截止)
+1/77
80后老木虫征婚
+1/56
招聘二维光电材料与器件方向博士后若干名
+1/32
安徽农业大学农林经济管理交叉课题组急招农林经济管理26级全日制博士
+1/32
双一流南京林业大学/化学工程学院-国家海外优青团队招青年师资、师资博后(自动转编)
+1/15
培养聚乙烯和聚丙烯球晶的温度
+1/9
合成
+1/8
知名外资仪器公厂家招技术销售
+1/4
南方科技大学物理系何佳清讲席教授课题组2026年招收热电或电镜方向博士生1名
+1/3
武汉科大周志辉教授招青年教师、博士后和实验助理(碳/高分子/富氧/PV/MOF/COF膜)
+1/1
香港城市大学校长奖学金/香港政府奖学金博士生招聘
+1/1
flyinuniverse58(金币+10):虽然和我的目标有些距离,但还是感谢你为我提供了思路(我在题目中说了y是2个不同的函数,不是x) 2010-11-19 14:50:25
flyinuniverse58(金币+10):没人回复了,都给你吧,谢谢 2010-12-23 08:10:11
flyinuniverse58(金币+10):没人回复了,都给你吧,谢谢 2010-12-23 08:10:11
|
请楼主 依葫芦画瓢,保证算的出来!! 1 设定目标函数. (M函数书写)% 可以是任意的 例如: function f=mydata(a,data) %y的值目标函数值 或者是第三维的,a=[a(1) ,a(2)] 列向量 x=data(1, ; %data 是一2维数组,x=x1y=data(2, ; %data 是一2维数组,x=x2 f=a(1)*x+a(2)*x.*y; %这里的a(1), a(2)为目标函数的系数值。 f的值相当于ydata的值 2 然后给出数据xdata和ydata的数据和拟合函数lsqcurvefit 例如: x1=[1.0500 1.0520 1.0530 1.0900 1.0990 1.1020 1.1240 1.1420... 1.1490 1.0500 1.0520 1.0530 1.0900 1.0990 1.1020 1.1240 1.1420 1.1490]; x2=[3.8500 1.6500 2.7500 5.5000 7.7000 3.3000 4.9500 8.2500 11.5500... 1.6500 2.7500 3.8500 7.7000 3.3000 5.5000 8.2500 11.5500 4.9500]; ydata=[56.2000 62.8000 62.2000 40.8000 61.4000 57.5000 44.5000 54.8000... 53.9000 64.2000 62.9000 64.1000 63.0000 62.2000 64.2000 63.6000... 52.5000 62.0000]; data=[x1;x2]; %类似于将x1 x2整合成一个2维数组。 a0= [-0.0014,0.07]; option=optimset('MaxFunEvals',5000); format long; [a,resnorm]=lsqcurvefit(@mydata,a0,data,ydata,[],[],option); yy=mydata(a,data); result=[ydata' yy' (yy-ydata)'] % a的值为拟合的目标函数的参数值 利用lsqcurvefit进行拟合的 它完整的语法形式是: % [x,resnorm,residual,exitflag,output,lambda,jacobian] =lsqcurvefit(fun,x0,xdata,ydata,lb,ub,options) |
2楼2010-11-19 10:19:10












回复此楼
; %data 是一2维数组,x=x1