信号发生器数字信号理课程设计报告.docVIP

中南大学 数字信号处理课程设计报告 专业班级: 指导老师： 姓 名: 学 号: 目 录 一、课程设计要求心得体会课程设计要求 二.设计过程 第一题: ⑴、信号发生器 用户根据测试需要，可任选以下两种方式之一生成测试信号： ①、直接输入（或从文件读取）测试序列； ②、输入由多个不同频率正弦信号叠加组合而成的模拟信号公式（如式 1-1 所示）、采 样频率（Hz）、采样点数，动态生成该信号的采样序列，作为测试信号。 100sin(2πf1t) +100sin(2πf2t ) +....+100sin(2πfnt ) ? ? ? ?（1-1） ⑵、源代码 t=0:1/20000:0.01; f1=2000,f2=5000,f3=8000; ft=100*sin(2*pi*f1*t)+100*sin(2*pi*f2*t)+100*sin(2*pi*f3*t); plot(t,ft); grid on; axis([0 0.01 1.1*min(ft) 1.1*max(ft)]); xlabel(t),ylabel(ft); ⑶、结果 ⑷分析 在给定的时域内，输入采样序列，生成了含有三个不同频率的混合信号波形。 第二题 ⑴、频谱分析 使用 FFT 对产生的测试信号进行频谱分析并展示其幅频特性与相频特性，指定需要滤除的频带，通过选择滤波器类型（IIR / FIR），确定对应的滤波器（低通、高通）技术指标。 ⑵、源代码 fs=20000;N=512;T=1/fs; t=0:1/fs:0.01; f1=2000,f2=5000,f3=8000; ft=100*sin(2*pi*f1*t)+100*sin(2*pi*f2*t)+100*sin(2*pi*f3*t); subplot(4,1,1); plot(t,ft);grid on; axis([0 0.01 1.1*min(ft) 1.1*max(ft)]); xlabel(t),ylabel(ft); title(抽样信号的连续形式); subplot(4,1,2); stem(t,ft);grid on; axis([0 0.01 1.1*min(ft) 1.1*max(ft)]); xlabel(t),ylabel(ft); title(实际抽样信号); k=0:N-1; Fw=fft(ft,N); subplot(4,1,3); plot(k,abs(Fw));grid on; axis([0 600 0 10000]); title(抽样信号幅度谱) subplot(4,1,4); plot(k,angle(Fw));grid on; axis([0 100 -4 4]); title(抽样信号相位谱) ⑶、结果 ⑷分析 在抽样信号的幅度谱中发现了三个频率的幅度谱，需要进行滤波。 第三题 ⑴、题目 滤波器设计 根据以上技术指标（通带截止频率、通带最大衰减、阻带截止频率、阻带最小衰减） ，设计数字滤波器，生成相应的滤波器系数，并画出对应的滤波器幅频特性与相频特性。 ①IIR DF 设计：可选择滤波器基型（巴特沃斯或切比雪夫型）； ②FIR DF 设计：使用窗口法（可选择窗口类型，并比较分析基于不同窗口、不同阶数 所设计数字滤波器的特点）。 ⑵、源代码 %IIR DF设计： %Step1:产生含有3个正弦分量的信号 Fs=20000; t=(1:100)/Fs; %抽样频率、时间轴 s1=100*sin(2*pi*t*2000);s2=100*sin(2*pi*t*5000);s3=sin(2*pi*t*8000); s=s1+s2+s3; subplot(511); %组成信号、指定图形位置 plot(t,s);title( 原始信号); xlabel(时间(s) ); ylabel(信号波形 ); %Step2:产生一个4阶IIR带通滤波器 %通带为100Hz到200Hz,并得出其幅频响应 [b,a]=ellip(2,0.5,20,[1000,2500]*2/Fs); %得到滤波器的系数矩阵 [H,w]=freqz (b,a,512); %H为滤波器的系统函数 subplot(512); plot(w*Fs/(2*pi),abs(H)); title(IIR带通滤波器幅频响应 ); xlabel( 频率(Hz) );ylabel( 幅度 ); axis([0 10000 0 1.5]); figure %FIR DF设计： As=70; ws=0