[討論] 想請問一維光子晶體的能帶模擬

看板MATLAB作者 (人生如戲)時間13年前 (2013/02/02 21:30), 編輯推噓6(6推 0噓 14→)
留言20則, 3人參與, 最新討論串1/1
我想用平面波展開法解一維光子晶體能帶 這是我的code:http://www.wretch.cc/blog/yang200016/13552369 我模擬出的結果覺得很奇怪，能隙gap的大小會和平面波的個數有關， 這是正常的嗎?(甚至第一個gap超小) 小弟還是matlab新手，code可能有些地方會化簡為繁 希望有人一起討論.... -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.116.22.151
推
02/03 15:37, , 1F
這個不是MATLAB問題吧
02/03 15:37, 1F
→
02/03 15:38, , 2F
做一下收斂性吧
02/03 15:38, 2F
推
02/04 09:14, , 3F
取樣數變多了計算的"頻率範圍"也變大了...感覺怪怪的
02/04 09:14, 3F
→
02/04 19:46, , 4F
取樣數變多了，所以fft的input參數變多，相對output
02/04 19:46, 4F
→
02/04 19:46, , 5F
output參數(即頻率)不是也會變多嗎??
02/04 19:46, 5F
→
02/04 20:14, , 6F
我好像已經解決了 好像輸入的數字要是奇數
02/04 20:14, 6F
→
02/04 20:15, , 7F
才能和書本化的一樣 (希望高手能糾錯)
02/04 20:15, 7F
→
02/04 20:16, , 8F
而且我code也改了，針對奇數用的
02/04 20:16, 8F
推
02/04 23:21, , 9F
不一定要奇數吧 要項次夠多 一班好像都是2N+1
02/04 23:21, 9F
→
02/04 23:22, , 10F
可以用寬度來做項數之間的關西~~來當收斂性分析
02/04 23:22, 10F
→
02/04 23:22, , 11F
可看出何時會收斂~
02/04 23:22, 11F
→
02/04 23:22, , 12F
一維的無限域計算應該不會太複雜~~展高次問題不大
02/04 23:22, 12F
推
02/04 23:24, , 13F
機械系 陳聯文也有做類似的東西 聲子光子~可以去問問
02/04 23:24, 13F
→
02/04 23:25, , 14F
台大應力 吳政忠也做蠻多類似的東西
02/04 23:25, 14F
→
02/04 23:25, , 15F
中央光電也有一位教授~~有出書~光子晶體
02/04 23:25, 15F
→
02/04 23:26, , 16F
參考一下吧~~我不是很熟就是了
02/04 23:26, 16F
推
02/04 23:29, , 17F
我看到很嚴重的程式錯誤
02/04 23:29, 17F
→
02/04 23:30, , 18F
你定義i為虛數 這個系統已經默認了 不需要特別寫
02/04 23:30, 18F
→
02/04 23:31, , 19F
for 迴圈又有用到i=1:n 這樣i就變成數字了 非虛數
02/04 23:31, 19F
推
02/04 23:35, , 20F
不過程式好像也沒用到
02/04 23:35, 20F
‣ 更多 janghang 的文章
文章代碼(AID): #1H3HJOXP (MATLAB)