数字信号分析理论与实践 (18).pdf

《数字信号分析理论与实践 (18).pdf》由会员分享，可在线阅读，更多相关《数字信号分析理论与实践 (18).pdf（20页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数字信号分析理论与实践理论知识1.信号函数和信号波形2.信号时域波形分析3.信号频域频谱分析4.信号时差域相关分析5.信号幅值域分析6.信号数字滤波方法7.信号时频域分析8.计算机测试系统结构大作业1.信号发生器和数字电子琴设计2.声音信号采集和频谱分析程序设计3.图像采集和图像直方图均衡程序设计4.MP3音乐信号滤波和数字均衡器设计5.选作：基于Arduino的计算机测控应用数字信号分析理论与实践项目4:音乐数字均衡器设计数字信号分析理论与实践3060Hz沉闷。60100Hz沉重。100200Hz丰满200500Hz力度。5001KHz明朗。1K2KHz透亮。2K4Kz尖锐。4K8Kz清脆。

2、8K16Kz纤细。不同信号频率的音感特征：频率通带1、FFT滤波滤波2、滤波器组、滤波器组3、FIR2/FIRLS数字信号分析理论与实践3060Hz沉闷。60100Hz沉重。100200Hz丰满200500Hz力度。5001KHz明朗。1K2KHz透亮。2K4Kz尖锐。4K8Kz清脆。8K16Kz纤细。V1V9归一化频率F/Fs4410030441006044100160001v1v9v数字信号分析理论与实践1、FFT滤波滤波信号FFTXRXI频域滤波HRHI镜像数字信号分析理论与实践2、滤波器组、滤波器组)(kx.1v2v3v9v数字信号分析理论与实践3、FIR2/FIRLS项目4:音乐数字

3、均衡器设计1.新建文件2.安放控件项目4:音乐数字均衡器设计3.修改控件属性1项目4:音乐数字均衡器设计3.为老版本Malab添加mp3文件读取功能在工程文件目录下建立mp3driver路径，为老版本Malab添加MP3文件读写功能。将驱动文件拷贝到其中，不需安装。从从“http:/labrosa.ee.columbia.edu/matlab/mp3read.html”下载下载项目4:音乐数字均衡器设计项目4:音乐数字均衡器设计项目4:音乐数字均衡器设计4.读mp3/wav文件按钮回调函数function pushbutton1_Callback(hObject,eventdata,handl

4、es)global abc;global len;global Fs;global x;m0=mfilename(fullpath);i=findstr(m0,);m1=m0(1:i(end)-1);m2=m1,mp3Driver;addpath(m2);FileName,PathName,filterindex=uigetfile(*.mp3,File Selector);set(handles.edit1,String,FileName);abc=fullfile(PathName,FileName);全局变量，文件名，大小，采样频率，数据全局变量，文件名，大小，采样频率，数据添加添加mp

5、3Driver工作路径工作路径选择选择mp3文件文件项目4:音乐数字均衡器设计SIZ=MP3READ(abc,size)len=SIZ(1);s1=sprintf(%f,len);set(handles.edit2,String,s1);x,Fs,nbits,OPTS =mp3read(abc,661500,1,1,1);dt=1.0/Fs;N=length(x);dt=1.0/Fs;T=N*dt;t=linspace(0,T,N);plot(handles.axes1,t,x);Amax=max(x);ylim(handles.axes1,-Amax,Amax);获取选中的获取选中的MP3文

6、件大小文件大小读读15秒的数据秒的数据x,Fs,nbits,OPTS =mp3read(文件名文件名,读入点数读入点数,单单/立体声立体声,选抽率选抽率,起始点起始点);mp3read格式：格式：项目4:音乐数字均衡器设计5.播放声音按钮function pushbutton2_Callback(hObject,eventdata,handles)global x;global Fs;wavplay(x,Fs);6.声音片选取滚动条function slider1_Callback(hObject,eventdata,handles)global abc;global len;global F

7、s;global x;p1=get(handles.slider1,Value);delay=floor(len-661500)*p1);if delay1 delay=1;endx,Fs,nbits,OPTS =mp3read(abc,661500,1,1,delay);dt=1.0/Fs;N=length(x);T=N*dt;t=linspace(0,T,N);plot(handles.axes1,t,x);Amax=max(x);ylim(handles.axes1,-Amax,Amax);按滚动条位置选声音片按滚动条位置选声音片项目4:音乐数字均衡器设计6.数字均衡按钮function

8、 pushbutton3_Callback(hObject,eventdata,handles)global Fs;global x;global y;v1=get(handles.slider2,Value)v2=get(handles.slider3,Value)v3=get(handles.slider4,Value)v4=get(handles.slider5,Value)v5=get(handles.slider6,Value)v6=get(handles.slider7,Value)v7=get(handles.slider8,Value)v8=get(handles.slider

9、9,Value)v9=get(handles.slider10,Value)f=0 0.0028 0.0057 0.0113 0.0227 0.0454 0.0907 0.1814 0.3628 0.7256 1;m=0 v1 v2 v3 v4 v5 v6 v7 v8 v9 0;b=fir2(100,f,m);y=filter(b,1,x);N=length(x);dt=1.0/Fs;T=N*dt;t=linspace(0,T,N);plot(handles.axes2,t,y);Amax=max(x);ylim(handles.axes2,-Amax,Amax);读取数字均衡器设定读取数字均

10、衡器设定数字均衡数字均衡/滤波滤波项目4:音乐数字均衡器设计7.滤波信号播放按钮function pushbutton4_Callback(hObject,eventdata,handles)global yglobal Fswavplay(y,Fs);8.数字均衡滤波器改进f=0 0.003 0.006 0.011 0.023 0.045 0.091 0.181 0.363 0.726 1;m=0 v1 v2 v3 v4 v5 v6 v7 v8 v9 0;b=fir2(100,f,m)y=filter(b,1,x);IIR滤波器组FIR滤波器组FFT频域滤波项目4:音乐数字均衡器设计运行效果：项目4:音乐数字均衡器设计项目练习:1.编程实现该项目,2.尝试用不同的数字均衡滤波算法.数字信号分析理论与实践