C语言——一些处理脚本(2)——(待完善)

最新推荐文章于 2024-06-13 20:48:31 发布

病树前头

最新推荐文章于 2024-06-13 20:48:31 发布

阅读量358

点赞数 11

分类专栏： C/C++学习与运用文章标签： c语言开发语言

本文链接：https://blog.csdn.net/yin1331102028yin/article/details/127028141

版权

C/C++学习与运用专栏收录该内容

25 篇文章 0 订阅

订阅专栏

# include <stdio.h>
# include <math.h>
# include <string.h>
# include <stdlib.h>

# define N   40

double th0=2.182;
double ph0=1.658;
double l=4.0;

int main(void)
{
  int i;
  double x[N], y[N], z[N];
  double xm, ym, zm;
  FILE *fp;

  ym=-cos(ph0)*sin(th0)*cos(th0/2)+sin(th0/2)-cos(th0)*sin(th0/2);
  ym=ym*l;

  xm=(l*cos(ph0)*sin(th0)+ym*cos(th0/2))/sin(th0/2);
  xm=-xm;

  zm=l*l-xm*xm-l*l*cos(th0/2)*cos(th0/2)-2*xm*l*cos(th0/2)-ym*ym-l*l*sin(th0/2)*sin(th0/2)+2*ym*l*sin(th0/2);
  zm=sqrt(zm);

  for(i=0;i<N;i++)
  {
      x[i]=0;
      y[i]=i*4;
      z[i]=-y[i]*(ym-2*l*sin(th0/2))/zm;
  }
  fp=fopen("snp0","w");
  for(i=0;i<N;i++)
  { 
    fprintf(fp,"%.3f    %.3f    %.3f\n",x[i],y[i],z[i]); 
  }
  fclose(fp);
  return 0;
}

# include <stdio.h>
# include <math.h>
# include <string.h>
# include <stdlib.h>

# define N   40

double th0=2.182;
double ph0=1.658;
double l=4.0;

int main(void)
{
  int i;
  double x[N], y[N], z[N];
  double xm, ym, zm;
  FILE *fp;

  ym=-cos(ph0)*sin(th0)*cos(th0/2)+sin(th0/2)-cos(th0)*sin(th0/2);
  ym=ym*l;

  xm=(l*cos(ph0)*sin(th0)+ym*cos(th0/2))/sin(th0/2);
  xm=-xm;

  zm=l*l-xm*xm-l*l*cos(th0/2)*cos(th0/2)-2*xm*l*cos(th0/2)-ym*ym-l*l*sin(th0/2)*sin(th0/2)+2*ym*l*sin(th0/2);
  zm=sqrt(zm);

  for(i=0;i<N;i++)
  {
      x[i]=0;
      y[i]=i*4;
      z[i]=y[i]/(ym-2*l*sin(th0/2))*zm;
  }
  fp=fopen("snp0","w");
  for(i=0;i<N;i++)
  { 
    fprintf(fp,"%.3f    %.3f    %.3f\n",x[i],y[i],z[i]); 
  }
  fclose(fp);
  return 0;
}