draw_weiyi.m

这是用于工程测量进行大坝变形分析的软件

for i=1:No                                     
     for j=1:size(u,1) 
        if j~=1 
           total_distchange_angel=alfa(u(j,i),v(j,i));%每个点的总位移的方位角
           project_distchange(j,i)=distchange(j,i)*cos(alfa(u(j,i),v(j,i))-str2num(get(edit1,'String'))/180*pi);%每个点在所求方向上的总位移
        end
    end
end  
for j=1:time
    fprintf(fit2,'观测日期：%5d年%3d月%3d日\n',myyear(j),mymonth(j),myday(j));
    for i=1:No
        fprintf(fit2,'z%1d   %5.2f 毫米\n',i,project_distchange(j,i)*1000);
    end
end
fig3=figure;
fig3_button1=uicontrol(fig3, 'Style',' push ', 'Position',[10 5 70 25],'String' ,'上翻','CallBack','fig3_pre');
fig3_button2=uicontrol(fig3, 'Style',' push ', 'Position',[110 5 70 25],'String' ,'下翻','CallBack','fig3_next');
begin_year_month=[str2num(get(edit2,'string'))  str2num(get(edit4,'string'))];
end_year_month=[str2num(get(edit3,'string'))  str2num(get(edit5,'string'))];
[begin_num_want_to_nihe,end_num_want_to_nihe]=get_nihe_num(begin_year_month,end_year_month,myyear,mymonth,time);
draw_weiyi_current_num=end_num_want_to_nihe;
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;