readdat.asv

一维FDTD

close all;clear all;clc;

dos('del Ex');
dos('del Hy');
T = 0;
for NSTEPS=
dos('a');
ex = ones(200);
hy = ones(200);
fp = fopen('Ex','r');
for k = 2:200
    ex(k) = fscanf(fp, '%f',1);
end
ex(1) = 0;
fclose(fp);

fp = fopen('Hy','r');
for k = 1:199
    hy(k) = fscanf(fp,'%f',1);
end
hy(200) = 0;
fclose(fp);

subplot(211);
plot(ex);
subplot(212);
plot(hy);