function varargout = cylinder(varargin)
% CYLINDER MATLAB code for cylinder.fig
% CYLINDER, by itself, creates a new CYLINDER or raises the existing
% singleton*.
%
% H = CYLINDER returns the handle to a new CYLINDER or the handle to
% the existing singleton*.
%
% CYLINDER('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in CYLINDER.M with the given input arguments.
%
% CYLINDER('Property','Value',...) creates a new CYLINDER or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before cylinder_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to cylinder_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help cylinder
% Last Modified by GUIDE v2.5 29-Apr-2023 01:11:39
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @cylinder_OpeningFcn, ...
'gui_OutputFcn', @cylinder_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
% --- Executes just before cylinder is made visible.
function cylinder_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to cylinder (see VARARGIN)
% Choose default command line output for cylinder
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes cylinder wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% --- Outputs from this function are returned to the command line.
function varargout = cylinder_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
% --- Executes on button press in radiobutton1.
function radiobutton1_Callback(hObject, eventdata, handles)%画△g主剖面
% hObject handle to radiobutton1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of radiobutton1
global D
global t
G=6.67*1e-2;
x=-200:10:200;
g=(2*G*t*D)./(x.^2+D^2);
plot(x,g,'*-','MarkerEdgeColor','r');
box on;
xlabel('X(m)');
ylabel('△g.u.');
function edit2_Callback(hObject, eventdata, handles)%对埋深D赋值
% hObject handle to edit2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit2 as text
% str2double(get(hObject,'String')) returns contents of edit2 as a double
global D
get(hObject,'String');
D=str2double(get(hObject,'String'));
handles.D=D;
guidata(hObject,handles);
% --- Executes during object creation, after setting all properties.
function edit2_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
function edit3_Callback(hObject, eventdata, handles)%对线密度t赋值
% hObject handle to edit3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit3 as text
% str2double(get(hObject,'String')) returns contents of edit3 as a double
global t
get(hObject,'String');
t=str2double(get(hObject,'String'));
handles.t=t;
guidata(hObject,handles);
% --- Executes during object creation, after setting all properties.
function edit3_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
% --- Executes on button press in radiobutton2.
function radiobutton2_Callback(hObject, eventdata, handles)%画重力异常平面等值线图
% hObject handle to radiobutton2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of radiobutton2
global D
global t
G=6.67*1e-2;
x=-200:10:200;
y=-200:10:200;
for i=1:size(x,2)
for j=1:size(x,2)
g(j,i)=(2*G*t*D)/(x(i)^2+D^2);
end
end
contourf(x,y,g,12);
colorbar;
box on;
xlabel('X(m)');
ylabel('Y(m)');
% --- Executes on button press in radiobutton3.
function radiobutton3_Callback(hObject, eventdata, handles)
% hObject handle to radiobutton3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of radiobutton3
global D
global t
global az
global el
G=6.67*1e-2;
x=-200:10:200;
y=-200:10:200;
for i=1:size(x,2)
for j=1:size(x,2)
g(j,i)=(2*G*t*D)/(x(i)^2+D^2);
end
end
mesh(x,y,g);
view([az,el]);
colorbar;
box on;
xlabel('X(m)');
ylabel('Y(m)');
zlabel('△g.u.');
% --- Executes on button press in radiobutton5.
function radiobutton5_Callback(hObject, eventdata, handles)%画Vxz主剖面
% hObject handle to radiobutton5 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of radiobutton5
global D
global t
G=6.67*1e-2;
x=-200:10:200;;
Vxz=((4*G*t*(-D)).*x)./((x.^2+D^2).^2);
plot(x,Vxz,'*-','MarkerEdgeColor','r');
box on;
xlabel('X(m)');
ylabel('Vxz');
% --- Executes on button press in radiobutton7.
function radiobutton7_Callback(hObject, eventdata, handles)%画Vxz平面图
% hObject handle to radiobutton7 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of radiobutton7
global D
global t
G=6.67*1e-2;
x=-200:10:200;
y=-200:10:200;
for i=1:size(x,2)
for j=1:size(x,2)
Vxz(j,i)=(4*G*t*(-D)*x(i))/((x(i)^2++D^2).^2);
end
end
contourf(x,y,Vxz,12);
colorbar;
box on;
xlabel('X(m)');
ylabel('Y(m)');
zlabel('Vxz');
% --- Executes on button press in radiobutton8.
function radiobutton8_Callback(hObject, eventdata, handles)%画Vxz三维图
% hObject handle to radiobutton8 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with
