基于python實現(xiàn)matlab filter函數(shù)過程詳解
matlab中的filter函數(shù):
y = filter(b,a,x)
python實現(xiàn)matlab中的filter函數(shù)
def filter_matlab(b,a,x):
y = []
y.append(b[0] * x[0])
for i in range(1,len(x)):
y.append(0)
for j in range(len(b)):
if i >= j :
y[i] = y[i] + b[j] * x[i - j ]
j += 1
for l in range(len(b)-1 ):
if i >l:
y[i] = (y[i] - a[l+1] * y[i -l-1])
l += 1
i += 1
return y
example:
取
b = [8,-3.5,0.5]
a = [1,-1.5,0.25]
x = [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]
y = filter_matlab(b,a,x)
函數(shù)的結(jié)果和matlab的filter函數(shù)結(jié)果一致,為
[8, 24.5, 52.25, 94.75, 156.5625, 243.65625, 363.84375, 527.3515625, 747.56640625, 1042.01171875, 1433.6259765625, 1952.43603515625, 2637.74755859375, 3541.0123291015625, 4729.581604003906, 6291.619323730469, 8342.533584594727, 11033.395545959473, 14561.959922790527, 19187.090997695923]
以上就是本文的全部內(nèi)容,希望對大家的學習有所幫助,也希望大家多多支持腳本之家。
相關(guān)文章
python3新特性函數(shù)注釋Function Annotations用法分析
這篇文章主要介紹了python3新特性函數(shù)注釋Function Annotations用法,結(jié)合實例形式分析了Python3函數(shù)注釋的定義方法與使用技巧,需要的朋友可以參考下2016-07-07
對Python進行數(shù)據(jù)分析_關(guān)于Package的安裝問題
下面小編就為大家?guī)硪黄獙ython進行數(shù)據(jù)分析_關(guān)于Package的安裝問題。小編覺得挺不錯的,現(xiàn)在就分享給大家,也給大家做個參考。一起跟隨小編過來看看吧2017-05-05

