MATLAB——求解特征方程

本代码主要利用MATLAB工具实现MATLAB——求解特征方程，简单明了，易于理解

function pdemodel
[pde_figax]=pdeinit;
pdetool（‘appl_cb‘1）;
set（ax‘DataAspectRatio‘[1 1.2 1]）;
set（ax‘PlotBoxAspectRatio‘[1.5 1 1]）;
set（ax‘XLim‘[-1.5 1.5]）;
set（ax‘YLim‘[-1 1.3999999999999999]）;
set（ax‘XTickMode‘‘auto‘）;
set（ax‘YTickMode‘‘auto‘）;

% Geometry description：
pderect（[-1 1 0 -1]‘R1‘）;
pderect（[0 1 1 -0]‘R2‘）;
set（findobj（get（pde_fig‘Children‘）‘Tag‘‘PDEEval‘）‘String‘‘R1+R2‘）

% Boundary conditions：
pdetool（‘changemode‘0）
pdesetbd（7...
‘dir‘...
1...
‘1‘...
‘0‘）
pdesetbd（6...
‘dir‘...
1...
‘1‘...
‘0‘）
pdesetbd（5...
‘dir‘...
1...
‘1‘...
‘0‘）
pdesetbd（4...
‘dir‘...
1...
‘1‘...
‘0‘）
pdesetbd（3...
‘dir‘...
1...
‘1‘...
‘0‘）
pdesetbd（2...
‘dir‘...
1...
‘1‘...
‘0‘）
pdesetbd（1...
‘dir‘...
1...
‘1‘...
‘0‘）

% Mesh generation：
setappdata（pde_fig‘Hgrad‘1.3）;
setappdata（pde_fig‘refinemethod‘‘regular‘）;
setappdata（pde_fig‘jiggle‘char（‘on‘‘mean‘‘‘））;
pdetool（‘initmesh‘）

% PDE coefficients：
pdeseteq（4...
‘1.0‘...
‘0.0‘...
‘10.0‘...
‘1.0‘...
‘0：10‘...
‘0.0‘...
‘0.0‘...
‘[0 100]‘）
setappdata（pde_fig‘currparam‘...
[‘1.0 ‘;...
‘0.0 ‘;...
‘10.0‘;...
‘1.0 ‘]）

% Solve parameters：
setappdata（pde_fig‘solveparam‘...
str2mat（‘0‘‘1000‘‘10‘‘pdeadworst‘...
‘0.5‘‘longest‘‘0‘‘1E-4‘‘‘‘fixed‘‘Inf‘））

% Plotflags and user data strings：
setappdata（pde_fig‘plotflags‘[1 1 1 1 1 1 1 1 0 0 0 1 1 0 0 0 0 1]）;
setappdata（pde_fig‘colstring‘‘‘）;
setappdata（pde_fig‘arrowstring‘‘‘）;
setappdata（pde_fig‘deformstring‘‘‘）;
setappdata（pde_fig‘heightstring‘‘‘）;

% Solve PDE：
pdetool（‘solve‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6148  2018-05-19 17:07  MATLAB——求解特征方程\.DS_Store

     文件       1665  2018-01-10 13:52  MATLAB——求解特征方程\pdemodel.m

     目录          0  2018-08-23 11:41  MATLAB——求解特征方程

----------- ---------  ---------- -----  ----

                 7813                    3