fig3_pre.m

来自「这是用于工程测量进行大坝变形分析的软件」· M 代码 · 共 21 行

if draw_weiyi_current_num~=begin_num_want_to_nihe
    draw_weiyi_current_num=draw_weiyi_current_num-1;
end
xmax=max(max(x));
xmin=min(min(x));
ymax=max(max(y));
ymin=min(min(y));
scale=(xmax-xmin)*50;
quiver([y(draw_weiyi_current_num,:) ymax+15 ymin],[x(draw_weiyi_current_num,:) xmin-10 xmin] ,[v(draw_weiyi_current_num,:) 0.001 0.01*sin(change_angel)].*scale,[u(draw_weiyi_current_num,:) 0 0.01*cos(change_angel)].*scale,0);
text(ymax+10,xmin-8,'一毫米');
text(ymin-10,xmin,'滑坡方向')
hold on
for i=1:No
    plot(y(draw_weiyi_current_num,i),x(draw_weiyi_current_num,i),'r*')
    text(1+y(draw_weiyi_current_num,i),x(draw_weiyi_current_num,i),[num2str(i) '号点'])
end
title(['这是' num2str(myyear(draw_weiyi_current_num)) '年' num2str(mymonth(draw_weiyi_current_num)) '月时期的位移图'])
xlabel('各点y坐标/米')
ylabel('各点x坐标/米')
axis equal
hold off;