draw_nihe.m

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

fprintf(fit2,['从' get(edit2,'string') '年' get(edit4,'string')  '月到' get(edit3,'string')  '年' get(edit5,'string') '月的拟合成果\n']);
fsize=[500  300];
figpos=[100 100 fsize(1:2)];
fig2=figure;
nihe_number=str2num(get(edit6,'string'));
if  get(radio2,'value')==1 & (length(nihe_number)~=1 | nihe_number==0)
	fig2_button1=uicontrol(fig2, 'Style',' push ', 'Position',[10 5 70 25],'String' ,'上翻','CallBack','fig2_pre');
	fig2_button2=uicontrol(fig2, 'Style',' push ', 'Position',[110 5 70 25],'String' ,'下翻','CallBack','fig2_next');
end
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);
x0=month_time([begin_num_want_to_nihe:end_num_want_to_nihe]);
for i=1:time                                     
	z_total_distchange_angel=alfa(z_u(i),z_v(i));%每个点的总位移的方位角
	z_distchange=sqrt(z_u(i)^2+z_v(i)^2);%每期整体总位移
	z_project_distchange(i)=z_distchange*cos(abs(alfa(z_u(i),z_v(i))-str2num(get(edit1,'String'))/180*pi));%整体在所求方向上的位移
end  
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  
if get(radio1,'value')==1
	y0=z_project_distchange([begin_num_want_to_nihe:end_num_want_to_nihe]); %把第i个点的位移量赋给y
    if get(radio1_1,'value')==1%用直线拟合====================================
        linefunction(x0,y0);
    elseif get(radio1_2,'value')==1%==========================用样条拟合==================== 
        before_cigema_squre=str2num(get(edit7,'string'));
        polyfunction(x0,y0,before_cigema_squre);
    end
elseif get(radio2,'value')==1%各点分别拟合
    draw_nihe_current_point=1;   
    if get(radio2_1,'value')==1%用直线拟合==========================
        if nihe_number==0
                y0=project_distchange([begin_num_want_to_nihe:end_num_want_to_nihe],draw_nihe_current_point); %把第i个点的位移量赋给y
                linefunction(x0,y0,draw_nihe_current_point);
        elseif length(nihe_number)==1
                y0=project_distchange([begin_num_want_to_nihe:end_num_want_to_nihe],nihe_number); %把第i个点的位移量赋给y
                linefunction(x0,y0,nihe_number(draw_nihe_current_point));          
        elseif length(nihe_number)>1
                y0=project_distchange([begin_num_want_to_nihe:end_num_want_to_nihe],nihe_number(draw_nihe_current_point)); 
                linefunction(x0,y0,nihe_number(draw_nihe_current_point));
        end
    elseif get(radio2_2,'value')==1%==========================用样条拟合======================================================
        before_cigema_squre=str2num(get(edit8,'string'));
        y0=project_distchange([begin_num_want_to_nihe:end_num_want_to_nihe],draw_nihe_current_point);%把第i个点的位移量赋给y
        if nihe_number==0
                polyfunction(x0,y0,before_cigema_squre,draw_nihe_current_point);
        elseif length(nihe_number)==1|length(nihe_number)>1
                polyfunction(x0,y0,before_cigema_squre,nihe_number(draw_nihe_current_point));
        end
    end
end