makemovie.m

非连续波导时间域有限差分法仿真（TD-FDM）

for i=1:MAX_Y
    for j=1:MAX_X
        cnt = nmbr(i,j);
        if cnt==0
            M(i,j,:) = 0;
        else
            M(i,j,:) = Hz(cnt,:);
        end
    end
end

for i=1:MAX_Y
    for j=1:MAX_X
        xx(i,j) = (j-1)*hx;
        yy(i,j) = (i-1)*hy;
    end
end

figure;

input('wanna start the movie?');

scale = [0 1.2];

for k=1:step
    contourf(xx,yy,M(:,:,k));
    caxis(scale);
    colorbar('ylim',scale);
%colorbar;
    axis equal;
    axis([xx(1,1) xx(1,MAX_X) yy(1,1) yy(MAX_Y,1)]);
    
    frame = getframe;
    X = frame2im(frame);
    dir = 'images/';
    number = num2str(k);
%    ext = 'png';
    filename = [dir, number];
    imwrite(X, eval('filename'), 'bmp');    
    
end