✅作者简介热爱科研的Matlab仿真开发者擅长毕业设计辅导、数学建模、数据处理、程序设计科研仿真。完整代码获取 定制创新 论文复现点击Matlab科研工作室 关注我领取海量matlab电子书和数学建模资料个人信条做科研博学之、审问之、慎思之、明辨之、笃行之是为博学慎思明辨笃行。 内容介绍一、引言质量 - 弹簧 - 阻尼器系统是物理学和工程学中常见的动力学模型广泛应用于机械振动、结构动力学等领域。熵作为热力学和信息论中的重要概念反映了系统的无序程度或不确定性。对质量 - 弹簧 - 阻尼器系统进行熵生成分析有助于深入理解系统的能量耗散和不可逆过程。Conformable 模型是一种新兴的非整数阶导数模型相较于传统整数阶导数它在描述复杂动力学系统时具有独特优势能够更精确地刻画系统的记忆和遗传特性。本文将基于 Conformable 模型对质量 - 弹簧 - 阻尼器系统进行熵生成分析揭示系统在非整数阶动力学下的熵演化规律。二、质量 - 弹簧 - 阻尼器系统模型⛳️ 运行结果 部分代码function dy mechsys_01(t,y,coefficients,definition)%% Input parameters:% ------------------% coefficients (4x1) - a_1*d2x a_2*d1x a_3*d0x 0, a_4 - gamma% definition (char) - {traditional|conformable}%% Output parameter:% ------------------% ODE system - d1y d2x, y(0),% d1x y, x(0),m coefficients(1);b coefficients(2);k coefficients(3);if numel(coefficients) 4g 1.0; % If g is not included in coefficients (default)elseg coefficients(4); % Read g from the last coefficientend% Note: if definition traditional, g is ignored% Find auxiliar paramters like natural frequency and damping factoromega_n sqrt(k/m);zeta b/(2*sqrt(k*m));% Select the model definition to use:switch definitioncase conformable % From Don Juans manuscriptdy(1,1) y(2);dy(2,1) -2*zeta*((omega_n).^g).*(t.^(g-1)).*y(2) - ...((omega_n).^(2*g)).*(t.^(2*(g-1))).*y(1);otherwise % Equivalent to definition traditional (default)dy(1,1) y(2);dy(2,1) -2*zeta*omega_n*y(2) - omega_n^2*y(1);endend 参考文献更多免费数学建模和仿真教程关注领取
