作者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