lempel ziv matlab,基于Python的LempelZiv算法的熵估計

此函數允許估計時間序列的熵。它基于Lempel-Ziv壓縮算法。對于長度為n的時間序列,熵估計為:

E=(1/n和L_i)^-1 ln(n)

式中,L逯i是從位置i開始的最短子串的長度,該子串之前沒有從位置1出現到i-1。當n接近無窮大時,估計的熵收斂到時間序列的實熵。在

我想用Python實現is,我是這樣做的:def contains(small, big):

for i in range(len(big)-len(small)+1):

if big[i:i+len(small)] == small:

return True

return False

def actual_entropy(l):

n = len(l)

sequence = [l[0]]

sum_gamma = 0

for i in range(1, n):

for j in range(i+1, n+1):

s = l[i:j]

if contains(s, sequence) != True:

sum_gamma += len(s)

sequence.append(l[i])

break

ae = 1 / (sum_gamma / n ) * math.log(n)

return ae

但是,我發現當數據越來越大時,它的計算速度太慢了。例如,我使用了一個包含23832個元素的列表作為輸入,所消耗的時間如下:(data can be foundhere)

^{pr2}$

我有成千上萬的這樣的名單要計算,這么長的時間是無法忍受的。我應該如何優化此功能并使其更快地工作?在

本文來自互聯網用戶投稿,該文觀點僅代表作者本人,不代表本站立場。本站僅提供信息存儲空間服務,不擁有所有權,不承擔相關法律責任。
如若轉載,請注明出處:http://www.pswp.cn/news/372600.shtml
繁體地址,請注明出處:http://hk.pswp.cn/news/372600.shtml
英文地址,請注明出處:http://en.pswp.cn/news/372600.shtml

如若內容造成侵權/違法違規/事實不符,請聯系多彩編程網進行投訴反饋email:809451989@qq.com,一經查實,立即刪除!

相關文章

Android使用繪圖Path總結

