附合导线计算程序设计

1、.5 附合导线计算程序设计#include<stdio.h>#define PI 3.1415926#include<math.h>#define LEN sizeof(struct node)#include<stdlib.h>struct zuobiaodouble x;double y;struct angleint degree;int minute;float second;struct nodedouble x;double y;struct angle ang1;struct angle fwj;double s;double dx;doubl

2、e dy;struct node *next;double jtod(struct angle jiao)double degree;degree=jiao.degree;degree+=(jiao.minute/60.0);degree+=(jiao.second/3600.0);return degree;struct angle dtoj(double jiaodu)struct angle fwj;fwj.degree=floor(jiaodu);fwj.minute =floor(jiaodu-floor(jiaodu)*60); fwj.second =(jiaodu-floor(

3、jiaodu)*60-fwj.minute)*60;.return fwj;double jtoh(struct angle jiaodu)double hu;hu=jtod(jiaodu);hu=hu*PI/180;return hu;main()FILE *fp;struct zuobiao zb0,zbN;struct node *head,*p1,*p2;struct angle angle1,f0,fB,fN;int n,i;double t,fBsecond,jiaodu;double fX=0,fY=0,K=0,S=0;fB.degree=0;fB.minute=0;fB.sec

4、ond=0;if(fp=fopen("d:prodata2.txt","r")=NULL)printf("cannot open this filen");fscanf(fp,"%d",&n);fscanf(fp,"%lf%lf",&zb0.x,&zb0.y);fscanf(fp,"%d%d%f",&f0.degree,&f0.minute,&f0.second);fscanf(fp,"%lf%lf",&am

5、p;zbN.x,&zbN.y);fscanf(fp,"%d%d%f",&fN.degree,&fN.minute,&fN.second);fB.degree=f0.degree;fB.minute=f0.minute;fB.second=f0.second;fBsecond=jtod(fB);for(i=0;i<n;i+)fscanf(fp,"%d%d%f",&angle1.degree,&angle1.minute,&angle1.second);fBsecond+=jtod(angle1)

6、;fBsecond-=180;if(fBsecond>360) fBsecond-=360;if(fBsecond<0) fBsecond+=360;fBsecond-=jtod(fN);fBsecond*=3600;printf("%lfn",fBsecond);fBsecond/=-n;rewind(fp);for(i=0;i<11;i+) fscanf(fp,"%lf",&t);.head=p1=p2=(struct node *)malloc(LEN); /fscanf(fp,"%d%d%f",&am

7、p;p1->fwj.degree ,&p1->fwj.minute ,&p1->fwj.second ); for(i=0;i<n;i+)fscanf(fp,"%d%d%f",&p1->ang1.degree ,&p1->ang1.minute ,&p1->ang1.second ); p1->ang1.second +=fBsecond;if(p1->ang1.second >=60)p1->ang1.second -=60;p1->ang1.minute +=

8、1;if(p1->ang1.minute >=60)p1->ang1.minute -=60;p1->ang1.degree +=1;if(p1->ang1.second <0)p1->ang1.second +=60;p1->ang1.minute -=1;if(p1->ang1.minute <0)p1->ang1.minute +=60;p1->ang1.degree -=1;p2->next =p1;p2=p1;p1=(struct node *)malloc(LEN);p2->next=NULL;p1

9、=head;i=0;dofscanf(fp,"%lf",&p1->s);printf("s%d=%.3lf mn",i+1,p1->s);S+=p1->s;i+;p1=p1->next ;while(p1->next!=NULL);p1->s=0;printf("S=%8.3lf mn",S);.p1=head;/p1->dx=p1->s*cos(jtoh(p1->fwj );/p1->dy=p1->s*sin(jtoh(p1->fwj );/p1=p1-

10、>next;jiaodu=jtod(f0);dojiaodu+=jtod(p1->ang1);jiaodu-=180;if(jiaodu>360) jiaodu-=360;if(jiaodu<0) jiaodu+=360;p1->fwj =dtoj(jiaodu);p1->dx=p1->s*cos(jtoh(p1->fwj);p1->dy=p1->s*sin(jtoh(p1->fwj);p1=p1->next ;while(p1->next !=NULL);p1=head;fX=zb0.x;fY=zb0.y;dofX

11、+=p1->dx;fY+=p1->dy;p1=p1->next ;while(p1->next !=NULL);fX-=zbN.x;fY-=zbN.y;K=sqrt(fX*fX+fY*fY);K/=S;K=1/K;p1=head;if(K>4000)do/ p2=p1->next ;p1->dx+=(-fX/S*p1->s);/ p2->x=p1->x+p2->dx;p1->dy+=(-fY/S*p1->s);/ p2->y=p1->y+p1->dy;p1=p1->next ;while(p1

12、->next !=NULL);p1=head;p1->x=zb0.x+p1->dx;p1->y=zb0.y+p1->dy;.dop2=p1->next;p2->x=p1->x+p2->dx;p2->y=p1->y+p2->dy;p1=p1->next ;while(p2->next !=NULL);else printf("the error is out of the limite!"),exit(0);printf("fX=%6.3lfmfY=%6.3lfmfS=%6.3lfm

13、K=1/ %dn",fX,fY ,sqrt(fX*fX+fY*fY),(int)K); p1=head;doprintf("%8.3lf %8.3lfn",p1->x,p1->y); p1=p1->next;while(p1->next !=NULL);fclose(fp);if(fp=fopen("d:procoordinate.dat","w")=NULL)printf("cannot open this file!n");fprintf(fp,"fx=%.3lf mn",fX);fprintf(fp,"fy=%.3lf mn",fY);fprintf(fp,"fs=%.3lf mn",sqrt(fX*fX+fY*fY);fprintf(fp,"K=1 / %dn",(int)K);fprintf(fp,"Th