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

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 用vmmare虚拟机用sentaurus仿真的时候,调用terminal程序,输入swb指令弹出这个,打不开workbench,桌面上面的sentaurus workbench也打不开
  • ¥75 使用winspool.drv的SetPrinter设置打印机失败
  • ¥15 simulink 硬件在环路仿真
  • ¥15 python动态规划:N根火柴摆出的最大数字
  • ¥20 (标签-excel)
  • ¥200 求idea10和MyEclipse7.1
  • ¥20 vb6.0截取当前窗体保存为jpg文件
  • ¥20 苹果手机不使用大疆sdk怎么获取遥控器控制信息或如何接入大疆sdk并且成功上架sdk
  • ¥20 woocommerce 注册按键重定向
  • ¥100 求书法图像文字切割代码