Scilab "for" 循环问题

Question

我正在尝试使用 scilab 显示多个图表，但它根本不起作用我的代码基于 python Python: 将 matplotlib.pyplot 导入为 plt 将 numpy 导入为 np

g = 9.81; v0 = 30; ax = 0; x0 = 0; az = -g; z0 = 0;
for i in np.linspace(0.5,9,20):
  theta = i*10*np.pi/180
  vx0 = v0*np.cos(theta)
  vz0 = v0*np.sin(theta)
  t=np.linspace(0,2*vz0/g,100)
  x =vx0*t+x0
  z =0.5*az*t**2+vz0*t+z0
  plt.plot(x,z)

这是 scilab 代码：

g = 9.81; v0 = 30; ax = 0; x0 = 0; az = -g; z0 = 0;
for i=0.5:9:20 t=linspace(0,2*v0*sin(i*10*%pi/180)/g,100) 
    x=v0*cos(i*10*%pi/180)*t+x0 
    z=0.5*az*t**2+v0*sin(i*10*%pi/180)*t+z0 
    plot(x,z) 
    ;end

我的目标是获得： https://imgur.com/dlDIKNx

Answer 1

Scilab 语法非常接近你python语法。

对您的一些修正 python 代码将为您提供合适的 scilab 代码：

g = 9.81; v0 = 30; ax = 0; x0 = 0; az = -g; z0 = 0;
colorindex=1;
for i=linspace(0.5,9,20)
      theta = i*10*%pi/180
  vx0 = v0*cos(theta)
  vz0 = v0*sin(theta)
  t=linspace(0,2*vz0/g,100)
  x =vx0*t+x0
  z =0.5*az*t**2+vz0*t+z0
  plot2d(x,z,colorindex)
  colorindex=colorindex+1;
end

'enveloppe'（或者你怎么称呼它）仍然没有什么问题，但我认为你现在有修复它的线索。

Answer 2

您可以通过在@Marvin 的代码末尾添加来获得'enveloppe'：

vx0 = v0
vz0 = 0
t=linspace(0,sqrt(-(-v0/az)**2+2*v0**2/az^2),100)
x =vx0*t+x0
z =0.5*az*t**2 + 0.5*az*(-v0/az)**2-v0*(v0/az)
plot2d(x,z,colorindex)

Scilab "for" 循环问题

Scilab "for" loop issue

python

graph

scilab