MATLAB——求指定区域上泊松方程的数值解

本代码主要利用MATLAB工具实现MATLAB——求指定区域上泊松方程的数值解，简单明了，易于理解

function pdemodel
[pde_figax]=pdeinit;
pdetool（‘appl_cb‘1）;
set（ax‘DataAspectRatio‘[1 1 1]）;
set（ax‘PlotBoxAspectRatio‘[921.59999999999991 614.39999999999998 614.39999999999998]）;
set（ax‘XLimMode‘‘auto‘）;
set（ax‘YLim‘[-1.5 1.5]）;
set（ax‘XTickMode‘‘auto‘）;
set（ax‘YTickMode‘‘auto‘）;
pdetool（‘gridon‘‘on‘）;
% Geometry description：
pdeellip（-0.00081433224755711464-0.00244299674267067780.994299674267101091.001628664495114...
0‘E1‘）;
set（findobj（get（pde_fig‘Children‘）‘Tag‘‘PDEEval‘）‘String‘‘E1‘）
% Boundary conditions：
pdetool（‘changemode‘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‘）
pdetool（‘refine‘）

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

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

% Plotflags and user data strings：
setappdata（pde_fig‘plotflags‘[4 1 1 1 1 1 1 1 0 0 0 1 1 0 0 0 0 1]）;
setappdata（pde_fig‘colstring‘‘u-（1-x.^2-y.^2）/4‘）;
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

     文件       1655  2018-01-10 13:52  MATLAB——求指定区域上泊松方程的数值解\pdemodel.m

     目录          0  2018-08-23 11:41  MATLAB——求指定区域上泊松方程的数值解

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

                 7803                    3