coord_for_elec.m

来自「MATLAB二维电阻抗断层成像算法!用于医学成像,里面包括有限元剖分正问题,及反」· M 代码 · 共 33 行

function [X,Y] = coord_for_elec(height,no,r)
% 
%
%

elecA = height*no;
p = 2*pi*r;
otherA = p - elecA;
wall = otherA/no;

elecangle = (height*pi)/(pi*r); 
ele = elecangle/3;
wallangle = (wall*pi)/(pi*r);

k = wallangle*ones(1,no*3+no);

k(1:4:no*3+no) = ele;
k(2:4:no*3+no) = ele;
k(3:4:no*3+no) = ele;

k = diag(k);

m = tril(ones(no*3+no));

l = m*k;
angle = sum(l,2);

[x,y] = pol2cart(angle,r);
s = size(x,1);
X = [x(s);x(1:s-1)];
s = size(y,1);
Y = [y(s);y(1:s-1)];