一. 简介

  离散傅立叶、离散余弦和离散小波变换是图像、音频信号常用基础操作,时域信号转换到不同变换域以后,会导致不同程度的能量集中,信息隐藏利用这个原理在变换域选择适当位置系数进行修改,嵌入信息,并确保图像、音频信号经处理后感官质量无明显变化。

二. 数学公式

一维离散傅立叶变换对定义

一维离散傅里叶变换:

一维离散傅里叶逆变换:

一维离散余弦变换对定义

一维离散余弦正变换:

一维离散余弦反变换:

一维连续小波变换对定义

一维连续小波变换,其总h(t)是小波母函数:

一维连续小波逆变换:

二维离散傅立叶变换对定义

二维离散傅里叶变换:

二维离散傅里叶逆变换:

二维离散余弦变换对定义

二维离散余弦正变换:

二维离散余弦反变换:

三. 代码实现

1. 用离散傅里叶变换分析合成音频和图像

1)分析合成音频文件

第一步:读取音频文件数据。

  uigetfile() 是文件对话框函数,提供图形界面供用户选择所需文件,返回目标的目录名和文件名。

  函数原型:y= wavread (FILE)

  功能:读取微软音频格式(wav)文件内容

  输入参数:file 表示音频文件名,字符串

  返回参数:y 表示音频样点,浮点型

第二步:一维离散傅立叶变换

  fft 函数对输入参数进行一维离散傅立叶变换并返回其系数,对应频率从 0 到 fs(采样频率),使用 fftshift 将零频对应系数移至中央。为了更好地观察频谱,可计算离散样点对应的频率值。

第三步:一维离散傅立叶逆变换

  ifft 函数对输入参数进行一维离散傅立叶逆变换并返回其系数。

第四步:观察结果

  figure(n)表示创建第 n 个图形窗。

  subplot 是子绘图函数,第一、二个参数指明子图像布局方式,例如,若参数为 2,3 则表示画面共分为 2 行,每行有 3 个子图像。第三个参数表明子图像序号,排序顺序为从左至右,从上至下。

  plot 是绘图函数,默认使用方式为 plot(y),参数 y 是要绘制的数据;如果需要指明图像横轴显示序列,则命令行为 plot(x, y),默认方式等同于 plot ([0..len-1], y),len为序列y的长度。

len=;
[fn,pn]=uigetfile('*.wav','请选择音频文件');
[x,fs]=wavread(strcat(pn,fn),len); %strcat()连接生成文件绝对路径 %一维离散傅里叶变换
xf=fft(x);
f1=[:len-]*fs/len;
xff=fftshift(xf);  % fftshift()将零频对应系数移至中央
h1=floor(len/);   %取不大于len/2的最大整数
f2=[-h1:h1]*fs/len; %一维离散傅里叶逆变换
xsync=ifft(xf);
figure;
subplot(,,);plot(x);title('原始图像');
subplot(,,);plot(xsync);title('synthesiaze audio');
subplot(,,);plot(f1,abs(xf));title('fft coef. of audio'); %abs()求整数的绝对值
subplot(,,);plot(f2(:len),abs(xff));title('fftshift coef. of audio');

 

(2)分析合成图像文件

第一步:读取图像文件数据

  函数原型:A = imread(filename,fmt)

  功能:读取 fmt 指定格式的图像文件内容

  输入参数:filename 表示图像文件名,字符串。Fmt 表示图像文件格式名,字符串,函数支持的图像格式包括:JPEG,TIFF,GIF,BMP 等等,当参数中不包括文件格式名时,函数尝试推断出文件格式。

  返回参数:A 表示图像数据内容,整型。

  rgb2gray 函数将 RGB 图像转换为灰度图。

第二步:二维离散傅立叶变换

  fft2 函数对输入参数进行二维离散傅立叶变换并返回其系数,使用 fftshift 将零频对应系数移至中央。

第三步:二维离散傅立叶逆变换

  ifft2 函数对输入参数进行二维离散傅立叶逆变换并返回其系数。

