信号处理基础
实验目的
了解音频和图像数据系数特点,掌握音频和图像文件的离散傅里叶,离散余弦和离散小波变换等基本条件
实验环境
(1)windows10
(2)Matlab2020
(3)bmp格式图像文件
(4)wav格式音频文件
1.用离散傅里叶变换分析合成音频和图像分析合成音频文件包括一下步骤:
(1)读取音频文件数据;
clc;
clear;
len=40000;
[fn,pn]=uigetfile('*.wav','请选择音频文件');
filename=strcat(pn,fn);
[x,fs]=audioread(filename,[1,len]);
注意:新版本的matlab中已经不再支持wavread()函数,替代函数为audioread(filename,N),其中N必须为[m,n]格式,如[2,100],且m,n均为正数
进行一维离散傅里叶变换:(其中fft()函数为一维离散傅里叶变换函数,fftshift函数将零频对应系数移至中央)
xf=fft(x);
f1=[0:len-1]*fs/len;
xff=fftshift(xf);
h1=floor(len/2);
f2=[-h1:h1]*fs/len;
进行一维离散傅里叶逆变换:
xsync=ifft(xf);
显示图像代码:
figure;
subplot(2,2,1);plot(x);title('original audio');
subplot(2,2,2);plot(xsync);title('synthesize audio');
subplot(2,2,3);plot(f1,abs(xf));title('fft coef. of audio');
subplot(2,2,4);plot(f2(1:len),abs(xff));title('fftshift coef. of auio');
结果: