www.gusucode.com > pde 案例源码 matlab代码程序 > pde/Convert2DMeshToPetFormExample.m

    %% Convert 2-D Mesh to [p,e,t] Form
% This example shows how to convert a mesh in object form to |[p,e,t]|
% form.
%%
% Create a 2-D PDE geometry and incorporate it into a model object. View
% the geometry.
model = createpde(1);

R1 = [3,4,-1,1,1,-1,-.4,-.4,.4,.4]';
C1 = [1,.5,0,.2]';
% Pad C1 with zeros to enable concatenation with R1
C1 = [C1;zeros(length(R1)-length(C1),1)];
geom = [R1,C1];
ns = (char('R1','C1'))';
sf = 'R1-C1';
gd = decsg(geom,sf,ns);

geometryFromEdges(model,gd);
pdegplot(model,'EdgeLabels','on')
xlim([-1.1 1.1])
axis equal
%%
% Create a mesh for the geometry. View the mesh.
generateMesh(model);
pdemesh(model)
axis equal
%%
% Convert the mesh to |[p,e,t]| form.
[p,e,t] = meshToPet(model.Mesh);
%%
% View the sizes of the |[p,e,t]| matrices.
size(p)
%%
size(e)
%%
size(t)