demo_comp.sci

来自「用来实现三维阻抗及光学断层成像重建的matlab程序」· SCI 代码 · 共 775 行 · 第 1/2 页

// 
tfac = 0.0000001;
// 
%v=size(sm1)
%v=size(sm1)
MR = [sm1,zeros(%v(1),%v(2));zeros(%v(1),%v(2)),sm1];
// 
 
//! mtlb(Jrr) can be replaced by Jrr() or Jrr whether Jrr is an m-file or not
 
//! mtlb(Jri) can be replaced by Jri() or Jri whether Jri is an m-file or not
 
//! mtlb(Jir) can be replaced by Jir() or Jir whether Jir is an m-file or not
 
//! mtlb(Jii) can be replaced by Jii() or Jii whether Jii is an m-file or not
MJ = [mtlb(Jrr),mtlb(Jri);mtlb(Jir),-mtlb(Jii)];
// 
dat = [dvrG;dviG];
// 
//load datacom MJ MR
// 
sol = (MJ'*MJ+tfac*MR)\MJ'*dat;
// 
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
sreal = sol(1:size(mtlb(%simp),1));
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
 
//! mtlb(end) can be replaced by end() or end whether end is an m-file or not
simag = sol(size(mtlb(%simp),1)+1:mtlb(end));
// 
mtlb_load('datacom','fc');
// 
xset('window',max(winsid()+1))
h1 = max(winsid());
 
//!! Scilab set function is a partial emulation of the Matlab one
set(h1,'NumberTitle','off');
 
//!! Scilab set function is a partial emulation of the Matlab one
set(h1,'Name','Simulated inhomogeneities');
%v$2 = int((1-1)/2);%v$1 = 1-1-2*%v$2
xsetech([%v$1/2,%v$2/1,1/2,1/1]);
 
//! mtlb(srf) can be replaced by srf() or srf whether srf is an m-file or not
 
//!! Unknown variable vtx ,the original calling sequence is used
 
//!! Unknown variable vtx ,the original calling sequence is used
 
//!! Unknown variable vtx ,the original calling sequence is used
 
//!! Unknown variable trimesh ,the original calling sequence is used
trimesh(mtlb(srf),vtx(:,1),vtx(:,2),vtx(:,3));
 
//!! Unknown variable axis ,the original calling sequence is used
axis('image');
 
//!! Scilab set function is a partial emulation of the Matlab one
set(xget('window'),'Colormap',[0,0,0]);
 
//!! Unknown variable hidden ,the original calling sequence is used
hidden('off');
mtlb_hold('on');
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
repaint_inho(real(mat),real(mat_ref),mtlb(vtx),mtlb(%simp));
xtitle('Original conductivity distribution',' ',' ');
 
//!! Unknown variable camlight ,the original calling sequence is used
camlight('left');
 
//!! Unknown variable lighting ,the original calling sequence is used
lighting('flat');
// 
%v$2 = int((2-1)/2);%v$1 = 2-1-2*%v$2
xsetech([%v$1/2,%v$2/1,1/2,1/1]);
 
//! mtlb(srf) can be replaced by srf() or srf whether srf is an m-file or not
 
//!! Unknown variable vtx ,the original calling sequence is used
 
//!! Unknown variable vtx ,the original calling sequence is used
 
//!! Unknown variable vtx ,the original calling sequence is used
 
//!! Unknown variable trimesh ,the original calling sequence is used
trimesh(mtlb(srf),vtx(:,1),vtx(:,2),vtx(:,3));
 
//!! Unknown variable axis ,the original calling sequence is used
axis('image');
 
//!! Scilab set function is a partial emulation of the Matlab one
set(xget('window'),'Colormap',[0,0,0]);
 
//!! Unknown variable hidden ,the original calling sequence is used
hidden('off');
mtlb_hold('on');
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
repaint_inho(imag(mat),imag(mat_ref),mtlb(vtx),mtlb(%simp));
xtitle('Original permittivity distribution',' ',' ');
 
//!! Unknown variable camlight ,the original calling sequence is used
camlight('left');
 
//!! Unknown variable lighting ,the original calling sequence is used
lighting('flat');
// 
// 
disp('And their reconstructed images');
xset('window',max(winsid()+1))
h2 = max(winsid());
 
//!! Scilab set function is a partial emulation of the Matlab one
set(h2,'NumberTitle','off');
 
//!! Scilab set function is a partial emulation of the Matlab one
set(h2,'Name','Reconstructed conductivity distribution');
%v$2 = int((1-1)/3);%v$1 = 1-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
 
//! mtlb(fc) can be replaced by fc() or fc whether fc is an m-file or not
fc = slicer_plot(2.63,sreal,mtlb(vtx),mtlb(%simp),mtlb(fc));
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=2.60',' ',' ');
%v$2 = int((2-1)/3);%v$1 = 2-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(2.1,sreal,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=2.10',' ',' ');
%v$2 = int((3-1)/3);%v$1 = 3-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(1.72,sreal,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=1.70',' ',' ');
%v$2 = int((4-1)/3);%v$1 = 4-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(1.1,sreal,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=1.10',' ',' ');
%v$2 = int((5-1)/3);%v$1 = 5-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(0.83,sreal,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=0.80',' ',' ');
%v$2 = int((6-1)/3);%v$1 = 6-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(0.1,sreal,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=0.10',' ',' ');
// 
xset('window',max(winsid()+1))
h3 = max(winsid());
 
//!! Scilab set function is a partial emulation of the Matlab one
set(h3,'NumberTitle','off');
 
//!! Scilab set function is a partial emulation of the Matlab one
set(h3,'Name','Reconstructed scaled permittivity distribution');
%v$2 = int((1-1)/3);%v$1 = 1-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(2.63,simag,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=2.60',' ',' ');
%v$2 = int((2-1)/3);%v$1 = 2-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(2.1,simag,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=2.10',' ',' ');
%v$2 = int((3-1)/3);%v$1 = 3-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(1.72,simag,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=1.70',' ',' ');
%v$2 = int((4-1)/3);%v$1 = 4-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(1.1,simag,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=1.10',' ',' ');
%v$2 = int((5-1)/3);%v$1 = 5-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(0.83,simag,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=0.80',' ',' ');
%v$2 = int((6-1)/3);%v$1 = 6-1-3*%v$2
xsetech([%v$1/3,%v$2/2,1/3,1/2]);
 
//! mtlb(vtx) can be replaced by vtx() or vtx whether vtx is an m-file or not
 
//! mtlb(%simp) can be replaced by %simp() or %simp whether %simp is an m-file or not
fc = slicer_plot(0.1,simag,mtlb(vtx),mtlb(%simp),fc);
 
//!! Unknown variable view ,the original calling sequence is used
view(2);
 
//!! Unknown variable grid ,the original calling sequence is used
grid();
 
//!! Unknown variable colorbar ,the original calling sequence is used
colorbar();
 
//!! Unknown variable axis ,the original calling sequence is used
axis('off');
xtitle('z=0.10',' ',' ');
// 
disp('Done');
// 
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
// This is part of the EIDORS suite.
// Copyright (c) N. Polydorides 2001
// Copying permitted under terms of GNU GPL
// See enclosed file gpl.html for details.
// EIDORS 3D version 1.0
// MATLAB version 5.3 R11
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%