| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 384 人关注过本帖
标题:程序运行结果没出数据,是不是因为没建接受数据的文件夹?
只看楼主 加入收藏
wzy88127
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2013-6-19
结帖率:0
收藏
已结贴  问题点数:20 回复次数:1 
程序运行结果没出数据,是不是因为没建接受数据的文件夹?
/*
  dA/dt-s*dA/dx=a*d2A/dx2+b*A-(c*|A|2-|A|4+d*|B|2)A
  dB/dt-s*dB/dx=a*d2B/dx2+b*B-(c*|B|2-|B|4+d*|A|2)B
  
  (subcritical)ai=bi=ci=di=ei=0
  
  */
  #include<stdio.h>
  #include<stdlib.h>
  #include<math.h>
  
  
  #define PI 3.14159265         /*shokisettei*/
  #define OMG 2.
  #define NL 200000
  #define NT 1000
  #define NS 100
  
  #define AR 1.0
  #define AI 0.0
  #define BR 1.0
  #define BI 0.0
  #define CR 1.0
  #define CI 2.0
  #define DR 1.0
  #define DI 0.0
  #define ER 1.0
  #define EI 0.0
  
  #define MR -0.06
  #define MI 0.0
  #define NR -0.06
  #define NI 0.0
  #define S 0.6
  #define DEE 2.0
  
  int M,s=200000;
  double  EPS;
  double u0,v0,w0,z0;
  double u[1000],v[1000],w[1000],z[1000];
  double u2[1000],v2[1000],w2[1000],z2[1000];
  double DX=2.0/10.0;
  double DT=1.0/100.0;
  double D2,D,K,DX2;
  
  #define L 20
  
  int i,l,m,n=0;
  char fname[6][20];
  FILE  *f[6],*fp;
  
  double keisan();
  double kyoukai();
  void   shokijouken();
  
  main (int argc,char*argv[])   /*input EPS finename*/
  {
  
    if((argc!=2)&&(argc!=3)){
      fprintf(stderr,"usage:%sEPS[filename]/n",argv[0]);
      exit(1);
    }
   
    EPS=atof(argv[1]);
    M=(L*5)+1;
    sprintf(fname[0],"g%2.2g1",EPS);
    sprintf(fname[1],"g%2.2g2",EPS);
    sprintf(fname[2],"f%2.2g1",EPS);
    sprintf(fname[3],"f%2.2g2",EPS);
    sprintf(fname[4],"f%2.2g",EPS);
    sprintf(fname[5],"s%2.2g",EPS);
   

    for(i=0;i<=4;i++){
        if((f[i]=fopen(fname[i],"w"))==NULL){  /*file open*/
           printf("Cannot open file %s¥n",fname[i]);
           exit(1);
               }
        }
    shokijouken(argc,argv[2]);
   
    for(n=1;n<=NL;n++){/*keisann 翻译成汉语是计算*/

    u2[0]=keisan(u0,u[0],u[1],v0,v[0],v[1],w[0],z[0],-1,-1);
     v2[0]=keisan(v0,v[0],v[1],u0,u[0],u[1],z[0],w[0],1,-1);
     w2[0]=keisan(w0,w[0],w[1],z0,z[0],z[1],u[0],v[0],-1,1);
     z2[0]=keisan(z0,z[0],z[1],w0,w[0],w[1],v[0],u[0],1,1);
    for (i=1;i<=M-1;i++){
   
     u2[i]=keisan(u[i-1],u[i],u[i+1],v[i-1],v[i],v[i+1],w[i],z[i],-1,-1);
     v2[i]=keisan(v[i-1],v[i],v[i+1],u[i-1],u[i],u[i+1],z[i],w[i],1,-1);
     w2[i]=keisan(w[i-1],w[i],w[i+1],z[i-1],z[i],z[i+1],u[i],v[i],-1,1);
     z2[i]=keisan(z[i-1],z[i],z[i+1],w[i-1],w[i],w[i+1],v[i],u[i],1,1);
    }
   
    for(i=0;i<=M-1;i++)
       u[i]=u2[i];
       v[i]=v2[i];
       w[i]=w2[i];
       z[i]=z2[i];
    }
    /*kyoukaijoukenn翻译成汉语是边界条件*/
    u[0]=kyoukai(u[1],u0,v[1],v0,w[1],w0,z[1],z0,-1);
    v[0]=kyoukai(v[1],v0,u[1],u0,z[1],z0,w[1],w0,1);
    w[0]=kyoukai(w[1],w0,z[1],z0,u[1],u0,v[1],v0,1);
    z[0]=kyoukai(z[1],z0,w[1],w0,v[1],v0,u[1],u0,-1);
   
    u[M-1]=kyoukai(u[M],u[M-2],v[M],v[M-2],w[M],w[M-2],z[M],z[M-2],-1);
    v[M-1]=kyoukai(v[M],v[M-2],u[M],u[M-2],z[M],z[M-2],w[M],w[M-2],1);
    w[M-1]=kyoukai(w[M],w[M-2],z[M],z[M-2],u[M],u[M-2],v[M],v[M-2],1);
    z[M-1]=kyoukai(z[M],z[M-2],w[M],w[M-2],v[M],v[M-2],u[M],u[M-2],-1);

    if(n%NS==0){/*file ni kakikomi翻译成汉语是在这个文件里写进去*/
              l=(M-1)/2;
              fprintf(f[0],"%d%e/n",n,sqrt(u[l]*u[l]+v[l]*v[l]));
              fprintf(f[1],"%d%e/n",n,sqrt(w[l]*w[l]+z[l]*z[l]));
    }
    if(n%NT==0)
       m=(M-1)/100;                                             
       K=2*PI/sqrt(2);
       for(i=0;i<=M-1;i+=m){         
          fprintf(f[2],"%f%d%e/n",i*DX,n,sqrt(u[i]*u[i]+v[i]*v[i]));
          fprintf(f[3],"%f%d%e/n",i*DX,n,sqrt(w[i]*w[i]+z[i]*z[i]));
          fprintf(f[4],"%f%d%e/n",i*DX,n,2*u[i]*cos(K*DX*i-OMG*DT*n)
                                           -2*v[i]*sin(K*DX*i-OMG*DT*n)
                                           +2*w[i]*cos(K*DX*i+OMG*DT*n)
                                           +2*z[i]*sin(K*DX*i+OMG*DT*n));
          {
          fprintf(f[2],"/n");
          fprintf(f[3],"/n");
          fprintf(f[4],"/n");
       }
    }
    for(i=0;i<=4;i++) fclose(f[i]);
   
   
    shokijouken(argc,argv[2]);
   
    for(n=1;n<=s;n++){/*keisann*/

       u2[0]=keisan(u0,u[0],u[1],v0,v[0],v[1],w[0],z[0],-1,-1);
       v2[0]=keisan(v0,v[0],v[1],u0,u[0],u[1],z[0],w[0],1,-1);
       w2[0]=keisan(w0,w[0],w[1],z0,z[0],z[1],u[0],v[0],-1,1);
       z2[0]=keisan(z0,z[0],z[1],w0,w[0],w[1],v[0],u[0],1,1);
   
   for (i=1;i<=M-1;i++){
      u2[i]=keisan(u[i-1],u[i],u[i+1],v[i-1],v[i],v[i+1],w[i],z[i],-1,-1);
      v2[i]=keisan(v[i-1],v[i],v[i+1],u[i-1],u[i],u[i+1],z[i],w[i],1,-1);
      w2[i]=keisan(w[i-1],w[i],w[i+1],z[i-1],z[i],z[i+1],u[i],v[i],-1,1);
      z2[i]=keisan(z[i-1],z[i],z[i+1],w[i-1],w[i],w[i+1],v[i],u[i],1,1);
   }
    for(i=0;i<=M-1;i++){
        u[i]=u2[i];
        v[i]=v2[i];
        w[i]=w2[i];
        z[i]=z2[i];
    }
  /*kyoukaijoukenn边界条件*/
    u[0]=kyoukai(u[1],u0,v[1],v0,w[1],w0,z[1],z0,-1);
    v[0]=kyoukai(v[1],v0,u[1],u0,z[1],z0,w[1],w0,1);
    w[0]=kyoukai(w[1],w0,z[1],z0,u[1],u0,v[1],v0,1);
    z[0]=kyoukai(z[1],z0,w[1],w0,v[1],v0,u[1],u0,-1);

     u[M-1]=kyoukai(u[M],u[M-2],v[M],v[M-2],w[M],w[M-2],z[M],z[M-2],-1);
     v[M-1]=kyoukai(v[M],v[M-2],u[M],u[M-2],z[M],z[M-2],w[M],w[M-2],1);
     w[M-1]=kyoukai(w[M],w[M-2],z[M],z[M-2],u[M],u[M-2],v[M],v[M-2],1);
     z[M-1]=kyoukai(z[M],z[M-2],w[M],w[M-2],v[M],v[M-2],u[M],u[M-2],-1);
}
  
   if((f[5]=fopen(fname[5],"wb"))==NULL){
     printf("Cannot open file %s/n",fname[5]);
     exit(1);
  }
  if((fwrite(&u0,sizeof(u0),1,f[5]))!=1){
     fprintf(stderr,"fwrite error%s/n",fname[5]);
     exit(1);
   }
  if(fwrite(u,sizeof(double),M+1,f[5])!=102){
     fprintf(stderr,"fwrite error%s/n",fname[5]);
     exit(1);
   }
  if((fwrite(&v0,sizeof(v0),1,f[5]))!=1){
     fprintf(stderr,"fwrite error%s/n",fname[5]);
     exit(1);
   }
  if((fwrite(v,sizeof(double),M+1,f[5]))!=102){
     fprintf(stderr,"fwrite error%s/n",fname[5]);
     exit(1);
   }
  if((fwrite(&w0,sizeof(w0),1,f[5]))!=1){
     fprintf(stderr,"fwrite error%s/n",fname[5]);
     exit(1);
   }
  if((fwrite(w,sizeof(w0),M+1,f[5]))!=102){
     fprintf(stderr,"fwrite error%s/n",fname[5]);
     exit(1);
   }
  if((fwrite(&z0,sizeof(z0),1,f[5]))!=1){
     fprintf(stderr,"fwrite error%s/n",fname[5]);
     exit(1);
   }
  if((fwrite(z,sizeof(double),M+1,f[5]))!=102){
     fprintf(stderr,"fwrite error%s/n",fname[5]);
     exit(1);
   }
  fclose(f[5]);
}
double u0;
void shokijouken(int argc,char*argv2)
{
 double u0;
 if (argc==2){
  
  u0=1.e-10*(1/(0.5*sqrt(PI)))*exp(-(1/(0.5*0.5))*(DX*(-1-M/2))*(DX*(-1-M/2)));
  

  v0=w0=z0=0;
  for(i=0;i<=M;i++)
    u[i]=1.e-10*(1/(0.5*sqrt(PI)))*exp(-(1/(0.5*0.5))*(DX*(i-M/2))*(DX*(i-M/2)));

  v[i]=w[i]=z[i]=0;
  }
else{/*argc==3*/
  if((fp=fopen(argv2,"rb"))==NULL){
    printf("Cannot open file %s/n",*argv2);
    exit(1);
  }
  if((fread(&u0,sizeof(u0),1,fp))!=1){
   fprintf(stderr,"fread error %s/n",argv2);
   exit(1);
  }
  if((fread(u,sizeof(double),M+1,fp))!=102){
    fprintf(stderr,"fread error %s/n",argv2);
    exit(1);
  }
  if((fread(&v0,sizeof(v0),1,fp))!=1){
   fprintf(stderr,"fread error %s/n",argv2);
   exit(1);
  }
  if((fread(v,sizeof(double),M+1,fp))!=102){
   fprintf(stderr,"fread error %s/n",argv2);
   exit(1);
  }
  if((fread(&w0,sizeof(w0),1,fp))!=1){
   fprintf(stderr,"fread error %s/n",argv2);
   exit(1);
  }
  if((fread(w,sizeof(double),M+1,fp))!=102){
   fprintf(stderr,"fread error %s/n",argv2);
   exit(1);
  }
  if((fread(&z0,sizeof(z0),1,fp))!=1){
   fprintf(stderr,"fread error %s/n",argv2);
   exit(1);
  }
   if((fread(z,sizeof(double),M+1,fp))!=102){
   fprintf(stderr,"fread error %s/n",argv2);
   exit(1);
  }
  fclose(fp);
 }
}
/*keisann计算*/
double keisan(double u_x,double u,double ux,double v_x,double v,double vx,double w,double z,int a, int b)
{
  D2=DT/(DX*DX);
  D=DT/(2*DX);

  return u+D2*(AR*(ux-2*u+u_x)+a*AI*(vx-2*v+v_x))+DT*(EPS*(BR*u+a*BI*v)+(CR*u+a*CI*v)*(u*u+v*v)-DEE*(DR*u+a*DI*v)*(w*w+z*z)-(ER*u+a*EI*v)*(u*u+v*v)*(u*u+v*v))+b*D*S*(ux-u_x);
}
 /*kyoukaijoukenn边界条件*/
double kyoukai(double u1,double u0,double v1,double v0,double w1,double w0,double z1,double z0,int a)
{
 DX2=2*DX;

 return (MR*(u1-u0)+a*MI*(v1-v0)+NR*(w1-w0)+a*NI*(z1-z0))/DX2;
 system ("pause");
}
搜索更多相关主题的帖子: include 文件夹 
2013-06-19 22:28
pauljames
Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19Rank: 19
等 级:千里冰封
威 望:9
帖 子:1555
专家分:10000
注 册:2011-5-8
收藏
得分:20 
带参数运行的,比如 程序名 a.txt

经常不在线不能及时回复短消息,如有c/单片机/运动控制/数据采集等方面的项目难题可加qq1921826084。
2013-06-22 05:54
快速回复:程序运行结果没出数据,是不是因为没建接受数据的文件夹?
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.020438 second(s), 9 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved