bmiafig11.m

来自「beamlet变化的工具箱」· M 代码 · 共 85 行

function BMIAfig11(OpenNewWindow,HaveTitle,LoadData)

if nargin<2
    HaveTitle=1;
end
if nargin<1
    OpenNewWindow=1;
end
if OpenNewWindow
    figure
end
% BMIAfig11:
%clc; close all;
fprintf('\n');
disp('Figure 11')
disp('The *scale-thickness graph* of three different *sets* at various')
disp('scales. The left column gives the original curve(s). The right column')
disp('gives the corresponding thickness--scale plots, which portray')
disp('the sum of squares of all thicknesses at each given scale.')


fprintf('\n');
disp('Warning: It takes about 6 minutes (n=512) to create ');
disp('a figure.');

fprintf('\n');
disp('You can change the image size by modifying function *fig11draw*.');
disp('(Changing the value of *L*.)');
fprintf('\n');

fprintf('\n');
%k = input('Would you like to use pre-computed results? (0-no/1-yes)');
fprintf('\n');

if LoadData==0, 
  L = 9; n = 2^L; scale=(3:L);
  fig11_circle
  fig11_rw
  fig11_snowflake
else%if k==1, 
  disp('You need to download the file "fig11data.mat" on the web.')
  disp('To download, visit "http://www.isye.gatech.edu/~beamlab/data/".')
  disp('Install the data file at [Your BeamLab Directory]/Applications/Regularity/.')
  disp('Now waiting...');
  pause;
  disp('Loading the data...');
  load fig11data;
  % figure-1
  %figure(1); clf; 
  subplot(321); AutoImage(-img1); axis image; title('Circle'); 
  subplot(322); plot(scale,betaThick1); hold on; plot(scale,betaThick1,'d'); 
  axis([2 10 10 60]); axis square; title('Beta Thickness'); drawnow;
  % figure-2
  %figure(2); clf; 
  subplot(323); AutoImage(-img2); 
  axis image; title('Random Walk'); subplot(324); plot(scale,betaThick2); 
  hold on; plot(scale,betaThick2,'d'); 
  axis square; title('Beta Thickness'); drawnow;
  % figure-3
  %figure(3); 
  subplot(325); AutoImage(-img3); axis image; title('Snow Flake'); 
  subplot(326); plot(scale,betaThick3); hold on; plot(scale,betaThick3,'d'); 
  axis square; title('Beta Thickness'); drawnow;
end;


fprintf('\n');

disp('For more information, visit subdirectory [Applications/Regularity].');
disp('*** END Figure 11 ***');

%
% Part of BeamLab Version:200
% Built:Friday,23-Aug-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%
%
% Part of BeamLab Version:200
% Built:Saturday,14-Sep-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%