標(biāo)題: C語言計算凸輪源程序 [打印本頁]
作者: qhb    時間: 2019-7-17 17:07
標(biāo)題: C語言計算凸輪源程序
  1. #include <stdio.h>
  2. #include <math.h>
  3. int main()
  4. {  int  r0=20,  r1=10;
  5.   double  a,b,c,m,d,x,y,x1,y1,s,h;

  7.   for( a=0;a>=0&&a<=360;a=a+5)
  8.   {  if(a>=0&&a<=120)
  9.     {  s=30*a/180;    m=(atan2(30,((r0+s)*3.14))*180)/3.14;
  10.        b=30*sin(3.14*a/180)/3.14+(r0+s)*cos(3.14*a/180);  c=30*cos(3.14*a/180)/3.14-(r0+s)*sin(3.14*a/180);
  11.     x=(r0+s)*sin(3.14*a/180);
  12.        y=(r0+s)*cos(3.14*a/180);
  13.         d=b/sqrt(b*b+c*c);  h=-c/sqrt(b*b+c*c);
  14.         x1=x-r1*h;     y1=y-r1*d;
  15.     printf("%f\t%f\t%f\t%f\t%f\t%f\n",a,x,y,x1,y1,m);}
  16.      if(a>120&&a<=180)
  17. {  s=20;     m=0;
  18.       b=(r0+s)*cos(3.14*a/180);    c=-(r0+s)*sin(3.14*a/180);
  19.   //b=-320*a/(180*3.14)*sin(3.14+3.14*a/180)+(r0+s)*cos(3.14+3.14*a/180)
  20.   x=(r0+s)*sin(3.14*a/180);
  21.         y=(r0+s)*cos(3.14*a/180);
  22.         d=b/sqrt(b*b+c*c);  h=-c/sqrt(b*b+c*c);
  23.          x1=x-r1*h;     y1=y-r1*d;
  24.   printf("%f\t%f\t%f\t%f\t%f\t%f\n",a,x,y,x1,y1,m);}
  25.   if(a>180&&a<=225)
  26. {  s=20-160*(a*3.14/180-3.14)/(3.14*3.14);   m=(atan2((320*(a*3.14-180*3.14)),(r0+s)*180*3.14*3.14))*180/3.14;
  27.       b=-320*(a-180)/(180*3.14)*sin(3.14*a/180)+(r0+s)*cos(3.14*a/180);
  28.          c=-320*(a-180)/(180*3.14)*cos(3.14*a/180)-(r0+s)*sin(3.14*a/180);
  29.    x=(r0+s)*sin(3.14*a/180);
  30.          y=(r0+s)*cos(3.14*a/180);
  31.          d=b/sqrt(b*b+c*c);  h=-c/sqrt(b*b+c*c);
  32.          x1=x-r1*h;     y1=y-r1*d;
  33.    printf("%f\t%f\t%f\t%f\t%f\t%f\n",a,x,y,x1,y1,m);}
  34.   if(a>225&&a<=270)
  35. {  s=40*((315-a)/90)*((315-a)/90);
  36.       m=(atan2(160*(315-a),90*3.14*(r0+s)))*180/3.14;
  37.    b=-160*(a-315)/(90*3.14)*sin(3.14*(a-45)/180)+(r0+s)*cos(3.14*(a-45)/180);
  38.    c=-160*(a-315)/(90*3.14)*cos(3.14*(a-45)/180)-(r0+s)*sin(3.14*(a-45)/180);
  39.    x=(r0+s)*sin(3.14*a/180);
  40.          y=(r0+s)*cos(3.14*a/180);
  41.          d=b/sqrt(b*b+c*c);  h=-c/sqrt(b*b+c*c);
  42.          x1=x-r1*h;     y1=y-r1*d;
  43.           printf("%f\t%f\t%f\t%f\t%f\t%f\n",a,x,y,x1,y1,m);}
  44.   if(a>270&&a<=360)
  45. {  s=0;  m=0;
  46.       b=(r0+s)*cos(3.14*a/180);
  47.    c=-(r0+s)*sin(3.14*a/180);
  48.    x=(r0+s)*sin(3.14*a/180);
  49.           y=(r0+s)*cos(3.14*a/180);
  50.          d=b/sqrt(b*b+c*c);  h=-c/sqrt(b*b+c*c);
  51.          x1=x-r1*h;     y1=y-r1*d;
  52.     printf("%f\t%f\t%f\t%f\t%f\t%f\n",a,x,y,x1,y1,m);}
  53.   }
  54.   return 0;
  55.   
  56. }

復(fù)制代碼






歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1