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 if nargin<2
0040 ah=gca;
0041 if nargin<1
0042 ax=1;
0043 end
0044 end
0045 axfield={'XLim' 'YLim' 'ZLim'; 'XTick' 'YTick' 'ZTick'; 'XMinorTick' 'YMinorTick' 'ZMinorTick'; 'XTickLabel' 'YTickLabel' 'ZTickLabel'; 'XScale' 'YScale' 'ZScale'};
0046 tryglobal=nargout>0;
0047 digith=1;
0048 digitw=0.5;
0049
0050 prefix={'y','z','a','f','p','n','�','m','','k','M','G','T','P','E','Z','Y'};
0051 marg=[2 0.5 0.25 0.25];
0052 ntreq=[3 2 2 1];
0053
0054 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]};
0055 ngrid=length(lgridtem);
0056 lgrid=cell(ngrid,1);
0057 agrid=zeros(ngrid,1);
0058
0059 for i=1:ngrid
0060 igridtem=[lgridtem{i} 100];
0061 igrid=zeros(1,sum(igridtem(1:2:end)));
0062 ntem=length(igridtem)/2;
0063 k=0;
0064 tick0=10;
0065 for j=1:ntem
0066 nstep=igridtem(2*j-1);
0067 igrid(k+1:k+nstep)=tick0+(0:nstep-1)*(igridtem(2*j)-tick0)/igridtem(2*j-1);
0068 k=k+nstep;
0069 tick0=igridtem(2*j);
0070 end
0071 agrid(i)=sum(log10([igrid(2:end) 100]./igrid).^2);
0072 lgrid{i}=igrid';
0073 end
0074 minsubsp=1;
0075 delcheck=[log10(2) log10(5) 2];
0076 delval=[1 2 5];
0077 dosubtick=0;
0078
0079 ngrid=length(lgrid);
0080 loggrid=cell(ngrid,1);
0081 for i=1:ngrid
0082 loggrid{i}=log10(lgrid{i})-1;
0083 end
0084
0085 getgca=get(ah);
0086 set(ah,'Units','points','FontUnits','points');
0087 getgcac=get(ah);
0088 set(ah,'Units',getgca.Units,'FontUnits',getgca.FontUnits);
0089 if ax==1
0090 widthc=getgcac.Position(3)/getgcac.FontSize;
0091 axdir=[1 0];
0092 else
0093 widthc=2*getgcac.Position(4)/getgcac.FontSize;
0094 axdir=[0 1];
0095 end
0096 axdir=max(abs(axdir),1e-10);
0097 a=getgca.(axfield{1,ax})(1);
0098 b=getgca.(axfield{1,ax})(2);
0099
0100 ntick=0;
0101 tickint=[];
0102 tickdp=[];
0103 ticksi=[];
0104 subtick=[];
0105 if strcmp(getgca.(axfield{5,ax}),'log')
0106 width10=widthc/log10(b/a);
0107 ai3=3*ceil(log10(a)/3);
0108 bi3=3*floor(log10(b)/3);
0109 if ai3>=-24 && bi3<=24
0110
0111 if tryglobal && a>=10^(bi3-1)
0112 gi=bi3;
0113 s=prefix{9+gi/3};
0114 globalsi=1;
0115 else
0116 gi=0;
0117 globalsi=0;
0118 s='';
0119
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 ticklab{j}=sprintf(sprintf('%%.%df%%s',tickdp(i)),tickint(i)*10^(-tickdp(i)),prefix{(ticksi(i)/3)*(tickint(i)~=0)+9});
0360 end
0361 end
0362 set(ah,axfield{2,ax},tps);
0363 set(ah,axfield{4,ax},ticklab);
0364 set(ah,axfield{3,ax},'on');
0365 end
0366 end