基于python實(shí)現(xiàn)matlab filter函數(shù)過程詳解
matlab中的filter函數(shù):
y = filter(b,a,x)
python實(shí)現(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)容,希望對大家的學(xué)習(xí)有所幫助,也希望大家多多支持好吧啦網(wǎng)。
相關(guān)文章:
1. Python+unittest+requests 接口自動化測試框架搭建教程2. Python的文本常量與字符串模板之string庫3. 利用CSS制作3D動畫4. 存儲于xml中需要的HTML轉(zhuǎn)義代碼5. 完美解決vue 中多個echarts圖表自適應(yīng)的問題6. jsp+servlet簡單實(shí)現(xiàn)上傳文件功能(保存目錄改進(jìn))7. 一款功能強(qiáng)大的markdown編輯器tui.editor使用示例詳解8. .Net加密神器Eazfuscator.NET?2023.2?最新版使用教程9. Java GZip 基于內(nèi)存實(shí)現(xiàn)壓縮和解壓的方法10. SpringBoot+TestNG單元測試的實(shí)現(xiàn)
