python可视化欧拉公式

欧拉公式图解

欧拉公式

$$e^{i\theta }=\cos \theta +i\sin \theta$$
对于每个实数$\theta$，都可以在复平面上的单位元找到对应的一点，如图：

图片来源见水印。

如果$\theta$是随着时间变化的，即$\theta (t)=wt$，其中$w$是角频率（角速度），那么公式变为：

$$e^{i\theta (t)}=\cos (\theta (t))+i\sin (\theta (t))=e^{iwt}=\cos (wt)+i\sin (wt)$$

相应的图形则是在复平面的实轴，虚轴以外再增加一个时间轴$t$。函数随时间变化的动态图如下：

可以看到，$e^{i\theta (t)}$沿着实轴投影，得到$i\sin \theta (t)$的变化图；沿着虚轴投影，得到$\cos \theta (t)$的变化图，沿着时间轴投影，则得到复平面上的单位圆：

增大$w$，对比如下：

可以看到$w$越大，单位时间内角度（相角）变化越大，转一圈所需的时间变少，即函数的周期（$T=\frac{2\pi }{w}$）变小。

python代码：

import matplotlib.pyplot as plt
import numpy as np

plt.rcParams['font.sans-serif'] = ['SimHei']  # 用于正常显示中文标签
plt.rcParams['axes.unicode_minus'] = False  # 用来正常显示负号

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
font = {'size': 18}

lim_value = 1.3
w_list = [2, 4]  # 角速度
t = np.arange(0, 10, 0.02)
width = 1.5
for i, w in enumerate(w_list):
    re = np.cos(w * t)
    im = np.sin(w * t)
    width = 1.5
    if i == 0:
        width = 4  # 为了对比明显
        # theta不随时间变化的情况
        ax.plot(re, im, zs=0, zdir='z', label=r'$\theta$不随时间变化', linewidth=width)
        ax.plot(np.zeros(np.size(im)), im, zs=0, zdir='z', color='black')
        ax.plot(re, np.zeros(np.size(re)), zs=0, zdir='z', color='black')
    # e^{jwt}
    ax.plot(re, im, t, label=r'$e^{i\theta},\theta=%st$' % w, linewidth=width)
    # cos(wt)
    ax.plot(re, t, zs=lim_value, zdir='y', label=r'$cos(\theta),\theta=%st$' % w, linewidth=width)
    ax.plot(np.zeros(np.size(re)), t, zs=lim_value, zdir='y', color='black')
    # sin(wt)
    ax.plot(im, t, zs=-lim_value, zdir='x', label=r'$sin(\theta),\theta=%st$' % w, linewidth=width)
    ax.plot(np.zeros(np.size(im)), t, zs=-lim_value, zdir='x', color='black')

    ax.text(0, -1.8, 0, s='微信公众号：二进制人工智能', zdir='x')

stick = np.arange(-1, 1.25, 0.25)
stick_list = ['{}i'.format(k) for k in stick]
plt.yticks(stick, stick_list)
ax.set_xlabel(r'Re（$cos(\theta)$）', font)
ax.set_ylabel(r'Im（$isin(\theta)$）', font)
ax.set_zlabel('t', font)
plt.xlim((-lim_value, lim_value))
plt.ylim((-lim_value, lim_value))
plt.legend()
plt.show()

角频率和角速度的区别[1]

1、定义不同

点沿着一个以弧度为单位的圆（一个圆周为$2\pi$,即：$360°=2\pi$)，在单位时间内所走的弧度即为角速度。角频率，也称圆频率，表示单位时间内变化的相角弧度值。

2、研究范围不一样

角频率是在任意的周期性运动中。角速度是在圆周运动中，或者至少是瞬时的圆周运动中。

3、物理意义不一样

角频率是单位时间内的振动次数（频率）与$2\pi$之积。点沿着一个以弧度为单位的圆，在单位时间内所走的弧度即为角速度。

4、计算方法不一样

角频率的计算公式是$w=2\pi /T$。角速度的计算公式为$w=v/r$，其中$v$为某时刻的线速度。

欧拉公式证明[2]

有许多方式可以证明欧拉公式，这里泰勒级数进行证明。

把函数$e^x$、$\cos x$和$\sin x$写成泰勒级数形式：
$$e^x=1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+\cdots$$

$$\cos x=1-\frac{x^2}{2!}+\frac{x^4}{4!}-\frac{x^6}{6!}+\cdots$$
$$\sin x=x-\frac{x^3}{3!}+\frac{x^5}{5!}-\frac{x^7}{71}+\cdots$$

将$x=i\theta$代入$e^x$，得到：
$$e^{i\theta }=1+i\theta +\frac{(i\theta {)}^2}{2!}+\frac{(i\theta {)}^3}{3!}+\frac{(i\theta {)}^4}{4!}+\frac{(i\theta {)}^5}{5!}+\frac{(i\theta {)}^6}{6!}+\frac{(i\theta {)}^7}{7!}+\frac{(i\theta {)}^8}{8!}+\cdots$$
$$=1+i\theta -\frac{{\theta }^2}{2!}-\frac{i{\theta }^3}{3!}+\frac{{\theta }^4}{4!}+\frac{i{\theta }^5}{5!}-\frac{{\theta }^6}{6!}-\frac{i{\theta }^7}{7!}+\frac{{\theta }^8}{8!}+\cdots$$
$$=(1-\frac{{\theta }^2}{2!}+\frac{{\theta }^4}{4!}-\frac{{\theta }^6}{6!}+\frac{{\theta }^8}{8!}-\cdots )+i(\theta -\frac{{\theta }^3}{3!}+\frac{{\theta }^5}{5!}-\frac{{\theta }^7}{7!}+\cdots )$$
$$=\cos \theta +i\sin \theta$$

我们还可以推出：

$$\sin \theta =\frac{e^{i\theta }-e^{-i\theta }}{2i}$$

$$\cos \theta =\frac{e^{i\theta }+e^{-i\theta }}{2}$$

参考：

[1 https://zhidao.baidu.com/question/327045010031888125.html

[2] https://blog.csdn.net/zoujiachi666/article/details/70943355

十分感谢三连支持！最靓的公式送给最靓的你：

$$e^{i\pi }+1=0$$