My Hacking Stuff

如何選擇Machine Learning的學習方法

2011-11-25T19:53:00.001-08:00

[image from Stanford Machine Learning Class (http://ml-class.org/)]

Machine Learning中的學習方法有很多，

而使用哪種方法也是我一直不懂的，

今天看到這份文件而有了大概。

簡而言之，

我們應該要先考慮「資料量的大小」和「特徵的數量」而使用不同Learning Method，

又說SVM幾乎沒有人自己寫，拿現成的用(e.g. libsvm)即可，

類神經網路則相對難一點，但是普遍有不錯的效果。

Life Timer

2011-11-20T09:10:00.001-08:00

I wrote a website which counts down the time to users' death.

http://stanford.edu/~heron93/LifeTimer/

By using JavaScript, html, you can easily grab user's input, and calculate to distance from now to the death time.

We use life expectancy data form Wikipedia (67.2) to determine the time of death.

If you are older the 67.2 years old, this site will always tell you that you'll die in one second, which means you are already older than average age, and please treasure the time you have.

Time and tide wait for no man.

初學Coding注意事項

2011-10-29T00:34:00.000-07:00

最近看了不少人的code，有一些心得，覺得可以先至少注意以下幾點：

每寫一段落就要測試是否正確才能往下寫。
每寫一件事情就寫上註解(comment)。
排版一定要正確，上下引號一起打，忘記分號很好笑。
Debug時，印出你關心的變數看看是否跟想的一樣。
想不通的時候一樣要在紙上拿筆想，可以畫圖或寫字。

為什麼要這麼做的原因不重要，重要的是初學時一定要這麼做！

Install Fedora(Linux) in Virtualbox under Windows

2011-06-06T10:38:00.000-07:00

好多人問過我這個，
所以錄製成影片希望有幫助。

LED Cube (444)

2011-06-03T10:36:00.000-07:00

製作這個Cube最早的動機乃希望找尋一個3D模型輸出方式，

而其實這的確可以算是一種低畫素、單一顏色的實體3D模型輸出(3D printing)。

製作方法：

利用LED現有的針腳，

焊接使得同一層LED共陰，

每一直條（垂直地面）則共陽，

於是總共可以得到：

4個陰極（4層）
16個陽級（16個構成一個平面）

使用Arduino控制16個陽級是否提高電壓（給電），

而陰極則輪流接地，

如此每一層會輪流亮，

而當速度很快、

快到每一週期小於視覺暫留時間，

整個Cube就可以任意發亮。

找個木板、鑽動以固定LED（倒插），也許覺得多餘，但是實際操作會發現沒有固定根本就無法焊接成形。（間距2.5cm）
同一層要把陰極焊接在一起。
注意絕對不能把正負極焊接在一起！
做完四層之後，將正極焊接（垂直）起來，Cube外型就完成了。
把線接出來。
連結到Arduino、寫程式控制之。

大功告成。

打地鼠遊戲機

2011-05-06T06:02:00.000-07:00

遊戲介紹

LED以亂序的方式亮起
按亮的LED之對應按鍵即可消去、另一個LED亮起，同時發出短聲響
有人性化的按鈕擺放安排

製作步驟

準備一張紙板當作基底
根據手指結構畫記按扭擺放位置
鑲上按鈕
焊接上單心線
每個按鍵上方裝上LED
把線製成排線
撰寫程式、連接上AVR單晶片、大功告成！

Wii Remote 紅外線裝置

2011-04-21T03:52:00.000-07:00

以Wii Bar做紅外線輸出，
Wii Remote前端CMOS感光元件做為輸入，
讀取Wii Remote前有紅外線出現的位置，
將運用來判斷相對位置。

Wii Bar的紅外線二極體(IR LED)左右各五個，
但是Wii Remote只有拍攝到兩個點乃因為設計上左右的五個LED角度安排導致，
如此的設計可以增加紅外線的強度，
維基上寫此紅外線可達5m，
而我以一顆最常見(1.5V)的IR LED也可以達到同樣效果，
但是距離只能約3.3m(大於2.5m後訊號不穩定)。

More shaking, bigger circle

2011-04-17T03:28:00.000-07:00

在Linux環境底下，
讀取Wii remote的三維加速度值，
取三維加速度值之平均對應至圓的半徑，
即當remote搖晃越劇烈，圓會越大！

使用：
● cwiid(python)
● pygame
● wii remote

在python中匯入cwiid,pygame模組，
接著讀入Wii remote(wm)，
即wm=cwiid.Wiimote()
而wm.state['acc']就有當前狀況下的remote加速度，
wm.get_acc_cal裡則有兩組常數，
視為scale wm.state['acc']至實際值所需的上下界(min,max)

Formula:

acc=(x-xmin)/(xmax-xmin)

其中

● x=wm.state['acc']

● xmin=wm.get_acc_cal[0]

● xmax=wm.get_acc_cal[1]

(以上皆為虛擬碼)

藍天對對樂

2011-04-14T07:31:00.000-07:00

藍天即代表藍天之子－－附中。

使用Visual Basic語言，
花費約一個晚上時間完成，
將對新生做學校介紹時使用。

善用物件導向中各個物件的屬性、發生事件等，
配合需求寫成，並不難。

推翻

2011-02-20T08:10:00.000-08:00

「生命中不能承受的輕」一書，
觀書有感阿！

書裡提到，
生命中有很多的岔路，
有不同的選擇，
但是我們無法判斷好壞，
因為沒有人可以回到從前，
選擇另一條路走走看。

所謂的好壞是比較而來，
既然沒有辦法比較，
就沒有好壞。

如此一來，
其實我們一輩子面臨的選擇，
都沒有好壞之分，
若是有，
不過是自己以為罷了。

這種討論，
即是科學上的實驗設計問題，
無法創造出對照組與實驗組，
所以後頭的什麼也解釋不了，
至於無法創造那兩組的原因，
不單單只是「生命不可重複」
更可以延伸為「時間不可倒流」

於是我又想，
任何建立在時間上的實驗，
是否因此都被推翻？

現有的知識裡，
我們不能使時間倒流，
因而，
我們不可能創造一樣環境的實驗組與對照組，
以已知的四維度解釋，
因為兩個「不同的東西」，
遠永不可能「相同」，則
前三維度相同則第四維度一定不同，
第四維度相同則前三維度一定不同，
在不能創造兩相同環境下，
所有的實驗將被推翻，
或是說它們只不過是「假設成立」而已。

ODCR using similarity 2 (使用相似度與特徵判別數字)

2011-02-06T23:30:00.000-08:00

這回不直接跟標準影像比較相似度，
而間接透過「特徵向量」比較。

方法：
1.運算標準影像(數字0~9)的特徵向量
2.運算運判斷之影像的特徵向量
3.計算上述兩者在多為空間中的距離（方根）
4.距離最小者及猜測為該數字

p.s.此程式擷取50個特徵

ODCR using similarity (使用相似度判別數字)

2011-02-04T21:02:00.000-08:00

這是一個簡單製作OCR的方法，
針對數字部份，
我們比較將要判斷之影像與標準影像的相似度，
找到最相近的則可以判斷數字。

此方法的前提：
1.有標準影像存在
2.要判斷之影像完整

Timer for Rubik's Cube

2010-08-07T04:58:00.000-07:00

p.s.這篇文章約一年前所撰

恩...

那是國二時候的事情了吧...

那時候開始學魔術方塊，

從看著公式轉到最後可以轉10幾秒，

也從班上3個人迷，

到最後有10個以上都會轉成功呢...！

我記得當時的我想要一個魔術方塊的計時器，

他的設計是這樣的：

手放在計時器上為準備，

離開之後開始計時，

再碰一次是停止計時。

不過那時網路上賣的一個要2000元，

父母親看到都嚇到了，

我還是吵著要買，

其實我也嚇到了ˊˋ...

最後，

當然沒有買成囉...

不過在今天暑假，

終於開始我的Microcontroller的旅程，

也寫出許許多多簡單的程式燒再裡面玩，

忽然想起一兩年前的樣子，

於是把現成的東西拼拼揍揍很快就變成那魔術方塊的計時器了。

扣除時間、知識、人力、分消耗之工具外，

一整套下來不會超過200元...

整整差了10倍...

以下是一些介紹，

可以選擇性的看：

影片中有四個主要部分：

1.Arduino板（核心）

2.四個七段顯示器

3.電池＝＝（先前用的是接到插座的變壓器）

4.開關

其實有做過晶片的人會知道很簡單，

用C語言把程式燒在AVR的單晶片上面，

在把八個針腳接到那七段顯示器上面，

在四個針腳決定是哪一個七段顯示器（有四個），

這樣就差不多了...

原來 1=0 (part2)

2010-06-29T07:27:00.000-07:00

_____________（接續第一集）_______________

過了一週，

可愛的老人笑笑的又對我說：「is 0＝ 1？」

呃...

被刺了一刀，

這好陰險...

才剛討論完0.999......是否等於1

而這「0＝1」卻是多麼的心機...

用了之前的結論，

「我們為了方便，定義0.999......＝1」，

然而此時兩邊同減去1，

再加負號，

便是0.000......1＝0，

接著：

竟然，

我們在無意之間也定義了「1＝0」...！

於是這結果使人慌張...

若1＝0，

兩邊分別乘上1、2、3、...等，

就可以很輕易的證明0＝1＝2＝3＝4＝.......

老人也慌了，

此時有個自稱布林的小姐走過來，

說：「在計算機裡，1 表示true，0 表示false」

（p.s.計算機泛指所有的計算元件，如電腦...等）

布小姐說的沒有錯，

1、0可能是一個開關的開或關，

或是一個電容的放電、充電等，

如今，若他們相等，

這理論可以讓世界上所有的電腦瓦解，

更可以讓科學界近乎所有的理論瓦解...

人類呀，

聰明的人類怎麼可能會做出這種事情？

竟然定義了一個毀掉所有科學的科學...

從定義到推廣，

我們必須要再回到定義，

最最根本的問題出在，

過去那些偉大的科學家們說：「無限大分之一為零」

而所謂的無限大是：沒有極限的...

沒有極限的，

然而我剛剛不是才在討論第無限項？...

其實，

當我們在討論最後一項（第無限項）發生的事情時，

他已經變成是「有限項」了。

無限是無法看到的，

是無法被討論的，

於是，

讓我們在從新整理一次思緒。

原來，

是人們定義了一個不能被討論的科學，

於是我們可以為那最後一項的變化隨需求而任意的定義。

既然，

科學無法完整的解釋這名詞，

再換個方式吧...

以天地永恆的時間作為分母，

而你存在的時間量放置分子，

在「無限大分之一為零」的定義下，

這個值是，是零。

換句話說，

你根本沒有存在過...

「沒有存在過...」老人小小聲的反覆我說的話，

「年輕人，謝謝你...我想我懂了...」天空有點雨絲...

要說感恩的是我，

讓我又了解了更多，

我說的並不是科學。

原來 1=0

2010-06-29T03:35:00.000-07:00

這大約是一年前的故事了，

有個滿臉鬍渣的老老人跑來與我討論：

「is 0.999......＝1？」

當下，

我不假思索的回答：「當然」。

然而，

這個問題卻讓我想起當初接觸到這東西的時候，

「當然」不應該說的這麼肯定，

明明不一樣的東西怎麼會相等呢？

...於是，和那老人聊到近半夜，

嘗試去找0.999......＝1的證明中的瑕疵，

卻發現其實近在眼前。

大家所聞的大部分會是以下的作法：

恩...

也許就這樣你便相信了。

後面的"......"其實隱藏的很多，

也訴說了為什麼0.999......＝1，

於是老人說：把...換成有"999"的方式出現吧！

注意到第無限位，

其實*10之後便會是空的。

原先的作法並沒有注意到這問題，

應該說是將第無限位和第無限-1位「視」為沒有差別的，

如此，

可以解決不少惱人的問題，

到這裡，

老人與我有了一個簡單的結論：

「0.999......＝1是被定義出來的」

而為何要這樣定義呢？

還記得高一的limit裡面的→（approach）？

因為我們說0.999.....＝1所以才可以解決無窮等比，

也因為他們等於所以才有微積分，

才有電腦螢幕上美麗的動畫。

恩...
0.999......＝1是被定義出來的

過了一週，

可愛的老人笑笑的又對我說：「is 0＝ 1？」

_____________（未完待續）______________

動態判讀

2010-06-23T02:54:00.000-07:00

看過「關鍵報告」嗎？

電影裡的人物，

不用持有任何的遙控器、

更不需鍵盤滑鼠等輸入裝置，

就可以用手勢控制大螢幕上的系統。

比起一般透過需要摸到的操控介面，

他是無須直接接觸到輸入裝置的----攝影機

隨著視訊的普及，

新一代的電腦，

網路攝影機是必備的裝置，

即使是舊的電腦，

也可以容易的買到市面上不超過500元的一般網路攝影機。

透過網路攝影機，

擷取一連串的影像，

就可以做到電影中神奇的事情了。

簡單來說，

分析「移動中」的物體，

也就是「與前一畫面不同」的畫素，

即是一種簡易的方式。

（左上為現在影像、右上為背景、中下為兩影像差）

這回，

我利用opencv+python+vpython，

讓程式啟動時，

先抓一張影像作為背景，

在讓之後的影像跟這一張背景比較（計算絕對值），

分析不同處。

知道哪裡不同之後就可以做讓電腦做反應了，

比如一有人經過，

電腦便發出聲音。

而我則是配合上vpython，

抓三個點，

若是有人從右走到左，

螢幕上的三個球，

就從右邊開始變大，

也就是說，

人移動到哪裡，

哪裡的球會變大。

import sys
#import
    #for opencv
from visual import *
from opencv import cv
from opencv import highgui
    #for music playing
"""
import pygame
pygame.init()
sound=pygame.mixer.Sound("beep.wav")
"""

if len(sys.argv)>=2 and sys.argv[1]=="demo_on":
    demo_on=True
else:
    demo_on=False

#arguments
sensitive=20

#window define
if demo_on:
    highgui.cvNamedWindow("Camera",highgui.CV_WINDOW_AUTOSIZE)
    highgui.cvNamedWindow("Background",highgui.CV_WINDOW_AUTOSIZE)
    highgui.cvNamedWindow("Foreground",highgui.CV_WINDOW_AUTOSIZE)

#set window position
if demo_on:
    highgui.cvMoveWindow("Camera",0,0)
    highgui.cvMoveWindow("Background",600,0)
    highgui.cvMoveWindow("Foreground",300,400)

#capture setting
capture=highgui.cvCreateCameraCapture(0)
highgui.cvSetCaptureProperty(capture,highgui.CV_CAP_PROP_FRAME_WIDTH, 160)
highgui.cvSetCaptureProperty(capture,highgui.CV_CAP_PROP_FRAME_HEIGHT,120)

#capture first capture
frame=highgui.cvQueryFrame(capture)
frame_size=cv.cvGetSize(frame)

#set background & foreground -> img and mat
background_img=cv.cvCreateImage(frame_size,cv.IPL_DEPTH_8U,1)
foreground_img=cv.cvCreateImage(frame_size,cv.IPL_DEPTH_8U,1)

background_mat=cv.cvCreateMat(frame.height,frame.width,cv.CV_32FC1)
foreground_mat=cv.cvCreateMat(frame.height,frame.width,cv.CV_32FC1)
frame_mat=cv.cvCreateMat(frame.height,frame.width,cv.CV_32FC1)

#convert img to mat (color: RGB -> Gray)
cv.cvCvtColor(frame,background_img,cv.CV_BGR2GRAY)
cv.cvCvtColor(frame,foreground_img,cv.CV_BGR2GRAY)

cv.cvConvert(foreground_img,frame_mat)
cv.cvConvert(foreground_img,background_img)
cv.cvConvert(foreground_img,foreground_img)

cv.cvConvert(background_img,background_mat)

#setting for visual python
#left and right are invert
ball=[
        visual.sphere(red=255,pos=(-5,0,0),color=color.red),
        visual.sphere(green=255,pos=(0,0,0),color=color.green),
        visual.sphere(blue=255,pos=(5,0,0),color=color.blue)
        ]

while 1:
    frame=highgui.cvQueryFrame(capture)
    if frame is None:
        break

    cv.cvCvtColor(frame,foreground_img,cv.CV_BGR2GRAY)
    cv.cvConvert(foreground_img,frame_mat)

    cv.cvAbsDiff(frame_mat,background_mat,foreground_mat)
    cv.cvThreshold(foreground_mat,foreground_img,100,255.0,cv.CV_THRESH_BINARY)

    #main stuff
    consider=[
            foreground_mat[frame.height/2][frame.width/4*3],
            foreground_mat[frame.height/2][frame.width/4*2],
            foreground_mat[frame.height/2][frame.width/4*1]
            ]

    #foreground_mat

    index=0
    for i in consider:
        print str(i)+ "  ",
        """
        print index,
        if  i >= sensitive:
            print " X |",
        else :
            print "   |",
        """
        ball[index].radius=2*(i/255.0)+1.0
        index+=1
    print

    if demo_on:
        highgui.cvShowImage ("Camera", frame)
        highgui.cvShowImage ("Background", background_img)
        highgui.cvShowImage ("Foreground", foreground_img)

    """
    k = highgui.cvWaitKey(10)
    if k == '\x1b':
        break
    """
    #reset background
    """
    if k == 's':
        cv.cvCvtColor(frame,background_img,cv.CV_BGR2GRAY)
        cv.cvConvert(foreground_img,background_img)
        cv.cvConvert(background_img,background_mat)
        print "background reset."
        """

if demo_on:
    highgui.cvDestroyWindow("Camera")
    highgui.cvDestroyWindow("Background")
    highgui.cvDestroyWindow("Foreground")

#end

運用這個技術，

還有更多東西值得挑戰，

陸續進行中。

常態分佈

2010-06-19T05:39:00.000-07:00

配合上一回的信心水準，

這回，

我利用python的亂數模擬一個事件（例如Heron帥不帥），

一百（若是太大則逼近實際值，則沒有意義）為一個單位，

並計算每一單位裡的支持率（例如，覺得Heron很帥的人數與總人數之比值），

重複這一單位的動作100,000次（此量值取決於電腦的處理量100*100,000），

裡用這100,000次得到的支持率使用vpython繪畫一張曲線圖，

而就可以看到常態分布（高斯分布）囉！！

Code:

import random
from visual.graph import *
from math import sqrt

#create empty array
percent=[]
for i in range(101):
    percent.append(0)

#create main data array
n_test=100000
for i in range(n_test):
    pos=0
    #use 100 for each random. sample to avoid float
    for j in range(100):
        if random.random()<0.5:
            pos+=1
    percent[pos]+=1


#picture drawing
##setting
graph=gdisplay(\
        x=0,y=0,\
        width=800,height=600,\
        title="normal distribution",\
        xtitle="percent",ytitle="times",\
        foreground=color.green,\
        background=color.black,\
        xmax=100,\
        xmin=0,ymin=0)
curve=gcurve(display=graph.display)
y_max=0
for i in range(101):
    print "[%d] : %d" % (i,percent[i])
    curve.plot(pos=(i,percent[i]),color=color.green)
    if percent[i]>y_max : y_max=percent[i]

#try to create 68-85-99.7
print "use real data to test testing data"

##real data
p=0.5
n=100
u=sqrt(p*(1-p)/n)

##bar drawing
bar_68=[
        gvbars(display=graph.display,color=color.red,delta=0.2).plot(pos=((p-u)*100,y_max+100)),
        gvbars(display=graph.display,color=color.red,delta=0.2).plot(pos=((p+u)*100,y_max+100))
        ]
bar_95=[
        gvbars(display=graph.display,color=color.yellow,delta=0.2).plot(pos=((p-2*u)*100,y_max+100)),
        gvbars(display=graph.display,color=color.yellow,delta=0.2).plot(pos=((p+2*u)*100,y_max+100))
        ]
bar_997=[
        gvbars(display=graph.display,color=color.blue,delta=0.2).plot(pos=((p-3*u)*100,y_max+100)),
        gvbars(display=graph.display,color=color.blue,delta=0.2).plot(pos=((p+3*u)*100,y_max+100))
        ]

##clac
area_1=0#68%
area_2=0#95%
area_3=0#99.7%
for i in range(101):
    consider=float(i)/100
    if consider>=p-u and consider<=p+u:   area_1+=percent[i]
    if consider>=p-2*u and consider<=p+2*u:   area_2+=percent[i]
    if consider>=p-3*u and consider<=p+3*u:   area_3+=percent[i]

print "area_1 n=" + str(area_1)
print "area_2 n=" + str(area_2)
print "area_3 n=" + str(area_3)

p_area_1=float(area_1)/n_test*100
p_area_2=float(area_2)/n_test*100
p_area_3=float(area_3)/n_test*100

print "68%%   area (red)   >>%f" % p_area_1
print "95%%   area (yellow) >>%f" % p_area_2
print "99.7%% area (blue)  >>%f" % p_area_3

print "Exit."

信心水準計算機

2010-06-19T04:17:00.000-07:00

統計學上，

表示當擁有一組樣本數據之後，

可以有68%的信心說真實的統計值在平均直的加減一個標準差之間，

95是兩個邊準差，

99.7是三個。

於是我簡易用python的做了一個計算機，

輸入欲統計的數值（例如贊成Heron聰不聰明），

之後運算出其信賴區間。

計畫要利用亂數配合上vpython繪畫亂數產生事件之機率，

看看是不是有常態分布的出現。

from math import sqrt

print
print "=================================================="
print "||Welcome to Heron's Confidence-level Calculator||"
print "=================================================="
print

pas=float(raw_input("Enter n+ >>"))
neg=float(raw_input("Enter n- >>"))

n=pas+neg
p=(pas)/n
u=sqrt((1-p)*p/n)

print "total n :%.0f" % n
print "p(n+):%f" % p
print "u :%f"%u

print

print "under 68%%   confidence :from %f to %f" % (p-u,p+u)
print "under 95%%   confidence :from %f to %f" % (p-2*u,p+2*u)
print "under 99.7%% confidence :from %f to %f" % (p-3*u,p+3*u)

print
print "Exit."
print

Vpython + Chemistry

2010-06-14T02:41:00.001-07:00

有鑑於學習高中選修化學第二章（分子間作用力和固體的結構）時，

老師教學上很難敘述三維的分子結構，

自己學習上雖能夠聽懂，

但是某些狀況下卻是有一種「看不到東西就不能說服我」的感覺，

於是花約2小時的時間，

把正在學習中的python程式能力運用在此。

程式碼下載(.py)

介紹影片：

Code :

#!/usr/bin/env ipython
from visual import *

#welcome messenge
print
print "================================================="
print "||Welcome to Heron's Little Chemical Program!! ||"
print "================================================="
print

#cheating or not?
cheat_on = (raw_input("cheat or not ?? (y/n) >> ") == "y")
print

#choose atom type
atom_type=["","ion","metal"]
print "choose atom type"
for i in range(1,len(atom_type)):
    print str(i)+"."+str(atom_type[i])+" ",
atom_type_chosen=str(atom_type[int(raw_input(">> "))])
print

#choose structure type
print "choose structure type"
if atom_type_chosen == "ion":
    structure_type=["","Body-centered cubic","Face-centered cubic","Zincblende"]
elif atom_type_chosen == "metal":
    structure_type=["","hcp","Face-centered cubic","Body-centered cubic"]
else:
    print "System Error, call Heron now...!"
    raise SystemExit
for i in range(1,len(structure_type)):
    print str(i)+"."+str(structure_type[i])+" ",
structure_type_chosen = str(structure_type[int(raw_input(">> "))])
print

#insert length
l=float(raw_input("l= "))
if atom_type_chosen == "ion":
    r_p=float(raw_input("(r+)= "))
    r_m=float(raw_input("(r-)= "))
elif atom_type_chosen == "metal":
    r=float(raw_input("r= "))
    r_p=r_m=r
else:
    print "System Error, call Heron now...!"
    raise SystemExit

l_h=l/2

#print out messenge
print "-------------------------------------------"
print "atom type = " + atom_type_chosen
print "structure type = " + structure_type_chosen
print "l=%f r_p=%f r_m=%f" % (l,r_p,r_m)
print "-------------------------------------------"

if structure_type_chosen =="hcp":
    print "hcp structure is not complete yet...\n"
    print "Exit..."
    raise SystemExit

elif structure_type_chosen == "Body-centered cubic" or\
     structure_type_chosen == "Face-centered cubic":

    #start
    print "drawing ..."

    #center
    if not(atom_type_chosen == "metal" and\
           structure_type_chosen == "Face-center cubic"):
        center=sphere(radius=r_p,color=color.red,opacity=0.6)
        if atom_type_chosen == "metal":
            center.color=color.green
        if cheat_on : center.material=materials.earth

    #corner
        #position
    corners=[
        sphere(pos=( l_h, l_h, l_h)),
        sphere(pos=( l_h, l_h,-l_h)),
        sphere(pos=( l_h,-l_h, l_h)),
        sphere(pos=(-l_h, l_h, l_h)),
        sphere(pos=(-l_h,-l_h, l_h)),
        sphere(pos=( l_h,-l_h,-l_h)),
        sphere(pos=(-l_h, l_h,-l_h)),
        sphere(pos=(-l_h,-l_h,-l_h))
        ]
        #settings
    for corner in corners:
        corner.color=color.green
        corner.opacity=0.6
        corner.radius=r_m
        if cheat_on:
            if atom_type_chosen == "ion":
                corner.material=materials.wood
            else:
                corner.material=materials.earth

    #face & edge for "Face-centered"
    if structure_type_chosen == "Face-centered cubic":
        #faces
            #position
        faces=[
                sphere(pos=( l_h,   0,   0)),
                sphere(pos=(   0, l_h,   0)),
                sphere(pos=(   0,   0, l_h)),
                sphere(pos=(-l_h,   0,   0)),
                sphere(pos=(   0,-l_h,   0)),
                sphere(pos=(   0,   0,-l_h))
                ]
            #setting
        for face in faces:
            face.color=corners[0].color
            face.radius=r_m
            face.opacity=0.5
            if cheat_on:
                if atom_type_chosen == "ion":
                    face.material=materials.marble
                else:
                    face.material=materials.earth

        #edge
            #position
        if atom_type_chosen != "metal":
            edges=[
                    sphere(pos=(   0, l_h, l_h)),
                    sphere(pos=(   0, l_h,-l_h)),
                    sphere(pos=(   0,-l_h, l_h)),
                    sphere(pos=(   0,-l_h,-l_h)),
                    sphere(pos=( l_h,   0, l_h)),
                    sphere(pos=( l_h,   0,-l_h)),
                    sphere(pos=(-l_h,   0, l_h)),
                    sphere(pos=(-l_h,   0,-l_h)),
                    sphere(pos=( l_h, l_h,   0)),
                    sphere(pos=( l_h,-l_h,   0)),
                    sphere(pos=(-l_h, l_h,   0)),
                    sphere(pos=(-l_h,-l_h,   0))
                    ]
                #setting
            for edge in edges:
                edge.color=color.red
                edge.opacity=0.5
                edge.radius=r_p

elif structure_type_chosen == "Zincblende":
    #center
        #position
    centers=[
            sphere(pos=( l/4, l/4, l/4)),
            sphere(pos=( l/4,-l/4,-l/4)),
            sphere(pos=(-l/4, l/4,-l/4)),
            sphere(pos=(-l/4,-l/4, l/4))
            ]
        #setting
    for center in centers:
        center.opacity=0.5
        center.radius=r_p
        if cheat_on:
            center.material=materials.earth
        else:
            center.color=color.blue

    #face
        #position
    faces=[
            sphere(pos=( l_h, l_h, l_h)),
            sphere(pos=( l_h, l_h,-l_h)),
            sphere(pos=( l_h,-l_h, l_h)),
            sphere(pos=(-l_h, l_h, l_h)),
            sphere(pos=(-l_h,-l_h, l_h)),
            sphere(pos=( l_h,-l_h,-l_h)),
            sphere(pos=(-l_h, l_h,-l_h)),
            sphere(pos=(-l_h,-l_h,-l_h))
            ]
        #setting
    for face in faces:
        face.color=color.yellow
        face.opacity=0.5
        face.radius=r_m
        if cheat_on:
            face.material=materials.wood

    #edge
        #position
    edges=[
            sphere(pos=( l_h,   0,   0)),
            sphere(pos=(   0, l_h,   0)),
            sphere(pos=(   0,   0, l_h)),
            sphere(pos=(-l_h,   0,   0)),
            sphere(pos=(   0,-l_h,   0)),
            sphere(pos=(   0,   0,-l_h)),
            ]
        #setting
    for edge in edges:
        edge.color=faces[0].color
        edge.optical=0.5
        edge.radius=r_m
        if cheat_on:
            edge.material=faces[0].material

#error messenge
else:
    print "System Error, call Heron now...!"
    raise SystemExit

print "done..."

注意：
windows底下，需安裝python之後安裝vpython方可在雙擊.py檔之後執行。
download

打招呼機器人

2010-03-06T01:16:00.000-08:00

這是我個人的第一台機器人，

完成於高二寒假。

他的唯一功能便是打招呼，

後來增加伸懶腰的進階功能。

詳細介紹

功能：

打招呼
伸懶腰

材料

Arduino 2009
servo motor
木條*1
人體紅外線感測模組
布偶*1
水銀開關
書架

運用技術

將撰寫完成的C程式燒入avr，插電後即可使用
水銀開關運用水銀液體的特性於玻璃容器中依循重力流動，並利用水銀的導電性開啟或是關閉開關 see more :

約10cm的木條鎖在servo motor上，並將布偶套上，如此可以模擬布偶敬禮的模樣
人體紅外線模組收人體紅外線的變化，上頭亦可以調整感應之後的延長時間

成品