第7章 分形维数应用与MATLAB实现.pptVIP

第七章 MATLAB优化算法案例分析与应用 第7章 分形维数应用与MATLAB实现 第七章 MATLAB优化算法案例分析与应用 被誉为大自然的几何学的分形（Fractal）理论，是现代数学的一个新分支，但其本质却是一种新的世界观和方法论。分形维数反映了复杂形体占有空间的有效性，它是复杂形体不规则性的量度。分形理论在现在图像处理和信号分析处理领域应用越来越广泛，本章主要借助分形维数理论，对二维图像和语音信号进行分析计算，让读者真正掌握分形盒维数的计算。 第七章 MATLAB优化算法案例分析与应用 7.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信息维数。 分形盒维数应用较广泛，在用数字图像盒维数法求得分维值时增大图像的大小可以降低分维值计算的误差。信号时域短时分形盒维数进行低信噪比的带噪信号的计算机仿真表明，该方法能较准确地检测低信噪比下的语音端点，并且其算法也相对简单。 第七章 MATLAB优化算法案例分析与应用 7.2 二维图像分形盒维数分析 图7- 1 大树分形计算 第七章 MATLAB优化算法案例分析与应用 Rmax=sqrt((M-1)^2+(N-1)^2); %求最大距离 Nr=zeros(1,floor(Rmax)); for k=1:floor(Rmax) for i=1:M for j=1:N for m=1:M for n=1:N if k==1 TotalGary=TotalGary+double(abs(I(i,j)-I(m,n))); end if k=sqrt((i-m)^2+(j-n)^2)sqrt((i-m)^2+(j-n)^2)(k+1) Nr(k)=Nr(k)+1; end end; end; end; end; end; k=[1:floor(Rmax)]; E=2.*TotalGary.*ones(1,floor(Rmax))./Nr(1,floor(Rmax)); [P ,s]=polyfit(log(k),log(E),1); D=3-P(1); 第七章 MATLAB优化算法案例分析与应用 7.3 基于短时分形维数的语音信号检测 图7- 2 信号为正方形网格所覆盖 第七章 MATLAB优化算法案例分析与应用 分形维数计算如式（7.2）所示。 短时分形盒维数 % 求一维信号分形维数 clc,clear,close all % 清屏、清工作区、关闭窗口 warning off % 消除警告 feature jit off % 加速代码执行 data = csvread(lod78.csv); samplerate = 365; % 采样率 freqsol = 400; % 频率分辨率 timesol = 8