数字信号处理实验八.docx
《数字信号处理实验八.docx》由会员分享,可在线阅读,更多相关《数字信号处理实验八.docx(15页珍藏版)》请在冰豆网上搜索。
数字信号处理实验八
数字信号处理实验报告
实验八音频频谱分析设计与实现
班级:
20100821
学号:
2010061312
姓名:
王晓浩
实验八音频频谱分析仪设计与实现
1、实验原理
MATLAB是一个数据信息和处理功能十分强大的工程实用软件,其数据采集工具箱为实现数据的输入和输出提供了十分方便的函数和命令。
本实验要求基于声卡与MATLAB实现音频信号频谱分析仪的设计原理与实现,功能包括:
(1)信号的输入,从声卡输入、从WAV文件输入、从标准信号发生器输入;
(2)信号的波形分析,包括幅值、频率、周期、相位的估计,以及统计量峰值、均值、均方值和方差的计算;
(3)信号频谱分析,频率、周期的估计,图形显示幅值谱、相位谱、实频谱、虚频谱和功率谱的曲线。
1、频率(周期)检测
对周期信号来说,可以用时域波形分析来确定信号的周期,也就是计算相邻的两个信号波峰的时间差或过零点的时间差。
这里采用过零点(ti)的时间差T(周期)。
频率即为f=1/T,由于能够求得多个T值(ti有多个),故采用它们的平均值作为周期的估计值。
2、幅值检测
在一个周期内,求出信号的最大值ymax与最小值ymin的差的一半,即A=(ymax-ymin)/2,除以第一个以外的A值的平均作为幅值的估计值。
3、相位检测
采用过零法,即通过判断与同频零相位信号过零点时刻,计算其时间差,然后换成相应的相位差。
Fi=2pi(1-ti/T),同样,以fi的平均值作为相位的估计值。
4、数字信号统计量估计
(1)峰值P的估计
在样本数据x中找出最大值与最小值,其差值为双峰值,双峰值的一半即为峰值。
(2)均值估计
式中,N为样本容量,下同。
(3)均方值估计
(4)方差估计
5、频谱分析原理
是与分析只能反映信号的幅值随时间的变化情况,除单位频率分量的简单波形外,很难明确提示信号的频率组成和个频率分量的大小,而频谱分析能很好的解决此问题。
(1)DFT与FFT
对于给定的时域信号y,可以通过Foutier变换得到频域信息Y。
Y可按下式计算
式中,N为样本容量,t=1/Fs为采样间隔。
采样信号的频谱是一个连续的频谱,不可能计算出所有的点的值,故采用离散傅立叶变换(DFT),即
式中,f=Fs/N。
但上式的计算效率很低,因为有大量的指数(等价于三角函数)运算,故实际中多采用快速傅立叶变换(FFT)。
其原理是将重复的三角函数计算的中间结果保存起来,以减少重复三角函数计算带来的时间浪费。
由于三角函数计算的重复量相当大,故FFT能极大地提高运算效率。
(2)频率、周期的估计
对于Y(kf),如果当kf=f时,Y(kf)取最大值,则f为频率的估计量,由于采样间隔的误差,f也存在误差,其误差最大值为f/2。
周期T=1/f。
从原理上可以看出,如果在标准信号中混有噪声,用上述方法仍能够精确地估计出原标准信号的频率和周期。
(3)频谱图
为了直观的表示信号的频率特性,工程上常常将傅立叶变换的结果用于图形的方式表示,即频谱图。
以频谱f为横坐标,|Y(f)|为纵坐标,可以得到幅值谱;
以频谱f为横坐标,argY(f)为纵坐标,可以得到相位谱;
以频谱f为横坐标,ReY(f)为纵坐标,可以得到实频谱;
以频谱f为横坐标,ImY(f)为纵坐标,可以得到虚频谱;
根据采样定理,只有频率不超过Fs/2的信号才能被正确采集,即傅立叶变换的结果中频率大于Fs/2的部分是不正确的部分,故不在频谱图中显示。
即横坐标f属于[0,Fs/2]。
2、实验内容
请参考书中所给的界面和代码,查阅相关资料,构建交互界面,设计一个音频频谱分析仪,实现相应的功能。
3、实验程序及结果
(1)
分析仪界面:
输入界面:
(2)声卡输入方式相关程序及结果:
functionsoundcard_Callback(hObject,eventdata,handles)
set(findobj('Tag','recordtime'),'enable','on');
h=findobj('Tag','filename');
set(h,'enable','off');
h=findobj('Tag','freq');
set(h,'enable','off');
h=findobj('Tag','amp');
set(h,'enable','off');
h=findobj('Tag','phase');
set(h,'enable','off');
set(handles.channel,'enable','off');
set(handles.fileopen,'enable','off');
set(handles.gensig,'enable','off');
set(handles.wavetype,'enable','off');
set(handles.add,'enable','off');
set(handles.startrecord,'enable','on');
functionstartrecord_Callback(hObject,eventdata,handles)
Fs=str2double(get(findobj('Tag','samplerate'),'String'));
handles.y=wavrecord(str2double(get(findobj('Tag','recordtime'),'String'))*Fs,Fs,'int16');
handles.inputtype=1;
guidata(hObject,handles);
plot(handles.time,handles.y);
title('WAVE');
ysize=size(handles.y)
set(handles.samplenum,'String',num2str(ysize
(1)));
滚动条程序:
functionscale_Callback(hObject,eventdata,handles)
val=get(hObject,'value');
val=10^(val*5+1);
x=get(handles.xmove,'Value')*str2double(get(handles.samplenum,'String'));
axis(handles.time,[xval+xmin(handles.y)max(handles.y)]);
functionscale_CreateFcn(hObject,eventdata,handles)
ifisequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor',[.9.9.9]);
end
functionxmove_Callback(hObject,eventdata,handles)
val=get(handles.scale,'value');
val=10^(val*5+1);
x=get(hObject,'Value')*str2double(get(handles.samplenum,'String'));
axis(handles.time,[xval+xmin(handles.y)max(handles.y)]);
functionxmove_CreateFcn(hObject,eventdata,handles)
ifisequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor',[.9.9.9]);
end
functionscale2_Callback(hObject,eventdata,handles)
val=get(hObject,'value');
val=10^(val*5+1);
xlim(handles.plot1,[0val]);
xlim(handles.plot2,[0val]);
xlim(handles.plot3,[0val]);
xlim(handles.plot4,[0val]);
functionscale2_CreateFcn(hObject,eventdata,handles)
ifisequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor',[.9.9.9]);
end
(3)WAV文件输入方式相关程序及结果:
functionWAVfile_Callback(hObject,eventdata,handles)
h=findobj('Tag','filename');
set(h,'enable','on');
h=findobj('Tag','freq');
set(h,'enable','off');
h=findobj('Tag','amp');
set(h,'enable','off');
h=findobj('Tag','phase');
set(h,'enable','off');
set(findobj('Tag','recordtime'),'enable','off');
set(handles.channel,'enable','on');
set(handles.fileopen,'enable','on');
set(handles.gensig,'enable','off');
set(handles.wavetype,'enable','off');
set(handles.add,'enable','off');
set(handles.startrecord,'enable','off');
functionfileopen_Callback(hObject,eventdata,handles)
temp=wavread(get(findobj('Tag','filename'),'String'));
channel=str2double(get(handles.channel,'String'));
handles.y=temp(:
channel);
handles.inputtype=2;
guidata(hObject,handles);
plot(handles.time,handles.y);
title('WAVE');
ysize=size(handles.y)
set(handles.samplenum,'String',num2str(ysize
(1)));
(4)信号发生器输入方式相关程序及结果:
functiongenerator_Callback(hObject,eventdata,handles)
h=findobj('Tag','filename');
set(h,'enable','off');
h=findobj('Tag','freq');
set(h,'enable','on');
h=findobj('Tag','amp');
set(h,'enable','on');
h=findobj('Tag','phase');
set(h,'enable','on');
set(findobj('Tag','recordtime'),'enable','off');
set(handles.channel,'enable','off');
set(handles.fileopen,'enable','off');
set(handles.gensig,'enable','on');
set(handles.wavetype,'enable','on');
set(handles.add,'enable','on');
set(handles.startrecord,'enable','off');
functiongensig_Callback(hObject,eventdata,handles)
Fs=str2double(get(findobj('Tag','samplerate'),'String'));
N=str2double(get(findobj('Tag','samplenum'),'String'));
x=linspace(0,N/Fs,N);
soundtype=get(handles.wavetype,'Value');
frequency=str2double(get(handles.freq,'String'));
amp=str2double(get(handles.amp,'String'));
phase=str2double(get(handles.phase,'String'));
switchsoundtype
case1
y=amp*sin(2*pi*x*frequency+phase);
case2
y=amp*sign(sin(2*pi*x*frequency+phase));
case3
y=amp*sawtooth(2*pi*x*frequency+phase,0.5);
case4
y=amp*sawtooth(2*pi*x*frequency+phase);
case5
y=amp*(2*rand(size(x))-1);
otherwise
errordlg('Illegalwavetype','Chooseerrer');
end
ifget(handles.add,'Value')==0.0
handles.y=y;
else
handles.y=handles.y+y;
end
handles.inputtype=3;
guidata(hObject,handles);
plot(handles.time,handles.y);
title('WAVE');
axis([0N-str2double(get(handles.amp,'String'))str2double(get(handles.amp,'String'))]);
正弦波:
方波:
三角波:
钜齿波:
白噪声:
分析界面:
(5)时域分析相关程序及结果:
functiontimeanalyse_Callback(hObject,eventdata,handles)
Fs=str2double(get(findobj('Tag','samplerate'),'String'));
N=str2double(get(findobj('Tag','samplenum'),'String'));
ifhandles.inputtype==0
msgbox('Nowaveexist!
Pleasechooseainputtype!
');
return;
end
n=1;
ymax=max([handles.y
(1)handles.y
(2)]);
ymin=min([handles.y
(1)handles.y
(2)]);
from=str2double(get(handles.pointfrom,'String'));
to=str2double(get(handles.pointto,'String'));
iffrom<1|to-from<5;
msgbox('Errorrange!
');
return;
end
fori=from+2:
to-1;
ifhandles.y(i-1)<0&handles.y(i-2)<0&handles.y(i)>=0&handles.y(i+1)>0
ifhandles.y(i)==0
ti(n)=i;
else
ti(n)=i-handles.y(i)/(handles.y(i)-handles.y(i-1));
end
amp(n)=(ymax-ymin)/2;
ymax=0;
ymin=0;
n=n+1;
else
ifymaxymax=handles.y(i);
end
ifymin>handles.y(i)
ymin=handles.y(i);
end
end
end
n=n-1;
%freqenceandperiodicity
fori=1:
n-1
T(i)=ti(i+1)-ti(i);
end
freq=Fs/mean(T);
set(handles.outt,'String',1/freq);
set(handles.outfreq,'String',num2str(freq));
set(handles.outamp,'String',num2str(mean(amp(2:
n-1))));
phase=2*pi*(1-(ti(1:
n-1)-1)./T+floor((ti(1:
n-1)-1)./T));
set(handles.outphase,'String',num2str(mean(phase)));
set(handles.outpeak,'String',(max(handles.y(from:
to))-min(handles.y(from:
to)))/2);
set(handles.outmean,'String',mean(handles.y(from:
to)));
set(handles.outmeansquare,'String',mean(handles.y(from:
to).^2));
set(handles.outs,'String',std(handles.y(from:
to))^2);
以正弦波分析为例:
(6)频域分析相关程序及结果:
functionfreqanalyse_Callback(hObject,eventdata,handles)
Fs=str2double(get(findobj('Tag','samplerate'),'String'));
N=str2double(get(findobj('Tag','samplenum'),'String'));
ifhandles.inputtype==0
msgbox('Nowaveexist!
Pleasechooseainputtype!
');
return;
end
from=str2double(get(handles.pointfrom,'String'));
to=str2double(get(handles.pointto,'String'));
sample=handles.y(from:
to);
f=linspace(0,Fs/2,(to-from+1)/2);
Y=fft(sample,to-from+1);
[C,I]=max(abs(Y));
set(handles.foutt,'String',1/f(I));
set(handles.foutfreq,'String',f(I));
Y=Y(1:
(to-from+1)/2);
plot(handles.plot1,f,2*sqrt(Y.*conj(Y)));
plot(handles.plot2,f,angle(Y));
plot(handles.plot3,f,real(Y));
plot(handles.plot4,f,imag(Y));
plot(handles.plot5,f,abs(Y).^2);
xlabel(handles.plot1,'freqency(Hz)');
xlabel(handles.plot2,'freqency(Hz)');
xlabel(handles.plot3,'freqency(Hz)');
xlabel(handles.plot4,'freqency(Hz)');
xlabel(handles.plot5,'freqency(Hz)');
ylabel(handles.plot1,'amplitude');
ylabel(handles.plot2,'phase(rad)');
ylabel(handles.plot3,'real');
ylabel(handles.plot4,'Imaginary');
ylabel(handles.plot5,'power');
以正弦波为例:
(7)播放键程序:
functionplay_Callback(hObject,eventdata,handles)
Fs=str2double(get(findobj('Tag','samplerate'),'String'));
wavplay(handles.y,Fs);