C语言实现matlab的butter函数

--------- 2019/04/18 更新 ---------

发现有网友反馈说不同的VS编译器的complex.h有些许不同, 故而自己写了一个用到复数的版本.

增加一个免complex.h的文件,  https://download.youkuaiyun.com/download/zone53/10813113

 

--------- 以下是原文 ---------

代码下载:  http://download.youkuaiyun.com/download/zone53/10029280

 

找了半天没在网上找到关于matlab中butter的C函数,不太清楚是别的开源库有还是怎么样,于是自己手翻了个,

 

只写了低通和带通部分,高通的带阻因为没用到,所以没写,有需要的也可以自己加下.不难.

 

由于对matlab不是特别了解,于是在翻的过程中,有部分matlab代码是省略掉了,但应该不影响结果,随机测试几组数据,结果都很好.

 

低通:

带通:

 

代码是在VS2015下面调的,因为程序内设计到了复数运算,故需包含头文件 #include "complex.h"

若没有这个头文件,那复数相关的部分,也可以自己写.

 

代码参考了 <徐士良常用算法程序集（C语言描述）第三版及程序>,此乃神书.

 

另,代码写的有点长,请到这里去下载; http://download.youkuaiyun.com/download/zone53/10029280

 

写的过程中卡在matlab的一个内置函数 eig,求矩阵的全部特征值,

矩阵的那部分早就还给老师了,但是手握神书而不知,最后得老师指点,直接用了神书内的代码.

但在查找相关资料的过程中,发现了C++的一些矩阵开源库.强推eigen

有关eigen库的一些基本使用方法:  http://blog.youkuaiyun.com/r1254/article/details/47418871    

Eigen:C++开源矩阵计算工具 -- Eigen的简单用法: http://blog.youkuaiyun.com/augusdi/article/details/12907341

C++矩阵库Eigen快速入门:  http://www.cnblogs.com/python27/p/EigenQuickRef.html

 

那就这样吧.

 

blog.youkuaiyun.com/zone53