第四步:观察结果

  imshow 是二维数据绘图函数,mesh 通过三维平面显示数据。

例:

[fn,pn]=uigetfile('*.bmp','请选择图像文件');
[x,map]=imread(strcat(pn,fn),'bmp');
I=rgb2gray(x);  %rgb2gray函数将 RGB 图像转换为灰度图 %二维离散傅立叶变换
xf=fft2(I);
xff=fftshift(xf); %二维离散傅立叶逆变换
xsync=ifft2(xf); %结果
figure;
subplot(,,);imshow(x);title('original image');
subplot(,,);imshow(uint8(abs(xsync)));title('synthesize image');
subplot(,,);mesh(abs(xf));title('fft coef. of image');
subplot(,,);mesh(abs(xff));title('fftshift coef. of image');

2. 用离散余弦变换分析合成音频和图像

1)分析合成音频文件

第一步:读取音频文件数据。

第二步:一维离散余弦变换

  dct 函数对输入参数进行一维离散余弦变换并返回其系数,对应频率从0到 fs(采样频率)。

第三步:一维离散余弦逆变换

  idct 函数对输入参数进行一维离散余弦逆变换并返回其系数。离散余弦变换常用于图像压缩,可以尝试只使用部分系数重构语言,通过观察可发现,原始音频和合成后音频两者差别不大。

第四步:观察结果

例:

len=;
[fn,pn]=uigetfile('*.wav','请选择音频文件');
[x,fs]=wavread(strcat(pn,fn),len); %一维离散余弦变换
xf=dct(x);
f1=[:len-]*fs/len;
h1=floor(len/);
f2=[-h1:h1]*fs/len; %一维离散余弦逆变换
xsync=idct(xf);
[row,col]=size(x);
xff=zeros(row,col);
xff(:row,:col)=xf(:row,:col);
y=idct(xff); %结果
figure;
subplot(,,);plot(x);title('原始图像');
subplot(,,);plot(xsync);title('synthesiaze audio');
subplot(,,);plot(f1,abs(xf));title('fft coef. of audio');
subplot(,,);plot(f2(:len),abs(xff));title('fftshift coef. of audio');

(2)分析合成图像文件

第一步:读取图像文件数据

第二步:二维离散余弦变换

  dct 函数对输入参数进行二维离散余弦变换并返回其系数。

第三步:二维离散余弦逆变换

  idct2 函数对输入参数进行二维离散余弦逆变换并返回其系数。可以尝试使用部分系数重构图像,例如使用系数矩阵中4/5的数据,其它部分置零。

  为了保证图像能正确显示,使用uint8 对重构图像原始数据进行了数据类型转换,确保其取值范围在 0 到 255 之间。

第四步:观察结果

请输入命令显示四个子图,分别是原始图像、使用全部系数恢复的图像,使用部分系数恢复的图像和用三维立体图方式显示系数。

例:

[fn,pn]=uigetfile('*.bmp','请选择图像文件');
[x,map]=imread(strcat(pn,fn),'bmp');
I=rgb2gray(x); %二维离散余弦变换
xf=dct2(I); %二维离散余弦逆变换
xsync=uint8(idct2(xf));
[row,col]=size(I);
lenr=round(row*/);
lenc=round(col*/);
xff=zeros(row,col);
xff(:lenr,:lenc)=xf(:lenr,:lenc);
y=uint8(idct2(xff)); %结果
figure;
subplot(,,);imshow(x);title('original image');
subplot(,,);imshow(uint8(abs(xsync)));title('synthesize image');
subplot(,,);imshow(uint8(abs(y)));title('part synthesize image');
subplot(,,);mesh(abs(xff));title('fftshift coef. of image');

3. 用离散小波变换分析合成音频和图像

1)分析合成音频文件

第一步:读取音频文件数据。

第二步:一维离散小波变换

  wavedec 函数对输入参数进行一维离散小波变换并返回其系数C 和各级系数长度L。第二个参数指明小波变换的级数,第三个参数指明小波变换使用的小波基名称。

