python实现matlab中的filter函数

东方耀

python实现matlab中的filter函数

1. clc;clear all;close all;
   
2. 
   
3. 
   
4. h = [8, -3.5, 0.5, -1.5, 0.2517, 18, 19, 20];
   
5. a = [1];
   
6. x = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16];
   
7. 
   
8. y = filter(h,a,x);
   
9. 
   
10. % y的值：[8,12.5000000000000,17.5000000000000,21,24.7517000000000,46.5034000000000,87.2551000000000,148.006800000000,208.758500000000,269.510200000000,330.261900000000,391.013600000000,451.765300000000,512.517000000000,573.268700000000,634.020400000000]
    
11. 
    

复制代码

1. 
   
2. import numpy as np
   
3. from scipy.signal import lfilter
   
4. 
   
5. 
   
6. def filter_matlab(b, a, x):
   
7.     y = []
   
8.     y.append(b[0] * x[0])
   
9.     for i in range(1, len(x)):
   
10.         y.append(0)
    
11.         for j in range(len(b)):
    
12.             if i >= j:
    
13.                 y[i] = y[i] + b[j] * x[i - j]
    
14.                 j += 1
    
15.         for l in range(len(b) - 1):
    
16.             if i > l:
    
17.                 y[i] = (y[i] - a[l + 1] * y[i - l - 1])
    
18.                 l += 1
    
19.         i += 1
    
20.     return y
    
21. 
    
22. 
    
23. h = [8, -3.5, 0.5, -1.5, 0.2517, 18, 19, 20]
    
24. a = [1]
    
25. x = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]
    
26. # y = filter_matlab(h, a, x)    # 这个基本对 但是有bug 以后再细看
    
27. # https://stackoverflow.com/questions/8922657/matlab-filter-with-scipy-lfilter
    
28. y = lfilter(h, a, x)
    
29. 
    
30. # According to their documentation for Matlab filter() and SciPy lfilter()
    
31. 
    
32. print("python的结果(个数=%d)：" % len(y), y)
    
33. 
    
34. # python的结果(个数=16)： [  8.      12.5     17.5     21.      24.7517  46.5034  87.2551 148.0068
    
35. #  208.7585 269.5102 330.2619 391.0136 451.7653 512.517  573.2687 634.0204]

复制代码

Dummy

让天下人人学会人工智能！人工智能的前景一片大好！

Dummy

请问怎样才有198的阅读权限？

Dummy

请问怎样才有198的阅读权限？

Dummy

请问怎样才有198的阅读权限？

东方耀

Dummy 发表于 2022-4-4 16:08
请问怎样才有198的阅读权限？

联系我微信：dfy_88888