readdat.m

一维FDTD

close all;clear all;clc;

dos('del Ex');
dos('del Hy');
dos('FDTD_1.1_change_twoSource.exe');
T = 1;
for steps = 1:400
    dos(['FDTD_1.1_change_twoSource.exe ' num2str(T) ' ' num2str(1)]);
    ex = ones(200,1);
    hy = ones(200,1);
    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);
    axis([1 200 0 1]);
    subplot(212);
    plot(hy);
    axis([1 200 -1 1]);
    T = T+1;
    pause(0.1);
end