Question

mathematica将数值求解的参数微分方程组画图

通过l = 10; \[Omega] = 5; R = 10; g = 9.8;
NDSolve[{x''[
t] == -((g*(z[t])*(x[t]) + ((x'[t])^2 + (y'[t])^2 + (z'[t])^2) x[
t])/l^2) + 2 \[Omega]*y[t] + \[Omega]^2*x[t],
y''[t] == -((g*(z[t])*(y[
t]) + ((x'[t])^2 + (y'[t])^2 + (z'[t])^2) y[t])/l^2) +
2 \[Omega]*x[t] + \[Omega]^2*y[t],
z''[t] == -((g*(z[t])^2 + ((x'[t])^2 + (y'[t])^2 + (z'[t])^2) z[t])/
l^2) - g, x[0] == R, y[0] == 0, z[0] == -l, x'[0] == 0,
y'[0] == -\[Omega]*R, z'[0] == -l}, {x, y, z}, {t, 0, 180}]
数值求解了参数方程组，现在希望把结果表示在o－xyz坐标系中，也就是消去参数画出函数图形，应该怎样写命令？

Answer 1

追问

如果我想把得到的（x,y,z）坐标通过乘以一个3*3矩阵进行坐标变化转化到另一个坐标系中再画图应当怎样做？The tensor is
{{Cos[\[Omega]t], -Sin[\[Omega]t], 0}, {Sin[\[Omega]t],
cos[\[Omega]t], 0}, [0, 0, 1]}

追答

你自己乘以变换矩阵得到新的参数方程不就可以了.