第一章Matalab简介及信号处理基础.ppt

第二步：二维离散傅立叶变换 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函数对输入参数进行二维离散余弦变换并返回其系数。 * 第一章 Matalab简介及信号处理基础 1.1 Matalab简介 【实验目的】 【实验环境】 【原理简介】 【实验步骤】 【思考题】 【实验目的】 了解Matalab软件使用环境，熟悉Matalab使用方法，掌握信息隐藏实验中的各种常用Matalab命令。 【实验环境】 (1)WindowsXP或Vista操作系统 (2)MATLAB7.1版本软件 【原理简介】 Matalab是广泛使用的一种可视化科学计算软件，它不但具有语法结构简单、数值计算高效、图形功能完备和图像处理方便的特点，而且是信号处理和信息隐藏处理中使用最多的软件。 【实验步骤】 常用命令 变量设置 整数操作 随机序列常用命令 矩阵常用操作命令 位操作 绘图操作 文件操作 常用命令 dir：列出当前目录下的所有文件 clc：清除命令窗 clear all：清除环境（从内存中清除所有变量） who：将内存中的当前变量以简单形式列出 close all: 关闭所有的Figure窗口 变量设置 变量命名 变量名以字母打头，后最多可跟19个字母或数字，如x,y,ae3和d3er45等都是合法的变量名，不能使用内部函数或命令名作为变量名。MATLAB中的变量名区分大小，ab与 Ab表示两个不同的变量。 变量设置 变量赋值 表达式赋值 矩阵赋值 通过引用特定的位置可以单独改变某个矩阵元素 可以引用已定义的矩阵，重新定义一个新矩阵 整数操作 fix(x)：截尾取整 floor(x)：不超过x 的最大整数(高斯取整) ceil(x) ：大于x 的最小整数 随机序列常用命令 rand：均匀分布随机矩阵 randn:正态分布随机矩阵 randsrc：产生均匀分布数组 矩阵常用操作命令 矩阵的输入 直接输入创建矩阵 以“[”和“ ]”作为首尾，同行的元素用“，”或空格隔开，不同行的元素用“；”或按Enter键来分隔；矩阵的元素可以是数字也可以是表达式，如果