我的scilab并没有展示任何关于我试图制造“振荡”的东西é犬齿联轴节é从python到scilab的问题,但它应该显示3个表示振荡的图形
下面是python:
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
m1=1;m2=1;k1=1;k2=1;k3=1;m=1;k=1;a=1;
def equadiff(var,t):
x1=var[0]
x2=var[1]
v1=var[2]
v2=var[3]
a1=( k2*x2-(k1+k2)*x1)/m1
a2=( k2*x1-(k2+k3)*x2)/m2
return ([v1,v2,a1,a2])
def Oscillateur(x10,x20,v10,v20,tmax):
t=np.linspace(0,tmax,300)
fig = plt.figure()
ax = fig.add_subplot(1,1,1,aspect="equal")
ax.spines["left"].set_position("zero")
ax.spines["right"].set_color("none")
ax.spines["bottom"].set_position("zero")
ax.spines["top"].set_color("none")
ax.xaxis.set_ticks_position("bottom")
ax.yaxis.set_ticks_position("left")
plt.axis("tight")
var0=[x10,x20,v10,v20]
solution=odeint(equadiff,var0,t)
x1=solution[:,0]
x2=solution[:,1]
v1=solution[:,2]
v2=solution[:,3]
plt.plot(t,x1,t,x2)
plt.axis("tight")
plt.show()
Oscillateur(a,a,0,0,10)
Oscillateur(a,-a,0,0,10)
Oscillateur(2*a,0,0,0,50)
这是我的scilab,问题在这段代码中,我猜应该是那些ax.yaxis.set\u ticks\u position(“左”),但我不知道如何解决它并显示我的图形:
m1=1;m2=1;k1=1;k2=1;k3=1;m=1;k=1;a=1;
function S=equadiff(var,t)
x1=var(0)
x2=var(1)
v1=var(2)
v2=var(3)
a1=(k2*x2-(k1+k2)*x1)/m1
a2=(k2*x1-(k2+k3)*x2)/m2
disp([v1,v2,a1,a2])
endfunction
function T=Oscillateur(x10,x20,v10,v20,tmax)
t=linspace(0,tmax,300)
figure()
ax=subplot(1,1,1,aspect="equal")
var0=[x10,x20,v10,v20]
solution=ode(var0,t0,t,equadiff)
x1=solution(:,0)
x2=solution(:,1)
v1=solution(:,2)
v2=solution(:,3)
plot(t,x1,t,x2)
a=gca()
a.x_location = "zero";
axis("tight");
endfunction
disp(Oscillateur(a,a,0,0,10))
disp(Oscillateur(a,-a,0,0,10))
disp(Oscillateur(2*a,0,0,0,50))
这里有很多错误:
equadiff
返回未声明的S
李>subplot
没有名为“aspect”的参数李>t0
未定义李>equadiff
参数的顺序必须是t,var
李>x
如果t是行向量,则参数必须是列向量李>Oscillateur
返回未声明的T
李>disp
无法返回此空返回李>x_location
属性不存在修复所有这些错误将为您提供:
相关问题 更多 >
编程相关推荐