matlab低通滤波器库函数代码_MATLAB低通滤波器的设计代码

单位冲激响应及其幅频响应及其代码

单位冲激响应的绘制思路比较简单，

就是将一条直线用

stem

函数绘制出来，

其图形必须为

23

个点。其代码如下：

n=0:22;

x=n./n;

stem(n,x); title('h(n)');

axis([0,25,0,1.3]);

单位冲激响应的幅频响应要用到专门的函数

m

文件。该函数文件可以在主

程序中调用多次，节省篇幅。函数

m

文件代码如下：

function xk=dft(xn,N)

n=[0:1:N-1];

k=n;

WN=exp(-j*2*pi/N);

nk=n'*k;

WNnk=WN.^nk;

xk=xn* WNnk;

不同窗函数低通幅频响应对比模块

窗函数不同，

会对设计出的低通滤波器的性能不同。

应尽量选取旁瓣小主瓣

窄的窗函数。

为了细致观察窗函数不同对设计的低通滤波器的不同影响，

本次试

验中设计了不同窗函数的低通滤波器的比较环节。

各个窗函数设计低通滤波器的

思路是相似的，只是其中的窗函数是不同的。

用矩形窗设计低通滤波器的代码：

D=str2num(get(handles.edit1,'string'));

wn=D*pi;

N=23;

b=fir1(N,wn,boxcar(N+1));