MATLAB——求解矩形域内的波动方程

本代码主要利用MATLAB工具实现MATLAB——求解矩形域内的波动方程，简单明了，易于理解

function pdemodel
[pde_figax]=pdeinit;
pdetool（‘appl_cb‘1）;
set（ax‘DataAspectRatio‘[1 1.5 1]）;
set（ax‘PlotBoxAspectRatio‘[1.5 1 1]）;
set（ax‘XLim‘[-1.5 1.5]）;
set（ax‘YLim‘[-1.5 1.5]）;
set（ax‘XTickMode‘‘auto‘）;
set（ax‘YTick‘[ -1.2...
 -1...
 -0.79999999999999993...
 -0.59999999999999987...
 -0.39999999999999991...
 -0.19999999999999996...
 0...
 0.19999999999999996...
 0.39999999999999991...
 0.59999999999999987...
 0.79999999999999993...
 1...
 1.2...
]）;

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

% Boundary conditions：
pdetool（‘changemode‘0）
pdesetbd（4...
‘dir‘...
1...
‘1‘...
‘0‘）
pdesetbd（3...
‘neu‘...
1...
‘0‘...
‘0‘）
pdesetbd（2...
‘dir‘...
1...
‘1‘...
‘0‘）
pdesetbd（1...
‘neu‘...
1...
‘0‘...
‘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（3...
‘1.0‘...
‘0.0‘...
‘0.0‘...
‘1.0‘...
‘linspace（0531）‘...
‘atan（cos（pi/2*x））‘...
‘3*sin（pi*x）.*exp（sin（pi/2*y））‘...
‘[0 100]‘）
setappdata（pde_fig‘currparam‘...
[‘1.0‘;...
‘0.0‘;...
‘0.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 1 0 1 31 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

     文件       1772  2018-01-10 13:52  MATLAB——求解矩形域内的波动方程\pdemodel.m

     目录          0  2018-08-23 11:41  MATLAB——求解矩形域内的波动方程

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

                 7920                    3