功能结果= mdof_simulation（m，c，k，f，fs）
输入：
M：质量矩阵（n * n）
C：阻尼矩阵（n * n）
K：刚度矩阵（n * n）
F：外力矩阵（n，n）
FS：采样频率
其中n是自由度的程度，n是动态力的数据点的长度
输出：
结果：是一个结构组成
结果。分离：位移（n * n）
结果。velocity：速度（n * n）
结果.Acceleration：加速度（n * n）
结果.Parameters.freq =自然频率（n * 1）
结果.Parameters.Dampratio =阻尼比率（n * 1）
结果.parameters.modeShape =模式形状矩阵（n * n）
参考：
Chopra，Anil K.“结构的动态。理论与应用。”地震工程（2017年）。
引用
Ayad Al-Rumaithi（2021）。多程度自由系统求解器（https://www.mathwands.com/matlabcentral/fileexchange/69509-multi-degree-of-freedom-system-solver），Matlab中央文件交换。检索到。
评论和评级（11.）
您还可以从以下列表中选择一个网站：
F的第一栏是在时间步骤1的力。第二列在时间步骤2 ......等等
步骤之间的时间差等于dt = 1 / fs
你好Ayad.
我的问题与下面的问题相同
您是否可以详细说明F：外力矩阵（N，N）的输入要求？如果我有3到3个矩阵，我想在3块质量上施加力量。第一列将是0.0,1000。F矩阵中的第二个和第三等等是什么
提前致谢
使用绘图功能
//www.tatmou.com/help/matlab/ref/plot.html.
喜欢：
绘图（t，结果。velocity（1，:)）;
其中t是时间向量
非常感谢。代码不会在图中绘制位移，速度和加速度，我无法将其运行。有人可以帮助这个数字的代码吗？
您是否可以详细说明F：外力矩阵（N，N）的输入要求？