这个怎么用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速度')
有帮助望题主给个采纳支持一下
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 1无用
悬赏问题
- ¥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 求书法图像文字切割代码