%局部坐标与全局坐标的转换函数
function [R]=CoordTransform(x,y,BarLength)
l=(x(2)-x(1))/BarLength; m=(y(2)-y(1))/BarLength;
R=[l m 0 0;0 0 l m];
end