【多重分形谱】基于MATLAB的图像多重分形谱分析仿真
FPGA教程目录
MATLAB教程目录
-----------------------------------------------------------------------
1.软件版本
matlab2017b
3.核心代码
function varargout = Tops(varargin)
% TOPS MATLAB code for Tops.fig
% TOPS, by itself, creates a new TOPS or raises the existing
% singleton*.
%
% H = TOPS returns the handle to a new TOPS or the handle to
% the existing singleton*.
%
% TOPS('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in TOPS.M with the given input arguments.
%
% TOPS('Property','Value',...) creates a new TOPS or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before Tops_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to Tops_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help Tops% Last Modified by GUIDE v2.5 04-Apr-2013 07:31:51% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...'gui_Singleton', gui_Singleton, ...'gui_OpeningFcn', @Tops_OpeningFcn, ...'gui_OutputFcn', @Tops_OutputFcn, ...'gui_LayoutFcn', [] , ...'gui_Callback', []);
if nargin && ischar(varargin{1})gui_State.gui_Callback = str2func(varargin{1});
endif nargout[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
elsegui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT% --- Executes just before Tops is made visible.
function Tops_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to Tops (see VARARGIN)% Choose default command line output for Tops
handles.output = hObject;% Update handles structure
guidata(hObject, handles);% UIWAIT makes Tops wait for user response (see UIRESUME)
% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line.
function varargout = Tops_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structure
varargout{1} = handles.output;% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
parameter;SIZE = str2num(get(handles.edit2,'string'));[filename,pathname]=uigetfile({'*.jpg';'*.bmp';'*.gif'},'选择图片');
if isequal(filename,0)disp('Users Selected Canceled');
else
str=[pathname filename];
I0 = imread(str);
I1 = imresize(I0,[SIZE,SIZE]);
axes(handles.axes1);
imshow(I1);
end;% --- Executes on button press in pushbutton2.
function pushbutton2_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
parameter;
I_gray = rgb2gray(I1);
axes(handles.axes2);
imshow(I_gray);
Inew = I_gray;% --- Executes on button press in pushbutton3.
function pushbutton3_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
parameter;
I_10 = im2bw(I1);
axes(handles.axes2);
imshow(I_10);
Inew = I_10;% --- Executes on button press in pushbutton4.
function pushbutton4_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton4 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
clc;
clear;
close all;function edit2_Callback(hObject, eventdata, handles)
% hObject handle to edit2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit2 as text
% str2double(get(hObject,'String')) returns contents of edit2 as a double% --- Executes during object creation, after setting all properties.
function edit2_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end% --- Executes on button press in pushbutton5.
function pushbutton5_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton5 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
parameter;wg = str2num(get(handles.edit3,'string'));
qs = str2num(get(handles.edit5,'string'));
qstep = str2num(get(handles.edit6,'string'));[Rs,Cs] = size(Inew);
npix = sum(sum(Inew));
%计算零即数量的像素在第i盒尺寸L
width = Rs;
p = log(width)/log(2);
MaxBox = power(Rs,2)/power(2,2);
NsL = double(zeros(MaxBox,p));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for Ig = p-1:-1:0siz = 2^(p-Ig);sizm1 = siz - 1;index = log2(siz);count = 0;for i=1:siz:(width-siz+1)for j=1:siz:(width-siz+1)count = count + 1;sums = sum(sum(Inew(i:i+sizm1,j:j+sizm1)));NsL(count,index) = sums;endend
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Qran = qs;
logl = zeros(p,1);
for i=1:plogl(i) = log(power(2,i));
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
PsL = double(zeros(MaxBox,p));
for j=1:pNo_boxes = power(Rs,2)/power(power(2,j),2);norm = sum(NsL(1:No_boxes,j));for i=1:No_boxesPsL(i,j) = NsL(i,j)/norm;end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%f(alpha),alpha
for kj=1:pcount = 0;No_boxes = power(Rs,2)/power(power(2,kj),2);for q = -Qran:qstep:Qran qsum = 0;for i=1:No_boxesif PsL(i,kj) ~= 0qsum = qsum + power(PsL(i,kj),q);endendfqnum = 0;aqnum = 0;smuiqL = 0;for i=1:No_boxesif PsL(i,kj) ~= 0muiqL = power(PsL(i,kj),q)/qsum;fqnum = fqnum + (muiqL * log(muiqL));aqnum = aqnum + (muiqL * log(PsL(i,kj)));smuiqL = smuiqL + muiqL;end endcount = count + 1;fql(kj,count) = fqnum;aql(kj,count) = aqnum;qval(count) = q;end
end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%alpha_q
for i=1:countline = polyfit(logl,aql(:,i),1);aq(i) = line(1);yfit = polyval(line,logl);sse = sum(power(aql(:,i)-yfit,2));sst = sum(power(aql(:,i)-mean(aql(:,i)),2));ar2(i) = 1-(sse/sst);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%f_q
for i=1:countline = polyfit(logl,fql(:,i),1);fq(i) = line(1);yfit = polyval(line,logl);sse = sum(power(fql(:,i)-yfit,2));sst = sum(power(fql(:,i)-mean(fql(:,i)),2));fr2(i) = 1-(sse/sst);
end
axes(handles.axes3);
plot(aq,fq,'b-o');
xlabel('alpha','FontSize',10);
ylabel('f(alpha)','FontSize',10);
pause(1);axes(handles.axes5);
plot(qval,fq,'b-o');
xlabel('q','FontSize',10);
ylabel('f(alpha)','FontSize',10);
pause(1);%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Qran = 20*qs;
for kj=1:pcount = 0;No_boxes = power(Rs,2)/power(power(2,kj),2);for q = -Qran:10*qstep:Qran qsum = 0;for i=1:No_boxesif PsL(i,kj) ~= 0qsum = qsum + power(PsL(i,kj),q);endendfqnum = 0;aqnum = 0;smuiqL = 0;for i=1:No_boxesif PsL(i,kj) ~= 0muiqL = power(PsL(i,kj),q)/qsum;fqnum = fqnum + (muiqL * log(muiqL));aqnum = aqnum + (muiqL * log(PsL(i,kj)));smuiqL = smuiqL + muiqL;end endcount = count + 1;fql(kj,count) = fqnum;aql(kj,count) = aqnum;qval(count) = q;end
end
for i=1:countline = polyfit(logl,aql(:,i),1);aq(i) = line(1);yfit = polyval(line,logl);sse = sum(power(aql(:,i)-yfit,2));sst = sum(power(aql(:,i)-mean(aql(:,i)),2));ar2(i) = 1-(sse/sst);
end
for i=1:countline = polyfit(logl,fql(:,i),1);fq(i) = line(1);yfit = polyval(line,logl);sse = sum(power(fql(:,i)-yfit,2));sst = sum(power(fql(:,i)-mean(fql(:,i)),2));fr2(i) = 1-(sse/sst);
end
for i=1:countDq(i) = (aq(i)*qval(i)-fq(i))/(qval(i)-1);
endaxes(handles.axes4);
plot(qval,Dq,'b-o');
xlabel('q','FontSize',10);
ylabel('D(q)','FontSize',10);
pause(1);axes(handles.axes6);
plot(aq,qval,'b-o');
xlabel('alpha','FontSize',10);
ylabel('q','FontSize',10);
pause(1);aqs = aq;
fqs = fq;
qvals = qval;
Dqs = Dq;function edit3_Callback(hObject, eventdata, handles)
% hObject handle to edit3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit3 as text
% str2double(get(hObject,'String')) returns contents of edit3 as a double% --- Executes during object creation, after setting all properties.
function edit3_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction edit5_Callback(hObject, eventdata, handles)
% hObject handle to edit5 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit5 as text
% str2double(get(hObject,'String')) returns contents of edit5 as a double% --- Executes during object creation, after setting all properties.
function edit5_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit5 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
endfunction edit6_Callback(hObject, eventdata, handles)
% hObject handle to edit6 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit6 as text
% str2double(get(hObject,'String')) returns contents of edit6 as a double% --- Executes during object creation, after setting all properties.
function edit6_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit6 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');
end% --- Executes on button press in pushbutton6.
function pushbutton6_Callback(hObject, eventdata, handles)
% hObject handle to pushbutton6 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
parameter;fid = fopen('save_alpha.txt','wt');
for i = 1:length(aqs)fprintf(fid,'%f\n',aqs(i));
end
fclose(fid);fid = fopen('save_falpha.txt','wt');
for i = 1:length(fqs)fprintf(fid,'%f\n',fqs(i));
end
fclose(fid);fid = fopen('save_Q.txt','wt');
for i = 1:length(qvals)fprintf(fid,'%f\n',qvals(i));
end
fclose(fid);fid = fopen('save_D.txt','wt');
for i = 1:length(Dqs)fprintf(fid,'%f\n',Dqs(i));
end
fclose(fid);
4.操作步骤与仿真结论
运行tops.m文件。会得到如下的结果:
我们选择任意一个图形进行测试:仿真效果如下所示:
A09-34
5.参考文献
[1]周孝华, 宋坤. 高频金融时间序列的异象特征分析及应用——基于多重分形谱及其参数的研究[J]. 财经研究, 2005, 31(7):10.
【多重分形谱】基于MATLAB的图像多重分形谱分析仿真相关推荐
- 【MATLAB教程案例27】基于matlab的图像配准算法的仿真与分析——sift,surf,kaze等
FPGA教程目录 MATLAB教程目录 目录 1.软件版本 2.图像配准理论概述 3.基于sift特征提取的图像配准
- 基于matlab的图像HOG特征提取算法仿真
目录 1.算法概述 2.仿真效果 3.MATLAB仿真源码 1.算法概述 方向梯度直方图(Histogram of Oriented Gradient, HOG)特征是一种在计算机视觉和图像处理中用来 ...
- 基于MATLAB的图像融合设计
基于MATLAB的图像融合设计 摘 要 图像融合能够将不同类型传感器获取的同一对象的图像数据进行空间配准.并且采用一定的算法将不同类型的传感器获取的同一对象的图像数据所含用的信息优势或互补性有机地结合 ...
- matlab用游程编码压缩图像,基于Matlab的图像压缩编码
开发与应用 计算机与信息技术 ·23· 基于 Matlab 的图像压缩编码 杨晓 李悦 (贵州大学 计算机与信息学院,贵州 贵阳 550025) 摘 要 本文描述了图像编码压缩方法的主要分类,介绍了每 ...
- 【图像修复】基于matlab损坏图像修复【含Matlab源码 731期】
一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[图像修复]基于matlab损坏图像修复[含Matlab源码 731期] 点击上面蓝色字体,直接付费下载,即可. 获取代码方式2: 付费专栏 ...
- 基于matlab的图像形状与分类毕业设计(含源文)
基于matlab的图像形状与分类 摘 要 数字图像处理是一门新兴技术,随着计算机硬件的发展,数字图像的实时处理已经成为可能,由于数字图像处理的各种算法的出现,使得其处理速度越来越快,能更好的为人们服务 ...
- matlab由图像得到复振幅程序,基于Matlab的空间滤波实验的计算机仿真
35 基于Matlab 的空间滤波实验的计算机仿真 张奇辉,王 洪,蓝发超 (华南理工大学物理科学与技术学院,广东 广州 510640) 摘 要:利用阿贝-波特实验装置和空间滤波系统,从改变频谱入手改 ...
- 基于MATLAB的图像去雾系统
基于MATLAB的图像去雾系统 (一)图像去雾基础原理 1.雾霭的形成机理 雾实际上是由悬浮颗粒在大气中的微小液滴构成的气溶胶,常呈现乳白色,其底部位于地球表面,所以也可以看作是接近地面的云.霭其实跟 ...
- 基于matlab的图像形状与分类的方法比较
基于matlab的图像形状与分类的方法比较 分类就是根据被识别对象的若干特征将其归入某一类别. 形状特征 (一)特点:各种基于形状特征的检索方法都可以比较有效地利用图像中感兴趣的目标来进行检索,但它们 ...
最新文章
- 第四周项目二-太乐了
- iPhone开发中的一些小技巧
- S3C2440的SPI控制器
- 微信小程序实战,用vue3实现每日浪漫情话推荐~
- Linux 网络配置 ifconfig、ip、netstat 命令的使用
- windows 截屏快捷键x220_电脑截屏快捷键是什么啊
- 苹果计算机系统是什么,苹果电脑系统和Win电脑系统有什么不同
- 使用大量数据提高算法性能的前提条件
- python监测网页变化_Python利用Last-Modified实现监控网页变化
- OpenCv之黑白上色小技巧,奥黛丽赫本神颜,超惊艳滴~
- 【C++】【记录自己的自学过程】用string简单的写一个拼字小游戏,只需要掌握string find方法即可。
- CS61C Spring 2021——Project 3: CS61CPU要求及实现思路
- .Net程序员的职业规划
- Verilog HDL的简单应用之74ls148功能的实现
- SOLIDWORKS 2022软件新增功能详细介绍
- E-NCAP增加黑暗场景下AEB测试,红外夜视或成主流配置
- 一男老师每日百词转载+连载(2)
- 消息被服务器拒绝访问,服务器上登录的Firebase访问被拒绝消息在哪里?
- 3D视觉——线激光测量胶水高度
- 插件推荐:json解析—Gson以及GsonFormat插件的运用
热门文章
- Ubuntu中C语言和C++的区别
- python-win32com-excel-分列,设置坐标轴格式
- 公司技术大咖分享会--后记
- MAC 系统下使用邮件客户端登录腾讯企业邮箱失败问题
- Java在歌唱比赛中,共有10位评委进行打分,在计算歌手得分时,去掉一个最高分,去掉一个最低分, 然后剩余的8位评委的分数进行平均,就是该选手的最终得分.输入每个评委的评分,求某选手的得分
- 3.19	外协加工(通过BOM体现加工物料总成本,非系统标准工序外协功能)
- Arch Linux 记录
- Cookie保存登录状态
- 基于SSM校园宿舍管理系统
- pandas数据分析之数据运算(逻辑运算、算术运算、统计运算、自定义运算)