今天我决定用Python编写简单的程序,只是为了在考试之前练习。首先,我想绘制罪和cos的情节,这并不是那么难。但后来,我决定挑战自己,画出正切的情节。
import pylab as p
x= p.arange(-1.0,1.0,0.1)
y= (p.sin(2*p.pi*x)) / (p.cos(2*p.pi*x))
p.plot(x,y,'g-',lw=1)
p.show()
它返回......呃......这个:
正如您所看到的,它看起来更像是ECK图而不是切线图。有谁知道什么是错的?
答案 0 :(得分:4)
如果您增加x
中的点数,
import pylab as p
import numpy as np
x = p.linspace(-1.0, 1.0, 1000)
y = (p.sin(2 * p.pi * x)) / (p.cos(2 * p.pi * x))
p.plot(x, y, 'g-', lw=1)
p.show()
你得到这样的东西:
注意y-range
有多大。 Matplotlib无法显示切线曲线中的大部分小值,因为范围太大。
可以通过忽略渐近线附近的极大值来改善绘图。使用Paul's workaround处理渐近线,
import pylab as p
import numpy as np
x = p.linspace(-1.0, 1.0, 1000)
y = (p.sin(2 * p.pi * x)) / (p.cos(2 * p.pi * x))
tol = 10
y[y > tol] = np.nan
y[y < -tol] = np.nan
p.plot(x, y, 'g-', lw=1)
p.show()
你得到了
答案 1 :(得分:3)
import pylab as p
x= p.arange(-1.0,1.0,0.01) # <- 0.01 step size.
y= (p.sin(2*p.pi*x)) / (p.cos(2*p.pi*x))
# y = p.tan(2*p.pi*x)
p.plot(x,y,'g-',lw=1)
p.ylim([-4, 4]) # <- Restricting the ylim so we don't see the ~inf values.
p.show()
如果你没有设置ylim,这将是结果。 (值接近无穷大)
设置ylim的结果。