作者granzi (乌木)
看板Statistics
标题[程式] Python写的Yate's Algorithm
时间Tue Aug 4 01:52:12 2009
------------------------------------------------------------------------
[软体程式类别]:
Python
[程式问题]:
实作Yate's Algorithm
[软体熟悉度]:
中
[问题叙述]:
实作Yate's Algorithm
输入2^k factorial design the 结果
输出Yates Algorithm 的结果 *[m
[程式范例]:
from array import array
import math
def YatesAlgorithm(dataArray):
dataLength = len(dataArray)
effectsNum = math.log(dataLength, 2)
dataNew = dataArray[:]
if ((effectsNum - int(effectsNum)) != 0.0):
print("The length of data array must be 2^k")
else:
for i in range(1, int(effectsNum)+1):
dataTmp = dataNew[:]
for j in range(1, int(dataLength/2)+1):
temp1 = dataNew[2*(j-1)] + dataNew[2*j-1]
temp2 = dataNew[2*j-1] - dataNew[2*(j-1)]
dataTmp[j-1] = temp1
dataTmp[int(dataLength/2)+j-1] = temp2
print(dataNew)
dataNew = dataTmp[:]
return(dataNew)
-----------------------------------------------------------------------------
谢谢指教
--
凡发生之事必合理
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 67.221.66.201