/* FDTD-1D-2.2.c 1D FDTD simulation of a lossy dielectric medium */
/* Simulation of a sinusoidal wave or a Gauss pulse transmiting in a lossy dielectric medium exsiting in free space */
/*8888888888888 Using a new formulation and using flux density 8888888888888*/
/* The Fourier Transform has been added */
/* 1D FDTD simulation of a frequency dependent material */
# include <math.h>
# include <stdlib.h>
# include <stdio.h>
# define KE 200 /* KE is the number of cells to be used */
# define NF 10 /* NF is the number of the frequence in fourier transform */
# define Pi 3.1415926
void main ()
{
double ex_low_m1, ex_low_m2,ex_high_m1,ex_high_m2;
double dx[KE],ix[KE],ex[KE],hy[KE],sx[KE];
double ga[KE],gb[KE],gc[KE];
int n, m,k, kc,NSTEPS;
int kstart_dielectric; /* 介质的左边界*/
int kend_dielectric; /* 介质的右边界*/
double T;
double t0, spread,pulse, pulse1, freq_in, pulse2; /*pulse1 is Gauss pulse, Pulse2 is a sinusoidal wave */
double epsilon_r,epsilon_0; /* 介质的相对介电常数 和 真空介电常数*/
double sigma; /* 介质的电导率 */
double ddx, dt; /* The FDTD Cell size and the time interval */
double freq[NF], arg[NF], ampn[NF][KE],phasen[NF][KE];
double real_part[NF][KE], imag_part[NF][KE];
double real_in[NF], imag_in[NF], amp_in[NF], phase_in[NF];
double mag[KE] ;
double tau, chil, del_exp;
FILE *fp ;
/* Initialize */
ddx=0.01; /* Set the cell size to 1cm */
dt=ddx/(2*3e8) ; /* calculate the time step */
epsilon_0=8.85419e-12;
kc=5; /* 激励源所引入的位置 */
t0= 50.0; /* Center of the incident pulse */
spread=10.0; /* Width of the incident pulse */
T=0;
NSTEPS=1;
ex_low_m1=0; /* 在n+1时存储左边界处前一个格点ex[1] 结果 */
ex_low_m2=0; /* 在n+2时存储左边界处前一个格点ex[1] 结果 */
ex_high_m1=0; /* 在n+1时存储右边界处前一个格点ex[KE-2] 结果 */
ex_high_m2=0; /* 在n+2时存储右边界处前一个格点ex[KE-2] 结果 */
for (k=0; k<KE; k++)
{ ga[k]=1.0; /* initialize to free space */
gb[k]=0.;
gc[k]=0; /* initialize to free space */
ex[k]=0.;
hy[k]=0.;
dx[k]=0.;
ix[k]=0;
sx[k]=0;
mag[k]=0;
for (m=0;m<=2;m++)
{ real_part[m][k]=0.; /* Real and imaginary parts of fourier transform */
imag_part[m][k]=0.;
ampn[m][k]=0.; /* Amplituds and phase of the fourier transform */
phasen[m][k]=0.;
}
}
for (m=0;m<=2;m++)
{ real_in[m]=0.; /* Fourier Trans. of input pulse */
imag_in[m]=0.;
}
/* Parameters for the Fourier Transform */
freq[0]=50.e6;
freq[1]=200.e6;
freq[2]=500.e6;
for (m=0;m<=2;m++)
{ arg[m]=2*Pi*freq[m]*dt;
printf("%2d %6.2f MHz %7.5f \n ", m, freq[m]*1e-6,arg[m]);
}
/* These parameters specify the input pulse */
// printf("input frequence (MHz) -->");
// scanf("%lf", &freq_in);
// freq_in=freq_in*1e6;
// printf("%8.0f \n",freq_in);
printf("Dielectric start at kstart_dielectric -->");
scanf("%d", &kstart_dielectric);
printf("Dielectric end at kend_dielectric -->");
scanf("%d", &kend_dielectric);
printf("Epsilon_r -->");
scanf("%lf", &epsilon_r);
printf("Conductivity -->");
scanf("%lf", &sigma);
printf("chil----->");
scanf("%lf", &chil);
tau=1000.; /* Make sure tau is >0. */
if (chil>0.0001)
{printf("tau (in microseconds)--->"); scanf("%lf",&tau); del_exp=exp(-dt/tau);}
printf ("%d %d %f %f %f %f \n", kstart_dielectric,kend_dielectric,epsilon_r, sigma,tau,chil);
tau=1.e-6*tau;
printf("del_exp= %f \n",del_exp) ;
for (k=kstart_dielectric; k<=kend_dielectric;k++)
{ ga[k]=1/(epsilon_r+sigma*dt/epsilon_0+chil*dt/tau);
gb[k]=sigma*dt/epsilon_0; /* initialize to dielectric medium */
gc[k]=chil*dt/tau;
}
for (k=0;k<KE;k++)
{ printf("%2d %4.2f %4.2f \n",k,ga[k],gb[k]); }
/* Main part of program */
while (NSTEPS>0)
{ printf("NSTEPS -->"); /* NSTEPS is the number of times that the Main loop has executed */
scanf("%d", &NSTEPS);
printf ("%d \n", NSTEPS);
printf("程序正在运行,请稍侯! \n");
for (n=1; n<=NSTEPS; n++)
{
T=T+1; /* T keeps track of the total number of times the main loop is executed */
/* Main FDTD Loop */
/* Calculate the Dx field */
for (k=1; k<KE; k++)
{ dx[k]=dx[k]+ 0.5*(hy[k-1]-hy[k]); }
/* Put a sinusoidal wave or the Gauss pulse at the cell kc */
pulse1=exp(-0.5*(pow((t0-T)/spread,2.0)));
// pulse2=sin(2*Pi*freq_in*dt*T);
pulse=pulse1; /* 脉冲的类型由pulse1和pulse2两个变量来选择 */
dx[kc]=dx[kc]+pulse;
printf("%5.1f %f %6.2f \n", T,pulse,dx[kc]);
/* Calculate the Ex from the Dx */
for (k=0;k<KE;k++)
{ ex[k]=ga[k]*(dx[k]-ix[k]-sx[k]);
ix[k]=ix[k]+gb[k]*ex[k];
sx[k]=del_exp*sx[k]+gc[k]*ex[k];
}
/* Calculate the fourier trasform of Ex */
for (k=0;k<KE;k++)
{ for (m=0;m<=2;m++)
{ real_part[m][k]=real_part[m][k]+cos(arg[m]*T)*ex[k];
imag_part[m][k]=imag_part[m][k]-sin(arg[m]*T)*ex[k];
}
}
/* Fourier Transform of the input pulse */
if (T<100)
{
for (m=0;m<=2;m++)
{ real_in[m]=real_in[m]+cos(arg[m]*T)*ex[10];
imag_in[m]=imag_in[m]-sin(arg[m]*T)*ex[10];
}
}
/* Absorbing boundary conditions */
ex[0]=ex_low_m2;
ex_low_m2=ex_low_m1;
ex_low_m1=ex[1];
ex[KE-1]=ex_high_m2;
ex_high_m2=ex_high_m1;
ex_high_m1=ex[KE-2];
/* Calculate the hy field */
for (k=0; k<KE-1;k++)
{ hy[k]=hy[k]+0.5*(ex[k]-ex[k+1]); }
}
/* End of the Main FDTD Loop */
/* At the end of the calculation, print out the Ex and Hy field */
for (k=0; k<KE; k++)
{ printf("%3d %6.2f %6.2f \n", k, ex[k], hy[k]);
}
/* Write the E field out to a file "Ex" */
fp=fopen("Ex.txt","w");
fprintf(fp, "T=%5.0f \n", T);
for (k=0;k<KE;k++)
{
fprintf (fp," %5d %6.2f %6.2f \n", k,dx[k],ex[k]);
}
fclose(fp);
/* Write the H field out to a file "Hy" */
fp=fopen("Hy.txt","w");
fprintf(fp, "T=%5.0f \n", T);
for (k=0;k<KE;k++)
{
fprintf (fp,"%5d %6.2f \n",k, hy[k]);
}
fclose(fp);
printf(" %5.0f \n", T);
/* calculate the amplitude and phase of each frequency */
/* Amplitude and the phase of the input pulse */
for (m=0; m<=2;m++)
{ amp_in[m]=sqrt(pow(imag_in[m],2.0)+pow(real_in[m],2.0));
phase_in[m]=atan2(imag_in[m],real_in[m]);
printf("%d Input pulse: %8.4f %8.4f %8.4f %7.2f \n", m,real_in[m],imag_in[m],amp_in[m],(180.0/Pi)*phase_in[m]);
for(k=0;k<KE;k++)
{ ampn[m][k]=(1.0/amp_in[m]*sqrt(pow(real_part[m][k],2.0)+pow(imag_part[m][k],2.0)));
phasen[m][k]=atan2(imag_part[m][k],real_part[m][k])-phase_in[m];
}
}
/* write the amplitude field out to a files "amp " */
fp=fopen ("amp0.txt","w");
for (k=0;k<KE;k++)
{ fprintf(fp,"%d %8.5f \n ", k,ampn[0][k] ); }
fclose(fp);
fp=fopen ("amp1.txt","w");
for (k=0;k<KE;k++)
{ fprintf(fp,"%d %8.5f \n ", k,ampn[1][k] ); }
fclose(fp);
fp=fopen ("amp2.txt","w");
for (k=0;k<KE;k++)
{ fprintf(fp,"%d %8.5f \n ", k,ampn[2][k] ); }
fclose(fp);
}
printf("程序运行完毕 \n");
}
运行之后出现如下错误
root@gallup-virtual-machine:/home/source# gcc -Wall li.c -o li
li.c:7:6: 警告: ‘main’的返回类型不是‘int’ [-Wmain]root@gallup-virtual-machine:/home/source#
解决方法:
1、 root@gallup-virtual-machine:/home/source# gcc li.c -lm -o li
文章浏览阅读3.2k次。本文研究全球与中国市场分布式光纤传感器的发展现状及未来发展趋势,分别从生产和消费的角度分析分布式光纤传感器的主要生产地区、主要消费地区以及主要的生产商。重点分析全球与中国市场的主要厂商产品特点、产品规格、不同规格产品的价格、产量、产值及全球和中国市场主要生产商的市场份额。主要生产商包括:FISO TechnologiesBrugg KabelSensor HighwayOmnisensAFL GlobalQinetiQ GroupLockheed MartinOSENSA Innovati_预计2026年中国分布式传感器市场规模有多大
文章浏览阅读1.1k次,点赞2次,收藏12次。常用组合逻辑电路结构——为IC设计的延时估计铺垫学习目的:估计模块间的delay,确保写的代码的timing 综合能给到多少HZ,以满足需求!_基4布斯算法代码
文章浏览阅读3.3k次,点赞3次,收藏5次。OpenAI Manager助手(基于SpringBoot和Vue)_chatgpt网页版
文章浏览阅读2.2k次。USACO自1992年举办,到目前为止已经举办了27届,目的是为了帮助美国信息学国家队选拔IOI的队员,目前逐渐发展为全球热门的线上赛事,成为美国大学申请条件下,含金量相当高的官方竞赛。USACO的比赛成绩可以助力计算机专业留学,越来越多的学生进入了康奈尔,麻省理工,普林斯顿,哈佛和耶鲁等大学,这些同学的共同点是他们都参加了美国计算机科学竞赛(USACO),并且取得过非常好的成绩。适合参赛人群USACO适合国内在读学生有意向申请美国大学的或者想锻炼自己编程能力的同学,高三学生也可以参加12月的第_usaco可以多次提交吗
文章浏览阅读394次。1.1 存储程序1.2 创建存储过程1.3 创建自定义函数1.3.1 示例1.4 自定义函数和存储过程的区别1.5 变量的使用1.6 定义条件和处理程序1.6.1 定义条件1.6.1.1 示例1.6.2 定义处理程序1.6.2.1 示例1.7 光标的使用1.7.1 声明光标1.7.2 打开光标1.7.3 使用光标1.7.4 关闭光标1.8 流程控制的使用1.8.1 IF语句1.8.2 CASE语句1.8.3 LOOP语句1.8.4 LEAVE语句1.8.5 ITERATE语句1.8.6 REPEAT语句。_mysql自定义函数和存储过程
文章浏览阅读188次。半导体二极管——集成电路最小组成单元。_本征半导体电流为0
文章浏览阅读2.8k次,点赞3次,收藏18次。游戏水面特效实现方式太多。咱们这边介绍的是一最简单的UV动画(无顶点位移),整个mesh由4个顶点构成。实现了水面效果(左图),不动代码稍微修改下参数和贴图可以实现岩浆效果(右图)。有要思路是1,uv按时间去做正弦波移动2,在1的基础上加个凹凸图混合uv3,在1、2的基础上加个水流方向4,加上对雾效的支持,如没必要请自行删除雾效代码(把包含fog的几行代码删除)S..._unity 岩浆shader
文章浏览阅读5k次。广义线性模型是线性模型的扩展,它通过连接函数建立响应变量的数学期望值与线性组合的预测变量之间的关系。广义线性模型拟合的形式为:其中g(μY)是条件均值的函数(称为连接函数)。另外,你可放松Y为正态分布的假设,改为Y 服从指数分布族中的一种分布即可。设定好连接函数和概率分布后,便可以通过最大似然估计的多次迭代推导出各参数值。在大部分情况下,线性模型就可以通过一系列连续型或类别型预测变量来预测正态分布的响应变量的工作。但是,有时候我们要进行非正态因变量的分析,例如:(1)类别型.._广义线性回归模型
文章浏览阅读69次。环境保护、 保护地球、 校园环保、垃圾分类、绿色家园、等网站的设计与制作。 总结了一些学生网页制作的经验:一般的网页需要融入以下知识点:div+css布局、浮动、定位、高级css、表格、表单及验证、js轮播图、音频 视频 Flash的应用、ul li、下拉导航栏、鼠标划过效果等知识点,网页的风格主题也很全面:如爱好、风景、校园、美食、动漫、游戏、咖啡、音乐、家乡、电影、名人、商城以及个人主页等主题,学生、新手可参考下方页面的布局和设计和HTML源码(有用点赞△) 一套A+的网_垃圾分类网页设计目标怎么写
文章浏览阅读614次,点赞7次,收藏11次。之前找到一个修改 exe 中 DLL地址 的方法, 不太好使,虽然能正确启动, 但无法改变 exe 的工作目录,这就影响了.Net 中很多获取 exe 执行目录来拼接的地址 ( 相对路径 ),比如 wwwroot 和 代码中相对目录还有一些复制到目录的普通文件 等等,它们的地址都会指向原来 exe 的目录, 而不是自定义的 “lib” 目录,根本原因就是没有修改 exe 的工作目录这次来搞一个启动程序,把 .net 的所有东西都放在一个文件夹,在文件夹同级的目录制作一个 exe._.net dll 全局目录
文章浏览阅读1.5k次。本文为转载,原博客地址:http://blog.csdn.net/hujingshuang/article/details/46910259简介 BRIEF是2010年的一篇名为《BRIEF:Binary Robust Independent Elementary Features》的文章中提出,BRIEF是对已检测到的特征点进行描述,它是一种二进制编码的描述子,摈弃了利用区域灰度..._breif description calculation 特征点
文章浏览阅读4.1k次,点赞21次,收藏79次。本文是《基于SpringBoot的房屋租赁管理系统》的配套原创说明文档,可以给应届毕业生提供格式撰写参考,也可以给开发类似系统的朋友们提供功能业务设计思路。_基于spring boot的房屋租赁系统论文