grafica.m

地震资料处理的matlab源码软件包

function grafica

global deltat;
global A;
global ngraf;
global grafArange;
global numero;
global tmax;
global Tmax;
global Uidens;
global Uiori;
global archivo;
global letras;
global rielaxes;
global grafaxes;
global textaxes;
global trange;
global Limiteizq;
global Limiteder;

color=['y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w' 'y' 'g' 'c' 'r' 'b' 'm' 'w'];
nplot=length(ngraf);

valden = get(Uidens,'Value');if valden < 0.1,valden=0.1,end;
valori = get(Uiori,'Value');if valori > 0.99,valori=0.99,end;

if Tmax==0.0;
   Tmax=tmax;
end;

xo     =  valori*Tmax;
xf     =  xo+Tmax*(valden)^2;if xf > Tmax;xf=Tmax;end;
trange =  [xo xf]; 

axes(rielaxes);
fill([0 0 xo xo 0], [0 1 1 0 0], 'w', [xo xo xf xf xo], [0 1 1 0 0], 'b', [xf xf Tmax Tmax xf], [0 1 1 0 0], 'w');
axis('off');

axes(grafaxes);
j=0;
plot([0 Tmax],[0 0],'w');hold on;
for i=ngraf;
    j=j+1;
    plot([0:deltat(i):deltat(i)*(numero(i+1)-numero(i)-1)],A([numero(i):numero(i+1)-1]),color(i));
end;  

j=0;
for i=ngraf;
    xtext = 0.05+j*0.25;
    text('Position',[xtext 1.1],'Units','normalized','Color',color(i),'String',archivo(i,1:letras(i)));
    j=j+1;
end;

if Limiteizq~=Limiteder;
   plot([Limiteizq Limiteizq], grafArange, 'r-');
   plot([Limiteder Limiteder], grafArange, 'r-');
end;

xlabel('tiempo (seg)');
if archivo(ngraf(1),letras(ngraf(1))-2:letras(ngraf(1)))=='vel';
   ylabel('velocidad (cm/seg)');
elseif archivo(ngraf(1),letras(ngraf(1))-2:letras(ngraf(1)))=='des';
   ylabel('desplazamiento (cm)');
else
   ylabel('aceleraci髇 (gales)');
end;

axis([trange grafArange]);
hold off;