matlab fmincon函数格式,求助 fmincon 函数调用格式说明

matlab fmincon函数格式,求助 fmincon 函数调用格式说明functionKineticsEstclearallclck0=[0.50.50.50.50.5];%参数初值lb=[00000];%参数下限ub=[+inf+inf+inf+inf+inf];%参数上限x0=[0.18830.25070.0467…

大家好,又见面了,我是你们的朋友全栈君。

function KineticsEst

clear all

clc

k0 = [0.5  0.5  0.5  0.5  0.5];         % 参数初值

lb = [0  0  0  0  0];                   % 参数下限

ub = [+inf  +inf  +inf  +inf  +inf];    % 参数上限

x0 = [0.1883  0.2507  0.0467  0.0899  0.1804  0.1394  0.1046];

KineticsData1;

yexp = ExpData(:,2:5);                  % yexp: 实验数据[x1        x4        x5        x6]

% 使用函数fmincon()进行参数估计

[k,fval,flag,output] = fmincon(@ObjFunc4Fmincon,k0,[],[],[],[],lb,ub,[],[],x0,yexp);

fprintf(‘\n使用函数fmincon()估计得到的参数值为:\n’)

fprintf(‘\tk1 = %.4f\n’,k(1))

fprintf(‘\tk2 = %.4f\n’,k(2))

fprintf(‘\tk3 = %.4f\n’,k(3))

fprintf(‘\tk4 = %.4f\n’,k(4))

fprintf(‘\tk5 = %.4f\n’,k(5))

fprintf(‘  The sum of the squares is: %.1e\n\n’,fval)

k_fmincon = k;

% ——————————————————————

function f = ObjFunc4Fmincon(k,x0,yexp)

tspan = [0.00 : 0.01 : 0.20];

[t x] = ode45(@KineticEqs,tspan,x0,[],k);

y(:,1) = x(:,1);

y(:,2:4) = x(:,4:6);

f = sum((y(:,1)-yexp(:,1)).^2) + sum((y(:,2)-yexp(:,2)).^2)   …

+ sum((y(:,3)-yexp(:,3)).^2) + sum((y(:,4)-yexp(:,4)).^2);

% ——————————————————————

function dxdt = KineticEqs(t,x,k)

q = 8.75 + k(5);

dxdt =  …

[ ( k(5)-q*x(1)- k(1)*x(1)*x(2)-k(4)*x(1)*x(6)*sqrt(0.9) )

( 7.0-q*x(2) – k(1)*x(1)*x(2)-2*k(2)*x(2)*x(3) )

( 1.75 -q*x(3) – k(2)*x(2)*x(3) )

( -q*x(4) + 2*k(1)*x(1)*x(2)-k(3)*x(4)*x(5) )

( -q*x(5) + 3*k(2)*x(2)*x(3)-k(3)*x(4)*x(5) )

( -q*x(6) + 2*k(3)*x(4)*x(5)-k(4)*x(1)*x(6)*sqrt(0.9) )

( -q*x(7) + 2*k(4)*x(1)*x(6)*sqrt(0.9) )

];

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/143047.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • Linux下对MySql数据库备份与恢复

    Linux下对MySql数据库备份与恢复MySQL命令行导入数据库:1,将要导入的.sql文件移至bin文件下,这样的路径比较方便2,同上面导出的第1步3,进入MySQL:mysql-u用户名-p如我输入的命令行:mysql-uroot-p (输入同样后会让你输入MySQL的密码)4,在MySQL-Front中新建你要建的数据库,这时是空数据库,如新建一个名为news的目标数据库5,输入:mysq

  • 地理加权分析_地理加权回归中的拟合度

    地理加权分析_地理加权回归中的拟合度地理加权回归分析完成之后,与OLS不同的是会默认生成一张可视化图,像下面这张一样的:这种图里面数值和颜色,主要是系数的标准误差。主要用来衡量每个系数估计值的可靠性。标准误差与实际系数值相比较小时,这些估计值的可信度会更高。较大标准误差可能表示局部多重共线性存在问题。根据官方的说法,需要检查超过2.5倍标准差的地方……这些地方可能会有问题。虽然在软件里面,默认只显示这样一张图,但

  • EM算法 实例讲解「建议收藏」

    EM算法 实例讲解「建议收藏」第一次接触EM算法,是在完成半隐马尔科夫算法大作业时。我先在网上下载了两份Baum-Welch算法的代码,通过复制粘贴,修修补补,用java实现了HMM算法(应用是韦小宝掷两种骰子的问题)。然后,参考有关半隐马尔科夫算法的论文,照着论文中的公式修改隐马尔科夫算法,完成了大作业。现在回想起来,就隐隐约约记得有一大堆公式。最近,我看到一篇很好的文章,对EM算法的计算有了进一步的了解,文章链接为http

  • C#简单串口调试助手【傻瓜教程】chanson_chapter01[通俗易懂]

    C#简单串口调试助手【傻瓜教程】chanson_chapter01[通俗易懂]简易串口调试助手制作教程【C#.Winform】chanson_chapter01[说明]:本教程基于win10x64位操作系统,软件环境为MicrosoftVisualStudioCommunity2015(以下简称VS)(新建工程文件)【1】打开VS,[文件]–>[新建]–>[项目],弹出窗体如下:(布局方法)【2】开始布局

  • mshta简介

    mshta简介mshta.exe,HTMLApplicationHOST,HTML应用程序主机,hta文件解释器(也能解释html),就如同IE是html的解释器,CMD是批处理文件的解释器,WScript.exe、CScript.exe是vbs/js脚本的解释器。mshta的命令行参数似乎并未公开,mshta/?等方式得不到任何帮助信息。大家如果感兴趣可以去csdn或微软的msdn查查。从网上我仅找到…

  • Java 反射 -超详细讲解(附源码)「建议收藏」

    Java 反射 -超详细讲解(附源码)「建议收藏」  之前也有学习过反射,可是很快就忘了,也不知道有什么用,怎么用,故特此写下此文,以加深对java反射的理解1:反射概述  JAVA反射机制是在运行状态中1,对于任意一个类,都

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号