使用ode45为弹簧摆时出错

“求解弹簧摆动系统的EOM，然后将EOM集成10秒，从静止时与径向位置静态。春天被拆开。然后，将EOM集成在休息时以各种初始开始于120秒以米为单位和程度的条件。“

目标：使用带有'x'标记的蓝线和使用黑色线与“x”标记的质量的最终位置绘制质量的初始位置。

步骤1：

                          清除变量
                         
                          CLC.
                         
                          m = 0.1416;％ 公斤
                         
                          k = 7.1;％n / m
                         
                          g = 9.81;％m / s ^ 2
                         
                          r_0 = 9;％ 厘米
                         
                          r = 29;％ 厘米
                         
                          θ= 30;％程度

第2步：

列表变量：

M：弥撒，

k：弹簧常数，

R：RADIUS，

R_0：未拉伸半径，

θ：初始角位移，

G：重力，

R_M：R方向上的质量线性位移，状态变量

R_MDDOT：R方向质量的加速度（EOMS）

theta_m： - 在-direction，状态变量中的质量位移

Theta_mddot： - 重定程中的群众的加速度（EOMS）

第3步：

                          Syms.m l r.
                         
                          Syms.Thetadot Thetaddot RDOT RDOT
                         
                          EQN（1）= m * thetaddot == m *（（l + r）* thetaddot + 2 * rdot * thetadot）;
                         
                          EQN（2）= M * RDDOT == M *（rddot-（l + r）* thetadot ^ 2）;

第四步：

找到未知的变量

Theta双点，以及R双点：分别在θ和r方向上的加速度

2等式，2个未知数

第5步：

                          未知数=解决（EQN，RDDOT，THETADDOT）
                         
                          Acceleoms =解决（EQN，RDDOT，THETADDOT）;
                         
                          Acceleoms = struct2cell（加速）;

第6步：

解决运动方程（EOM）

                          [t，s] = ode45（@（t，s）EOM（m，k，r，r_0，theta，g），linspace（0,10,100），[0,0]）;
                         

你现在关注这个问题

3评论
表演隐藏＆nbsp2旧的评论

答案（0）

也可以看看

类别

标签

下载188bet金宝搏

释放

社区宝藏狩猎

标签

下载188bet金宝搏

释放

也可以看看

社区宝藏狩猎

你现在关注这个问题

使用ode45为弹簧摆时出错

3评论 表演隐藏＆nbsp2旧的评论

答案（0）

也可以看看

类别

标签

下载188bet金宝搏

释放

社区宝藏狩猎

标签

下载188bet金宝搏

释放

也可以看看

社区宝藏狩猎

发生错误

3评论
表演隐藏＆nbsp2旧的评论