location.m

我写的一个地震单台定位的程序

clf
global data;
%initinal the windows and kinds of parameters
set(gcf,'unit','normalized','position',[0.1 0.4 0.85 0.35]);
set(gcf,'defaultuicontrolunits','normalized');
set(gcf,'defaultuicontrolfontsize',11);
set(gcf,'defaultuicontrolfontname','隶书');
set(gcf,'menubar','none');
hmenu = uimenu('Label','Open');
hmenu1 = uimenu('Label','Method','Visible','off');
hsub1 = uimenu(hmenu1,'Label','Super_cycle');
hsub2 = uimenu(hmenu1,'Label','Azimath');
hsub3 = uimenu(hmenu1,'Label','Mkmod');
set(gcf,'defaultuicontrolhorizontal','left');
str = 'using two methods to location a seismology';
set(gcf,'name',str,'numbertitle','off');
haxes = axes('position',[0.05 0.15 0.45 0.60],'visible','on');
%place control components
htext = uicontrol(gcf,'Style','text','Position',[0.52 0.87 0.26 0.1],...
    'String','parameters input');
htext1 = uicontrol(gcf,'Style','text','Position',[.52 .77 .10 .1],'String','Start:');
hedit1 = uicontrol(gcf,'Style','edit','Position',[0.62 0.77 0.16 0.1]);
htext2 = uicontrol(gcf,'Style','text','Position',[.52 .67 .10 .1],'String','Standing:');
hedit2 = uicontrol(gcf,'Style','edit','Position',[0.62 0.67 0.16 0.1]);
htext3 = uicontrol(gcf,'Style','text','Position',[.52 .57 .10 .1],'String','Sample:');
hedit3 = uicontrol(gcf,'Style','edit','Position',[0.62 0.57 0.16 0.1]);
htext4 = uicontrol(gcf,'Style','text','Position',[.52 .47 .10 .1],'String','Alfa:');
hedit4 = uicontrol(gcf,'Style','edit','Position',[0.62 0.47 0.16 0.1]);
htext5 = uicontrol(gcf,'Style','text','Position',[.52 .37 .26 .1],'String','Parameter Output');
htext6 = uicontrol(gcf,'Style','text','Position',[.52 .27 .10 .1],'String','Epicenter:');
hedit6 = uicontrol(gcf,'Style','edit','Position',[0.62 0.27 0.16 0.1]);
htext7 = uicontrol(gcf,'Style','text','Position',[.52 .17 .10 .1],'String','Magnitude:');
hedit7 = uicontrol(gcf,'Style','edit','Position',[0.62 0.17 0.16 0.1]);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
htext8 = uicontrol(gcf,'Style','text','Position',[.52 .07 .1 .1],'String','Sita:');
hedit8 = uicontrol(gcf,'Style','edit','Position',[.62 .07 .16 .10]);
htext9 = uicontrol(gcf,'Style','text','Position',[.8 .85 .18 .10],'String','plot channel');
hpopup = uicontrol(gcf,'Style','popup','Position',[0.8 0.75 0.18 0.12],...
    'String','UP|EW|NS');
houtlatlon = uicontrol(gcf,'Style','text','Position',[.8 .65 .20 .10],'String','long & lat Out');
%hlist = uicontrol(gcf,'Style','list','Position',[0.8 0.23 0.28 0.37],...
%    'String','Grid on|Box on|Hidden off|Axis off');
%hpush = uicontrol(gcf,'Style','push','Position',[0.8 0.05 0.18 0.15],...
%    'String','Apply');
%set(hedit,'CallBack','CallEdit(hedit,hpopup,hlist)');
hlontex = uicontrol(gcf,'Style','text','Position',[0.8 0.53 0.12 0.1],'String','longitude:');
hlonedi = uicontrol(gcf,'Style','edit','Position',[0.8 0.43 0.16 0.1]);
hlattex = uicontrol(gcf,'Style','text','Position',[0.8 0.33 0.12 0.1],'String','latitude:');
hlatedi = uicontrol(gcf,'Style','edit','Position',[0.8 0.23 0.16 0.1]);
htext9 = uicontrol(gcf,'Style','text','Position',[0.00 0.85 0.50 0.1],...
    'String','First input the lon&lat of the station');
htext10 = uicontrol(gcf,'style','text','Position',[0.05 0.75 0.05 0.1],'String','lon:');
hinlonedi = uicontrol(gcf,'Style','edit','Position',[0.1 0.75 0.2 0.1]);
htext11 = uicontrol(gcf,'Style','text','Position',[0.25 0.75 0.05 0.1],'String','lat:');
hinlatedi = uicontrol(gcf,'Style','edit','Position',[0.3 0.75 0.2 0.1]);
%set(hpopup,'CallBack','CallEdit(hpopup,hlist)');
%set(hlist,'CallBack','CallEdit(hpopup,hlist)');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
set(hsub1,'Callback','super_cycle();');
set(hsub2,'Callback','azimath;');
set(hsub3,'Callback','mkmod');
set(gcf,'UserData',[hedit1,hedit2,hedit3,hedit4,hedit6,hedit7,hedit8,hlonedi,hlatedi...
    hinlonedi hinlatedi]);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
set(hmenu,'CallBack',[...
    'set(hmenu1,''Visible'',''on'');' ...
    '[name,path] = uigetfile(''*.*'',''please choose a file'',''please open a file'');' ...
    'data = load(strcat(path,name));' ...
    'plot(haxes,data(:,1))']);