作业帮 > 综合 > 作业

matlab里的一段程序 尤其是程序的前面5行,

来源:学生作业帮 编辑:拍题作业网作业帮 分类:综合作业 时间:2024/04/28 02:46:47
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)
%频域中上抽样的效果
%使用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') %求离散系统频响特性.
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
里面有丰富的例子。应该会对你有所帮助。