实训报告信号.docVIP

无锡太湖学院 实训报告 实训课程名称：数字信号处理课程设计 学 院：工学院 专 业： 指导教师： 完成日期： 2016年月日 滤波器源程序： ws=0.125*pi;wp=0.075*pi;Rp=0.25;As=50; [N,wc]=buttord(wp/pi,ws/pi,Rp,As); [b,a]=butter(N,wc); [h,w]=freqz(b,a); plot(w/pi,abs(h)); title(滤波器的频率响应); xlabel(w/π); ylabel(幅度/dB); 滤波器的频率响应如下： 3.滤波 源程序如下： [y,fs,nbits]=wavread(F:\song\zhaohaiyan1.wav); ws=0.125*pi;wp=0.075*pi;Rp=0.25;As=50; [N,wc]=buttord(wp,ws,Rp,As); [b,a]=butter(N,wc); [Rows,Cols]=size(y); y(:,1)=y(:,1)-mean(y(:,1)); Y=abs(fft(y)+eps);%为防止出现0值 Y(1)=Y(2); L=Rows; w=0:(L-1); w=2*pi/Rows*w; x=filter(b,a,y); X=abs(fft(x)+eps);%为防止出现0值 X(1)=X(2); figure(1); subplot(2,2,1);plot(y);title(滤波前信号波形); subplot(2,2,2);plot(w/pi,Y(1:L));title(滤波前信号频谱);xlabel(w(单位pi));ylabel(幅值); subplot(2,2,3);plot(x);title(滤波后信号波形); subplot(2,2,4);plot(w/pi,X(1:L));title(低通滤波后信号频谱);xlabel(w(单位pi));ylabel(幅值); wavwrite(x,F:\song\zhaohaiyan2.wav); sound(x); 运行结果： 4.GUI界面 源程序： function varargout = untitled1(varargin) % UNTITLED1 M-file for untitled1.fig % UNTITLED1, by itself, creates a new UNTITLED1 or raises the existing % singleton*. % % H = UNTITLED1 returns the handle to a new UNTITLED1 or the handle to % the existing singleton*. % % UNTITLED1(CALLBACK,hObject,eventData,handles,...) calls the local % function named CALLBACK in UNTITLED1.M with the given input arguments. % % UNTITLED1(Property,Value,...) creates a new UNTITLED1 or raises the % existing singleton*. Starting from the left, property value pairs are % applied to the GUI before untitled1_OpeningFunction gets called. An % unrecognized property name or invalid value makes property application % stop. All inputs are passed to untitled1_OpeningFcn via varargin. % % *See GUI Options on GUIDEs Tools menu. Choose GUI allows only one % instance to run (singleton). % % See also: GUIDE, GUIDATA, GUIHANDLES % Copyright 2002-2003 The MathWorks, Inc. % Edit the a