极坐标放样参数计算.rar_Matlab极坐标放样参数计算_放样_放样参数计算_极坐标放样

function[] = lofting() %% parameter setting n=4; % 'n' is the number of lofting points. If n has been changed, line 10 should also been changed. sheet=1; % The number of sheet which saves lofting points. Angle = cell(4,1); RP6 = [313744.881,499447.504]; % given point one RP7 = [313737.15,499474.883]; % given point two %% read lofting points [data,ID,all] = xlsread('lofting points.xlsx',sheet); % Take care of lofting points' position. % Both .xls and .xlsx can be used, so long as the filename extension is correctly written. %% judge reference point and back vision point [St,BV] = station_or_backvision(RP6,RP7,data); if St == RP6 xlswrite('lofting points.xlsx',{'Station is Point 6'},sheet,'F1'); else xlswrite('lofting points.xlsx',{'Station is Point 7'},sheet,'F1'); end %% input known ditection direction_0 = atan2(BV(:,2)-St(:,2),BV(:,1)-St(:,1)); %% calculate distances for k=1:n D_k = distance(St,data(k,:)); %eval(['D_',num2str(k),'=','D_k',';']); %eval([['D_',num2str(k)],'=sqrt(D2_k)']); %% calculate angles direction_k = atan2(data(k,2)-St(:,2),data(k,1)-St(:,1)); rad_k = direction_k-direction_0; if (rad_k<0) rad_k = rad_k+2*pi; end degree_k = {rad2degree(rad_k)}; %% write results into lofting points.xls Distance(k,1) = D_k; Angle(k,1) = degree_k; end xlswrite('lofting points.xlsx',{'放样距离/m','放样角度'},sheet,'D1'); xlswrite('lofting points.xlsx',Distance,sheet,'D2'); xlswrite('lofting points.xlsx',Angle,sheet,'E2'); %% Display the position of points figure plot(St(:,2),St(:,1),'r*'... ,'MarkerSize',10); hold on plot(BV(:,2),BV(:,1),'r+'... ,'MarkerSize',5) for k=1:n plot(data(k,2),data(k,1),'go'... ,'MarkerSize',5) end hold off legend('Station','Back Vision'); %% calculate sides i=1; for k=1:n-1 for p= k+1:n length(i,:) = distance(data(k,:),data(p,:)); side(i,:) = strcat(ID(k+1,1),',',ID(p+1,1)); i = i+1; end end %% write sides into lofting points.xls xlswrite('lofting points.xlsx',side,sheet,'A7'); xlswrite('lofting points.xlsx',length,sheet,'B7'); xlswrite('lofting points.xlsx',{'边','边长/m'},sheet,'A6');