弹簧质量采用数值

9的观点(30天)
リッディチャワン
リッディチャワン 2021年11月18日
编辑: VBBV 约10小时前
我写这个程序求解弹簧质量方程采用数值,但它显示了一个错误的数值
M = 1.0;
K = 10;
C = 1;
l = 0.5;
g = 9.8;
x0 = [0;
0);
tspan = [0, 10];
(t, x) =数值(@func tspan, x0);
次要情节(2,1,1)
情节(t) x (: 1));
网格(“上”);
%轴([0 10 1 1]);
包含(“时间[s]”);
ylabel (“x [m]”);
次要情节(2,1,2)
情节(t) x (:, 2));
网格(“上”);
%轴([0 10 1 1]);
包含(“时间[s]”);
ylabel (“美元\点{x} $ (m / s) ');
%導関数の定義
函数xdot = func (t, x)
M = 1.0;
K = 10;
C = 1;
l = 0.5;
g = 9.8;
xdot = [x (2);
- k * (xl) * (1) / mc * x (2) / M + g] * dt;
结束

答案(1)

艾伦·史蒂文斯
艾伦·史蒂文斯 2021年11月18日
像这样
x0 = [0;
0);
tspan = [0, 10];
(t, x) =数值(@func tspan, x0);
次要情节(2,1,1)
情节(t) x (: 1));
网格(“上”);
%轴([0 10 1 1]);
包含(“时间[s]”);
ylabel (“x [m]”);
次要情节(2,1,2)
情节(t) x (:, 2));
网格;
%轴([0 10 1 1]);
包含(“时间[s]”);
ylabel (“dx / dt (m / s)”);
函数xdot = func (x ~)
M = 1.0;
K = 10;
C = 1;
l = 0.5;
g = 9.8;
xdot = [x (2);
- k * (x (1) - l) * (1) / mc * x (2) / M + g);%(1)不仅仅是x。不需要dt
结束

类别

找到更多的在编程帮助中心文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!