1、三、编程题。(每题 30 分,共 60 分) 1、利用计算机分析信号x(t) sin20t 的频谱,试利用 Matlab 编写相关的程序,实 20t 西南大学网络与继续教育学院课程考试试题卷西南大学网络与继续教育学院课程考试试题卷 类别:网教专业:电气工程及其自动化 2018 年 6 月 课程名称【编号】 :数字信号处理【1077】 A 卷 大作业满分:100 分 一、简答题。 (20 分) 1、试画出输入倒位序、输出自然顺序,基-2 的快速傅里叶算法的基本的蝶形单元的 流图。 答: 现信号的抽样、频谱的计算、频谱的显示等,并对每条语句加以必要的注释。 答: t= -1:0.01:1;%以 0
2、.01 秒周期进行抽样,并加矩形窗截断,满足抽样定理 x=sinc(20*t);%得到点数为 201 的有限长序列 h=fft(x,1024);%做 1024 点的快速傅里叶变换,满足频域抽样定理 ff=1000*(0:511)/1024;%将数字频率转换为模拟频率,单位为Hz plot(ff,abs(h(1:512);%显示信号的幅度谱,由于对称性,只显示一半 - 1 - 二、计算题。 (20 分) 1、某频谱分析仪在谱分析时,一次抽样的点数为 2 的整数次幂,要求频率分辨力 10Hz,信号最高频率4KHz,试确定最小的记录长度,最小抽样频率及一个记录中最少点 数 N。 答: 2、设计一个线
3、性相位 FIR 低通滤波器, 已知抽样频率为 s 21.5103(rad /sec), 通带截止频率为 p 21.5102(rad /sec),阻带起始频率为 st 23102(rad /sec), 阻带衰减不小于-50dB。要求利用窗函数法进行设计,试给出设计过程和相应的 Matlab 程 序。 答: 程序如下: N=33;wc=.3; wn=hamming(N); n=0:N-1;a=(N-1)/2; hd=sin(wc*pi*(n-a)./(pi*(n-a); hd(a+1)=wc; hn=hd.*wn; h,w=freqz(hn,1); plot(w/pi,20*log10(abs(h); grid; xlabel(归一化频率) ylabel(幅度/dB) - 2 -