matlab ode45求解齿轮动力学,Matlab拟合动力学参数遇到问题(ode45)
function [rho2,F] = rho2_F(k,r,ss,Ne,Np)
sy = sum(r.^2);
rho2 = 1 - ss./sy
F = (sy - ss).*(Ne-Np)./(Np.*ss) % F:F比
fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n')
fprintf('\tk1 = %.4f ± %.4f\n',k(1),ci(1,2)-k(1))
fprintf('\tE1 = %.4f ± %.4f\n',k(2),ci(2,2)-k(2))
fprintf('\tk2 = %.4f ± %.4f\n',k(3),ci(3,2)-k(3))
fprintf('\tE2 = %.4f ± %.4f\n',k(4),ci(4,2)-k(4))
fprintf('\tk3 = %.4f ± %.4f\n',k(5),ci(5,2)-k(5))
fprintf('\tE3 = %.4f ± %.4f\n',k(6),ci(6,2)-k(6))
fprintf('\tk4 = %.4f ± %.4f\n',k(7),ci(7,2)-k(7))
fprintf('\tE4 = %.4f ± %.4f\n',k(8),ci(8,2)-k(8))
fprintf('\tk5 = %.4f ± %.4f\n',k(9),ci(9,2)-k(9))
fprintf('\tE5 = %.4f ± %.4f\n',k(10),ci(10,2)-k(10))
fprintf('\tk6 = %.4f ± %.4f\n',k(11),ci(11,2)-k(11))
fprintf('\tE6 = %.4f ± %.4f\n',k(12),ci(12,2)-k(12))
fprintf('\tk7 = %.4f ± %.4f\n',k(13),ci(13,2)-k(13))
fprintf('\tE7 = %.4f ± %.4f\n',k(14),ci(14,2)-k(14))
fprintf('\tk8 = %.4f ± %.4f\n',k(15),ci(15,2)-k(15))
fprintf('\tE8 = %.4f ± %.4f\n',k(16),ci(16,2)-k(16))
fprintf('\tk9 = %.4f ± %.4f\n',k(17),ci(17,2)-k(17))
fprintf('\tE9 = %.4f ± %.4f\n',k(18),ci(18,2)-k(18))
function f = fy(k,x0,y0,ts,tt)
for i=1:42,%42组试验点
[t Xsim] = ode45(@g,[0,ts(i)],x0(i,:),[],k,tt(i));
i
ysim(i,1) = Xsim(end,1);
ysim(i,2) = Xsim(end,2);
ysim(i,3) = Xsim(end,3);
ysim(i,4) = Xsim(end,4);
ysim(i,5) = Xsim(end,5);
ysim(i,6) = Xsim(end,6);
ysim(i,7) = Xsim(end,7);
ysim(i,8) = Xsim(end,8);
ysim(i,9) = Xsim(end,9);
ysim(i,10) = Xsim(end,10);
ysim(i,11) = Xsim(end,11);
ysim(i,12) = Xsim(end,12);
%12个组分
end;
f = [ysim(:,:)-y0(:,:)];
function dCdt =g(t,C,k,tt) % ODE模型
r1=k(1).*exp(-k(2)./8.314.*(1/tt-1/653)).*C(1).*C(2)-k(1).*exp(-k(2)./8.314.*(1/tt-1/653))./(0.001*tt*tt-1.514*tt+514.6).*C(4);
r2=k(3).*exp(-k(4)./8.314.*(1/tt-1/653)).*C(2).*C(2)-k(3).*exp(-k(4)./8.314.*(1/tt-1/653))./(0.00002*tt*tt-0.033*tt+11.44).*C(5);
r4=k(5).*exp(-k(6)./8.314.*(1/tt-1/653)).*C(2).*C(3)-k(5).*exp(-k(6)./8.314.*(1/tt-1/653))./(0.00003*tt*tt-0.046*tt+15.94).*C(6);
r5=k(7).*exp(-k(8)./8.314.*(1/tt-1/653)).*C(2).*C(4)-k(7).*exp(-k(8)./8.314.*(1/tt-1/653))./(0.00003*tt*tt-0.041*tt+14.14).*C(7);
r6=k(9).*exp(-k(10)./8.314.*(1/tt-1/653)).*C(3).*C(4)-k(9).*exp(-k(10)./8.314.*(1/tt-1/653))./(0.00003*tt*tt-0.034*tt+11.86).*C(8);
r8=k(11).*exp(-k(12)./8.314.*(1/tt-1/653)).*C(6);
r9=k(13).*exp(-k(14)./8.314.*(1/tt-1/653)).*C(7);
r10=k(15).*exp(-k(16)./8.314.*(1/tt-1/653)).*C(8) ;
r11=k(17).*exp(-k(18)./8.314.*(1/tt-1/653)).*C(2) ;
dC1dt=-r1;
dC2dt=-r1-2*r2-r4-r5-r11;
dC3dt=-r4-r6;
dC4dt=r1-r5-r6;
dC5dt=r2;
dC6dt=r4-r8;
dC7dt=r5-r9;
dC8dt=r6-r10;
dC9dt=r8;
dC10dt=r9;
dC11dt=r10;
dC12dt=r11;
dCdt = [dC1dt;dC2dt;dC3dt;dC4dt;dC5dt;dC6dt;dC7dt;dC8dt;dC9dt;dC10dt;dC11dt;dC12dt];
matlab ode45求解齿轮动力学,Matlab拟合动力学参数遇到问题(ode45)相关推荐
- matlab ode45求解齿轮动力学,ode45求解多自由度动力学方程实例.doc
ode45求解多自由度动力学方程实例 Ode45函数调用形式如下:[T,Y]=ode45(odefun,tspan,y0)相关参数介绍如下:参数名称参数说明odefun用于存放待求解的方程的m文件名, ...
- matlab ode45求解齿轮动力学,[转载]Matlab中解常微分方程的ode45 【转载】
ode是专门用于解微分方程的功能函数,他有ode23,ode45,ode23s等等,采用的是Runge-Kutta算法.ode45表示采用四阶,五阶runge-kutta单步算法,截断误差为(Δx)^ ...
- matlab ode45求解齿轮动力学,使用ODE45求解齿轮系统动力学方程后结果发散
EDA365欢迎您登录! 您需要 登录 才可以下载或查看,没有帐号?注册 x 1.1参数定义及动力学方程降阶 : g) c6 b" p' | p7 _8 Gfunction [dx,ff1 ...
- Matlab遗传算法求解TSO,遗传算法matlab程序实例.doc
遗传算法matlab程序实例.doc --------------------------------------------------------------------------------- ...
- matlab 变分不等式,求解变分不等式的matlab程序
满意答案 kide754861 2013.07.11 采纳率:46% 等级:11 已帮助:5599人 function x=Porjection() clc sigama=0.5; gama=0 ...
- matlab如何求解代数方程,用Matlab求解复杂代数方程
编写一个function,保存为KMN.m function [K1,M1,N1,K2,M2,N2,K3,M3,N3]=KMN(x,y,z,R,r,L,l) K1=(power(x,2)+power( ...
- matlab中常微分方法,MATLAB解常微分方程组的解法(好东西要共享)
1:问题 常微分方程的初值问题的标准数学表述为:y'=f(t,y),a<=t<=b,y(a)=y(0) :我们要求解的任何高阶常微分方程都可以用替换法化为上式所示的一阶形式,其中y为向量, ...
- matlab求微分方程的初值问题,Matlab 解常微分方程的初值问题.doc
Matlab 解常微分方程的初值问题 题目:Matlab 解常微分方程的初值问题 设计目的: 1.熟练掌握Matlab的基本编程方法,及其编程风格. 2.熟练掌握Matlab常用函数的使用. 3.与本 ...
- 齿轮-转子-轴承系统动力学matlab程序代码
齿轮-转子-轴承系统动力学matlab程序代码 详询TB机械动力学超市 感谢关注 期待你的到来
最新文章
- 2022-2028年中国抗菌不锈钢行业研究及前瞻分析报告
- UniApp H5 浏览器下载图片 兼容iOS、Android
- 「shell」替代rm,放入回收站
- 直播预告 | 第四范式2021发布会技术分论坛报名开启,6月23日线上见
- 【数据结构与算法】之深入解析“下一个更大元素I”的求解思路与算法示例
- HTML与CSS布局技巧总结
- android管理activity顺序,android activity 六个主要函数以及一些调用顺序
- 基于Cairngorm的Silverlight开发 - part3
- php枚举mysql,小技巧:取得MYSQL中ENUM(枚举)列的全部可能值。-PHP教程,PHP应用...
- 每隔一段时间执行php_用php脚本,你如何定时更新商品列表
- homework7_ZhankunLuo
- c语言按键长按双击,51单片机实现单按键单击、双击、长按功能
- 【软件需求工程】北理的恶龙们——软件需求规格说明
- win10 32位和64位的区别,win10 32位和64位有什么区别
- 百度搜索移动端流量词热度统计方法
- 计算机控制电机启动接线图,电机控制线路图大全(上下)
- php小写转大写,php怎么实现英文小写转大写
- 一个IT公司技术总监的职责
- VirusTotal 为 Chrome 和 Firefox 发布 VT4Browsers 扩展
- python证书有用吗_python有证书考吗
热门文章
- 测试两个主机之间的连通性_常用的光缆4种测试方法介绍
- 明显调用的表达式前的括号必须具有函数类型_Chisel(二) Scala语法 变量与函数...
- 【Servlet】getInitParameter()发生空指针导致HTTP500的解决方案
- 【算法分析与设计】埃氏筛素数算法
- 怎样成为一个好的技术领导者
- 从初创型到独角兽企业,监控架构演进的那些事儿
- React Suspense提供Redux的替代方案
- 在Exchange Server 2007中为出站邮件添加免责声明
- 分布式系统测试那些事儿——信心的毁灭与重建
- strlensizeof