-
Notifications
You must be signed in to change notification settings - Fork 2
/
pplane2.m
391 lines (294 loc) · 15.4 KB
/
pplane2.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
function varargout = pplane2(varargin)
% PPLANE2 MATLAB code for pplane2.fig
% PPLANE2, by itself, creates a new PPLANE2 or raises the existing
% singleton*.
%
% H = PPLANE2 returns the handle to a new PPLANE2 or the handle to
% the existing singleton*.
%
% PPLANE2('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in PPLANE2.M with the given input arguments.
%
% PPLANE2('Property','Value',...) creates a new PPLANE2 or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before pplane2_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to pplane2_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 pplane2
% Last Modified by GUIDE v2.5 31-Mar-2017 14:43:42
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @pplane2_OpeningFcn, ...
'gui_OutputFcn', @pplane2_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 pplane2 is made visible.
function pplane2_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 pplane2 (see VARARGIN)
% Choose default command line output for pplane2
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes pplane2 wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% --- Outputs from this function are returned to the command line.
function varargout = pplane2_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;
function txt_xdot_Callback(hObject, eventdata, handles)
% hObject handle to txt_xdot (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 txt_xdot as text
% str2double(get(hObject,'String')) returns contents of txt_xdot as a double
% --- Executes during object creation, after setting all properties.
function txt_xdot_CreateFcn(hObject, eventdata, handles)
% hObject handle to txt_xdot (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 txt_ydot_Callback(hObject, eventdata, handles)
% hObject handle to txt_ydot (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 txt_ydot as text
% str2double(get(hObject,'String')) returns contents of txt_ydot as a double
% --- Executes during object creation, after setting all properties.
function txt_ydot_CreateFcn(hObject, eventdata, handles)
% hObject handle to txt_ydot (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 txt_minx_Callback(hObject, eventdata, handles)
% hObject handle to txt_minx (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 txt_minx as text
% str2double(get(hObject,'String')) returns contents of txt_minx as a double
% --- Executes during object creation, after setting all properties.
function txt_minx_CreateFcn(hObject, eventdata, handles)
% hObject handle to txt_minx (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 txt_maxx_Callback(hObject, eventdata, handles)
% hObject handle to txt_maxx (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 txt_maxx as text
% str2double(get(hObject,'String')) returns contents of txt_maxx as a double
% --- Executes during object creation, after setting all properties.
function txt_maxx_CreateFcn(hObject, eventdata, handles)
% hObject handle to txt_maxx (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 txt_miny_Callback(hObject, eventdata, handles)
% hObject handle to txt_miny (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 txt_miny as text
% str2double(get(hObject,'String')) returns contents of txt_miny as a double
% --- Executes during object creation, after setting all properties.
function txt_miny_CreateFcn(hObject, eventdata, handles)
% hObject handle to txt_miny (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 txt_maxy_Callback(hObject, eventdata, handles)
% hObject handle to txt_maxy (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 txt_maxy as text
% str2double(get(hObject,'String')) returns contents of txt_maxy as a double
% --- Executes during object creation, after setting all properties.
function txt_maxy_CreateFcn(hObject, eventdata, handles)
% hObject handle to txt_maxy (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 btn_graph_phase_plane.
function btn_graph_phase_plane_Callback(hObject, eventdata, handles)
% hObject handle to btn_graph_phase_plane (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
global equ_xdot;
global equ_ydot;
equ_xdot_=get(handles.txt_xdot, 'String');
equ_ydot_=get(handles.txt_ydot, 'String');
equ_xdot = strrep(equ_xdot_, 'x', 'x(1)');
equ_xdot = strrep(equ_xdot, 'y', 'x(2)');
equ_ydot = strrep(equ_ydot_, 'x', 'x(1)');
equ_ydot = strrep(equ_ydot, 'y', 'x(2)');
xmin = str2double(get(handles.txt_minx, 'String'));
xmax = str2double(get(handles.txt_maxx, 'String'));
ymin = str2double(get(handles.txt_miny, 'String'));
ymax = str2double(get(handles.txt_maxy, 'String'));
simtime = str2double(get(handles.txt_simtime, 'String'));
resolution = [str2double(get(handles.txt_resolutionh, 'String')) ...
str2double(get(handles.txt_resolutionv, 'String'))];
phase_plot_2_interactive(@f, [xmin xmax; ymin ymax], simtime, ...
{strcat('$\dot{x}=', equ_xdot_, '$') strcat('$\dot{y}=', equ_ydot_, '$')}, ...
resolution);
function txt_simtime_Callback(hObject, eventdata, handles)
% hObject handle to txt_simtime (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 txt_simtime as text
% str2double(get(hObject,'String')) returns contents of txt_simtime as a double
% --- Executes during object creation, after setting all properties.
function txt_simtime_CreateFcn(hObject, eventdata, handles)
% hObject handle to txt_simtime (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
% --- Callback function of the ODE
function xdot = f(t, x)
global equ_xdot;
global equ_ydot;
xdot(1,1) = eval(equ_xdot);
xdot(2,1) = eval(equ_ydot);
function txt_resolutionh_Callback(hObject, eventdata, handles)
% hObject handle to txt_resolutionh (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 txt_resolutionh as text
% str2double(get(hObject,'String')) returns contents of txt_resolutionh as a double
% --- Executes during object creation, after setting all properties.
function txt_resolutionh_CreateFcn(hObject, eventdata, handles)
% hObject handle to txt_resolutionh (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 txt_resolutionv_Callback(hObject, eventdata, handles)
% hObject handle to txt_resolutionv (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 txt_resolutionv as text
% str2double(get(hObject,'String')) returns contents of txt_resolutionv as a double
% --- Executes during object creation, after setting all properties.
function txt_resolutionv_CreateFcn(hObject, eventdata, handles)
% hObject handle to txt_resolutionv (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 mnu_file_Callback(hObject, eventdata, handles)
% hObject handle to mnu_file (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function mnu_load_Callback(hObject, eventdata, handles)
% hObject handle to mnu_load (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
[file, path] = uigetfile('*.dat', 'Select file');
if file ~= 0
T = readtable(strcat(path, file));
data = table2cell(T);
set(handles.txt_xdot, 'String', data{1, :});
set(handles.txt_ydot, 'String', data{2, :});
set(handles.txt_minx, 'String', data{3, :});
set(handles.txt_maxx, 'String', data{4, :});
set(handles.txt_miny, 'String', data{5, :});
set(handles.txt_maxy, 'String', data{6, :});
set(handles.txt_simtime, 'String', data{7, :});
set(handles.txt_resolutionh, 'String', data{8, :})
set(handles.txt_resolutionv, 'String', data{9, :});
end
% --------------------------------------------------------------------
function mnu_save_Callback(hObject, eventdata, handles)
% hObject handle to mnu_save (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
ode2ndorder = {get(handles.txt_xdot, 'String'); ...
get(handles.txt_ydot, 'String'); ...
get(handles.txt_minx, 'String'); ...
get(handles.txt_maxx, 'String'); ...
get(handles.txt_miny, 'String'); ...
get(handles.txt_maxy, 'String'); ...
get(handles.txt_simtime, 'String'); ...
get(handles.txt_resolutionh, 'String');...
get(handles.txt_resolutionv, 'String')};
[file, path] = uiputfile('myode.dat', 'Save file name');
if file ~= 0
T = cell2table(ode2ndorder);
writetable(T, strcat(path, file));
end
% --------------------------------------------------------------------
function mnu_help_Callback(hObject, eventdata, handles)
% hObject handle to mnu_help (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function mnu_instruction_Callback(hObject, eventdata, handles)
% hObject handle to mnu_instruction (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
msgbox({'Please try the provided examples (.dat files).'}, 'Instruction');
% --------------------------------------------------------------------
function mnu_about_Callback(hObject, eventdata, handles)
% hObject handle to mnu_about (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
msgbox({'manurunga@yandex.com'; 'http://kataauralius.com/'}, 'About', 'modal');