求助！！！重建ECT 算法

冷凌川L

根据ECT数据重建图像。实验对象是一个6cm柱形容器内靠近管壁相对放置的两个有机玻璃棒，玻璃棒直径分别为2.5cm和1.5cm。使用12电极ECT传感器测量。给出灵敏度矩阵S，满场数据P2，空场数据P0，测量数据P1。另外nodeeit.m文件为重建图像剖分网格函数。可生成方形剖分网格数812个。要求是：根据所给数据，用一种图像重建算法重建ECT图像，使重建图像质量最好。数据见群共享压缩文件“ECT数据.rar”。


S: [66x812 double]
    P0: [66x1 double]
    P1: [66x1 double]
    P2: [66x1 double]
function [newnode,nodenum]=nodeeit();
aa=32;
x=(-1+1/aa:2/aa:1-1/aa);
y=(-1+1/aa:2/aa:1-1/aa);
for i=1:aa
    for j=1:aa
        node(aa*(i-1)+j,1)=x(i);
        node(aa*(i-1)+j,2)=y(j);
    end
end
j=1;
for i=1:aa*aa
    if node(i,1)^2+node(i,2)^2<=1
        newnode(j,1)=node(i,1);
        newnode(j,2)=node(i,2);
        nodenum(j)=i;
        j=j+1;
    end
end
就这些资料。。。