Question

matlab里的一段程序 求大神解释一下每句的含义。尤其是程序的前面5行，谢啦！！！

%频域中上抽样的效果
%使用fir2构造一个带限输入序列
clf
freq=[0 0.45 0.5 1]
mag=[0 1 0 0]
x=fir2(99,freq,mag)
%求取并画出输入谱
[Xz,w]=freqz(x,1,512,'whole')
subplot(2,1,1)
plot(w/pi,abs(Xz))
axis([0 1 0 1])
grid
xlabel('\ omega/ \pi')
ylabel('幅度')
title('输入谱')
subplot(2,1,2)

Answer 1

%频域中上抽样的效果
%使用fir2构造一个带限输入序列
clf %清除当前图形窗口图形。
freq=[0 0.45 0.5 1] %定义一个1*4的矩阵freq。
mag=[0 1 0 0] %定义一个1*4的矩阵mag。
x=fir2(99,freq,mag) %设置滤波器。

%fir2(n,f,m)函数fir2用来设计多通带任意响应FIR滤波器，该滤波器的幅频特性由向量对f和m确定，f为归一化频率向量，m为对应频率点上的幅度。当设计的滤波器在频率为π的幅度响应不是0时，滤波器的阶数n为偶数；
%上面是fir2的简介，来自百度百科。

%求取并画出输入谱
[Xz,w]=freqz(x,1,512,'whole') %求离散系统频响特性。

%关于freqz，可以看http://blog.sina.com.cn/s/blog_4c5154ea0100t41r.html

subplot(2,1,1) %把当前窗口分为两部分。在第一部分画图。
plot(w/pi,abs(Xz)) %画图。
axis([0 1 0 1]) %指定坐标轴范围。
grid %打开网格。
xlabel('\ omega/ \pi') %设定x轴说明。
ylabel('幅度') %设定x轴说明。
title('输入谱') %设定标题。
subplot(2,1,2) %在当前窗口第二个部分画图。

有问题再问。建议多多百度。

追问

(1)[H,w]=freqz(B,A,N)
(2)[H,w]=freqz(B,A,N,’whole’)
看了你给的那个网址对freqz的解释，但还是没看出来（1）和（2）的区别额。。。能说下吗

追答

[H,w]=freqz(B,A,N,’whole’)
意思是返回围绕整个单位圆N个采样点的频率响应。术语是双边频率特性？抱歉这方面的专业知识我不是太懂。
可以参看matlab官方的帮助文档：
http://cn.mathworks.com/help/signal/ref/freqz.html?searchHighlight=freqz
里面有丰富的例子。应该会对你有所帮助。

Answer 2

%频域中上抽样的效果
%使用fir2构造一个带限输入序列
clf  %清当前图形窗口
freq=[0 0.45 0.5 1];  %归一化频率
mag=[0 1 0 0];        %幅度
x=fir2(99,freq,mag);  %根据幅频特性(freq,mag确定)构造线性相位滤波器
%求取并画出输入谱
[Xz,w]=freqz(x,1,512,'whole') ; %求滤波器的频率特性
 %画出频率特性曲线
subplot(2,1,1);    
plot(w/pi,abs(Xz));     %归一化频率
axis([0 1 0 1])
grid
xlabel('\omega/ \pi');