一、简介
被誉为大自然的几何学的分形（Fractal）理论，是现代数学的一个新分支，但其本质却是一种新的世界观和方法论。分形维数反映了复杂形体占有空间的有效性，它是复杂形体不规则性的量度。分形理论在现在图像处理和信号分析处理领域应用越来越广泛，本章主要借助分形维数理论，对二维图像和语音信号进行分析计算，让读者真正掌握分形盒维数的计算。
1 分形盒维数概述
分形盒维数分为以下几类。
（1）相似维数：若某图形是由把全体缩小成1/a的b个相似形所组成，由于 ，则有 。
（2）Kolmogorov容量维数：用半径为 的d维球包覆某d维图形集合时，假定 是球的个数的最小值。容量维数 可用下式来定义：
（3）盒子维数（box-counting dimension）：在双对数坐标纸上画出 对 的曲线，其直线部分的斜率就是此分形对象的盒子孙维数 。 是小盒子的边长， 为盒子数。
（4）信息维数（information dimension）：把小盒子编号，如果知道分形中的点落入第i只盒子的概率是 ，定义“信息维数” ， 。
（5）关联维数（correlation dimension）：如果把在空间随机分布的某量坐标X处中的密度记为 ，则关联函数 ， 表示平均。可以是全体平均，也可以是空间平均。1983年，P. Grassberger和J. Procassia给出了关联维数的定义 。
（6）广义维数：H. G. E. Hentschel等提出了广义维数的概念，其定义是： ，其中 是 阶Renyi信息， 叫作q阶广义维数，有时也叫Renyi信息维数。
分形盒维数应用较广泛，在用数字图像盒维数法求得分维值时增大图像的大小可以降低分维值计算的误差。信号时域短时分形盒维数进行低信噪比的带噪信号的计算机仿真表明，该方法能较准确地检测低信噪比下的语音端点，并且其算法也相对简单。
% 求一维信号分形维数
clc % 清屏
clear all; % 删除workplace变量
close all; % 关掉显示图形窗口
data = csvread('lod78.csv');
samplerate = 365; % 采样率
freqsol = 400; % 频率分辨率
timesol = 800; % 时间分辨率
df = dbox(data,samplerate); % 分形维数
disp(['分形维数为: ',num2str(df)])
%%
lc;clear;close all
[FileName,PathName] = uigetfile({'*.wav'},'Load Wav File');
[x,fs] = wavread([PathName '/' FileName]);
handlesx = x;
handlesfs = fs;
time = 0:1/fs:(length(handlesx)-1)/fs;
figure,plot(time,handlesx(:,1));
figure,plot(time,handlesx(:,2));
title('Original Signal');
figure,specgram(handlesx(:,1), 1024, handlesfs);
figure,specgram(handlesx(:,2), 1024, handlesfs);
title('Spectrogram of Original Signal');
clc % 清屏
clear all; % 删除workplace变量
close all; % 关掉显示图形窗口
% Read the image
[sFileName,sPath]=uigetfile('*.*','Load Image');
sFullFileName=[sPath,sFileName];
im = imread(sFullFileName);
im1 =imresize(im,[32 32]); % 图像压缩到尺寸到32*32，尺寸太大，PC运算太慢
function df=dbox(s,Fs)
% df: 网格维数，
% s: 一维信号
% Fs: 采样频率
d=length(s);
da1=0;
da2=0;
if d/2==0
s=[s,s(d)]; % 补成奇数个,方便第二个for循环
d=d+1;
end
二、源代码
三、运行结果
四、备注
完整代码或者代写添加QQ1575304183