/* 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
文章浏览阅读331次。第一部分:准备工作1 安装虚拟机2 安装centos73 安装JDK以上三步是准备工作,至此已经完成一台已安装JDK的主机第二部分:准备3台虚拟机以下所有工作最好都在root权限下操作1 克隆上面已经有一台虚拟机了,现在对master进行克隆,克隆出另外2台子机;1.1 进行克隆21.2 下一步1.3 下一步1.4 下一步1.5 根据子机需要,命名和安装路径1.6 ..._创建一个hadoop项目
文章浏览阅读1.7k次。心脏滴血漏洞HeartBleed CVE-2014-0160 是由heartbeat功能引入的,本文从深入码层面的分析该漏洞产生的原因_heartbleed代码分析
文章浏览阅读1.4k次。前言ofd是国家文档标准,其对标的文档格式是pdf。ofd文档是容器格式文件,ofd其实就是压缩包。将ofd文件后缀改为.zip,解压后可看到文件包含的内容。ofd文件分析工具下载:点我下载。ofd文件解压后,可以看到如下内容: 对于xml文件,可以用文本工具查看。但是对于印章文件(Seal.esl)、签名文件(SignedValue.dat)就无法查看其内容了。本人开发一款ofd内容查看器,..._signedvalue.dat
文章浏览阅读1.8w次,点赞29次,收藏313次。整体系统设计本设计主要是对ADC和DAC的使用,主要实现功能流程为:首先通过串口向FPGA发送控制信号,控制DAC芯片tlv5618进行DA装换,转换的数据存在ROM中,转换开始时读取ROM中数据进行读取转换。其次用按键控制adc128s052进行模数转换100次,模数转换数据存储到FIFO中,再从FIFO中读取数据通过串口输出显示在pc上。其整体系统框图如下:图1:FPGA数据采集系统框图从图中可以看出,该系统主要包括9个模块:串口接收模块、按键消抖模块、按键控制模块、ROM模块、D.._基于fpga的信息采集
文章浏览阅读2.5w次。1.背景错误信息:-- [http-nio-9904-exec-5] o.s.c.n.z.filters.post.SendErrorFilter : Error during filteringcom.netflix.zuul.exception.ZuulException: Forwarding error at org.springframework.cloud..._com.netflix.zuul.exception.zuulexception
文章浏览阅读358次。1.介绍图的相关概念 图是由顶点的有穷非空集和一个描述顶点之间关系-边(或者弧)的集合组成。通常,图中的数据元素被称为顶点,顶点间的关系用边表示,图通常用字母G表示,图的顶点通常用字母V表示,所以图可以定义为: G=(V,E)其中,V(G)是图中顶点的有穷非空集合,E(G)是V(G)中顶点的边的有穷集合1.1 无向图:图中任意两个顶点构成的边是没有方向的1.2 有向图:图中..._给定一个邻接矩阵未必能够造出一个图
文章浏览阅读321次。(十二)、WDS服务器安装通过前面的测试我们会发现,每次安装的时候需要加域光盘映像,这是一个比较麻烦的事情,试想一个上万个的公司,你天天带着一个光盘与光驱去给别人装系统,这将是一个多么痛苦的事情啊,有什么方法可以解决这个问题了?答案是肯定的,下面我们就来简单说一下。WDS服务器,它是Windows自带的一个免费的基于系统本身角色的一个功能,它主要提供一种简单、安全的通过网络快速、远程将Window..._doc server2012上通过wds+mdt无人值守部署win11系统.doc
文章浏览阅读219次。python–xlrd/xlwt/xlutilsxlrd只能读取,不能改,支持 xlsx和xls 格式xlwt只能改,不能读xlwt只能保存为.xls格式xlutils能将xlrd.Book转为xlwt.Workbook,从而得以在现有xls的基础上修改数据,并创建一个新的xls,实现修改xlrd打开文件import xlrdexcel=xlrd.open_workbook('E:/test.xlsx') 返回值为xlrd.book.Book对象,不能修改获取sheett_xlutils模块可以读xlsx吗
文章浏览阅读8.2w次,点赞267次,收藏656次。运行Selenium出现'WebDriver' object has no attribute 'find_element_by_id'或AttributeError: 'WebDriver' object has no attribute 'find_element_by_xpath'等定位元素代码错误,是因为selenium更新到了新的版本,以前的一些语法经过改动。..............._unresolved attribute reference 'find_element_by_id' for class 'webdriver
文章浏览阅读198次。一:模态窗口//父页面JSwindow.showModalDialog(ifrmehref, window, 'dialogWidth:550px;dialogHeight:150px;help:no;resizable:no;status:no');//子页面获取父页面DOM对象//window.showModalDialog的DOM对象var v=parentWin..._jquery获取父window下的dom对象
文章浏览阅读1.7w次,点赞15次,收藏129次。算法(algorithm)是解决一系列问题的清晰指令,也就是,能对一定规范的输入,在有限的时间内获得所要求的输出。 简单来说,算法就是解决一个问题的具体方法和步骤。算法是程序的灵 魂。二、算法的特征1.可行性 算法中执行的任何计算步骤都可以分解为基本可执行的操作步,即每个计算步都可以在有限时间里完成(也称之为有效性) 算法的每一步都要有确切的意义,不能有二义性。例如“增加x的值”,并没有说增加多少,计算机就无法执行明确的运算。 _算法
文章浏览阅读1.5k次,点赞18次,收藏26次。网络安全的标准和规范是网络安全领域的重要组成部分。它们为网络安全提供了技术依据,规定了网络安全的技术要求和操作方式,帮助我们构建安全的网络环境。下面,我们将详细介绍一些主要的网络安全标准和规范,以及它们在实际操作中的应用。_网络安全标准规范