| 查看: 579 | 回复: 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函数时发现好像只能对一个因变量进行拟合。 敬请赐教！谢谢 |
回复此楼» 猜你喜欢
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
国自然评审意见
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有13人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
-
有没有仍没收到信息的
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有8人回复
-
基金不中，共勉
已经有11人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教matlab反应动力学参数估计遇到的问题，谢谢
已经有15人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- matlab拟合求值
已经有17人回复
- 在matlab中拟合曲面的表达式
已经有29人回复
- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- matlab拟合拟均相动力学参数
已经有6人回复
- 需要用matlab拟合数据的外行人急需帮助
已经有10人回复
- matlab非线性参数拟合问题
已经有7人回复
- 【求助】matlab曲面拟合表达式
已经有10人回复
- 【求助】使用Matlab拟合反应动力学方程问题
已经有7人回复
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
30