plotenergies.m

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

function plotenergies(en,meshname)
% Plots temporal evolution of field energies
%
% meshname -> string for indentification of mesh
% en -> trace of electric/magnetic energies during timestepping, N x 5 - array, 
%        en(:,1) = time, 
%        en(:,2) = electric energy of nodal solution
%        en(:,3) = magnetic energy of nodal solution
%        en(:,4) = electric energy of edge element solution
%        en(:,5) = magnetic energy of edge element solution     
%

if (size(en,2) ~= 5), error('Energy array must have 5 columns'); end;
if (nargin < 2)
  meshname = [];
else
  meshname = [meshname , ': '];
end

figure; clf;
subplot(1,2,1);
plot(en(:,1),en(:,2),'r-',en(:,1),en(:,4),'b-');
axis([0 3.1 0 ceil(max(en(:,4)')/50)*50]);
legend('Nodal scheme','Edge elements');
title([meshname , 'Electric energy']);
xlabel('time');
subplot(1,2,2);
plot(en(:,1),en(:,3),'r-',en(:,1),en(:,5),'b-');
axis([0 3.1 0 ceil(max(en(:,5)')/50)*50]);
legend('Nodal scheme','Edge elements');
title([meshname , 'Magnetic energy']);
xlabel('time');
drawnow;