如果不想使用面片,可以使用椭圆的参数方程: int x = a(t);y = v + B sin(t)
import numpy as np
from matplotlib import pyplot as plt
from math import pi
u=1. #x-position of the center
v=0.5 #y-position of the center
a=2. #radius on the x-axis
b=1.5 #radius on the y-axis
t = np.linspace(0, 2*pi, 100)
plt.plot( u+a*np.cos(t) , v+b*np.sin(t) )
plt.grid(color='lightgray',linestyle='--')
plt.show()
其给出:
由于2D旋转矩阵,椭圆可以旋转:
import numpy as np
from matplotlib import pyplot as plt
from math import pi, cos, sin
u=1. #x-position of the center
v=0.5 #y-position of the center
a=2. #radius on the x-axis
b=1.5 #radius on the y-axis
t_rot=pi/4 #rotation angle
t = np.linspace(0, 2*pi, 100)
Ell = np.array([a*np.cos(t) , b*np.sin(t)])
#u,v removed to keep the same center location
R_rot = np.array([[cos(t_rot) , -sin(t_rot)],[sin(t_rot) , cos(t_rot)]])
#2-D rotation matrix
Ell_rot = np.zeros((2,Ell.shape[1]))
for i in range(Ell.shape[1]):
Ell_rot[:,i] = np.dot(R_rot,Ell[:,i])
plt.plot( u+Ell[0,:] , v+Ell[1,:] ) #initial ellipse
plt.plot( u+Ell_rot[0,:] , v+Ell_rot[1,:],'darkorange' ) #rotated ellipse
plt.grid(color='lightgray',linestyle='--')
plt.show()
3条答案
按热度按时间olhwl3o21#
如果不想使用面片,可以使用椭圆的参数方程:
int x = a(t);y = v + B sin(t)
其给出:
由于2D旋转矩阵,椭圆可以旋转:
返回:
h6my8fg22#
matplotlib椭圆演示很不错。但是我不能在没有for循环的代码中实现它。我得到了一个轴的数字错误。这是我所做的,当然xy中心是我自己的坐标,具有相应的宽度和高度,基于我绘制椭圆的图像。
此代码部分基于this page上的第一个代码框。请参阅上面Chris的回复,以获得
matplotlib.patches.Ellipse
的链接。35g0bw713#
你看过matplotlib ellipse demo吗?这里使用
matplotlib.patches.Ellipse
。