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