基于python实现matlab filter函数过程详解
作者：chila 发布时间：2023-06-04 04:54:55
标签：python,matlab,filter,函数
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]
来源：https://www.cnblogs.com/chilalaa/p/13043716.html
0
投稿猜你喜欢
在默认情况下，大多数浏览器都会将有序列表中的数字序列的与其列表文字内容显示为相同的字体。这篇快速教程将教你如何使用有序列表(ol)和段落(p- 可能有些地方翻译得不好，请见谅在这个冠冕堂皇的标题之下，我想回答一个土方，有人在博客提出很久了。该土方很好的描述了很多图标设计师所面临的典型
- 之前在一个web系统的设计中，和另一个设计师讨论，“保存”和“取消”按钮该怎么设计。我的观点是，保存是比取消更常用的按钮，也是用户的主要目的
- 1、说明关键词传递以“形参变量名=实参”的形式参与实参关联，根据形参的名称进行参数传递，使实参和形参的顺序不一致。不用担心定义函数时参数的顺
- Python入门教程FROM:http://www.cnblogs.com/taowen/articles/11239.aspx作者：tao
- 1.php in_array方法说明PHP查找数组元素是否存在，一般会使用in_array方法。bool in_array ( mixed
- 今早打开 腾讯ISD的博客 ，看到一篇新的文章，《迷你屋视觉规范简介》，赶紧看了来学习。不过给我抓到问题咯，臭鱼不介意我在这说下吧：这套规范
- 前言在写程序时，我们会经常碰到程序出现异常，这时候我们就不得不处理这些异常，以保证程序的健壮性。处理异常的版本有以下几种，你通常的做法是哪种
- 自己从工艺品设计到平面设计到网络设计，虽然设计原则不离其宗，但经验下来的心得告诉自己，设计媒介的变化带来很多媒介自身的特殊性，下面总结下网站
- 我准备在ASP中连接MYSQL了，请问如何做？首先要正确安装MYSQLX，装好之后，可调用以下程序即可正常访问MYSQL：<％@&nb
- 作为设计师，我们都知道，一个极简的设计可以实现漂亮的效果。然而，很多设计师在实现上有些麻烦：要么是没有时间让使用如此少的元素制作的页面看起来
- 一、软件准备：以下均为截止2005-4-20的最新正式版本PHP（5.1.2）：http://www.php.netMySQL（5.0.19
- 本文实例为大家分享了Python曲线拟合的最小二乘法，供大家参考，具体内容如下模块导入import numpy as npimport ga
- numpy 中 的random模块有多个函数用于生成不同类型的随机数，常见的有 uniform、rand、random、randint、ra
- ImageField的使用笔记今天完善作业写的订单系统，主要是给每一个菜品增加图片，看起来美观一些，但是没想到这个小小的需求花了我一天时间，
- 如何制作一个WAP手机的WML网页？代码如下：＜％@Language=VBScriptMaxNoAds = 10'
- 看过数据库的备份与还原。大多数都是用组件来完成的。其实可通过sql语句来完成。 由于时间关系，未对参数进行验证和界面美化。代码
- 内容概要：print() 是一个常用函数。那么，您是否注意过，print() 会在显示当前语句后换行。如果遇到需要连续显示、不换行的情况，比
- 问：如何让内联框架里面的网页背景透明？（用iframe嵌套一网页，怎么能够设置其背景为透明以便能显示父框架中网页的背景图？）答：这是需要 I
- 译注：前两天看到一篇不错的英文文章，叫做 How browsers work，该文概要的介绍了浏览器从头到尾的工作机制，包括HTML等的解析
