Question

matlab程序问题

M=21
wc=pi/4
n=0:M-1
r=(M-1)/2
nr=n-r+eps*((n-r)==0)
hdn=sin(wc*nr)/pi./nr;
if rem(M,2)~=0,hdn(r+1)=wc/pi;end
wn1=boxcar(M)
hn1=hdn.*wn1'
subplot(211)
stem(n,hn1,'.')
line([0,20],[0,0])
hw1=fft(hn1,512)
w1=2*[0:511]/512
subplot(212)
plot(w1,20*log10(abs(hw1)))

谁能对这段程序解释下，尤其是第五行的那句究竟是什么意思？
哦，这样啊，那么第七行呢？ 好像即使删除了，也不影响图像呀~~

为什么这句和没有这句，程序绘出的图像好像没啥区别啊？

line([0,20],[0,0])这条语句又是什么意思呢？ 还有我如何才能将一个信号通过这个滤波器呢？用什么函数？是filter吗？是的话怎么用呢？
我再加点分吧，拜托了啊~~~~~

Answer 1

eps 叫做机器的浮点运算误差限
nr=n-r+eps*((n-r)==0),这句括号里面（n-r)==0）表示判断n-r如果为0，结果就是1，如果n-r不为0，结果就是0；这句总的来说让n-r为0的元素加上一个特别小的数，以便hdn=sin(wc*nr)/pi./nr;因为分母不能为0；
这个程序就是一个Sa信号（sinc）画出样本图，和求其幅频特性。

rem函数是取余数函数rem(M,2)~=0表示M处以2看看余数是不是等于0，这个题目中看似不起作用，是因为hdn=sin(wc*nr)/pi./nr;这步刚好算出来的是hdn(11)=0.25,和那个if里面的刚好一样，所以不起作用。
谢谢刚才采纳我的了，希望会对你学习有帮助(⊙o⊙)

line([0,20],[0,0]); 这是画一条直线，就是在横轴是0-20，纵轴是0-0，就是说y=0；
你可以试试这个line([0,20],[0,1]);就明白了。
是的，我一般就是filter的，用法一般是y=filter(a,b,x),a和b是滤波器的系数矩阵，x通过滤波器得到y；具体的用法要看你用的什么样的滤波器了

Answer 2

第七行的意思是：如果 M/2的余数不等于0.就把wc/pi的值赋给hdn（r+1）；结束。其中rem（）是取余函数。

Answer 3