plotitmod.m

来自「解决麦克斯韦的matlab源文件」· M 代码 · 共 24 行

function F = plotitmod(mesh,ev,nv,t,figno)
% Plot the modules of the leapfrog iterates
%
% mesh -> 2D triangulation
% ne -> vector of edge dofs of length #of active edges
% nv -> vector of nodal dofs of length #of active vertices

evf = zeros(mesh.Ne,1);
nvf = zeros(2*mesh.Nv,1);

evf(Eactidx(mesh)) = ev;
nvf(Nactidx(mesh)) = nv;
en = edge2vertex(mesh,evf);
s = sprintf('time = %f',t);

figure(figno); clf;
subplot(1,2,1);
plotmod(mesh,en,['Edge elements, ', s]);
view([-30 60]);
subplot(1,2,2);
plotmod(mesh,nvf,['Nodal elements, ', s]);
view([-30 60]);
drawnow;