数字信号课程设计精要.docxVIP

课程设计报 告 音频信号在MATLAB的生成及处理 班级：通信工程电子121班 姓名：周源思 学号：1207300502 指导老师：文珺 2014年 6 月 实验14 一、实验目的 加深理解采用数字信号处理方法对模拟信号处理的过程、掌握使用MATLAB处理的方法。 二、实验内容和步骤 1）采用MATLAB设计一个类似例4.3.1中的有多个不同频率的音频信号，长度3秒，自行选择合适的采样频率Fs；例4.3.1中的音频信号为钢琴音频，前三个音频为262,330,392，后面音频范围为493-2500Hz; 2）对产生的音频信号进行处理 ?设计一个预滤波器，对输入的信号进行预滤波，以消除超过1/2Fs頻率的信号，防止采样时出现信号混叠。设计参数如下： 通带截止频率fp=2500; 阻带截止频率fs=3000; 通带最大衰减 rp=1; 阻带最小衰减rs=100; ?设计一个IIR数字滤波器，对预滤波后的信号进行处理。要求提取出前三个低频基本音，并通过衰减第三个音频幅度，使得前三个钢琴的声基本平衡。根据要求可知滤波器在0~0.282Π频带内幅度平稳，在0.282Π~0.329Π幅度单调下降。设计参数如下： 通带截止频率fp=360; 阻带截止频率fs=420; 通带最大衰减 rp=1; 阻带最小衰减rs=-15; 3）画出处理前后的时域频域图，要求： ?音频信号处理前后时域波形图要求横坐标是时间，单位s，频域频谱图横坐标要求是模拟频率（范围0-Fs/2），单位Hz； ?滤波器幅度特性画两幅图，分别是纵坐标是dB、横坐标归一化数字频率（范围0-1）和纵坐标是幅度、横坐标是模拟频率（范围0-Fs/2），单位Hz； ?画出相应相位图； 回放音频信号并对处理前后的音频信号频谱进行对比、分析； 设计GUI； 三 、具体算法 滤波器都是使用巴特沃思数字低通滤波器，数字低通滤波器参数可通过双线性变换法，由模拟巴特沃思低通滤波器参数得到。给出滤波器参数 fp，fs,rp,rs 通带波纹幅度参数 阻带波纹幅度参数 过渡比参数 偏离参数 模拟巴特沃思低通滤波器系统函数由阶数N和3dB截止频率决定 四 、GUI设计 利用 GUI 图形用户界面设计工具制作数字滤波器处理面板，框架如下： 留有三个音频频率输入端口，主页面有六个窗口，分别可以画出原始信号，滤波后信号的时域，频域，和相位图。旁边有“预滤波器”按钮，执行如下程序 function pushbutton6_Callback(hObject, eventdata, handles) % hObject handle to pushbutton6 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global time global fx global piano2 figure(2) subplot(3,1,1); plot(time,piano2);%滤波后语音信号的时域图形 axis ([0 3 -2.2 2.2]); title(预滤波后音频信号时域图形);xlabel(时间(s)); ylabel(幅度); subplot(3,1,2); y2=fft(piano2,2048); %对信号做2048个FFT变换 plot(fx,abs(y2(1:1024)),b); title(预滤波后语音信号的幅频特性（模拟频率）范围0-Fs/2）);xlabel(频率（Hz）); ylabel(幅度); subplot(3,1,3); ph2=angle(y2(1:1024));%Fs/2的相位 plot(fx,ph2); title(预滤波后语音信号的相频图（模拟频率）范围0-Fs/2）);xlabel(频率（Hz）); ylabel(幅度); 将画出预滤波后信号的时域，频域，和相位图。 “数字滤波器”按钮，执行如下程序： function pushbutton5_Callback(hObject, eventdata, handles) % hObject handle to pushbutton5 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global Bz global Az global f