1

模块 OneDMaps： 

def LogisticMap(a,nIts,x):
    for n in xrange(0,nIts):
        return 4.*a*x*(1.-x)

实际程序： 

# Import plotting routines
from pylab import *
import OneDMaps 

def OneDMap(a,N,x,f):
    return x.append(f(a,N,x))

# Simulation parameters
# Control parameter of the map: A period-3 cycle
a = 0.98
# Set up an array of iterates and set the initital condition
x = [0.1]
# The number of iterations to generate
N = 100

#function name in OneDMaps module
func = LogisticMap

# Setup the plot
xlabel('Time step n') # set x-axis label
ylabel('x(n)') # set y-axis label
title(str(func) + ' at r= ' + str(a)) # set plot title

# Plot the time series: once with circles, once with lines
plot(OneDMap(a,N,x,func), 'ro', OneDMap(a,N,x,func) , 'b')

该程序应该从模块 OneDMaps.py 调用一个函数，然后根据它的迭代绘制它。我收到错误“无法将序列乘以浮点类型的非整数”，并且我尝试使用 LogisticMap(float(a)...) 但这不起作用。此外，我希望函数名称显示在绘图的标题中，但我得到“在 r=0.98 处，而不是在 r=0.98 处显示 LogisticMap。

4

1 回答 1

4

你设置一个list这样的：

x = [0.1]

然后你将它乘以一个float：

return 4.*a*x*(1.-x)

你不能那样做。也许你想x成为一个array而不是一个list？

x = array([0.1])

（这将按元素进行乘法）

请注意，添加列表会连接：

[0] + [1] == [0, 1]

乘以一个整数与多次连接相同：

[0] * 3 == [0, 0, 0]

但这对花车没有意义， 

[0] * 3.0 #raises TypeError as you've seen

（例如，如果乘以 3.5，您应该得到什么？）

于 2012-12-09T04:02:41.883 回答