求解三阶机电暂态模型微分方程matlab代码

您好，您可以使用下面的Matlab代码来求解三阶机电暂态模型微分方程：

function [x,t]=mechanical_oscillationm=1;
k=1;
b=0.5;
f=1;
w=sqrt(k/m);

% 设置初始位置和初速度
x0=0.1;
v0=0;

% 定义微分方程
ode=@(t,x) [x(2);-(b/m)*x(2)-(k/m)*x(1)+f*sin(w*t)];

% 设置时间范围
tspan=[0,10];

% 设置初始条件
y0=[x0;v0];

% 使用ode45求解微分方程
[t,x]=ode45(ode,tspan,y0);

% 绘