因为教研室项目要求，近期做了关于巴特沃斯滤波器部分，采用的是C#与Matlab混合编程的方式，由于是第一次写博客，还有许多不足的地方。

教研室用的VS版本为2017版，而MatlabR2012a和MatlabR2014b似乎并不支持VS2017版，经过几番折腾，确定Matlab为2016版。
Matlab2016a安装步骤及破解详见以下地址：
[http://jingyan.baidu.com/article/870c6fc300c2fab03ee4be70.html]
安装完成后，在安装目录下的toolbox\compiler\deploy\win64找到MCRInstaller点击安装，完成MCR的安装

1、安装及破解工作完成后打开Matlab软件，界面如下：

2、点击新建—>函数

function Y = b_filter(X,W1,W2,W3,W4,M)
    if M==0   %低通滤波器
        [n,Wn] = buttord(W1,W2,3,60);
        [b,a] = butter(n,Wn,‘low‘);
    elseif M==1 %高通滤波器
        [n,Wn] = buttord(W1,W2,3,60);
        [b,a] = butter(n,Wn,‘high‘);
    elseif M==2 %带通滤波器
        [n,Wn] = buttord([W1 W2],[W3 W4],3,60);
        [b,a] = butter(n,Wn);
    elseif M==3 %带阻
        [n,Wn] = buttord([W1 W2],[W3 W4],3,60);
        [b,a] = butter(n,Wn,‘stop‘);
    else%默认为低通滤波器
        [n,Wn] = buttord(W1,W2,3,60);
        [b,a] = butter(n,Wn,‘low‘);
    end
Y = filter(b,a,X);
end12345678910111213141516171819

在本次巴特沃斯低通滤波器的设计中，将通带最大衰减R_p取值为3dB，阻带最小衰减A_s取值为60dB，使用Matlab工具中的自带函数buttord得到阶次n和截止频率W_n，即：
[n, W_n] = buttord(W1,W2,3,60)，其中W1和W2为C#部分中用户输入的截止频率w1和w2经过计算得到的数字频率；再用自带函数butter得到滤波器系数，即[b,a] = butter(n, W_n,’low’)，其中b和a为分子和分母的系数；最后使用Matlab工具中的自带函数filter完成滤波，即filter（b,a,X),其中X为输入的原始数据。

3、接下来就是打包生成.dll文件，在命令窗口输入deploytool，如图：

点击回车，出现以下界面：

选择第三个，Library Complier