readdat.m

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

close all;clear all;clc;

dos('del output');
dos('FDTD_1D.exe');
for steps = 1:400
    dos('FDTD_1D.exe output');
    fp = fopen('output','r');
    
    Ex = ones(200,1);
    Hy = ones(200,1);
    
    for k = 1:11
        fgetl(fp);
    end
    for k = 2:200
        Ex(k) = fscanf(fp, '%f',1);
    end
    Ex(1) = 0;

    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]);
    pause(0.1);
end