0001 function s=v_xyzticksi(ax,ah)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
0035
0036
0037
0038
0039
0040 if nargin<2
0041 ah=gca;
0042 if nargin<1
0043 ax=1;
0044 end
0045 end
0046 axfield={'XLim' 'YLim' 'ZLim'; 'XTick' 'YTick' 'ZTick'; 'XMinorTick' 'YMinorTick' 'ZMinorTick'; 'XTickLabel' 'YTickLabel' 'ZTickLabel'; 'XScale' 'YScale' 'ZScale'};
0047 tryglobal=nargout>0;
0048 digith=1;
0049 digitw=0.5;
0050
0051 prefix={'y','z','a','f','p','n','�','m','','k','M','G','T','P','E','Z','Y'};
0052 marg=[2 0.5 0.25 0.25];
0053 ntreq=[3 2 2 1];
0054
0055 lgridtem={1; [1 20 1 50 1]; [1 20 4]; 9; [2 20 8]; [5 20 2 30 7]; [10 20 5 30 4 50 5]};
0056 ngrid=length(lgridtem);
0057 lgrid=cell(ngrid,1);
0058 agrid=zeros(ngrid,1);
0059
0060 for i=1:ngrid
0061 igridtem=[lgridtem{i} 100];
0062 igrid=zeros(1,sum(igridtem(1:2:end)));
0063 ntem=length(igridtem)/2;
0064 k=0;
0065 tick0=10;
0066 for j=1:ntem
0067 nstep=igridtem(2*j-1);
0068 igrid(k+1:k+nstep)=tick0+(0:nstep-1)*(igridtem(2*j)-tick0)/igridtem(2*j-1);
0069 k=k+nstep;
0070 tick0=igridtem(2*j);
0071 end
0072 agrid(i)=sum(log10([igrid(2:end) 100]./igrid).^2);
0073 lgrid{i}=igrid';
0074 end
0075 minsubsp=1;
0076 delcheck=[log10(2) log10(5) 2];
0077 delval=[1 2 5];
0078 dosubtick=0;
0079
0080 ngrid=length(lgrid);
0081 loggrid=cell(ngrid,1);
0082 for i=1:ngrid
0083 loggrid{i}=log10(lgrid{i})-1;
0084 end
0085
0086 getgca=get(ah);
0087 set(ah,'Units','points','FontUnits','points');
0088 getgcac=get(ah);
0089 set(ah,'Units',getgca.Units,'FontUnits',getgca.FontUnits);
0090 if ax==1
0091 widthc=getgcac.Position(3)/getgcac.FontSize;
0092 axdir=[1 0];
0093 else
0094 widthc=2*getgcac.Position(4)/getgcac.FontSize;
0095 axdir=[0 1];
0096 end
0097 axdir=max(abs(axdir),1e-10);
0098 a=getgca.(axfield{1,ax})(1);
0099 b=getgca.(axfield{1,ax})(2);
0100
0101 ntick=0;
0102 tickint=[];
0103 tickdp=[];
0104 ticksi=[];
0105 subtick=[];
0106 if strcmp(getgca.(axfield{5,ax}),'log')
0107 width10=widthc/log10(b/a);
0108 ai3=3*ceil(log10(a)/3);
0109 bi3=3*floor(log10(b)/3);
0110 if ai3>=-24 && bi3<=24
0111
0112 if tryglobal && a>=10^(bi3-1)
0113 gi=bi3;
0114 s=prefix{9+gi/3};
0115 globalsi=1;
0116 else
0117 gi=0;
0118 globalsi=0;
0119 s='';
0120 end
0121 g=10^gi;
0122 ag=a/g;
0123 bg=b/g;
0124 al=log10(ag);
0125 bl=log10(bg);
0126 ai=ceil(al);
0127 bi=floor(bl);
0128 ai3=3*ceil(ai/3);
0129 bi3=3*floor(bi/3);
0130 for ipan=1:4
0131
0132 margin=marg(ipan);
0133 incsi=3*ceil(min((2*digitw+margin)/axdir(1),(digith+margin)/axdir(2))/(3*width10));
0134 switch ipan
0135 case {1,2}
0136 ticksi=incsi*ceil(ai/incsi):incsi:incsi*floor(bi/incsi);
0137 case {3,4}
0138 ticksi=ai3:incsi:bi3;
0139 end
0140 ntick=length(ticksi);
0141 tickint=ones(1,ntick);
0142 tickdp=zeros(1,ntick);
0143 if width10>0.25
0144 ticki=ai:bi;
0145 subtick=10.^(ticki(ticki~=3*fix(ticki/3)));
0146 end
0147 if incsi==3
0148 ci=floor(al);
0149 cibi=ci:bi;
0150 ndec=bi-ci+1;
0151 if globalsi
0152 siq0=zeros(1,ndec);
0153 else
0154 siq0=3*floor((cibi)/3);
0155 end
0156 siw0=siq0~=0;
0157 dpq0=max(siq0-cibi+1,1);
0158 for jgrid=1:ngrid
0159 igrid=jgrid-(ipan<=2)*(2*jgrid-ngrid-1);
0160 lgridi=lgrid{igrid};
0161 ngridi=length(lgridi);
0162 intq=reshape(repmat(lgridi,1,ndec).*repmat(10.^(cibi+dpq0-siq0-1),ngridi,1),1,[]);
0163 dpq=reshape(repmat(dpq0,ngridi,1),1,[]);
0164 msk=dpq>0 & rem(intq,10)==0;
0165 intq(msk)=intq(msk)/10;
0166 dpq(msk)=dpq(msk)-1;
0167 widq=1+floor(log10(intq));
0168 widq=digitw*(widq+(dpq>0).*max(1,dpq+2-widq)+reshape(repmat(siw0,ngridi,1),1,[]));
0169 logvq=reshape(repmat(loggrid{igrid},1,ndec)+repmat(ci:ndec+ci-1,ngridi,1),1,[]);
0170
0171 msk=logvq>=al & logvq<=bl;
0172 widq=widq(msk);
0173 logvq=logvq(msk);
0174
0175 if numel(widq)>=ntreq(ipan) && all(min((widq(1:end-1)+widq(2:end)+2*margin)/axdir(1),2*(digith+margin)/axdir(2))<=2*width10*(logvq(2:end)-logvq(1:end-1)))
0176
0177 ntick=numel(widq);
0178 tickint=intq(msk);
0179 tickdp=dpq(msk);
0180 siq=reshape(repmat(siq0,ngridi,1),1,[]);
0181 ticksi=siq(msk);
0182 subtick=[];
0183 dosubtick=igrid>1;
0184 break;
0185 end
0186 end
0187
0188 if ntick<5
0189 ldeltamin=log10(bg- bg*10^(-min((digitw+margin)/axdir(1),(digith+margin)/axdir(2))/width10));
0190 ildelta=floor(ldeltamin);
0191 ix=find(ldeltamin-ildelta<=delcheck,1);
0192 jx=ildelta*3+ix;
0193 while 1
0194 deltax=floor(jx/3);
0195 deltav=delval(jx-3*deltax+1);
0196 delta=deltav*10^deltax;
0197 multq=ceil(ag/delta):floor(bg/delta);
0198 ntickq=numel(multq);
0199 if ntickq<=ntick || ntickq<ntreq(ipan)
0200 break;
0201 end
0202 intq=deltav*multq;
0203 lintq=floor(log10(intq));
0204 siq=3*floor((lintq+deltax)/3);
0205 dpq=siq-deltax;
0206 msk=dpq<0;
0207 intq(msk)=intq(msk).*10.^(-dpq(msk));
0208 dpq(msk)=0;
0209 msk=rem(intq,10)==0 & dpq>0;
0210 while any(msk)
0211 dpq(msk)=dpq(msk)-1;
0212 intq(msk)=intq(msk)/10;
0213 msk=rem(intq,10)==0 & dpq>0;
0214 end
0215 widq=1+floor(log10(intq));
0216 widq=digitw*(widq+(dpq>0).*max(1,dpq+2-widq)+(siq~=0));
0217 logvq=log10(multq)+log10(deltav)+deltax;
0218
0219 if all(min((widq(1:end-1)+widq(2:end)+2*margin)/axdir(1),2*(digith+margin)/axdir(2))<=2*width10*(logvq(2:end)-logvq(1:end-1)))
0220 ntick=ntickq;
0221 tickint=intq;
0222 tickdp=dpq;
0223 ticksi=siq;
0224 dosubtick=1;
0225 break
0226 end
0227 jx=jx+1;
0228 end
0229 end
0230 end
0231 if ntick>=ntreq(ipan)
0232 break
0233 end
0234 end
0235 end
0236 if ntick
0237 tickint(tickint==0)=0;
0238 tickexp=gi+ticksi-tickdp;
0239 tickpos=tickint .* 10.^tickexp;
0240 ratthresh=10^(minsubsp/width10);
0241 if dosubtick
0242 subtick=[];
0243 if ntick>1
0244 stepexp=min(tickexp(1:end-1),tickexp(2:end))-1;
0245 stepint=round((tickpos(2:end)-tickpos(1:end-1)).*10.^(-stepexp));
0246 stepleft=tickint(1:end-1).*10.^(tickexp(1:end-1)-stepexp);
0247 subbase=10.^ceil(log10(stepint)-1);
0248 substep=[-1 -3 5]*((1+[1; 2; 5]*(subbase./stepleft))>ratthresh);
0249 substep(stepint~=10*substep)=max(2-substep(stepint~=10*substep),0);
0250 substep=substep.*subbase;
0251 for i=1:ntick-1
0252 ss=substep(i);
0253 sl=stepleft(i);
0254 if ss
0255 subtick=[subtick (sl+(ss:ss:stepint(i)-ss))*10^stepexp(i)];
0256 if i==1 && sl/(sl-ss)>ratthresh
0257 subtick=[subtick (sl-(ss:ss:floor((tickpos(1)-a)/(ss*10^stepexp(i)))*ss))*10^stepexp(i)];
0258 elseif i==ntick-1 && (1+ss/(sl+stepint(1)))>ratthresh
0259 subtick=[subtick (sl+stepint(i)+(ss:ss:floor((b-tickpos(end))/(ss*10^stepexp(i)))*ss))*10^stepexp(i)];
0260 end
0261 end
0262 end
0263 end
0264 end
0265 [tps,ix]=sort([tickpos subtick]);
0266 nticks=length(tps);
0267 ticklab=cell(nticks,1);
0268 for j=1:nticks
0269 i=ix(j);
0270 if i>ntick
0271 ticklab{j}='';
0272 else
0273 ticklab{j}=sprintf(sprintf('%%.%df%%s',tickdp(i)),tickint(i)*10^(-tickdp(i)),prefix{ticksi(i)/3+9});
0274 end
0275 end
0276 if width10<2.5
0277 set(ah,axfield{3,ax},'off');
0278 end
0279 set(ah,axfield{2,ax},tps);
0280 set(ah,axfield{4,ax},ticklab);
0281 end
0282
0283 else
0284 for ipan=1:4
0285 margin=marg(ipan);
0286
0287 if tryglobal
0288 gi=3*floor(log10(max(abs(a),abs(b)))/3);
0289 s=prefix{9+gi/3};
0290
0291 else
0292 gi=0;
0293 s='';
0294 end
0295 g=10^gi;
0296 ag=a/g;
0297 bg=b/g;
0298 width1=widthc/(bg-ag);
0299 ldeltamin=log10(min((digitw+margin)/axdir(1),(digith+margin)/axdir(2))/width1);
0300 ildelta=floor(ldeltamin);
0301 ix=find(ldeltamin-ildelta<=delcheck,1);
0302 jx=ildelta*3+ix;
0303 while 1
0304 deltax=floor(jx/3);
0305 deltav=delval(jx-3*deltax+1);
0306 delta=deltav*10^deltax;
0307 multq=ceil(ag/delta):floor(bg/delta);
0308 ntickq=numel(multq);
0309 if ntickq<ntreq(ipan)
0310 break;
0311 end
0312 intq=deltav*multq;
0313 lintq=floor(log10(abs(intq)+(intq==0)));
0314 siq=3*floor((lintq+deltax)/3)*~tryglobal;
0315 dpq=siq-deltax;
0316 msk=dpq<0;
0317 intq(msk)=intq(msk).*10.^(-dpq(msk));
0318 dpq(msk)=0;
0319 msk=rem(intq,10)==0 & dpq>0;
0320 while any(msk)
0321 dpq(msk)=dpq(msk)-1;
0322 intq(msk)=intq(msk)/10;
0323 msk=rem(intq,10)==0 & dpq>0;
0324 end
0325 widq=1+floor(log10(abs(intq)+(intq==0)));
0326 widq=digitw*(widq+(dpq>0).*max(1,dpq+2-widq)+(siq~=0).*(intq~=0)+(intq<0));
0327
0328 if all((widq(1:end-1)+widq(2:end)+2*margin)<=2*axdir(1)*width1*delta) || (digith+margin<=axdir(2)*width1*delta);
0329 ntick=ntickq;
0330 tickint=intq;
0331 tickdp=dpq;
0332 ticksi=siq;
0333 if deltav>1 && width1*delta>0.5*deltav
0334 mults=ceil(ag*deltav/delta):floor(bg*deltav/delta);
0335 subtick=(mults(deltav*fix(mults/deltav)~=mults))*delta/deltav;
0336 else
0337 subtick=[];
0338 end
0339 break
0340 end
0341 jx=jx+1;
0342 end
0343 if ntick>=ntreq(ipan)
0344 break
0345 end
0346 end
0347 if ntick
0348 tickint(tickint==0)=0;
0349 tickexp=gi+ticksi-tickdp;
0350 tickpos=tickint .* 10.^tickexp;
0351 [tps,ix]=sort([tickpos subtick*10^gi]);
0352 nticks=length(tps);
0353 ticklab=cell(nticks,1);
0354 for j=1:nticks
0355 i=ix(j);
0356 if i>ntick
0357 ticklab{j}='';
0358 else
0359 if tickint(i)==0
0360 ticklab{j}=sprintf(sprintf('%%.%df%%s',tickdp(i)),0,prefix{(ticksi(i)/3)*(tickint(i)~=0)+9});
0361 else
0362 ticklab{j}=sprintf(sprintf('%%.%df%%s',tickdp(i)),tickint(i)*10^(-tickdp(i)),prefix{(ticksi(i)/3)*(tickint(i)~=0)+9});
0363 end
0364 end
0365 end
0366 set(ah,axfield{2,ax},tps);
0367 set(ah,axfield{4,ax},ticklab);
0368 set(ah,axfield{3,ax},'on');
0369 end
0370 end