我的家庭作业让我写一个代码,模拟泰勒教科书在经典力学中的数字。如果有人有兴趣知道,请拨打和。分岔图没有绘图/剧情只是没有出现
我能够重现他们中的一个是下面的代码(这可能是我确实遇到问题的代码很好的参考):
import nympy as np
import matplotlib.pyplot as plt
# We need to calculate the first fixed point
r1=np.array(np.arange(0,4,0.09))
x1 = np.zeros((len(r1),1))
# Now calculating the second fixed point
r2=np.array(np.arange(1,4,0.1))
x2 = (r2 -1)/r2
# Now finding when the fixed points split up again
r3=np.array(np.arange(3,4,0.1))
y1 = (((r3**2 - 2*r3 - 3)**0.5) + 1 + r3)/(2*r3)
y2 = ((-(r3**2 - 2*r3 - 3)**0.5) + 1 + r3)/(2*r3)
# Now finding the experimental values for 1/2 of a split
x3 = []
for r in np.arange(0,4,0.09):
x = 0.666
for i in range(100):
x = (r**2) * x * (1.0 -x) - (r**3) * (x**2)*((1-x)**2)
x3.append(x)
# Doing the same as above second 1/2
x4 = []
for r in np.arange(0,4,0.09):
x = 0.8
for i in range(100):
x = (r**2) * x * (1.0 -x) - (r**3) * (x**2)*((1-x)**2)
x4.append(x)
plt.plot(r1,x3,'bo', label='Experimental')
plt.plot(r1,x4,'bo')
plt.plot(r3,y2,'k-')
plt.plot(r3,y1,'k-')
plt.plot(r1,x1,'k-', label='Theoretical')
plt.plot(r2,x2,'k-')
plt.legend(loc=2)
plt.show()
,这里是第二图像的代码这似乎并不奏效。我不知道为什么。任何帮助,将不胜感激。这个数字只是没有绘制,我不知道为什么。
import numpy as np
import matplotlib.pyplot as plt
for r in n.arange(2.8,4,0.01):
x = 0.5
for i in range(150):
x = r*x*(1-x)
if i >= 125:
plt.plot(r,x,'k')
plt.xlim (2.8,4)
plt.show()
欢迎所以堆栈溢出。请阅读[如何提出一个好问题](http://stackoverflow.com/help/how-to-ask)。另外,请粘贴您的代码,而不是插入图片。在帖子中有一个“{}”按钮,可以将所有内容缩进四个空格,并显示为代码 –
您是否在使用IPython笔记本? – Xevaquor