我想使用MATLAB中的PDE工具箱解决一个简单的热方程,使用以下代码:
clear model;
width = 1e-5;
height = 1e-5;
gdm = [3 4 0 width width 0 0 0 height height]';
g = decsg(gdm, 'S1', ('S1')');
numberOfPDE = 1;
model = createpde(numberOfPDE);
geometryFromEdges(model, g);
% figure;
% pdegplot(model,'EdgeLabels','on');
% axis([-.1*width 1.1*height -.1*width 1.1*height]);
% title 'Geometry With Edge Labels Displayed';
hmax = .1 * width;
msh = generateMesh(model, 'Hmax', hmax);
p = msh.Nodes;
% figure;
% pdeplot(model);
% axis equal
% title 'Plate With Triangular Element Mesh'
% xlabel 'X-coordinate, meters'
% ylabel 'Y-coordinate, meters'
pulseDuration = 1e-13;