plotiterate.m

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

function F = plotiterate(mesh,ev,nv,t,figno)
% Plots the current iterate during the leapfrog iteration
%
% mesh -> 2D triangulation
% ne -> vector of edge dofs of length #of active edges
% nv -> vector of nodal dofs of length #of active vertices
% t -> time (for title)

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(2,2,1);
plotfield(mesh,en,'Edge elements');
%xlabel('Edge elements');
subplot(2,2,3);
plotmod(mesh,en,s);
view([-30 70]);
%xlabel('Edge elements');
subplot(2,2,2);
plotfield(mesh,nvf,'Nodal elements');
%xlabel('Nodal elements');
subplot(2,2,4);
plotmod(mesh,nvf,s);
view([-30 70]);
%xlabel('Nodal elements');
drawnow;