最近看到个视频,打算绘制个心型动态的曲线。
网上有很多关于心形曲线的资料,各种各样的形状,但是我比较倾心于桃心形。先在Matlab里简单看一下是什么样的。
桃心形曲线的方程如下
Matlab绘制效果如下:
Matlab程序如下。
clc
clear
t = -10:0.01:10
x = 16*(sin(t)).^3
y = 13*cos(t)-5*cos(2*t)-2*cos(3*t)-cos(4*t)
%ylim([-5,2])
%xlim([-3,3])
grid on;
hold on
plot(x,y,'linewidth',1)
光是一张图片,属实是太乏味了,不如直接找网图。为了让这个心形曲线更有意思点,我使用了matplotlib
中的animation
模块来绘制一个动态的心形。
首先需要进行初始化x,y
列表、窗口和曲线。
fig, ax = plt.subplots()
#初始化窗口和坐标轴
xdata, ydata = [], []
#初始化x,y列表
ln, = ax.plot([], [], 'r-', animated=False)
#初始化绘制曲线的参数。
然后使用两个核心函数init()
和updata()
计算曲线经过的点。
#init()函数初始化x,y轴范围
def init():
ax.set_xlim(-20,20)
ax.set_ylim(-20,15)
return ln,
#迭代器,根据传入的frame更新x,y值
def update(frame):
xdata.append(np.sin(frame)^3)
ydata.append(13*np.cos(frame)-5*cos(2*frame)-2*cos(3*frame)-cos(4*frame))
ln.set_data(xdata, ydata)
return ln,
最后调用一下函数,设置一下显示速度(interval),是否重复绘制(repeat)等参数即可。
ani = animation.FuncAnimation(fig, update, frames=np.linspace(0,10,1000),init_func=init,interval=5, repeat = False,blit=True)
ani.save('heart.gif',writer='imagemagick',fps=30)
%matplotlib qt5