《高等数学》同济大学出版:三角函数
编写 test_sin_cos.py 如下
# -*- coding: utf-8 -*-
""" 正弦函数 y=sin(x) 和 余弦函数 y=cos(x) 的曲线 """
import numpy as np
from matplotlib import pyplot as plt
# 用于正常显示中文标题,负号
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False
x = np.arange(-7.0, 7.01, 0.01)
y = np.sin(x)
y1 = np.cos(x)
print(len(x))
print('sin(x):', y[0], y[-1])
print('cos(x):', y1[0],y1[-1])
# 可视化
fig = plt.figure()
axes = fig.add_subplot(111)
axes.plot(x, y, label='sin(x)') # 画曲线
axes.plot(x, y1, label='cos(x)')
#axes.axis('scaled') # 用缩尺制图
axes.axis('equal')
plt.title('正弦函数 y=sin(x) 和 余弦函数 y=cos(x) 的曲线')
plt.xlabel('x')
plt.ylabel('y')
plt.xticks(range(-7,8,1), range(-7,8,1))
plt.yticks(range(-2,3,1), range(-2,3,1))
plt.legend()
plt.grid()
plt.show()
运行 python test_sin_cos.py
当0<=x<=2 时,y=sin(x) 与 y=cos(x) 的交点坐标是 (π/4, √2/2),
大约是 (0.785398, 0.70710678)。
编写 test_tg_tan.py 如下
# -*- coding: utf-8 -*-
""" 正切函数 y=tan(x) 的曲线 """
import numpy as np
from matplotlib import pyplot as plt
# 用于正常显示中文标题,负号
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False
pi = np.pi
print('pi/2:', pi/2)
x = np.arange(-1.56, 1.57, 0.01)
y = np.tan(x)
print(len(x))
print('tan(x):', y[0], y[-1])
# 可视化
fig = plt.figure()
axes = fig.add_subplot(111)
axes.plot(x, y, label='tan(x)') # 画曲线
#axes.axis('scaled') # 用缩尺制图
#axes.axis('equal')
plt.title('正切函数 y=tan(x) 的曲线')
plt.xlabel('x')
plt.ylabel('y')
plt.legend()
plt.grid()
plt.show()
运行 python test_tg_tan.py
看起来不像《高等数学》书上印的,尽信书?可否
编写 test_cotan.py 如下
# -*- coding: utf-8 -*-
""" 余切函数 y=cotan(x) 的曲线 """
import numpy as np
from matplotlib import pyplot as plt
# 用于正常显示中文标题,负号
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False
pi = np.pi
print('pi/2:', pi/2)
x = np.arange(0.01, pi, 0.01)
# cotan(x) = 1/tan(x) 求倒数
y = np.reciprocal(np.tan(x))
print(len(x))
print('cotan(x):', y[0], y[-1])
# 可视化
fig = plt.figure()
axes = fig.add_subplot(111)
axes.plot(x, y, label='cotan(x)') # 画曲线
#axes.axis('scaled') # 用缩尺制图
#axes.axis('equal')
plt.title('余切函数 y=cotan(x) 的曲线')
plt.xlabel('x')
plt.ylabel('y')
plt.legend()
plt.grid()
plt.show()
运行 python test_cotan.py