小波分析实验报告.docxVIP

/N-A.八十了A

指导老师：··革

姓名：·剑

班级：电信0906

学号：U200913743

小波分析实验目录

小波分析实验一

1.随机离散信号的傅里叶频谱分析

2.二维图像的傅里叶频谱分析

3.一维采样信号的傅里叶频谱分析

4.小波的分解与重构

(a)随机离散信号的分解与重构

(b)二维图像的分解与完全重构

(c)二维图像的分解与修改系数后的有损重构

小波分析实验二

1.基于小波的图像压缩

1.1图像压缩算法

(a)低频保留压缩

(b)部分细节压缩

(c)全局阈值压缩

(d)基于Birge-Massart策略中经验系数的分层阈值压缩

(e)基于小波包变换的图像压缩

(f)基于JPEG2K标准的图像压缩

1.2JPEG与JPEG2K的综合比较

1.3基于小波图像压缩技术的算法研究

(a)小波基的选择

(b)小波变换系数分析

2.基于小波的盲数字水印

2.1数字盲水印算法

2.2DCT与DWT盲水印综合比较

小波分析实验一

1、对随机生成的离散的一维信号，实现傅里叶变换，并显示原信号

及傅里叶频谱。

(1)实验图像

Figure1.(1)原始信号及其Fourier频谱(连续表示)

Figure1.(2)原始信号及其Fourier频谱(离散表示)

(2)实验代码

%对随机生成的离散的一维信号，实现傅里叶变换，并显示原信号及傅里叶频谱

clc

clear

8确定随机数最大值为1

MAXN=1;

S确定随机数总数

NUM=32;

s=rand(1,NUM)*MAXN;

subplot(2,1,1),

plot(s);

axis([ONUMOMAXN*1.5]);

xlabel(discreten);ylabel(magnituden);

title(原始信号(连续表示));

subplot(2,1,2),

Sf=abs(fft(s));

plot(Sf);

MAXY=max(Sf);

axis([ONUMOMAXY*1.5]);

xlabel(discreten);ylabel(magnitude);

title(Fourier频谱(连续表示));

figure,

subplot(2,1,1),

stem(s,*);

axis([ONUMOMAXN*1.5]);

xlabel(discreten);ylabel(magnituden);

title(原始信号(离散表示));

subplot(2,1,2),

stem(Sf,*);

MAXY=max(Sf)

axis([0ONUMOMAXY*1.5]);

xlabel(discreten);ylabel(magnitude);

title(Fourier频谱(离散表示));

(3)实验分析

代码实现对长度为32的随机一维行矩阵进行傅里叶变换，并显

示原始信号和频谱图。

2、对一幅离散二维图像(如512*512像素尺度的黑色背景叠加一个20*40像素尺寸的白色矩形，也可选择其他图像，如fingerprint等，可使用imread函数读取图像),实现傅里叶变换，并显示原图

及傅里叶频谱。

(1)实验图像

原始图像傅立叶频谱图像

原始图像

居中的傅立叶频谱图像使用对数变换的傅立叶频谱图像

Figure2.(1)原始图像及其Fourier频谱(一)

原始图像傅立叶频谱图像

原始图像

居中的傅立叶频谱图像使用对数变换的傅立叶频谱图像

Figure2.(2)原始图像及其Fourier频谱(二)

(2)实验代码

号对一幅离散二维图像实现傅里叶变换，并显示原图及傅里叶频谱

clc

clear

8fingerprint原图及频谱图

f=imread(fingerprint.bmp);

figure,

subplot(2,1,1),

imshow(f);title(原始图像!);

F=fft2(f):

S=abs(F):

S=gscale(S);

subplot(2,2,2).

imshow(S);title(傅里叶频谱图像!);

Fc=fftshift(F)

S=abs(Fc):

S=gscale(S):

subplot(2,2,3),

imshow(S);title(居中的傅里叶频谱图像);

S=fftshift(log(1+abs(F));

S=gscale(S);

subplot(2,2,4),

imsho