1、 现代信号处理课程设计报告现代信号处理课程设计报告 专业班级专业班级: :电子信息班电子信息班 目目 录录 1、课程设计要求课程设计要求 2、设计过程、设计过程 （1）设计题目）设计题目 （2）设计源代码）设计源代码 （3）设计结果）设计结果 （4）结果分析）结果分析 3、设计、设计总结与总结与心得体会心得体会 4、课程设计指导书课程设计指导书 5、鸣谢、鸣谢 一、一、课程设计要求课程设计要求 使用 MATLAB（或其它开发工具）编程实现上述内 容，写出课程设计报告。 二二. .设计过程设计过程 第一题第一题: : (1).(1).请设计实验对采样定理进行验证。 请画出其傅立叶变换图； 请画出
2、信号重建图（可以采用不同的内插方法） (2).源代码源代码 t=(0:0.00001:1); n=0:100; f1=50; y=sin(2*pi*f1*t); f=input(please f=); T=1/f; x=sin(2*pi*f1*n*T); m=fft(x); h=abs(m); figure(1); subplot(321) plot(t,y); subplot(322) stem(n,x,.); title(xulitu); subplot(323) plot(n,h); title(fupintu); subplot(324) xi=interp1(n,x,t*f1,lin
3、ear); plot(t,xi); title(chongjiantu); （3）结果）结果 （4）分析：分析： 采样原理： 对模拟信号进行采样可以看作是一个模拟信号通过一个电 子开关 S。设电子开关每隔周期 T 合上一次，每次合上的时间为， 在电子开关输出端得到其采样信号，一般很小, 越小，采样输出 脉冲的幅度越接近输入信号在离散时间点上的瞬时值。 第二题第二题 （1） 题目题目 请设计实验对序列傅立叶变换的频移特性 进行验证。 （2） 源代码源代码 n=0:100; x=sin(n).*exp(j*0*n); m=sin(n).*exp(j*1*n); figure(1); y=fft(x
4、); z=fft(m); figure(1); subplot(211) plot(n,y); title(0yixiangwei); subplot(212) plot(n,z); title(yidonghou); （3） 结果结果 （4）分析：分析：当时域附加一个相位移后,相应的原信号在频域的频谱会 发生相应的频率移位 第三题第三题 （1） 题目题目 请设计实验比较 DFT 和 FFT，比较它们 的计算结果和计算速度。 （4） 源代码源代码 N=32; x=ones(1,5),zeros(1,N-5); n=0:N-1; tic; Xk=dft(x,N); toc; tic; Zk=ff
5、t(x); toc; nagX=abs(Zk); magX=abs(Xk); k=(0:length(magX)-1)*N/length(magX); subplot(311); stem(n,x,.); subplot(312); stem(k,magX,.); subplot(313); stem(k,nagX,.); axis(0,32,0,5); functionXk=dft(xn,N) n=0:1:N-1; k=0:1:N-1; WN=exp(-j*2*pi/N); nk=n*k; WNnk=WN.nk; Xk=xn*WNnk; （3）结果）结果 Elapsed time is 0.
6、756333 seconds. Elapsed time is 0.000633 seconds. (4) 分析分析 快速傅里叶变换运算速度要比离散傅快速傅里叶变换运算速度要比离散傅里叶变换速度快里叶变换速度快 第四题第四题 （1） 题目题目 滤波器设计根据输入的数字滤波器的技术指 标，包括通带截止频率,通带最大衰减,阻带截止 频率,阻带最小衰减， 设计滤波器， 生成相应的滤 波器系数， 并画出对应的滤波器幅频、 相频特性。 IIR DF 设计：可选择滤波器基型（巴特沃斯或切比雪 夫型） ； （2） 源代码源代码 i=input(please input i(choose fuction)=); switch fix(i) case 1%lpdf wp=input(please input wp=); ws=input(please input ws=); ap=input(please input ap=
