转载
https://blog.csdn.net/zone53/article/details/78280901
C语言实现matlab的butter函数
发现有网友反馈说不同的VS编译器的complex.h有些许不同, 故而自己写了一个用到复数的版本.
增加一个免complex.h的文件, https://download.csdn.net/download/zone53/10813113
--------- 以下是原文 ---------
代码下载: http://download.csdn.net/download/zone53/10029280
找了半天没在网上找到关于matlab中butter的C函数,不太清楚是别的开源库有还是怎么样,于是自己手翻了个,
只写了低通和带通部分,高通的带阻因为没用到,所以没写,有需要的也可以自己加下.不难.
由于对matlab不是特别了解,于是在翻的过程中,有部分matlab代码是省略掉了,但应该不影响结果,随机测试几组数据,结果都很好.
低通:
带通:
代码是在VS2015下面调的,因为程序内设计到了复数运算,故需包含头文件 #include “complex.h”
若没有这个头文件,那复数相关的部分,也可以自己写.
代码参考了 <徐士良常用算法程序集(C语言描述)第三版及程序>,此乃神书.
另,代码写的有点长,请到这里去下载; http://download.csdn.net/download/zone53/10029280
写的过程中卡在matlab的一个内置函数 eig,求矩阵的全部特征值,
矩阵的那部分早就还给老师了,但是手握神书而不知,最后得老师指点,直接用了神书内的代码.
但在查找相关资料的过程中,发现了C++的一些矩阵开源库.强推eigen
有关eigen库的一些基本使用方法: http://blog.csdn.net/r1254/article/details/47418871
Eigen:C++开源矩阵计算工具 – Eigen的简单用法: http://blog.csdn.net/augusdi/article/details/12907341
C++矩阵库Eigen快速入门: http://www.cnblogs.com/python27/p/EigenQuickRef.html
那就这样吧.