第一章节_Matalab简介及信号处理基础幻灯片.ppt

分析合成图像文件包括以下步骤： 读取图像文件数据 二维离散傅立叶变换 二维离散傅立叶逆变换 观察结果 读取fmt指定格式的图像文件内容 [fn, pn] = uigetfile(*.bmp, 请选择图像文件); [x, map] = imread(strcat(pn, fn), bmp); I = rgb2gray(x); 第二步：二维离散傅立叶变换 xf = fft2(I); xff = fftshift(xf); fft2函数对输入参数进行二维离散傅立叶变换并返回其系数，使用fftshift将零频对应系数移至中央。 第三步：二维离散傅立叶逆变换 xsync = ifft2(xf); ifft2函数对输入参数进行二维离散傅立叶逆变换并返回其系数。 第四步：观察结果 figure; subplot(2, 2, 1); imshow(x);title(original image); subplot(2, 2, 2); imshow(uint8(abs(xsync)));title(synthesize image); subplot(2, 2, 3); mesh(abs(xf));title(fft coef. of image); subplot(2, 2, 4); mesh(abs(xff));title(fftshift coef. of image); imshow是二维数据绘图函数，mesh通过三维平面显示数据。 用离散余弦变换分析合成音频和图像 分析合成音频文件包括以下步骤： 读取音频文件数据 一维离散余弦变换 一维离散余弦逆变换 观察结果 第一步：一维离散余弦变换 xf = dct(x); dct函数对输入参数进行一维离散余弦变换并返回其系数，对应频率从0到fs（采样频率）。 第二步：一维离散余弦逆变换 xsync = idct(xf); [row,col]=size(x); xff=zeros(row,col); xff(1:row,1:col)=xf(1:row,1:col); y=idct(xff); idct函数对输入参数进行一维离散余弦逆变换并返回其系数。离散余弦变换常用于图像压缩，可以尝试只使用部分系数重构语言，通过观察可发现，原始音频和合成后音频两者差别不大。 第三步：观察结果 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); 分析合成图像文件包括以下步骤： 读取图像文件数据 二维离散余弦变换 二维离散余弦逆变换 观察结果 第一步：二维离散余弦变换 xf = dct2(I); dct2函数对输入参数进行二维离散余弦变换并返回其系数。 第二步：二维离散余弦逆变换 xsync = uint8(idct2(xf)); [row, col] = size(I); lenr = round(row * 4 / 5); lenc = round(col * 4 / 5); xff = zeros(row, col); xff(1:lenr, 1:lenc) = xf(1:lenr, 1:lenc); y = uint8(idct2(xff)); 第三步：观察结果 输入命令显示四个子图，分别是原始图像、使用全部系数恢复的图像，使用部分系数恢复的图像和用三维立体图方式显示系数。 figure; subplot(2, 2, 1); imshow(x);title(original image); subplot(2, 2, 2); imshow(uint8(abs(xsync)));title(synthesize image); subplot(2, 2, 3); imshow(uint8(abs(y)));title(part synthesize image); subplot(2, 2, 4); mesh(abs(xff));title(fftshift coef. of image); 用离散小波变换分析合成音频和图像 分析合成音频文件包括以下步骤： 读取音频文件数据 一维离散小波变换 一维离散小波逆变换 观察结果 第一步：一维离散小波变换 [C, L] = wavedec(x,2,db4); wavedec函数对输入参数进行一维离