jzh686 2021-10-13 20:32 采纳率: 66.7%
浏览 110
已结题

MATLAB的ode45函数解常微分方程的用法!!

img


img

这个怎么用MATLAB的ode函数来弄呀,我看其他教程都是只有一个方程,;两个方程怎么弄呀,请教一下!

  • 写回答

1条回答 默认 最新

  • joel_1993 2021-10-13 20:50
    关注

    你好,两个方程一样求解:
    代码供参考,该注释的都注释了

    
    % 设置物块质量和弹簧刚度
    m1 = 1;
    m2 = 1;
    k1 = 1;
    k2 = 1;
    k3 = 1;
    % 设x(1)对应物块1的位移;x(2)对应物块1的速度
    % 设x(3)对应物块2的位移;x(4)对应物块2的速度
    odefun = @(t, x) [ x(2); % 物块1的速度对时间积分得到位移
        (k2*x(3) - (k1+k2)*x(1)) / m1; % 物块1的加速度对时间积分得到速度
        x(4); % 物块2的速度对时间积分得到位移
        ( -(k2+k3)*x(3) + k2*x(1)) / m2]; % 物块2的加速度对时间积分得到速度
    x0 = [0.1; 0; 0; 1]; % 物块1初始位移0.1,初始速度0,物块2初始位移0,初始速度1
    tspan = [0, 10]; % 运行10s
    [t,x] = ode45(odefun, tspan, x0);%ode45求解
    % 画图
    plot(t,x(:,1),'r-');hold on
    plot(t,x(:,2),'b-')
    plot(t,x(:,3),'r--')
    hold on
    plot(t,x(:,4),'b--')
    legend('物块1位移','物块1速度','物块2位移','物块2速度')
    

    img

    有帮助望题主给个采纳支持一下

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 系统已结题 10月21日
  • 已采纳回答 10月13日
  • 修改了问题 10月13日
  • 创建了问题 10月13日

悬赏问题

  • ¥15 nrf52810-c三个a 程序
  • ¥15 lego-loam跑出来的roll误差很大
  • ¥50 求一个半透明没有锯齿的圆角窗体的实现例子
  • ¥15 STM32cubeMX里的FreeRTOS无法释放内存
  • ¥15 CATIA有些零件打开直接单机确定终止
  • ¥15 请问有会的吗,用MATLAB做
  • ¥15 phython如何实现以下功能?查找同一用户名的消费金额合并—
  • ¥15 ARIMA模型时间序列预测用pathon解决
  • ¥15 孟德尔随机化怎样画共定位分析图
  • ¥18 模拟电路问题解答有偿速度