plotinvsolnode.m

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

function [fh,fp]=Plotinvsolnode(sol,g,H);

%Plotinvsolnode Plots the solution of 2D EIT problem in linear basis 
% Function []=Plotinvsolnode(sol,g,H);
% plots the solution of 2D EIT problem. Conductivity in linear basis.
%
% INPUT
%
% sol = the solution, values in the nodes
% g = node coordinates (x,y) (N x 2 matrix)
% H = connectivity matrix




[solN,solM]=size(sol);
[gN,gM]=size(g);
[HN,HM]=size(H);

for ii=1:solM
  view(2);
  hold on
  for ij=1:HN
    X=g(H(ij,:),1);
    Y=g(H(ij,:),2);
    Z=sol(H(ij,:),ii);
    fp=patch(X,Y,Z,Z);
    set(fp,'edgecolor','none');
    hold on
  end    
end
axis('square')