readdat.m

来自「本程序是使用有限时域差分法计算模拟电磁波在光学薄膜中的传输」· M 代码 · 共 30 行

close all;clear all;clc;

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

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

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

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