readmodel.m

一个参照葛德彪书编写的3D FDTD源文件

md=textread('feiji.txt');
md_l=length(md);
for im=2:1+md(1,1)
    mmd(md(im,1):md(im,4),md(im,2),md(im,3))=1;%md(im,5);
    %plot3(md(im,1),md(im,2),md(im,3),'.');
    %plot3(md(im,4),md(im,2),md(im,3),'.');
    %hold on;
end
%for im=2+md(1,1):1+md(1,1)+md(1,2)
%    mmd(md(im,1),md(im,2):md(im,4),md(im,3))=md(im,5);
%    plot3(md(im,1),md(im,2),md(im,3),'.');
%    plot3(md(im,1),md(im,4),md(im,3),'.');
%end
%for im=2+md(1,1)+md(1,2):md_l
%    mmd(md(im,1),md(im,2),md(im,3):md(im,4))=md(im,5);
%    plot3(md(im,1),md(im,2),md(im,3),'.');
%    plot3(md(im,1),md(im,2),md(im,4),'.');
%end
%[mmdi,mmdj,mmdk]=size(mmd);
%xlabel('x');
%ylabel('y');
%zlabel('z');
%axis([1  mmdi+3  1  mmdj+3  1  mmdk+3])