下载文档原格式
电 子 科 技 大 学实 验 报 告学生姓名：Shrimp 学 号： 指导教师：一、实验室名称：数字信号处理实验室 二、实验项目名称：FFT 的实现 三、实验原理：一．FFT 算法思想：1．DFT 的定义：对于有限长离散数字信号{x[n]}，0 ≤ n ≤ N-1，其离散谱{x[k]}可以由离散付氏变换（DFT ）求得。
DFT 的定义为：21[][]N jnk Nn X k x n eπ--==∑，k=0,1,…N-1通常令2jNN eW π-=，称为旋转因子。
2．直接计算DFT 的问题及FFT 的基本思想：由DFT 的定义可以看出，在x[n]为复数序列的情况下，完全直接运算N 点DFT 需要（N-1）2次复数乘法和N （N-1）次加法。
因此，对于一些相当大的N 值（如1024）来说，直接计算它的DFT 所作的计算量是很大的。
FFT 的基本思想在于，将原有的N 点序列分成两个较短的序列，这些序列的DFT 可以很简单的组合起来得到原序列的DFT 。
例如，若N 为偶数，将原有的N 点序列分成两个（N/2）点序列，那么计算N 点DFT 将只需要约[(N/2)2 ·2]=N 2/2次复数乘法。
即比直接计算少作一半乘法。
因子（N/2）2表示直接计算（N/2）点DFT 所需要的乘法次数，而乘数2代表必须完成两个DFT 。
上述处理方法可以反复使用，即（N/2）点的DFT 计算也可以化成两个（N/4）点的DFT （假定N/2为偶数），从而又少作一半的乘法。
这样一级一级的划分下去一直到最后就划分成两点的FFT 运算的情况。
3．基2按时间抽取（DIT ）的FFT 算法思想：设序列长度为2L N =，L 为整数（如果序列长度不满足此条件，通过在后面补零让其满足）。
将长度为2L N =的序列[](0,1,...,1)x n n N =-，先按n 的奇偶分成两组：12[2][][21][]x r x r x r x r =+=，r=0,1,…,N/2-1DFT 化为：1/21/212(21)0/21/21221200/21/211/22/2[]{[]}[][2][21][][][][]N N N nk rk r kNNNn r r N N rk k rk NNN r r N N rk k rk N NN r r X k DFT x n x n Wx r Wx r W x r W Wx r W x r WWx r W ---+===--==--=====++=+=+∑∑∑∑∑∑∑上式中利用了旋转因子的可约性，即：2/2rkrkNN W W =。
一．选择与填空：1．若一系统为因果系统，则该系统的单位脉冲响应h(n)应满足0,0)(<≡n n h ，该系统的系统函数应满足 其收敛域为∞≤<-z R x ；若该系统为稳定系统，则该系统的单位脉冲响应h(n)应满足∞<∑∞-∞=n n h )(1=。
2．对下列序列，试求其零极点，Z 变换和收敛域。
a) )()(10n R n x = b) )(5sin )(n u n n x ⎪⎭⎫⎝⎛=π解:a) )1(111)(910110--=--=--z z z z z z X 其收敛域为∞≤<z 0,其零极点分布为: 零点:9个一阶零点9,,1 102 ==k ez k j k π;极点:一个9阶极点;0=zb) ⎪⎪⎭⎫ ⎝⎛-==-∞=--∞=∞=-∑∑∑n n n j n n n j n nz e z e j z nz X 0505021)5sin()(πππ ))(()2.0sin()2.0cos(21)2.0sin(111121552111515πππππππj j j j e z e z z zz z z e z e j ---------=+-=⎪⎪⎪⎭⎫⎝⎛---= 其收敛域为∞≤<z 1,其零极点分布为: 零点:1个一阶零点0=z ;极点:2个一阶极点;2.01πj e z =;2.02πj e z -=3．若对话音信号进行离散时间处理，已知话音信号为实信号，其频率范围为300~3400Hz ，则抽样频率为 >6800Hz 。
4．对连续时间非周期信号的谱分析工具是 傅立叶变换 ，对连续时间周期信号的谱分析工具是 傅立叶级数 ，对周期序列的谱分析工具是 离散傅立叶级数 ，对非周期序列的谱分析工具是 离散傅立叶变换 。
5．若某系统的群延时为常数，则该系统具有的相位特性为 线性相位 。
6．脉冲响应不变法存在频谱混叠的现象，是因为 z s →映射为多对一映射 ；双线性变换法不存在频谱混叠的现象，是因为 z s →映射为一对一映射 ；因此，对于同样指标设计的滤波器，双线性变换法的设计结果 优于 （优于/劣于）脉冲响应不变法的设计结果。
第二次编程作业加载信号(EEGdata.txt)；分析信号的幅度谱,确定低通FIR数字滤波器的指标；分别利用各种窗函数（Rectangular, Hamming, Kaiser）设计此FIR滤波器；对信号加随机噪声，并用设计的FIR滤波器对含噪声信号进行滤波。
要求：1）给出程序，画出滤波器幅度谱及其增益图(dB)；分析设计的滤波器是否达到指标。
2）利用设计的三种滤波器对加载的信号分别进行滤波，对比分析滤波结果。
程序：data = load('EEGdata.txt');n=1:length(data);xn=0.1*rand(length(data),1)-0.05;x=xn+data;[h0,w0]=freqz(data);[h1,w1]=freqz(x);figure(1)plot(w0/pi,abs(h0),w1/pi,abs(h1),'r');legend('Ô-dataÐźÅƵÆ×ͼ','¼ÓÔëÐźÅƵÆ×ͼ');wp = 0.15*pi;ws = 0.2*pi;DB = ws-wp;N = ceil(2*3.11*pi/DB+1);wc = (wp+ws)/2;wc = wc/pi;win = rectwin(N+1);b = fir1(N, wc, 'low', win);figure(2)freqz(b,1);wp = 0.15*pi;ws = 0.2*pi;DB = ws-wp;N = ceil(2*3.11*pi/DB+1);wc = (wp+ws)/2;wc = wc/pi;win = hamming(N+1);b1 = fir1(N, wc, 'low', win);figure(3)freqz(b1,1);fpts = [0.15, 0.2];mag = [1, 0];dev = [0.01, 0.01];[N, Wn, beta, ftype] = kaiserord(fpts, mag, dev)Kwin = kaiser(N+1, beta);b2 = fir1(N, Wn, Kwin);[h,omega] = freqz(b2,1,512);figure(4)freqz(b2,1);y=filter(b, 1, x);figure(5)subplot(2,1,1);plot(n,data,n,x,'g',n,y,'r');title('¾ØÐδ°Â˲¨·ù¶ÈÆ×Ч¹û¶Ô±Èͼ');lege nd('Ô-ÐźÅdata','¼ÓÔëÐźÅx','Â˲¨ºóÐźÅy');[H,W]=freqz(y);subplot(2,1,2);plot(w0/pi,abs(h0),W/pi,abs(H),'r');title('¾ØÐδ°Â˲¨ÆµÆ×Ч¹û¶Ô±Èͼ') ;legend('Ô-ÐźÅdata','Â˲¨ºóÐźÅy');y1=filter(b1, 1, x);figure(6)subplot(2,1,1);plot(n,data,n,x,'g',n,y1,'r');title('Hamming´°Â˲¨·ù¶ÈÆ×Ч¹û¶Ô±Èͼ'); legend('Ô-ÐźÅdata','¼ÓÔëÐźÅx','Â˲¨ºóÐźÅy1');[H,W]=freqz(y1);subplot(2,1,2);plot(w0/pi,abs(h0),W/pi,abs(H),'r');title('Hamming´°Â˲¨ÆµÆ×Ч¹û¶Ô±Èͼ');legend('Ô-ÐźÅdata','Â˲¨ºóÐźÅy1');y2=filtfilt(b2, 1, x); figure(7) subplot(2,1,1);plot(n,data,n,x,'g',n,y2,'r');title('Kaiser´°Â˲¨·ù¶ÈÆ×Ч¹û¶Ô±Èͼ');l egend('Ô-ÐźÅdata','¼ÓÔëÐźÅx','Â˲¨ºóÐźÅy2'); [H,W]=freqz(y2); subplot(2,1,2);plot(w0/pi,abs(h0),W/pi,abs(H),'r');title('Kaiser´°Â˲¨ÆµÆ×Ч¹û¶Ô±Èͼ');legend('Ô-ÐźÅdata','Â˲¨ºóÐźÅy2');figure(8)plot(n,data,n,y,'c',n,y1,'r',n,y2,'g');title('ÈýÖÖ´°º¯ÊýÂ˲¨Ð§¹û¶Ô±Èͼ');legend('Ô-dataÐźÅ','¾ØÐδ°Â˲¨','Hamming´°Â˲¨','Kaiser´°Â˲¨');0.10.20.30.40.50.60.70.80.9102468101214161800.10.20.30.40.50.60.70.80.91-3000-2000-1000Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )00.10.20.30.40.50.60.70.80.91-100-5050Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )00.10.20.30.40.50.60.70.80.91-3000-2000-1000Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )0.10.20.30.40.50.60.70.80.91-150-100-50050Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )0.10.20.30.40.50.60.70.80.91-2000-1500-1000-500Normalized Frequency (⨯π rad/sample)P h a s e (d e g r e e s )0.10.20.30.40.50.60.70.80.91-150-100-50050Normalized Frequency (⨯π rad/sample)M a g n i t u d e (d B )-0.2-0.100.10.2矩形窗滤波幅度谱效果对比图0.10.20.30.40.50.60.70.80.91051015矩形窗滤波频谱效果对比图-0.2-0.100.10.2Hamming 窗滤波幅度谱效果对比图0.10.20.30.40.50.60.70.80.91051015Hamming 窗滤波频谱效果对比图-0.2-0.100.10.2Kaiser 窗滤波幅度谱效果对比图0.10.20.30.40.50.60.70.80.9105101520Kaiser 窗滤波频谱效果对比图结果分析：滤波效果皆不尽如人意，原因是加载的噪声信号是随机信号，各频率皆有，滤波器通带部分的噪声无法滤掉。
１自测题一及答案一．选择与填空：1．数字信号是______(3)___。
（1）时间上和幅度上均连续的理想信号 （2）时间上离散而幅度上连续的信号（3）时间上离散、幅度上量化的非理想信号 （4）时间上连续而幅度上量化的非理想信号2．对于时间上离散的周期信号，可以___(4)___________。
（1）用傅立叶变换（FT ）求其频谱 （2）用离散傅立叶变换（DFT ）求其频谱 （3）可展成傅立叶级数（FS ） （4）可展成离散傅立叶级数（DFS ）3．已知x(n)为16点序列，其DFT 应为___(5)___，IDFT 应为____(6)___。
（1） ∑==15016)(161)(n nK W n x K X （2） nKn W n x K X 1616)()(∑==（3） ∑==1516)()(n nKWn x K X （4） ∑==1515)()(n nKW n x K X（5） ∑==1516)()(K nK WK X n x （6） ∑=-=15016)(161)(K nKW K X n x （7） ∑==15016)(161)(K nK W K X n x （8） ∑=-=15015)(161)(K nKW K X n x 4．数字信号处理中，时域上表征系统的瞬时输出经常用_____(2)_____。
（1）系统的单位脉冲响应 h(n) （2）离散卷积和（3）系统的传递函数H(Z) （4）系统的频响H(e j ω)5．如果对数字滤波器的相频特性不作要求，仅从滤波效率来看，IIR 数字滤波器的效率比FIR 数字滤波 器的效率要___高___，即在相同的通带、阻带滤波性能指标下，__IIR_____数字滤波器结构要简单些， 其滤波器的阶数N 要__低_____。
6．已知某线性相位FIR 数字滤波器的其中一个零点40πjeZ =，根据零点分布规律，其系统函数：H （Z ）=)1)(1(44ππjjeeA ---。
1.下面关于IIR滤波器设计说法正确的是（）。
A.双线性变换法的优点是数字频率和模拟频率成线性关系B.冲激响应不变法无频率混叠现象C.冲激响应不变法不适合设计高通滤波器D.双线性变换法只适合设计低通、带通滤波器【参考答案】: C2.题目及选项如下图所示A.AB.BC.CD.D【参考答案】: A3.下列关于窗函数设计法的说法中错误的是（）。
A.窗函数的长度增加，则主瓣宽度减小，旁瓣宽度减小B.窗函数的旁瓣相对幅度取决于窗函数的形状，与窗函数的长度无关C.为减小旁瓣相对幅度而改变窗函数的形状，通常主瓣的宽度会增加D.对于长度固定的窗，只要选择合适的窗函数就可以使主瓣宽度足够窄，旁瓣幅度足够小【参考答案】: B4.一个线性移不变系统稳定的充分必要条件是其系统函数的收敛域包括（）。
A.原点B.单位圆C.虚轴D.实轴【参考答案】: B5.题目及选项如下图所示A.AB.BC.CD.D【参考答案】: A6.已知一连续时间信号的最高截止频率是4kHz，则要从该连续时间信号的均匀离散采样值无失真地恢复原信号，则采样频率为（）。
A.2kHzB.4kHzC.6kHzD.8kHz【参考答案】: D7.题目及选项如下图所示A.AB.BC.CD.D【参考答案】: A8.题目及选项如下图所示A.AB.BC.CD.D【参考答案】: A9.对连续时间周期信号的谱分析工具是（）。
A.傅里叶变换B.傅里叶级数C.离散傅里叶变换D.离散傅里叶级数【参考答案】: B10.下面说法中正确的是（）。
A.连续时间非周期信号的频谱为连续周期函数B.连续时间周期信号的频谱为连续周期函数C.离散时间非周期信号的频谱为离散周期函数D.离散时间周期信号的频谱为离散周期函数【参考答案】: D11.题目及选项如下图所示A.AB.BC.CD.D【参考答案】: A12.要处理一个连续时间信号，对其进行采样的频率为6kHz，要无失真的恢复该连续信号，则该连续信号的最高频率可能是（）。
电子科技大学生命科学与技术学院标准实验报告（实验）课程名称数字信号处理2012-2013-第2学期电子科技大学教务处制表电子科技大学实验报告学生姓名：张华博学号：2011091010004 指导教师：李永杰实验地点：清水河校区实验时间：2013年月日一、实验室名称：科二 504 机房二、实验名称：FIR滤波器的特性及应用三、实验学时：2学时四、实验原理：五、实验目的：（详细填写）1. 熟悉相位信息的特性2. 测试FIR滤波器的滤波特性六、实验内容：（详细填写）1. 相位信息的作用2. 测试线性相位FIR滤波器的滤波特性七、实验器材（设备、元器件）：八、实验步骤：九、实验数据及结果分析：（详细填写）（包括程序、图、结果等）1、n=0:99;m=2/100:2/100:2;y=sin(0.1*pi*n)+2*cos(0.2*pi*n);XXX = fft(y);plot(m,abs(XXX));title('信号频谱图');figure(2);plot(m,unwrap(angle(XXX))); title('信号相位谱');x1=cos(0.85*pi*n);x2=0.2*n;Pha1=unwrap(angle(XXX))+x1;Pha2=unwrap(angle(XXX))+x2;figure(3)plot(m,Pha1);title('加随机相位的相位谱');figure(4)plot(m,Pha2);title('加线性相位的相位谱');cc1 = abs(XXX).*exp(i*Pha1);dd1 = ifft(cc1, 'symmetric');cc2 = abs(XXX).*exp(i*Pha2);dd2 = ifft(cc2, 'symmetric');figure(5)plot(n,y,'b',n,dd1,'r',n,dd2,': b');title('fft反变换与原信号对比图'); legend('原信号','加随机相位的信号','加线性相位的信号');0.20.40.60.811.21.41.61.820102030405060708090100信号频谱图00.20.40.60.81 1.2 1.4 1.6 1.82510152025信号相位谱00.20.40.60.81 1.2 1.4 1.6 1.82510152025加随机相位后的相位谱00.20.40.60.81 1.2 1.4 1.6 1.8251015202530354045加线性相位后的相位谱102030405060708090100-3-2-1123fft 反变换与原信号对比图2、h=[-1 2 -3 6 -3 2 -1]; %×Ô¼º¹¹Ôìh(n)ÐòÁÐ a=1; %µ±ÎªFIRÂ˲¨Æ÷ʱ£¬È¡a=1; zplane(h,a); M = 100;[H,w] = freqz(h,a,M); figure(2); plot(w,abs(H)); xlabel('\omega'); ylabel('Amplitude'); figure(3);plot(w,unwrap(angle(H))); xlabel('\omega');ylabel('Phase');h1=[-1 2 -3 -3 2 -1]; %×Ô¼º¹¹Ôìh(n)ÐòÁÐ a1=1; %µ±ÎªFIRÂ˲¨Æ÷ʱ£¬È¡a=1; figure(4); zplane(h1,a1); M = 100;[H1,w1] = freqz(h1,a1,M); figure(5); plot(w1,abs(H1));xlabel('\omega');ylabel('Amplitude');figure(6);plot(w1,unwrap(angle(H1)));xlabel('\omega');ylabel('Phase');h2=[-1 2 -3 0 3 -2 1]; %×Ô¼º¹¹Ôìh(n)ÐòÁÐa2=1; %µ±ÎªFIRÂ˲¨Æ÷ʱ£¬È¡a=1;figure(7);zplane(h2,a2);M = 100;[H2,w2] = freqz(h2,a2,M);figure(8);plot(w2,abs(H2));xlabel('\omega');ylabel('Amplitude');figure(9);plot(w2,unwrap(angle(H2)));xlabel('\omega');ylabel('Phase');h3=[-1 2 -3 3 -2 1]; %×Ô¼º¹¹Ôìh(n)ÐòÁÐa3=1; %µ±ÎªFIRÂ˲¨Æ÷ʱ£¬È¡a=1;figure(10);zplane(h3,a3);M = 100;[H3,w3] = freqz(h3,a3,M);figure(11);plot(w3,abs(H3));xlabel('\omega');ylabel('Amplitude');figure(12);plot(w3,unwrap(angle(H3)));xlabel('\omega');ylabel('Phase');-1.5-1-0.50.511.522.5-1.5-1-0.50.511.5Real PartI m a g i n a r y P a r t00.51 1.52 2.53 3.524681012141618A m p l i t u d e0.511.522.533.5-10-9-8-7-6-5-4-3-2-1P h a s e-2-1.5-1-0.500.511.52-1.5-1-0.500.511.5Real PartI m a g i n a r y P a r t00.51 1.52 2.53 3.5123456789ωA m p l i t u d e00.51 1.52 2.53 3.5-5-4-3-2-101234ωP h a s e-2-1.5-1-0.500.511.52-1.5-1-0.50.511.5Real PartI m a g i n a r y P a r t0.511.522.533.5012345678910A m p l i t u d e00.51 1.52 2.53 3.5-12-10-8-6-4-2P h a s e-1-0.500.51-1-0.8-0.6-0.4-0.200.20.40.60.81Real PartI m a g i n a r y P a r t00.51 1.52 2.53 3.524681012ωA m p l i t u d e00.51 1.52 2.53 3.5-5-4-3-2-11ωP h a s eh=[-1 2 -3 6 -3 2 -1]; a=1;coef = sum(abs(h))/sum(a);h = h/coef; zplane(h,a); M = 50;[H,w] = freqz(h,a,M); figure(2); plot(w/pi,abs(H)); figure(3);plot(w/pi,unwrap(angle(H))); n=0:49;xs=sin(0.95*pi*n); xn=sin(0.5*pi*n); x=xs+xn;y = filter(h, a, x); [h2,w2]=freqz(y) figure(4)plot(n,xs,n,y,'r');legend('Ô-ÐźÅ','Â˲¨ºó'); figure(5)plot(w1/pi,abs(h1),w2/pi,abs(h2),'r');title('ƵÆ×ͼ¶Ô±È');legend('Ô-ÐźÅ','Â˲¨ºó');-1.5-1-0.50.511.522.5-1.5-1-0.50.511.5Real PartI m a g i n a r y P a r t00.10.20.30.40.50.60.70.80.910.10.20.30.40.50.60.70.80.910.10.20.30.40.50.60.70.80.91-10-9-8-7-6-5-4-3-2-15101520253035404550-1.5-1-0.50.511.50.10.20.30.40.50.60.70.80.910510152025频谱图对比十、实验结论：（详细填写）1、给一信号相位谱加一线性相位后，做fft 反变换所得图像为原图平移一定相位；若加一随机相位，做反变换后，不仅相位发生平移，而且幅度也发生变化。
M2.22.30(b)n=[-10:10];>> x=sin(0.8*pi*n+0.8*pi); >> stem(n,x)title('2.30（b）的序列')>> xlabel('n')>> ylabel('x')2.30(c)n=[-10:10];>> x=cos(pi*n/5);>> stem(n,x)>> title('the real part')>> xlabel('n')>> ylabel('re(x[n])')>> title('the real part of 2.30(c)')x=sin(pi*n/10);>> stem(n,x)>> title('the imaginary part of 2.30(c)')>> xlabel('n')>> ylabel('im(x[n])')2.30(d)n=[-20:20];>> x=3*cos(1.3*pi*n)-4*sin(0.5*pi*n+0.5*pi);>> stem(n,x)>> title('2.30(d)的序列')>> xlabel('n')>> ylabel('x[n]')2.30(e)n=[-20:20];>> x=5*cos(1.5*pi*n+0.75*pi)+4*cos(0.6*pi*n)-sin(0.5*pi*n); >> stem(n,x)>> title('2.30(e)的序列')>> xlabel('n')>> ylabel('x[n]')M2.7编写的程序如下% Program M2.7% Computation the discreet time system of exercise 2.40alpha=input('alpha=');y0=1;y1=0.5*(y0+(alpha/y0));while abs(y1-y0)>0.0001y2=0.5*(y1+(alpha/y1));y0=y1;y1=y2;enddisp('y[n]收敛于根号alpha=');y2运行结果如下：alpha=9y[n]收敛于根号alpha=y2 =3.0000alpha=30y[n]收敛于根号alpha=y2 =5.4772经两组数据对比，知已用仿真证明，n趋于无穷时，Y[n]收敛于根号下alpha。
学院 姓名 学号 任课老师 选课号……………密……………封……………线……………以……………内……………答……………题……………无……………效…………电子科技大学二零零七至二零零 八学年第 1 学期期 期末 考试《数字信号处理》 课程考试题 A 卷 （ 120 分钟） 考试形式： 开卷 考试日期 2008年1月 24 日课程成绩构成：平时 20分， 期中 20分， 实验 20 分， 期末 40分 Test of Digital Signal Processing（Give the answer of the following problems in English as possible as you can ）Note: Problem 1 to 5 must be answered. You can choose to solute either one or both of problem6 and 7, if both, the higher score will be accounted .1. Which one of the following comments is not correct ? ( )(a) The DFT can be thought of the sampled version of the DTFT . (b) The FFT is the efficient computation algorithm for the DFT .(c) A linear-phase digital filter can be realized by either FIR or IIR systems .(d) If H(z) is the transfer function of an LPF filter , then H(-z) can be thought of an HPF filter in general .2. Two systems have impulse responses as follows ：]1[][][1−+=n n n h δδ ,]1[][][2−−=n n n h δδ, respectively. （1）Determine the DTFTs of h1[n] and h2[n].（2）Determine and sketch their magnitude responses . （3）Which type of linear-phase system are they belong to?（4）If input x[n] = 2cos(πn) ，determine the output sequences of these two systems respectively.3. Consider the system transfer function:23()3148z H z z z +=−+(1) Sketch the pole-zero plot of H(z);(2) Discuss the type of inverse Z-transform (left-sided, right-sided, and two-sided sequences) associated with each of all possible ROCs of H(z). (3) Determine the expression of causal impulse response h[n].4. Figure 4（a ） is identical to Figure 4.2 of our textbook, which illustrates the discrete-time processing of a continuous-time signal. The ideal sampler and ideal interpolator work at the same sampling frequency T F . Assume the discrete-time processing block represents an idealbandpass digital filter with frequency response)(ωj e H illustrated in Figure 4（b ）.Let )1502cos(5.1)1252sin()1002cos()(t t t t x a ×+×+×=πππ. (1) If Hz F T 100=, determine []x n , []y n and ()a y t . (2) If Hz F T150=, determine []x n , []y n and ()a y t .Figure 4(a))Figure 4(b)ω544−5−5. Let X[k] be an 8-point DFT of a length-8 real sequence x[n]. The 5 samples of X[k] are given by X[0]=0, X[1]=0, X[2]=4, X[3]=0, X[4]=0.(1) Determine and sketch X[k].(2) Determine the value of ∑=702][n n x .(3) If ]2[][8>+<=k X k G , sketch G[k] and determine the expression of g[n]which is the 8-point IDFT of G[k].ing windowed Fourier series method, design a causal bandpass FIR digital filter withfollowing specifications:center frequency 5kHzpass band edges 3.5 kHz and 6.5 kHzstop band edges 1.5 kHz and 8.5 kHzstop band attenuation 50 dB,passband ripple can be ignored.Sampling frequency F T=20kHz.Please determine the expression of the impulse response of the FIR digital filter designed.7. The normalized transfer function of order-2 lowpass Butterworth analog filter is as follows: 121)(2++=s s s H an .By means of bilinear transform, design an order-2 low-pass IIR digital filter with -3dB cutoff frequency at 1000Hz and sampling frequency at 4000Hz. You are expected to determine the expression of transfer function in closed form and sketch an implement block-diagram in any cannonic form.。
答案很详细，考试前或者平时作业的时候可以好好研究，祝各位考试成功！！电子科技大学微电子与固体电子学陈钢教授著数字信号处理课后答案1.2 教材第一章习题解答1. 用单位脉冲序列()n δ及其加权和表示题1图所示的序列。
解：()(4)2(2)(1)2()(1)2(2)4(3)0.5(4)2(6)x n n n n n n n n n n δδδδδδδδδ=+++-+++-+-+-+-+-2. 给定信号：25,41()6,040,n n x n n +-≤≤-⎧⎪=≤≤⎨⎪⎩其它（1）画出()x n 序列的波形，标上各序列的值；（2）试用延迟单位脉冲序列及其加权和表示()x n 序列； （3）令1()2(2)x n x n =-，试画出1()x n 波形； （4）令2()2(2)x n x n =+，试画出2()x n 波形； （5）令3()2(2)x n x n =-，试画出3()x n 波形。
解：（1）x(n)的波形如题2解图（一）所示。
（2）()3(4)(3)(2)3(1)6()6(1)6(2)6(3)6(4)x n n n n n n n n n n δδδδδδδδδ=-+-+++++++-+-+-+-（3）1()x n 的波形是x(n)的波形右移2位，在乘以2，画出图形如题2解图（二）所示。
（4）2()x n 的波形是x(n)的波形左移2位，在乘以2，画出图形如题2解图（三）所示。
（5）画3()x n 时，先画x(-n)的波形，然后再右移2位，3()x n 波形如 5. 设系统分别用下面的差分方程描述，()x n 与()y n 分别表示系统输入和输出，判断系统是否是线性非时变的。
（1）()()2(1)3(2)y n x n x n x n =+-+-； （3）0()()y n x n n =-，0n 为整常数； （5）2()()y n x n =； （7）0()()nm y n x m ==∑。
1.计算序列x(n)的256点DFT，需要的复数乘法次数为 ( ) 。
A.256B.256×256C.256×255D.128×8【参考答案】: B2.下列关于FIR滤波器的说法中正确的是（）A.FIR滤波器容易设计成线性相位特性B.FIR滤波器的脉冲响应长度是无限的C.FIR滤波器的脉冲响应长度是确定的D.对于相同的幅频特性要求，用FIR滤波器实现要比用IIR滤波器实现阶数低【参考答案】: C3.实序列傅里叶变换的实部和虚部分别为( )。
A.偶函数和奇函数B.奇函数和偶函数C.奇函数和奇函数D.偶函数和偶函数【参考答案】: A4.序列x(n)=u(n)的能量为（）。
A.1B.9C.11D.∞【参考答案】: D5.要处理一个连续时间信号，对其进行采样的频率为3kHz，要不失真的恢复该连续信号，则该连续信号的最高频率可能是为( )。
A.6kHzB.1、5kHzC.3kHzD.2kHz【参考答案】: B6.以下对有限长单位冲激响应（FIR）滤波器特点的论述中错误的是( )。
A.FIR滤波器容易设计成线性相位特性B.FIR滤波器的单位冲激抽样响应h(n)在有限个n值处不为零 C.系统函数H(z)的极点都在z=0处 D.实现结构只能是非递归结构【参考答案】: D7.下面关于IIR滤波器设计说法正确的是( )。
A.双线性变换法的优点是数字频率和模拟频率成线性关系B.冲激响应不变法无频率混叠现象 C.冲激响应不变法不适合设计高通滤波器 D.双线性变换法只适合设计低通、带通滤波器【参考答案】: C8.下列关于用冲激响应不变法设计IIR滤波器的说法中错误的是( )、A.数字频率与模拟频率之间呈线性关系B.能将稳定的模拟滤波器映射为一个稳定的数字滤波器 C.使用的变换是s平面到z平面的多值映射 D.可以用于设计低通、高通和带阻等各类滤波器【参考答案】: D9.在窗函数法的FIR滤波器设计法中，滤波器过渡带宽度近似等于窗函数幅度谱的主瓣宽度的（）倍。
18秋《数字信号处理》在线作业1
题目及选项如下图所示
A.A
B.B
C.C
D.D
正确答案:A
下列关于因果稳定系统说法错误的是（）。
A.极点可以在单位圆外
B.系统函数的Z变换收敛域包括单位圆
C.因果稳定系统的单位冲激响应是因果序列
D.因果稳定系统的系统函数的Z变换收敛域包括
正确答案:A
题目及选项如下图所示
A.A
B.B
C.C
D.D
正确答案:C
一个线性移不变系统稳定的充分必要条件是其系统函数的收敛域包括（）。
A.原点
B.单位圆
C.虚轴
D.实轴
正确答案:B
题目及选项如下图所示
A.A
B.B
C.C
D.D
正确答案:A。
第四章 快速傅立叶变换运算需要多少时间。
计算需要多少时间，用，问直拉点的，用它来计算每次复加速度为平均每次复乘需如果一台通用计算机的FFT DFT[x(n)]512s 5 s 50.1μμ 解: 解: ⑴ 直接计算:复乘所需时间: 复加所需时间:⑵用FFT 计算:复乘所需时间:复加所需时间:s N T N 01152.0 512log 105 log 105 2251262261=⨯⨯⨯=⨯⨯=--sT T T sN N T 013824.0 002304.0 512log 512105.0 log 105.0 2126262=+=∴=⨯⨯⨯=⨯⨯⨯=--sT T T sN N T 441536.1 130816.0 )1512(512105.0 )1(105.0 21662=+=∴=-⨯⨯⨯=-⨯⨯⨯=--s N T 31072.1 512105 105 26261=⨯⨯=⨯⨯=--运算一次完成。
点试用一个为了提高运算效率值求今需要从值的点实序列是两个已知IFFT N n y n x k Y k X DFT n y n x N k Y k X ,,)(),()(),(,)(),()(),(.2值的过程。
）（），（完成计算点）可用一次（）（）（综上所述，构造序列）（）（）（）（可得：）（）（）（再根据都是实序列，）（），（由原题可知：）（）（）（）（（）（）（性质：又根据可得序列点作对取序列依据题意解 ]Im[ ]Re[ ][][ ][ ).()( )()()( )()();()( ::n y n x IFFT N k jY k X k Z n z n y n z n x n jy n x n z n y n x n jy n x k Y jIDFT k X IDFT k jY k X IDFT DFT n z IFFT N k Z k jY k X k Z k Y n y k X n x +===+=+=+=++=⇔⇔。
