给出一组数据_如何同matlab中提供的工具箱拟合一条曲线,[MATLAB工具箱] 怎样用一组参数同时拟合两个曲线..._Mike昊的博客-程序员宅基地

怎样用一组参数同时拟合两个曲线

如题，小弟想拟合曲线r1和r2，这两条曲线都用的同一组参数A1 A2 A3 A4 A5 A6 A7 Aw E1 E2 E3 E4 E5 E6 E7 Ew进行非线性最小二乘拟合，请问怎样用matlab实现一组参数同时拟合两个方程？(如果分开拟合的话每个参数必然有两个值)。下面是相关的数据及方程，其中y1和y2分别是r1和r2的实验值，r1和r2的自变量值相同，均为x1.请各位大侠能鼎力相助，最好能给出代码或者计算方法，谢谢

x1=[0.51333,0.36301,0.34194,0.28502,0.26485,0.17771];

y1=[0.50536,0.61684,0.63566,0.68357,0.70947,0.72182];

y2=[2.5363,1.16376,0.75531,0.4164,0.27812,0.1259];

syms A1 A2 A3 A4 A5 A6 A7 Aw E1 E2 E3 E4 E5 E6 E7 Ew x;

k1=A1*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E1)/4098.802);

k2=A2*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E2)/4098.802);

k3=A3*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E3)/4098.802);

k4=A4*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E4)/4098.802);

k5=A5*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E5)/4098.802);

k6=A6*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E6)/4098.802);

k7=A7*exp ((-Aw*exp (-Ew/4098.802)*((0.133+0.173*x)/(0.483-0.346*x))-E7)/4098.802);

r1=-(x*k1*k2*k3*(-k5-k7))/(-(x*k1*k2-k3*(-x*k1-k2-x*k4-k6))*(-k5-k7)-x*k3*k4*(-x*k1+k7));

r2=-(x^2)*k1*k3*k4*k5/((x*k1*k2-k3*(-x*k1-k2-x*k4-k6))*(-k5-k7)-x*k3*k4*(-x*k1+k7));

100 <= A1 <= 10000;

10^11 <= A2 <= 10^14;

10^12 <= A3 <= 10^14;

100 <= A4 <= 10000;

10^12 <= A5 <= 10^14;

10^12 <= A6 <= 10^14;

10^12 <= A7 <= 10^14;

100 <= Aw <= 10000;

0 <= E1 <= 100;

50 <= E2 <= 140;

50 <= E3 <= 140;

0 <= E4 <= 100;

50 <= E5 <= 150;

100 <= E6 <= 150;

0 <= E7 <= 100;

0 <= Ew <= 100;

本文链接：https://blog.csdn.net/weixin_35645230/article/details/116106794

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

TP6 Auth权限使用_tp6 auth权限管理_人海中的海盗的博客-程序员宅基地

ThinkPHP6 Auth权限管理

linux 下多线程错误 undefined reference to `sem_init'_xqhrs232的博客-程序员宅基地

原文地址::http://www.cnblogs.com/stli/archive/2010/03/11/1683823.htmlundefined reference to `sem_init'undefined reference to `sem_post'undefined reference to `sem_wait'编译选项需要加入一个多线程g

【已解决】Android Studio无法识别虚拟机，虚拟机无法运行_Biligle的博客-程序员宅基地

本人很是费解，为什么AS能够检查出我的虚拟机，但是虚拟机不能运行，我检查后台进程之后，猛然发现后台居然没有进程！！！解决：然后我把电脑里的豌豆荚、金山手机助手、360手机助手卸载掉，把安全管家金山杀毒全部关掉，再次插入手机，虚拟机顺利运行了！！！这是为什么呢？也许是因为以上的软件，占用了我手机adb的进程，让我的虚拟机无法运行了。

'vue' 不是内部或外部命令，也不是可运行的程序或批处理文件_拦路雨g的博客-程序员宅基地

当各位小伙伴准备使用vue cli搭建第一个webpack脚手架时，是否会出现：'vue' 不是内部或外部命令，也不是可运行的程序或批处理文件。直接来解决方案：主要原因是当前环境无法找到vue-cli的路径；1. 使用下列命令查找prefix配置路径：npm config listprefix配置的路径就是npm全局下载的路径；我们使用npm install -g v...

IOS XCode 控制台打印中文方法_ios控制台打印汉子框架_qq910441958的博客-程序员宅基地

①首先定义一个打印方法：//打印输出#ifdef DEBUG#define DDLog(format, ...) printf("class: &lt;%p %s:(%d) &gt; method: %s \n%s\n", self, [[[NSString stringWithUTF8String:__FILE__] lastPathComponent] UTF8String], ...

java arff文件_怎么用java将sql查询结果存成weka的标准.arff文件？_weixin_39574943的博客-程序员宅基地

importweka.core.Instances;importweka.experiment.InstanceQuery;...InstanceQueryquery=newInstanceQuery();query.setDatabaseURL("jdbc_url");query.setUsername("the_user");query.setPassword("the_passw...