Path作為Android中一種相對復雜的繪圖方式,官方文檔中的有些解釋并不是很好理解,這里作一個相對全面一些的總結,供日后查看,也分享給大家,共同進步。 1.基本繪圖方法 addArc(RectF oval, float startAngle, float swee…

2017.3.23下午

下午通過對OSPF基本原理進一步的學習,對上午學習的內容進行了復習。 轉載于:https://www.cnblogs.com/bgd140206206/p/6606192.html

編寫Eclipse插件教程–第1部分

Eclipse是三個最受歡迎的Java開發IDE之一。 其成功的原因之一是其可擴展性。 對于任何知道該怎么做并且已經做到的人來說,編寫eclipse插件都可以非常輕松快捷。 不幸的是,第一次在Eclipse中進行操作可能會非常耗時且令人沮喪。 Eclipse框架非常龐大&…

簡單Window下 Android Studio的安裝

(1)首先安裝JDK 下載JDK 本人覺得官方網站下JDK比較慢,可以直接百度JDK,(如果是64位 百度搜索記得64位) 類似于這樣的下載 安裝可以看下教程,包括環境變量的配置 如何安裝JDK (2&…

日期處理一之NSLalendar的使用

一、日期和時間模式 日期和時間格式由日期和時間模式字符串組成,在日期和時間模式字符串中未加引號的A到‘Z’和a到‘z’被解釋為模式字母,用來表示日期或時間。字符串元素,文本可以使用單引號(‘’)引起來使用。定義以…

java的使用Pair要導入什么包,第三方jar包的使用

被導入的外部類所在源文件通常要打包成jar包,java中的jar文件裝的是 .class 文件。它是一種壓縮格式和zip兼容,被稱為jar包。JDK提供的許多類,也是以jar包的形式提供的。在用的時候呢,你的文件里有很多個類,把這些類和…

十大最受歡迎的新Eclipse插件

Eclipse Marketplace仍然是發現有趣且相關的Eclipse插件的地方。 通過Eclipse Marketplace客戶端,每月成功安裝100,000多個基于Eclipse的產品。 我們提供了過去30天 以來所有時間最受歡迎的插件列表。 我認為看看過去12個月中最受歡迎的新插件會很有趣。 以下列出了…

在桌面顯示我電腦

打開Windows PowerShell(一個像是命令提示符的東西[藍底白字],但不是命令提示符),在Windows PowerShell內輸入cmd回車,當返回如下信息: Microsoft Windows [版本 6.2.9200](c) 2012 Microsoft Corporation。…

《Java技術》第二次作業計科1501趙健宇

(一)學習總結 1.使用Eclipse關聯jdk源代碼,查看String類的equals()方法 equals()方法截圖 “”比較的是地址。equals方法他同樣使用號進行內存地址的比較。但是equals方法重寫如果號比較不相等,…

注射php,UPDATE注射(mysqlphp)的兩個模式

一.測試環境:OS:Windowsxpsp2php:php4.3.10(mysql4.1.9apache1.3.33二.測試數據庫結構:-----start-----數據庫:test----------------------------------------------------------------表的結構userinfo--CREATETABLEuserinfo(groudidvarchar(12)NOTNULL…

JUnit的內置Hamcrest Core Matcher支持

在用JUnit和Hamcrest改進assertEquals的文章中,我簡要討論了Hamcrest “ 核心 ”匹配器與JUnit的現代版本“結合”在一起的情況。 在那篇文章中,我特別關注了JUnit的assertThat(T,Matcher)靜態方法與Hamcrest核心is()匹…

FormsAuthentication使用指南,實現登錄

一般情況下,在我們做訪問權限管理的時候,會把用戶的正確登錄后的基本信息保存在Session中,以后用戶每次請求頁面或接口數據的時候,拿到Session中存儲的用戶基本信息,查看比較他有沒有登錄和能否訪問當前頁面。Session的…

Java開發高性能網站需要關注的事

轉自:http://www.javabloger.com/java-development-concern-those-things/ 近期各家IT媒體舉辦的業內技術大會讓很多網站都在披露自己的技術內幕與同行們分享,大到facebook,百度,小到剛起步的網站。facebook,百度之類的…

mysql到oracle數據遷移,mysql數據遷移到oracle

mysql: utf-8字符集oracle: Al32UTF8步驟:1. 從mysql導出csv格式數據,導出時設置default_character_setgbk,否則會亂碼2.上傳至linux服務器(oracle服務器),用以下命令將windows的/r/n換成/r:sed -i "s//r//g" file3.將gbk編碼轉為u…

測試抽象類和模板方法模式

摘自Wikipedia,“模板方法定義了算法的程序框架。 子類可以覆蓋一個或多個算法步驟,以允許不同的行為,同時確保仍然遵循總體算法。” 通常,此模式由兩個或多個類組成,一個是提供模板方法(非抽象&#xff09…

LoadRunner性能測試結果計數器指標說明

LoadRunner性能測試結果計數器指標說明 轉載2015-09-23 09:57:13標簽:loadrunner計數器針對性能測試結果分析過程中,面對大量的測試數據,反而感覺無從下手分析。今天我們就Windows操作系統計數器中的個別被監控對象進行簡單的說明。 Memory: …

Elasticsearch之需要注意的問題(es和jdk版本)

(1)在使用java代碼操作es集群的時候 要保證本地使用的es的版本和集群上es的版本保持一致。 (2)保證集群中每個節點的JDK版本和es基本配置一致 這個很簡單,不多說。 (3)es集群中jdk版本不一致 …

oracle視圖查詢數據慢,8i查詢DBA_FREE_SPACE視圖極慢的問題

還是那套古老的8.1.7.4,在該系統上檢查表空間使用情況的SQL運行緩慢,其SQL如下:SELECT D.TABLESPACE_NAME,SPACE "SUM_SPACE(M)",SPACE - NVL(FREE_SPACE, 0) "USED_SPACE(M)",ROUND((1 - NVL(FREE_SPACE, 0) / SPACE) * 100, 2) &q…

bzoj 4332:JSOI2012 分零食

描述 這里是歡樂的進香河,這里是歡樂的幼兒園。 今天是2月14日,星期二。在這個特殊的日子里,老師帶著同學們歡樂地跳著,笑著。校長從幼兒園旁邊的小吃店買了大量的零食決定分給同學們。聽到這個消息,所有同學都安安靜靜…

并發–順序線程和原始線程

我不久前參與了一個項目,該項目的報告流程如下: 用戶會要求舉報 報告要求將被翻譯成較小的部分 每個零件的報告將基于零件/節的類型由報告生成器生成 組成報告的各個部分將重新組合成最終報告,并返回給用戶 我的目標是展示如何從錯誤的實…