明霞山资源网 Design By www.htccd.com

matlab中的filter函数:

y = filter(b,a,x)

python实现matlab中的filter函数

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)

函数的结果和matlab的filter函数结果一致,为

[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]

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。

标签:
python,matlab,filter函数

明霞山资源网 Design By www.htccd.com
广告合作:本站广告合作请联系QQ:858582 申请时备注:广告合作(否则不回)
免责声明:本站文章均来自网站采集或用户投稿,网站不提供任何软件下载或自行开发的软件! 如有用户或公司发现本站内容信息存在侵权行为,请邮件告知! 858582#qq.com
明霞山资源网 Design By www.htccd.com

稳了!魔兽国服回归的3条重磅消息!官宣时间再确认!

昨天有一位朋友在大神群里分享,自己亚服账号被封号之后居然弹出了国服的封号信息对话框。

这里面让他访问的是一个国服的战网网址,com.cn和后面的zh都非常明白地表明这就是国服战网。

而他在复制这个网址并且进行登录之后,确实是网易的网址,也就是我们熟悉的停服之后国服发布的暴雪游戏产品运营到期开放退款的说明。这是一件比较奇怪的事情,因为以前都没有出现这样的情况,现在突然提示跳转到国服战网的网址,是不是说明了简体中文客户端已经开始进行更新了呢?