哈尔滨工程大学数字信号处理实验八.docVIP

实验原理 1、频率（周期）检测 对周期信号来说，可以用时域波形分析来确定信号的周期，也就是计算相邻的两个信号波峰的时间差或过零点的时间差。这里采用过零点（ti）的时间差T（周期）。频率即为f=1/T,由于能够求得多个T的值（ti有多个），故采用它们的平均值作为周期的估计值。 2、幅值检测 在一个周期内，求出信号最大值ymax与最小值ymin的差的一半，即A=（ymax-ymin）/2,同样，也会求出多个A值，但第一个A值对应的ymax和ymin不是在一个周期内搜索得到的，故以除以第一个以外的A值的平均作为幅值的估计值。 3、相位检测 采用过零法，即通过判断与同频零相位信号过零点时刻，计算其时间差，然后换成相应的相位差。φ=2п（１－ｔｉ／Ｔ），同样，以φ的平均值作为相位的估计值。 频率、幅值和相位估计的流程图见实验讲义64页图3.13。 4、数字信号统计量估计 （1）峰值p的估计 在样本数据x中找出最大值与最小值，其差值为双峰值，双峰值的一半即为峰值。 （2）均值估计 （3）均方值估计 (4) 差估估计 5、频谱分析原理 时域分析只能反映信号的幅值随时间的变化情况，除单频谱分量的简单波形外，很难明确提示信号的频率组成的各频率分量的大小，而频域分析能很好的解决此问题。 DFT与FFT 对于给定的时域信号y，可以通过Fourier变换得到频域信息Y。采样信号是一个不连续的频谱，不能算出所有点的值，所以采用离散Fourier变换DFT，但其计算效率较低，故实际中常用快速傅立叶变换。 频率、周期的估计 对于Y(k△f)，如果当k△f=f~时，Y(k△f)取最大值，则f为频率的估计值，由于采样间隔的误差，f~也存在误差，其误差最大为△f/2。周期T=1/f，从原理上可以看出，如果在标准信号中混有噪声，用上述方法仍能够精确地估计出原标准信号的频率和周期。 频谱图 为了直观地表示信号的频率特性，常常将傅立叶变换的结果用图形的方式表示，及频谱图。以频率f为横坐标，|Y(f)|为纵坐标，可以得到幅值谱；以频率f为横坐标，argY(f)为纵坐标，可以得到相位谱；以频率f为横坐标，ReY(f)为纵坐标，可以得到实频谱；以频率f为横坐标，ImY(f)为纵坐标，可以得到虚频谱。根据采样定理，只有频率不超过Fs/2的信号才能被正确采集，即傅立叶变换结果中频率大于Fs/2的部分是不正确的 ，故不在频谱中显示。 模块划分 模块化就是把程序划分成独立命名且可独立访问的模块，每个模块完成一个子功能，把这些模块集成起来构成一个整体，可以完成指定的功能满足用户需求。模块划分时应遵循如下规则：改进软件结构提高模块独立性；模块规模应该适中；深度，宽度，扇出都应适当；模块的作用域应该在控制域之内；力争降低接口的复杂程度；设计单入单出的模块口德模块；模块功能可以预测。对软件进行如讲义66页3.14图所示的模块划分。 参考代码如下： 声卡输入：“开始录音”按钮的回调函数如下 %首先获得设定的Fs值 Fs=str2double(get(findobj(Tag,samplerate),String)); %根据设定的录音实诚进行录音，将其存入handles.y中 handles.y=wavrecord(str2double(get(findobj(Tag,recordtime),String))*Fs,Fs,double); handles.inputtype=1; %保存handles结构体，使得handle.y在别的函数中也能使用 guidata(hObject,handles); %在波形显示区绘出波形 plot(handles.time,handles.y); title(WAVE); %将所采到的点的数量输出在“采样点数中” ysize=size(handles.y) set(handles.samplenum,String,num2str(ysize(1))); WAV文件输入：“打开文件”按钮的回调函数的部分内容 %从WAVw文件中读取的声音信息并临时存放到temp变量中 temp=wavread(get(findobj(Tag,filename),String)); %获得所选的声道 channel=str2double(get(handles.channel,String)); %将指定的声道的信息存放到handles.y中 handles.y=temp(:,channel); handles.inputtype=2; guidata(hObject,handles); plot(handles.time,handles.y); title(WAVE); ysize=size(handles