第三步:一维离散小波逆变换

  waverec 函数对输入参数进行一维离散小波逆变换并返回其系数。

  appcoef 函数返回小波系数近似分量,第一个参数 C、第二个参数 L 是 wavedec 的返回参数,为各级小波系数和其长度,第三个参数指明小波基名称,第四个参数指明级数。

  detcoef 函数返回小波系数细节分量,第一个参数 C、第二个参数 L 是 wavedec 的返回参数,为各级小波系数和其长度,第三个参数指明级数。

第四步:观察结果

例:

len=;
[fn,pn]=uigetfile('*.wav','请选择音频文件');
[x,fs]=wavread(strcat(pn,fn),len); %一维离散小波变换
[C,L]=wavedec(x,,'db4'); %一维离散小波逆变换
xsync=waverec(C,L,'db4');
cA2=appcoef(C,L,'db4',);
cD2=detcoef(C,L,);
cD1=detcoef(C,L,); %结果
figure;
subplot(,,);plot(x);title('原始图像');
subplot(,,);plot(xsync);title('synthesiaze audio');
subplot(,,);plot(cA2);title('app coef. of audio');
subplot(,,);plot(cD2);title('det coef. of audio');
subplot(,,);plot(cD1);title('det coef. of audio');

2)分析合成图像文件

第一步:读取图像文件数据

第二步:二维离散小波变换

  dwt2函数对输入参数进行二维一级离散小波变换并返回近似分量,水平细节分量,垂直细节分量和对角线细节分量。如果要对图像进行多级小波分解,使用wavedec2函数。

第三步:二维离散小波逆变换

  idwt2 函数对输入参数进行二维离散小波逆变换并返回其系数。可以尝试仅使用近似分量、水平细节分量、垂直细节分量或对角线细节分量重构图像。

第四步:观察结果

  输入命令显示六个子图,分别是原始图像、使用全部系数恢复的图像、小波系数近似分量、水平细节分量、垂直细节分量和对角线细节分量。

例:

[fn,pn]=uigetfile('*.bmp','请选择bmp格式图像文件');
[x,map]=imread(strcat(pn,fn),'bmp');
I=rgb2gray(x); %二维离散小波变换
sx=size(I);
[cA1,cH1,cV1,cD1]=dwt2(I,'bior3.7'); %二维离散小波逆变换
xsync=uint8(idwt2(cA1,cH1,cV1,cD1,'bior3.7',sx));
A1=uint8(idwt2(cA1,[],[],[],'bior3.7',sx));
H1=uint8(idwt2([],cH1,[],[],'bior3.7',sx));
V1=uint8(idwt2([],[],cV1,[],'bior3.7',sx));
D1=uint8(idwt2([],[],[],cD1,'bior3.7',sx)); %结果
figure;
subplot(,,);imshow(x);title('original image');
subplot(,,);imshow(uint8(abs(xsync)));title('synthesize image');
subplot(,,);mesh(A1);title('app coef. of image');
subplot(,,);mesh(H1);title('hor coef. of image');
subplot(,,);mesh(V1);title('ver coef. of image');
subplot(,,);mesh(D1);title('dia coef. of image');

 

