1148阅读
2回复

[求助]关于FDTD一维的一个C程序大家帮忙修改一下不能运行

UID ：95999

注册：2012-06-22
登录：2014-05-10
发帖：46
等级：仿真新人

0楼发表于: 2013-07-04 15:04:35

#include<math.h>
#include<stdlib.h>
#include<stdio.h>
#define KE 200
main()
{
    float ex[KE],hy[KE];
    int n,k,kc,ke,NSTEPS;
    float T;
    float t0,spread,pulse;
    FILE *fp,*fopen();
    for(k=1;k<KE;k++)
    {ex[k]=0;
    hy[k]=0;}
    kc=KE/2;
    t0=40.0;
    spread=12;
    T=0;
    NSTEPS=1;
    while(NSTEPS>0)
    {
        printf("NSTEPS-->");
        scanf("%d",&NSTEPS);
        printf("%d\n",NSTEPS);
        n=0;
        for(n=1;n<=NSTEPS;n++)
        {T=T+1;
        for(k=1;k<KE;k++)
        {ex[k]=ex[k]+0.5(hy[k-1]-hy[k]);}
        pulse=exp(-0.5*(pow((t0-T)/spread,2.0)));
        ex[kc]=pulse;
        printf("%5.1f  %6.2f\n",t0-T,ex[kc]);}
    for(k=0;k<KE-1;k++)
    {hy[k]=hy[k]+0.5*(ex[k]-ex[k+1]);}
    }
        for(k=1;k<KE;k++)
        {printf("%3d   %6.2f   %6.2f\n",k,ex[k],hy[k]);}
        fp=fopen("Ex","w");
            for(k=1;k<KE;k++)
            {fprintf(fp,"  %6.2f\n",ex[k]);}
            fclose(fp);
            fp=fopen("Hy","w");
                for(k=1;k<KE;k++)
                {fprintf(f ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线nqx2012

UID ：110358

注册：2013-08-04
登录：2015-01-25
发帖：24
等级：仿真新人

1楼发表于: 2013-09-03 08:06:37

你是哪里出的问题？

离线yaomei18

UID ：1577

注册：2007-04-04
登录：2025-04-28
发帖：221
等级：仿真三级

2楼发表于: 2013-11-16 10:22:07

我上传了一个C++一维FDTD你可以看看

发帖回复