如何使用数值求解二阶微分方程随机变量吗?
4视图(30天)
显示旧的评论
在季度汽车模型方程,(见附加文件以上skematic模型)
在哪里
=簧载质量
=非簧载质量
=非线性弹簧stifffness
=线性弹簧刚度
=由于阻尼器的阻尼系数
——Sinusodial函数振幅
和频率
是独立与正态分布的随机变量。
我怎么能把这个随机变量的方程采用数值吗?
我做了一个样本,任何模拟100个样本。
dt = 0.01;
tf = 30;
tspan = dt: dt: tf;
y0 = [0; 100; 0; 10];
[t、y] =数值(“eqsystem”tspan, y0);
情节(y (: 1))
网格
包含(“X-Displacement”)
ylabel (“y-Displacement”)
标题(“X和Y位移”)
持有在;
函数dydt = eqsystem (t, y)
ks = 1000;% N / m ^ 3 -高斯
ku = 1000;% N / m -高斯
女士= 10;% kg -高斯
μ= 20;% kg -高斯
c = 300;% Ns / m -高斯
一个= 0.10;% m -
ω= 2 *π;% rad / s -
dydt = 0 (4,1);
:dydt (1) = y (2);
:dydt (2) = ((- k / ms) * (y (1) - y (3)) ^ 3 - (c / ms) * (y (2) - y (4)));
:dydt (3) = y (4);
dydt (4) = ((ks /μ)* (y (1) - y (3)) ^ 3 + (y (2) - y (4)) + ku * (A * sin(ω* t) - y (3)));
结束
接受的答案
1月
2021年6月6日
dt = 0.01;
tf = 30;
tspan = dt: dt: tf;
y0 = [0; 100; 0; 10];
为k = 1:10 0
女士= 10 + randn;
μ= 20 + randn * 2;
…
fcn = @ (t, y) eqsystem (t, y,女士,亩,…)
[t、y] =数值(“eqsystem”tspan, y0);
结束
函数dydt = eqsystem (t y女士,亩,…)
dydt = 0 (4,1);
dydt (1) = y (2);
dydt (2) = (- k / ms) * (y y (1) - (3)) ^ 3 - (c / ms) * (y y (2) - (4));
dydt (3) = y (4);
dydt (4) = (ks /μ)* (y y (1) - (3)) ^ 3 + (y y (2) - (4) +…
ku *(*罪(ω* t) - y (3));
结束