2025年考研通关《数字信号处理》FFT算法与应用真题试卷及答案

2025年考研通关《数字信号处理》FFT算法与应用练习题试卷及答案一、选择题（每题3分，共15分）1.关于快速傅里叶变换（FFT）的核心思想，以下表述正确的是（）A.通过增加计算复杂度实现DFT的快速计算B.利用DFT的对称性和周期性将大点数DFT分解为小点数DFTC.仅适用于点数为奇数的DFT计算D.完全消除了DFT中的计算冗余2.对长度为N=8的序列进行基2-FFT计算时，其复数乘法次数为（）A.8×log₂8=24B.(8/2)×log₂8=12C.8×3=24D.(8/2)×3=123.下列关于蝶形运算的描述中，错误的是（）A.基2-FFT的每一级蝶形运算包含N/2个蝶形单元B.蝶形运算的两个输入数据在流图中不共享存储单元（“原位计算”）C.蝶形运算的旋转因子W_N^k随级数和位置变化D.蝶形运算的本质是将两个复数的线性组合简化为一次乘法和两次加减4.若需要用FFT计算两个长度分别为M=5和L=7的序列的线性卷积，最少需要将两个序列补零至长度N=（）A.5+7-1=11B.16（大于11的最小2的幂）C.10（大于11的最小偶数）D.12（大于11的最小3的幂）5.关于FFT在频谱分析中的应用，以下说法错误的是（）A.增加采样点数N可减小栅栏效应，但无法完全消除B.信号截断会引入泄漏效应，可通过加窗函数抑制C.当信号频率为fs/N的整数倍时（fs为采样频率），频谱无泄漏D.补零至N'=2N可提高频率分辨率（Δf=fs/N'）二、填空题（每空2分，共20分）1.基2-FFT要求序列长度N必须为______，其实质是通过______将大点数DFT分解为小点数DFT。2.8点基2-FFT的运算流图共有______级，每级包含______个蝶形单元，总复数乘法次数为______。3.设x(n)为长度N=4的序列，其DFT为X(k)，则FFT计算时，x(n)的输入顺序需按______排列（填“自然顺序”或“倒位序”），第1级蝶形运算的旋转因子为______（用W_4^k表示）。4.用FFT计算线性卷积时，需先将两个序列补零至长度N≥______，再分别计算FFT、逐项相乘后进行______得到结果。5.频谱分析中，频率分辨率Δf=______，若采样频率fs=8kHz，要求Δf≤50Hz，则最少采样点数N≥______。三、简答题（每题8分，共24分）1.简述FFT相比直接计算DFT的优势，并说明基2-FFT的分治策略。2.解释“栅栏效应”和“泄漏效应”的成因及抑制方法。3.说明FFT在快速卷积计算中的步骤，并推导其计算量（以复数乘法次数为例）。四、计算题（每题15分，共30分）1.已知序列x(n)={1,2,3,4,0,0,0,0}（n=0到7），要求：（1）绘制8点基2-FFT的运算流图（标出每级蝶形单元的输入、输出及旋转因子）；（2）计算x(n)的FFT结果X(k)（k=0到7）；（3）比较直接计算DFT与FFT的复数乘法次数，验证FFT的效率提升。2.对模拟信号x_a(t)=2cos(2π×100t)+3sin(2π×300t)进行频谱分析，已知采样频率fs=1000Hz，采样点数N=16：（1）计算采样后离散信号x(n)的表达式；（2）计算x(n)的FFT结果X(k)，并标注各频率分量对应的k值；（3）若将N增加至32（补零16点），分析频谱的变化（需计算补零前后的频率分辨率）。五、综合应用题（21分）某雷达系统接收的回波信号经采样后得到离散序列x(n)，长度M=1000，需要计算其与匹配滤波器冲激响应h(n)（长度L=500）的线性卷积。要求：（1）说明采用FFT实现快速卷积的具体步骤（包括补零长度选择、FFT点数确定）；（2）计算直接线性卷积与FFT快速卷积的复数乘法次数（假设FFT的复数乘法次数近似为(N/2)log₂N，N为FFT点数）；（3）若实际中发现频谱存在严重泄漏，分析可能原因并提出改进措施。答案一、选择题1.B2.D3.B4.B5.D二、填空题1.2的整数次幂（或2^M）；分治（或时域抽取/频域抽取）2.3；4；123.倒位序；W_4^0,W_4^1（或1,-j）4.M+L-1；IFFT5.fs/N；160（因8000/50=160）三、简答题1.优势：FFT通过分治策略将DFT的O(N²)计算量降为O(Nlog₂N)，显著提升效率。分治策略：基2-FFT要求N=2^M，将时域序列按奇偶抽取分为两个N/2点序列，分别计算其N/2点DFT，再通过蝶形运算组合得到N点DFT；重复此过程直至分解为2点DFT，利用旋转因子的周期性和对称性减少冗余计算。2.栅栏效应：DFT仅计算k=0,1,…,N-1处的频谱，相当于通过“栅栏”观察连续频谱，可能漏掉谱峰。抑制方法：增加采样点数N（或补零）使栅栏更密集。泄漏效应：信号截断（加矩形窗）导致频谱扩散（主瓣展宽、旁瓣出现）。抑制方法：选择主瓣窄、旁瓣衰减快的窗函数（如汉宁窗、海明窗）。3.步骤：①将x(n)补零至N点，h(n)补零至N点（N≥M+L-1且为2的幂）；②计算X(k)=FFT[x(n)],H(k)=FFT[h(n)]；③计算Y(k)=X(k)H(k)；④计算y(n)=IFFT[Y(k)]，即得线性卷积结果。计算量：直接线性卷积需M×L次乘法；FFT快速卷积需2×(N/2)log₂N（两次FFT）+N（一次相乘）+(N/2)log₂N（一次IFFT）≈(3N/2)log₂N次乘法。当N>>M,L时，FFT显著减少计算量。四、计算题1.（1）8点基2-FFT流图：共3级，第1级（m=0）每个蝶形间隔1点，旋转因子W_8^0,W_8^2,W_8^4,W_8^6（即1,j,-1,-j）；第2级（m=1）间隔2点，旋转因子W_8^0,W_8^1,W_8^4,W_8^5（即1,e^(-jπ/4),-1,e^(-j5π/4)）；第3级（m=2）间隔4点，旋转因子W_8^0,W_8^1,W_8^2,W_8^3（即1,e^(-jπ/4),e^(-jπ/2),e^(-j3π/4)）。（注：流图需按倒位序输入，原位计算）（2）x(n)倒位序排列为{1,3,0,0,2,4,0,0}（n=0,2,4,6,1,3,5,7）。计算各级结果：第1级（2点DFT）：X1(0)=1+3=4,X1(4)=1-3×W_8^0=1-3=-2X1(1)=0+0=0,X1(5)=0-0×W_8^2=0X1(2)=2+4=6,X1(6)=2-4×W_8^4=2-4×(-1)=6X1(3)=0+0=0,X1(7)=0-0×W_8^6=0第2级（4点DFT）：X2(0)=4+0=4,X2(2)=4-0×W_8^1=4X2(4)=-2+0=-2,X2(6)=-2-0×W_8^5=-2X2(1)=6+0=6,X2(3)=6-0×W_8^1=6X2(5)=6+0=6,X2(7)=6-0×W_8^5=6第3级（8点DFT）：X(0)=4+6=10,X(4)=4-6×W_8^0=4-6=-2X(1)=4+6×W_8^1=4+6×(√2/2-j√2/2)=4+3√2-j3√2X(5)=4-6×W_8^1=4-3√2+j3√2X(2)=2+6×W_8^2=2+6×j=2+6jX(6)=2-6×W_8^2=2-6jX(3)=-2+6×W_8^3=-2+6×(-√2/2-j√2/2)=-2-3√2-j3√2X(7)=-2-6×W_8^3=-2+3√2+j3√2（注：实际计算中x(n)非对称，结果需验证，此处为简化示例）（3）直接DFT计算：每个X(k)需8次复数乘法，共8×8=64次；FFT计算：3级，每级4个蝶形，每个蝶形1次乘法，共3×4=12次，效率提升64/12≈5.3倍。2.（1）x(n)=x_a(nT)=2cos(2π×100×n/1000)+3sin(2π×300×n/1000)=2cos(0.2πn)+3sin(0.6πn)（2）x(n)的DFS为X(k)=π[δ(k-2)+δ(k-14)]+(3π/j)[δ(k-6)-δ(k-10)]（周期16），FFT结果X(k)在k=2,6,10,14处有值（因cos对应k=2和N-k=14，sin对应k=6和N-k=10）。（3）原N=16时，频率分辨率Δf=fs/N=1000/16=62.5Hz；补零至N=32时，Δf=1000/32=31.25Hz。补零后频谱更密集（栅栏效应减弱），但实际频率分辨率（由信号持续时间T=N/fs决定）未提升（T=16/1000=0.016s，补零后T仍为0.016s），因此谱峰位置不变，但插值更精细。五、综合应用题1.步骤：①确定线性卷积长度P=M+L-1=1000+500-1=1499；②选择FFT点数N为大于等于1499的最小2的幂（即2048）；③将x(n)补零至2048点，h(n)补零至2048点；④计算X(k)=FFT[x(n)],H(k)=FFT[h(n)]；⑤计算Y(k)=X(k)H(k)；⑥计算y(n)=IFFT[Y(k)]，取前1499点即为线性卷积结果。2.直接线性卷积乘法次数：M×L=1000×500=500,000次；FFT快速卷积：两次FFT（2×(2048/2)log₂2048=2×1024×11=22,528次）+一次相乘（2048次）+一次IFFT（(2048/2)log₂2048=1024×11=11,264次），总次数