上机练习题6:数学物理方程的计算机求解和可视化.pptVIP

上机练习题6:数学物理方程的计算机求解和可视化.ppt

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
上机练习题6:数学物理方程的计算机求解和可视化

矿物三维周期结构——? 蛋白石 Page * Page * 光子晶体光纤 Page * 光子晶体的周期性结构 (a) 一维光子晶体 (b) 二维光子晶体 (c) 三维光子晶体 光子晶体——具有周期性结构的人工合成(或天然)材料。 光子晶体的电磁场本征值问题 Page * 根据Maxwell方程,光子晶体中电场满足: 平面波展开法:将E(r)展开为一系列平面波,解上述本征方程,本征值 为(ω/c)2, 为Hermitian算符,得到频率ω与波矢k之间的色散关系。 Page * 得到频谱空间中的本征方程: 一维光子晶体TE模式的本征方程: 对 和 作傅里叶变换: 取值范围 ,为z方向波矢,绘制出 与 之间的色散曲线。 Hermitian算符 本征值 一维光子晶体TE模式本证值问题 ?1. 求本征方程: 令 ,特征值为 ,特征矢量为 。 Fourier系数: 2.令 、 、 ,代入 中, 利用一个循环求出Q= 再使用:omega_c=eig(Q); omega_c=sort(sqrt(omega_c)); 3.求出 与 的,作出程序及图像。 Page * 求解思路 为 的函数 clear all d=2; %air gap a=10; %total periodicity eps_r=12.25; d_over_a=d/a; Max=50; %自定义 Q=zeros(2*Max+1) freq=[] for kz=-pi/a:pi/(10*a):pi/a for m=1:2*Max+1 for n=1:2*Max+1 M=m-Max; N=n-Max; kn=(1-1/eps_r)*d_over_a.*sinc(pi*(M-N)*d_over_a)+((M-N)==0)*1/eps_r; Q(m,n)=(2*pi*N/a+kz).^2.*kn; end end Page * 一维光子晶体TE模式本证值问题 omega_c=eig(Q); omega_c=sort(sqrt(omega_c)); freq=[freq;omega_c.]; end kz=-pi/a:pi/(10*a):pi/a; plot(kz,freq(:,1:10),.-) xlabel(Kz) ylabel(\omega/c) title(PBG of 1D PC with d=2,a=10,\epsilon_r=12.25) Page * 一维光子晶体TE模式本证值问题 背景为空气,d/a=0.2,周期性介质层的介电常数为12.25。 Page * 一维光子晶体TE模式能带曲线 考核内容及注意事项 注:要求交纸质版答卷和电子版答卷。纸质版答卷包括每题的解题思路和程序,12月23日(星期五)上午3、4节课到逸夫楼406教研室交纸质版答卷。电子版答卷于12月23日前发到邮箱397379519@。 Page * 考核内容和比重: 完成课后作业 40% 上机情况 10% 期末考试试题 50% * * * 数学物理建模与计算机辅助设计 数学物理建模与计算机辅助设计 上机练习专题6:数学物理方程的计算机求解和可视化 Page * 一维无界弦自由振动(无强迫力)齐次定解问题为: 直接利用达朗贝尔公式求解: ( ) 利用行波法求解数学物理定解问题 达朗贝尔解的物理意义 :正向波(或右行波) :反向波(或左行波) 传播速度为 a 实例1:一维无界弦振动定解问题: 根据达朗贝尔公式,位移为 : clear u0=2; u(1:560)=0; x1=-1;x2=1; x=linspace(-2,2,560); u(141:280)=2*u0*(x(141:280)-x1)/(x2-x1); u(281:420)=2*u0*(x2-x(281:420))/(x2-x1); uu=u; h=plot(x,u,linewidth,3); axis([2*x1,2*x2,-4,4]); set(h,EraseMode,xor) for at=2:140 lu(1:560)=0; ru

文档评论(0)

qiwqpu54 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档