Python 的多項式曲線擬合 

import numpy as np
import pylab as plt

# 隨機產生資料點來Fit
x = np.random.rand(7)
y = np.random.rand(7)

# 最高次項為4次方的多項式
line4 = polyfit(x, y, 4) 
pn4 = poly1d(line4) 

line3 = polyfit(x, y, 3) # 最高次項為3次方的多項式
pn3 = poly1d(line3)

xaxis = linspace(0, 1, 100)
plt.plot(x, y, 'o', xaxis, pn4(xaxis),'-g', xaxis, pn3(xaxis),'-b')
plt.legend(['data to fit', '4th degree poly', '3th degree poly'])
plt.axis([0,1,0,1])
plt.show()
再提高次方數可以擬合的更好，點到即止。

留言

張貼留言

這個網誌中的熱門文章

VBA EXCEL 工作表變化 馬上執行 的作法 Worksheet_Change

圖片
沒想到還會用到 Excel 中的 VBA，當初寫VBA就是一段窮極無聊的日子，用Excel寫星象表… 現在居然要幫一位老大哥寫這樣的程式…當真所有的事，當下是很難知道它發生的意義是什麼… 癈話不多說 先開啟EXCEL，存成副檔名為xlsm，這樣可以使用巨集功能。 然後看一上面，標籤列，有無開發人員選項。 如果沒有，從 檔案->選項->自定功能區->開發人員(打勾) 如果有，點擊 visual basic，就可以開始寫VBA 我們的目的是這樣的，如果 worksheet中有任何的更動，都會執行一段程式。 所以我們先點擊左邊的工作表1，然後右邊的程式編輯欄中上面的選擇列中，選worksheet。要記得這個設定是 對應到個別的worksheet ，所以如果有其它工作表，要個別寫對應的程式。如下 接下來選右邊的change 選完就會自動產生兩個程式名。 Worksheet_Change與 Worksheet_SelectionChange，如下。Worksheet_Change是指任何在worksheet 上的變動，就會執行這個程式。Worksheet_SelectionChange是指，選定EXCEL中的格有變化，如本來選A1改點擊了B1，這個動作就會觸發Worksheet_SelectionChange。 我們在兩個程式中寫一個MsgBox, 如下 各位可以試試效果如何。只要點擊左上方，excel的小圖，就會回到excel的worksheet中。 我們在 Worksheet_Change中增加三行程式，會在A1, A2, A3格中填入a, b, c。 各位是否有推論，這樣的寫法，會產生無窮迴圈？ 這個很重要！在寫 Worksheet_Change程式時，要很小心，如果是 工作表1自己的 Worksheet_Change要對工作表1自己做改變，就會進入無窮迴圈。 所以我們要能自由的開關觸發功能。用   Application.EnableEvents  我們加入  Application.EnableEvents ，再試試看。...
閱讀完整內容

python 找圖自動點擊

需要 import 這幾個 module import cv2 import numpy as np import pyautogui import time 滑鼠點擊 def click_image(image,pos,  action, timestamp,offset=5):     img = cv2.imread(image)     height, width, channels = img.shape     pyautogui.moveTo(pos[0] + offset, pos[1] + ,offset, timestamp)     pyautogui.click(button=action) 圖片尋找 def imagesearch(image, precision=0.8):     im = pyautogui.screenshot()     img_rgb = np.array(im)     img_gray = cv2.cvtColor(img_rgb, cv2.COLOR_BGR2GRAY)     template = cv2.imread(image, 0)     template.shape[::-1]     res = cv2.matchTemplate(img_gray, template, cv2.TM_CCOEFF_NORMED)     min_val, max_val, min_loc, max_loc = cv2.minMaxLoc(res)     if max_val < precision:         return [-1,-1]     return max_loc #返回圖片座標 使用方式 if __name__ == '__main__':     pp =im...
閱讀完整內容

Python pyserial 抓取系統內的 COM PORT

讀取COM1的資料，再把讀到的資料寫入一個 the_com1.txt的檔案中。 import serial def connectCom(): ser = serial.Serial("COM1", 19200) print ser.port while 1: f = open( 'the_com1.txt','a') sline = ser.readline() time1=datetime.datetime.now().strftime('%Y-%m-%d %H:%M:%S') out10=time1,"COM1",sline.strip() f.writelines(str(out10)) f.writelines("\n") f.close() 如果我們不知道到底有那一些COM PORT可以讀呢？可以藉由以下的片段來得知目前有哪一些COM PORT存在。 import serial def serial_ports(): p = ['COM%s' % (i + 1) for i in range(256)] rlt = [] for port in p: try: s = serial.Serial(port) s.close() rlt.append(port) except (OSError, serial.SerialException): pass return rlt
閱讀完整內容