科氏力与地球自转产生的惯性离心力差别
Posted
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了科氏力与地球自转产生的惯性离心力差别相关的知识，希望对你有一定的参考价值。
参考技术A 区别：1.只有运动物体才会受到科里奥利力,而惯性离心力不管物体运动还是静止都有.2.方向不同.惯性离心力方向是从物体所在位置自转的小圆圆心指向物体所在位置,方向一般有垂直于地表向上的分量；科里奥利力方向是垂直...运动学基于matlab GUI地球自转模拟含Matlab源码 1115期
一、简介
基于matlab GUI地球自转模拟
二、源代码
function varargout = globegui(varargin)
% GLOBEGUI M-file for globegui.fig created with GUIDE
%
% Creates a GUI for viewing a spinning Earth globe
%
% GLOBEGUI, by itself, creates a new GLOBEGUI or raises the existing
% one.
%
% H = GLOBEGUI returns the handle to a new GLOBEGUI or the handle to
% the existing one.
%
% GLOBEGUI('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in GLOBEGUI.M with the given input arguments.
%
% GLOBEGUI('Property','Value',...) creates a new GLOBEGUI or raises the
% existing one. Starting from the left, property value pairs are
% applied to the GUI before globegui_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to globegui_OpeningFcn via varargin.
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Last Modified by GUIDE v2.5 10-Mar-2021 09:50:13
%
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @globegui_OpeningFcn, ...
'gui_OutputFcn', @globegui_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 globegui is made visible.
function globegui_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 globegui (see VARARGIN)
% Check whether initialization has been done in case this GUI
% is a singleton and has already been opened.
% if isfield(handles,'running')
% disp('No initialization this time.')
% else
% handles.running = true;
% guidata(hObject,handles)
% disp('Initializing...')
% end
% Choose default command line output for globegui
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes globegui wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% --- Outputs from this function are returned to the command line.
function varargout = globegui_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 spinstopbutton.
function spinstopbutton_Callback(hObject, eventdata, handles)
% Spins the globe in the axes or stops it, renaming the button
% from "Spin" to "Stop" and back again. When its label is "Stop"
% this callback is executing an endless loop and is re-entered;
% thus it must have properties set as follows (which is the
% default behavior of a GUIDE GUI):
% Interruptible: 'on', BusyAction: 'queue'
%
% hObject handle to spinstopbutton (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
str = get(hObject,'String'); % get the current pushbutton string
% Determine which string the label matches
state = find(strcmp(str,handles.Strings));
% Toggle the button label to other string
set(hObject,'String',handles.Strings{3-state});
% If the index when entering was 1, start to spin the object
if (state == 1)
% globe = struct;
filming = handles.movie;
az = handles.azimuth;
hgrotate = handles.tform;
% Spin globe as long as the figure exists or until user
% interrupts by pressing the button a second time
while ishandle(handles.axes1)
% If button label changed since last iteration, stop now
if find(strcmp(get(hObject,'String'),...
handles.Strings)) == 1
% Save rotation state to restart at this orientation
handles.azimuth = az;
guidata(hObject,handles);
break
end
az = az + 0.01745329252; % Increment azimuth (in radians)
% to rotate east one degree
% Modify the hgtransform controling the two surface objects
set(hgrotate,'Matrix',makehgtform('zrotate',az));
drawnow % Refresh the screen
% If the Make movie button is checked, save frames
% but don't store more than one revolution
% NOTE: filming slows down the animation
% Need to test whether axes exists because user can quit
% during filming, destroying axes and figure
if ishandle(handles.axes1) && filming > 0 && filming < 361
globeframes(filming) = getframe(handles.axes1);
filming = filming + 1;
end
end
% Write captured frames to MAT-file if in movie mode
if (filming)
filename = sprintf('globe%i.mat',filming-1);
disp(['Writing movie to file ' filename]);
save (filename, 'globeframes')
end
end
% --- Executes on button press in quitbutton.
function quitbutton_Callback(hObject, eventdata, handles)
% hObject handle to quitbutton (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
%handles = guidata(hObject);
% Get the figure's handle, then destroy it
fig = handles.figure1;
close(fig)
% --- Executes during object creation, after setting all properties.
function spinstopbutton_CreateFcn(hObject, eventdata, handles)
% hObject handle to spinstopbutton (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - not created until after all CreateFcns called
% Creates the handles structure and places into it label strings
% to test button's current name against
handles.Strings = {'Spin';'Stop'};
% Commit the new struct element to appdata
guidata(hObject, handles);
% --- Executes during object creation, after setting all properties.
function axes1_CreateFcn(hObject, eventdata, handles)
% hObject handle to axes1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - not created until after all CreateFcns called
% Hint: place code in OpeningFcn to populate axes1
% Generate a colormap appropriate for terrain display
cmap = [0 0 0.2000; 0 0 0.2471;...
0 0 0.2941; 0 0 0.3412;...
0 0 0.3882; 0 0 0.4353;...
0 0 0.4824; 0 0 0.5294;...
0 0 0.5765; 0 0 0.6235;...
0 0 0.6706; 0 0 0.7176;...
0 0 0.7647; 0 0 0.8118;...
0 0 0.8588; 0 0 0.9059;...
0 0 0.9529; 0 0 1.0000;...
0 0.0556 1.0000; 0 0.1111 1.0000;...
0 0.1667 1.0000; 0 0.2222 1.0000;...
0 0.2778 1.0000; 0 0.3333 1.0000;...
0 0.3889 1.0000; 0 0.4444 1.0000;...
0 0.5000 1.0000; 0 0.5556 1.0000;...
0 0.6111 1.0000; 0 0.6667 1.0000;...
0 0.7222 1.0000; 0 0.7778 1.0000;...
0 0.8333 1.0000; 0 0.8889 1.0000;...
0 0.9444 1.0000; 0 1.0000 1.0000;...
0 0.4000 0.2000; 0.0253 0.4429 0.2043;...
0.0555 0.4857 0.2092; 0.0906 0.5286 0.2157;...
0.1306 0.5714 0.2251; 0.1755 0.6143 0.2382;...
0.2253 0.6571 0.2562; 0.2800 0.7000 0.2800;...
0.3684 0.7429 0.3396; 0.4586 0.7857 0.4041;...
0.5496 0.8286 0.4735; 0.6402 0.8714 0.5478;...
0.7296 0.9143 0.6269; 0.8165 0.9571 0.7110;...
0.9000 1.0000 0.8000; 0.8499 0.9538 0.7044;...
0.8099 0.9077 0.6144; 0.7787 0.8615 0.5302;...
0.7548 0.8154 0.4516; 0.7367 0.7692 0.3787;...
0.7231 0.7231 0.3115; 0.6769 0.6413 0.2499;...
0.6308 0.5580 0.1941; 0.5846 0.4744 0.1439;...
0.5385 0.3921 0.0994; 0.4923 0.3124 0.0606;...
0.4462 0.2368 0.0275; 0.4000 0.1667 0];
load topo % Get 1x1 degree terrain grid
% Make axes a slightly oversized unit box centered on 0,0,0
set(hObject,'xlim',[-1.02 1.02],...
'ylim',[-1.02 1.02],...
'zlim',[-1.02 1.02]);
% Create a spherical structure
[x,y,z] = sphere(50);
hgttilt = hgtransform;
hgrotate = hgtransform('parent',hgttilt);
% Set display properties
props.FaceColor= 'texture';
props.EdgeColor = 'none';
props.FaceLighting = 'gouraud';
props.Cdata = topo; % Use topo grid as a texturemap
props.Parent = hgrotate; % Make hgtransform surface parent
hsurf = surface(x,y,z,props); % Draw 3-D view
colormap(cmap) % Use special terrain colormap defined above
% Rotate the surface by 23.44 deg (0.4091 radians) around x-axis;
% this is the earth's tilt from normal to the ecliptic.
% To learn about geometric operations, type "doc hgtransform".
set(hgttilt,'Matrix',makehgtform('xrotate',0.4091));
% Create another mesh to be the graticule
[gx,gy,gz] = sphere(15);
% Decimate every other row to make mesh elements square
for j = 2:9
gx(j,:) = [];
gy(j,:) = [];
gz(j,:) = [];
end
三、运行结果
四、备注
版本：2014a
以上是关于科氏力与地球自转产生的惯性离心力差别的主要内容，如果未能解决你的问题，请参考以下文章