这微分方程组能解吗?求大佬
mathematica吧
全部回复
仅看楼主
level 11
圣翔天玥 楼主
在书上没见过微分方程组的求解的,但这次遇到一个问题是两个微分方程联立的,按照方程组的解法打进去解不出来,所以求各位大佬帮我解答疑难,二楼贴代码
2017年11月23日 13点11分 1
level 11
圣翔天玥 楼主
\[Rho] = 1000, g = 9.8, h2 = 1.2, h1 = 0.5, v = 4*10^-4, s =
0.1^2*Pi, p0 = 1.01*10^5, dh = 0.00001;
S = DSolve[{\[Rho] g h2 - \[CapitalDelta]p1 -
p0 (v - x2[t] s + x1[t] s)/v == \[Rho] dh x2''[t] &&
p0 (v - x2[t] s + x1[t] s)/
v - \[Rho] g h .1 - \[CapitalDelta]p2 == \[Rho] dh x1''[t],
x1[0] == 0, x2[0] == 0}, x1[t], x2[t], t]
Clear["Global`*"];
2017年11月24日 08点11分 7
level 11
圣翔天玥 楼主
2017年11月24日 08点11分 8
不知道为什么每次发代码都会被删,只有贴图了
2017年11月24日 08点11分
@圣翔天玥 x1[t],x2[t] 那要括起来吧?
2017年11月24日 08点11分
level 11
圣翔天玥 楼主
新的是这个样子的
2017年11月24日 09点11分 10
NDSolve 自变量需要指定数值范围,这一点和 DSolve 不一样
2017年11月24日 10点11分
@Alexander0620 这样啊,谢谢
2017年11月24日 10点11分
吧务
level 15
……要分清楚分号和逗号。
{\[Rho] = 1000, g = 9.8, h2 = 1.2, h1 = 0.5, v = 4*10^-4, s = 0.1^2*Pi, p0 = 1.01*10^5,
dh = 0.00001};
S = DSolve[{\[Rho] g h2 - \[CapitalDelta]p1 -
p0 (v - x2[t] s + x1[t] s)/v == \[Rho] dh x2''[t] &&
p0 (v - x2[t] s + x1[t] s)/v - \[Rho] g h .1 - \[CapitalDelta]p2 == \[Rho] dh x1''[
t], x1[0] == 0, x2[0] == 0}, {x1[t], x2[t]}, t]
2017年12月02日 08点12分 11
1