现代信息处理技术实验报告.docxVIP

目 录 实验一 短时傅里叶变化与小波变换 ……………………………………………1 实验目的 ………………………………………………………………1 实验内容 ………………………………………………………………1 短时傅里叶变换 …………………………………………………1 小波变换 …………………………………………………………14 实验二 图像处理 …………………………………………………………………27 实验目的 ………………………………………………………………27 实验内容 ………………………………………………………………27 PAGE 36 实验一 短时傅里叶变换与小波变换 实验目的 熟悉并掌握短时傅里叶变换的性质、参数以及不同信号的短时傅里叶变换； 熟悉并掌握小波变换的性质、参数以及不同信号的小波变换。 实验内容 短时傅里叶变换 Matlab中的短时傅里叶变换函数spectrogram S = spectrogram(x) S = spectrogram(x,window) S = spectrogram(x,window,noverlap) S = spectrogram(x,window,noverlap,nfft) S = spectrogram(x,window,noverlap,nfft,fs) 调用及参数描述： window is a Hamming window of length nfft. noverlap is the number of samples that each segment overlaps. The default. value is the number producing 50% overlap between segments. nfft is the FFT length and is the maximum of 256 or the next power of 2 greater than the length of each segment of x. Instead of nfft, you can specify a vector of frequencies, F. See below for more information. fs is the sampling frequency, which defaults to normalized frequency 短时傅里叶变换 正弦信号 生成信号长度为1s、采样频率为1kHz、周期分别为0.1s、1s和10s的正弦信号s，并画出这些正弦信号。 MATLAB程序如下： 运行结果如下： 用spectrogram画出这些正弦信号的短时傅里叶变换。 spectrogram(s,hamming(256),255,256,1000); Matlab程序如下： 运行结果如下： 窗口的影响 针对周期为0.1秒的正弦函数，分别调整hamming窗口大小为32、64、128、256，并画出该正弦信号的短时傅里叶变换。 MATLAB程序如下： 运行结果如下： 分析 短时傅立叶变换基本思想是给信号加滑动时间窗，并对窗内信号做傅立叶变换，得到信号的时变频谱。在短时傅里叶变换过程中，窗的长度决定频谱图的时间分辨率和频率分辨率，窗宽越大，截取的信号越长，信号越长，傅里叶变换后频率分辨率越高（能看到频谱的快变化），时间分辨率越差。也即在实际变换中，时间分辨率和频率分辨率之间不能兼得。 针对周期为0.1秒的正弦函数，窗口大小为128，分别调整窗口类型为hamming、rectwin和blackman，并画出该正弦信号的短时傅里叶变换。 MATLAB程序如下： 运行结果如下： 分析 一个窗是否合适：窗谱主瓣宽度就尽可能的窄，且能量集中在主瓣内，以获得较陡的过渡带；窗谱旁瓣与主瓣相比应尽可能的小，旁瓣能量衰减要快，以利于增加阻带衰耗。 Hamming窗在频率范围中的分辨率较高，旁瓣衰减较大，主瓣峰值与第一个旁瓣峰值衰减可达40db，频谱泄露少。频谱中高频分量弱、波动小，得到较平滑的谱。 Rectwin窗导致变换中带进了高频干扰和泄漏，甚至出现负谱现象。频率识别精度最高，幅值识别精度最低。 Blackman窗主瓣宽，旁瓣小，频率识别精度最低，但幅值识别精度最高。 不同信号的短时傅里叶变换 例四中的离散信号，其前500点是