作者janghang (人生如戏)
看板MATLAB
标题[讨论] 想请问一维光子晶体的能带模拟
时间Sat Feb 2 21:29:57 2013
我想用平面波展开法解一维光子晶体能带
这是我的code:
http://www.wretch.cc/blog/yang200016/13552369
我模拟出的结果觉得很奇怪，能隙gap的大小会和平面波的个数有关，
这是正常的吗?(甚至第一个gap超小)
小弟还是matlab新手，code可能有些地方会化简为繁
希望有人一起讨论....
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.116.22.151
1F：推 alamabarry:这个不是MATLAB问题吧 02/03 15:37
2F：→ alamabarry:做一下收敛性吧 02/03 15:38
3F：推 Jeffch:取样数变多了计算的"频率范围"也变大了...感觉怪怪的 02/04 09:14
4F：→ janghang:取样数变多了，所以fft的input参数变多，相对output 02/04 19:46
5F：→ janghang:output参数(即频率)不是也会变多吗?? 02/04 19:46
6F：→ janghang:我好像已经解决了 好像输入的数字要是奇数 02/04 20:14
7F：→ janghang:才能和书本化的一样 (希望高手能纠错) 02/04 20:15
8F：→ janghang:而且我code也改了，针对奇数用的 02/04 20:16
9F：推 alamabarry:不一定要奇数吧 要项次够多 一班好像都是2N+1 02/04 23:21
10F：→ alamabarry:可以用宽度来做项数之间的关西~~来当收敛性分析 02/04 23:22
11F：→ alamabarry:可看出何时会收敛~ 02/04 23:22
12F：→ alamabarry:一维的无限域计算应该不会太复杂~~展高次问题不大 02/04 23:22
13F：推 alamabarry:机械系 陈联文也有做类似的东西 声子光子~可以去问问 02/04 23:24
14F：→ alamabarry:台大应力 吴政忠也做蛮多类似的东西 02/04 23:25
15F：→ alamabarry:中央光电也有一位教授~~有出书~光子晶体 02/04 23:25
16F：→ alamabarry:参考一下吧~~我不是很熟就是了 02/04 23:26
17F：推 alamabarry:我看到很严重的程式错误 02/04 23:29
18F：→ alamabarry:你定义i为虚数 这个系统已经默认了 不需要特别写 02/04 23:30
19F：→ alamabarry:for 回圈又有用到i=1:n 这样i就变成数字了 非虚数 02/04 23:31
20F：推 alamabarry:不过程式好像也没用到 02/04 23:35