Matlab信号处理基础的更多相关文章

  1. Matlab编程基础

    平台:Win7 64 bit,Matlab R2014a(8.3) “Matlab”是“Matrix Laboratory” 的缩写,中文“矩阵实验室”,是强大的数学工具.本文侧重于Matlab的编程 ...

  2. MATLAB图像处理基础

    MATLAB图像处理基础 2.2.1 图像文件格式及图像类型 1.MATLAB支持的几种图像文件格式: ⑴JPEG(Joint Photogyaphic Expeyts Group):一种称为联合图像 ...

  3. Matlab图像处理基础知识

    Matlab图像处理基础知识 Matlab的图片以矩阵的形式存储,矩阵的行列值为图片的行列的色彩值. 1图像表达方式: 像素索引 图像被视为离散单元.如使用I(2,2)可以获取第二行第二列的像素值 空 ...

  4. Matlab入门基础

    matlab入门篇,一些基础用法记一下~ M语言是解释型语言 ​ who:查看当前变量 ​ whoes:查看当前变量及其维数.所占字节数等. ​ clear: 清除所有变量 ​ clear + 变量名 ...

  5. matlab绘图基础

    matlab绘制条形图并分组显示: a =[1 2 3] b =[4 5 6] >> d=[a;b] d = 1 2 3 4 5 6 >> bar(d,'group') 修改横 ...

  6. matlab语言基础

    (一)常用命令 清空命令窗口:clc 清空工作空间:clear 打开editor:edit 查看函数的帮助信息:help 函数名  让命令窗口信息紧凑显示:format compact (二)变量与常 ...

  7. matlab图像基础知识

    1.MATLAB支持的几种图像文件格式: ⑴JPEG(Joint Photogyaphic Expeyts Group):一种称为联合图像专家组的图像压缩格式. ⑵BMP(Windows Bitmap ...

  8. Matlab绘图基础——图形修饰处理(入门)

    引入--标题.色条.坐标轴.图例等 例一: set(groot,'defaultAxesLineStyleOrder','remove','defaultAxesColorOrder','remove ...

  9. Matlab绘图基础——用print函数保存图片(Print figure or save to file)

    print(figure_handle,'formats','-rnumber','filename')  %将图形保存为png格式,分辨率为number的(默认为72),最好指定的分辨率大一点,否则 ...

随机推荐

  1. Django之auth组件

    Django自带的用户认证 我们在开发一个网站的时候,无可避免的需要设计实现网站的用户系统.此时我们需要实现包括用户注册.用户登录.用户认证.注销.修改密码等功能,这还真是个麻烦的事情呢. Djang ...

  2. BZOJ.1535.[POI2005]SZA-Template(KMP DP)

    BZOJ 洛谷 \(Description\) 给定一个字符串\(s\),求一个最短的字符串\(t\)满足,将\(t\)拼接多次后,可以得到\(s\).拼接是指,可以将\(t\)放在当前串的任意位置, ...

  3. [HNOI/AHOI2018]道路

    Description: W 国的交通呈一棵树的形状.W 国一共有\(n - 1\)个城市和\(n\)个乡村,其中城市从\(1\)到\(n - 1\) 编号,乡村从\(1\)到\(n\)编号,且\(1 ...

  4. Java基础-对象的内存分配与初始化(一定要明白的干货)

    首先,什么是类的加载?类的加载由类加载器执行.该步骤将查找字节码(classpath指定目录),并从这些字节码中创建一个Class对象.Java虚拟机为每种类型管理一个独一无二的Class对象.也就是 ...

  5. Math.round()

    在 JAVA 中四舍五入采用 Math.round(T a) 函数,函数返回的是一个 long 类型的长整型,参数 a 可以是 double 也可以是 float. 查看 JDK 源码: public ...

  6. js间隔一段时间打印数据库中的值

    https://blog.csdn.net/qq_33599109/article/details/78417474

  7. Jmeter接口测试实例2-获取所有学生信息

    Jmeter实例2:获取所有学生信息 添加http协议—添加IP.路径.方法,添加信息头管理器,察看结果树,运行 如下图所示,响应结果中获取到所有学生信息

  8. Idea-每次修改JS文件都需要重启Idea才能生效解决方法

    最近开始使用Idea,有些地方的确比eclipse方便.但是我发现工程每次修改JS或者是JSP页面后,并没有生效,每次修改都需要重启一次Tomcat这样的确不方便.我想Idea肯定有设置的方法,不可能 ...

  9. C_关于递归算法的几个例子

    1.递归算法的定义: 2.递归与迭代的优劣 eg1:斐波那契数列:斐波那契数列(Fibonacci sequence),又称黄金分割数列.因数学家列昂纳多·斐波那契(Leonardoda Fibona ...

  10. Zepto tap 穿透bug、解决移动端点击穿透问题

    当两个层重叠在一起时,或是有个弹窗,使用Zepto的tap事件时,点击上面的一层时会触发下面一层的事件,特别是底层如果是input框时,必“穿 透”,“google”说原因是“tap事件实际上是在冒泡 ...