readdat.m

本程序是使用加入边界条件的有限时域差分法计算模拟电磁波在光学薄膜中的传输

close all;clear all;clc;

 dos('del output');
 dos('FDTD_init.exe output');%pause
for steps = 1:50
    dos('FDTD_1D.exe output');
    fp = fopen('output','r');

    Ex = ones(200,1);
    Hy = ones(200,1);

    for k = 1:13
        fgetl(fp);
    end
    for k = 1:201
        Ex(k) = fscanf(fp, '%f',1);
    end

    for k = 1:200
        Hy(k) = fscanf(fp,'%f',1);
    end
    fclose(fp);

    subplot(211);
    plot(Ex);
    axis([1 200 -1 1]);
    subplot(212);
    plot(Hy);
    axis([1 200 -1 1]);
    pause(0.1);
end