亚洲免费在线-亚洲免费在线播放-亚洲免费在线观看-亚洲免费在线观看视频-亚洲免费在线看-亚洲免费在线视频

Python求離散序列導數的示例

系統 2572 0

有一組4096長度的數據,需要找到一階導數從正到負的點,和三階導數從負到正的點,截取了一小段。

394.0
388.0
389.0
388.0
388.0
392.0
393.0
395.0
395.0
394.0
394.0
390.0
392.0

按照之前所了解的,對離散值求導其實就是求差分,例如第i點的導數(差分)為:

即在一個寬度為2m+1的窗口內通過計算前后m個值加權后的和得到。但是在實際使用過程中效果不是很好。于是想到了同樣在一個寬度為2k+1的窗口內,將這2k+1個點擬合成一個函數,然后求導就可以得到任意階數的導數值。

首先是函數擬合,使用from scipy.optimize import leastsq即最小二乘擬合

            
from scipy.optimize import leastsq
class search(object):
  def __init__(self, filename):
    self.filename = filename

  def func(self, x, p):
    f = np.poly1d(p)
    return f(x)

  def residuals(self, p, x, y, reg):
    regularization = 0.1 # 正則化系數lambda
    ret = y - self.func(x, p)
    if reg == 1:
      ret = np.append(ret, np.sqrt(regularization) * p)
    return ret

  def LeastSquare(self, data, k=100, order=4, reg=1, show=1): # k為求導窗口寬度,order為多項式階數,reg為是否正則化
    l = self.len
    step = 2 * k + 1
    p = [1] * order
    for i in range(0, l, step):
      if i + step < l:
        y = data[i:i + step]
        x = np.arange(i, i + step)
      else:
        y = data[i:]
        x = np.arange(i, l)
      try: 
        r = leastsq(self.residuals, p, args=(x, y, reg))
      except:
        print("Error - curve_fit failed")
      fun = np.poly1d(r[0]) # 返回擬合方程系數
      df_1 = np.poly1d.deriv(fun) # 求得導函數
      df_2 = np.poly1d.deriv(df_1)
      df_3 = np.poly1d.deriv(df_2)
      df_value = df_1(x)
      df3_value = df_3(x)
          

fun = np.poly1d(r[0]),fun返回的是一個 polynomial class,具體使用可以見官方文檔numpy.poly1d
polynomial對象可以使用deriv方法求導數,求得的依然是 polynomial對象。 df_value = df_1(x)所得到的就是x這個幾個點求得的導數值。

看似大功告成,但是求導的結果并不是很好,如下圖,實際最高點在100左右,但是擬合出來的曲線最高點在120左右,而原因在于使用多項式擬合很難準確擬合曲線。

Python求離散序列導數的示例_第1張圖片

于是想用高斯函數來實現對曲線的擬合,在matlab中試了下,三階高斯擬合可以很好的擬合曲線,

Python求離散序列導數的示例_第2張圖片

但是numpy以及sicpy中沒有找到類似poly1d這種對象,雖然可以自己定義高斯函數,如下

            
  def gaussian(self, x, *param):
    fun = param[0]*np.exp(-np.power(x - param[2], 2.) / (2 * np.power(param[4],    2.)))+param[1]*np.exp(-np.power(x - param[3], 2.) / (2 * np.power(param[5], 2.)))
    return fun
          

但是,在通過最小二乘擬合得到函數參數后只能得到擬合后的點,無法直接求導數..所以并不適合。

所以還是只能回到多項式擬合,如果4階多項式不能表征的話,更高階的呢

Python求離散序列導數的示例_第3張圖片

總體來說,效果還是可以接受的。

如果下階段找到好的高斯函數擬合方法,會繼續更新。

以上這篇Python求離散序列導數的示例就是小編分享給大家的全部內容了,希望能給大家一個參考,也希望大家多多支持腳本之家。


更多文章、技術交流、商務合作、聯系博主

微信掃碼或搜索:z360901061

微信掃一掃加我為好友

QQ號聯系: 360901061

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。

【本文對您有幫助就好】

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描上面二維碼支持博主2元、5元、10元、自定義金額等您想捐的金額吧,站長會非常 感謝您的哦!!!

發表我的評論
最新評論 總共0條評論
主站蜘蛛池模板: 兽皇在线观看 | 国产精品视频全国免费观看 | 国产 福利 在线 | 久久精品综合一区二区三区 | 毛片在线免费视频 | 欧美另类黑人巨大videos | 婷综合| 国产成人精品视频频 | 在线欧美亚洲 | 尤物视频在线 | 欧美乱大交xxxxxx喷潮免费 | 青青爽国产手机在线观看免费 | 久久国产精品广西柳州门 | 二级毛片免费观看全程 | 精品国产品国语在线不卡丶 | 久久精品资源站 | 日韩免费高清 | 欧美与黑人午夜性猛交久久久 | 精品日韩在线观看 | 日韩精品国产自在久久现线拍 | 色偷偷在线刺激免费视频 | 九九热精品视频在线播放 | 久久久久久久久久免免费精品 | 青娱乐国产在线视频 | 久久这里只有精品18 | 国产午夜亚洲精品国产 | 中国精品久久精品三级 | 国产三级不卡 | 久热re在线视频精品免费 | 黄视频在线观看www 黄视频在线观看网站 | 久久久久依人综合影院 | 狠狠大日本亚洲香蕉亚洲 | 日日碰日日摸日日澡视频播放 | 亚洲精品国产字幕久久vr | 日本三级一区 | 久久综合九色综合亚洲小说 | 欧美日韩三区 | 老司机午夜永久在线观看 | 精品哟哟哟国产在线观看不卡 | 福利视频在线播放 | 国产一级淫片a视频免费观看 |