我是靠谱客的博主 隐形火龙果,最近开发中收集的这篇文章主要介绍Cellular Automaton-Python for Programmer : UC Berkeley SelfPace Course,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

好久没上CSDN了噢,今天觉得无论如何得学一下python了,捣鼓了一晚,好弱啊~~ 在这贴个代码,以备忘。

----------------------------------
题外话:这学期要设计基于FPGA的CPU,我组设计很大程度参考MIPS架构,发现UCB 61C的CPU设计部分的课件(注:itunes U是有视频的,不过看视频比较耗时。。)还是很给力的哈哈。进一步发掘资料,发现还有很多好资源哈哈。其中SelfPace Course就是非常好的资源,建议CSDN的朋友们多看看。这学期还有在做编译器实现(系里采用的是stanford 的cool 语言方案 CS143)和B/S的数据库应用系统的专题实验,有在做相关东西的同学大家多交流。
---------------------------------

UCB 课程目录       http://www-inst.eecs.berkeley.edu/classes-eecs.html
UCB SelfPace课程   http://www-inst.eecs.berkeley.edu/~selfpace/cs9honline/P1/
UCB Selfpace course: python for  
                             Project 1: Cellular Automaton

author: JRH @XJTU
2012-3-25

s= ''
times= 50
ruleNumber= 30

def  D2Binary8bit(i):
   def  Denary2Binary(n):
    bStr= ''
     if n== 0: return  '0'
     while n> 0:
      bStr= str(n% 2)+bStr
      n=n>> 1
     return bStr

   def  plusZero(s):
     while  8len(s)> 0:
      s= '0'+s 
     return s
  s=Denary2Binary(i)
  strTemp=plusZero(s)
   #print i
   return strTemp

ruleBinary=D2Binary8bit(ruleNumber)
#print "ruleBinary is"
#print ruleBinary

def  rule(bit1,bit2,bit3):
   if bit1== 1  and bit2== 1  and bit3== 1: return ruleBinary[ 0]
   if bit1== 1  and bit2== 1  and bit3== 0: return ruleBinary[ 1]
   if bit1== 1  and bit2== 0  and bit3== 1: return ruleBinary[ 2]
   if bit1== 1  and bit2== 0  and bit3== 0: return ruleBinary[ 3]
   if bit1== 0  and bit2== 1  and bit3== 1: return ruleBinary[ 4]
   if bit1== 0  and bit2== 1  and bit3== 0: return ruleBinary[ 5]
   if bit1== 0  and bit2== 0  and bit3== 1: return ruleBinary[ 6]
   if bit1== 0  and bit2== 0  and bit3== 0: return ruleBinary[ 7]

def  update():
  temp= ''
  temp=temp+ '1'
   global times
   global s
  each= 2
   while each<times* 2+ 1:
    temp = temp+ str(rule( int(s[each- 2]), int(s[each- 1]), int(s[each])))
    each+= 1
  temp=temp+ '1'
  s=temp
   return temp



def  ini() : 
  t= 0
   global s
   while t< times* 2+ 1:
     if t==times:
      s+= '1'
      t+= 1
     else:
      s+= '0'
      t+= 1
   print s


def  main():
   print  "P1", times* 2 + 1,times+ 1
  ini()
  index= 0
   while index < times:
     print update()
    index+= 1

if __name__== '__main__':main()


-----------------------------------

对times和ruleNumber变量进行修改,可以得到不同规则(ruleNumber)下不同迭代次数(times)的结果:

如将times=3 , ruleNumber=30,本程序输出为:

 P1 7 4
 0001000
 1011101
 1010001
 1011011

程序输出为这样的格式,只要利用linux 的 > 到bpm文件(Portable Bitmap),即可制作一副由该Cellular Automaton 产生的图像了。下面简要介绍一下bpm:

 Portable Bitmap Format   例如,一个7*3 的bpm图像


第一行内容:  P1  宽  长

P1 7 3
0001000
0011100
0110010

---

下面将程序的输出重定向到 rule30.bpm文件中:

$ python project1.py > rule30.pbm

成功后即可用一般的图像查看软件查看该图像了。 还可有修改times和ruleNumber,产生更多图案。更多Cellular Automaton 内容,请参考http://atlas.wolfram.com/TOC/




最后

以上就是隐形火龙果为你收集整理的Cellular Automaton-Python for Programmer : UC Berkeley SelfPace Course的全部内容,希望文章能够帮你解决Cellular Automaton-Python for Programmer : UC Berkeley SelfPace Course所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(43)

评论列表共有 0 条评论

立即
投稿
返回
顶部