MATLAB在量子力學(xué)中的應(yīng)用

MATLAB MATLAB 語言語言 課程論文課程論文 MATLABMATLAB 在量子力學(xué)中的應(yīng)用在量子力學(xué)中的應(yīng)用 姓名 魏祎姓名 魏祎 學(xué)號 學(xué)號 1201124398912011243989 專業(yè) 通信工程專業(yè) 通信工程 班級 班級 20112011 級通信工程級通信工程 指導(dǎo)老師 湯全武指導(dǎo)老師 湯全武 學(xué)院 物理電氣信息學(xué)院學(xué)院 物理電氣信息學(xué)院 完成日期 完成日期 2012 12 112012 12 11 1 MATLABMATLAB 在量子力學(xué)中的應(yīng)用在量子力學(xué)中的應(yīng)用 魏祎 12011243989 2011 級通信班 摘要摘要 量子力學(xué)的應(yīng)用和成就是多方面的 迄今仍保持有旺盛的生命力 碩果頗傳 雖然 大學(xué)物理 中介紹的量子力學(xué)只是一些最基本的概念 但之中涉及了許多復(fù)雜的數(shù)值計 算問題 解微分方程的問題 圖像顯示問題 例如一維無限深勢阱問題 一維運動粒子的 波函數(shù)曲線問題 對其手工求解較為復(fù)雜 而 MATLAB 語言正是處理這些復(fù)雜問題的很好工 具 既能進行數(shù)值求解 又能繪制有關(guān)曲線 非常方便實用 另外利用其可以減少工作量 節(jié)約時間 加深理解對量子力學(xué)的理解 同時可以培養(yǎng)應(yīng)用知識的能力 關(guān)鍵詞關(guān)鍵詞 量子力學(xué) MATLAB 語言 一維無限深方勢阱 波函數(shù) 概率密度 一 問題的提出一 問題的提出 MATLAB 語言是當(dāng)今國際上科學(xué)界 尤其是自動控制理論 最具影響力 也 是最有活力的軟件 它提供了強大的科學(xué)運算 靈活的程序設(shè)計流程 高質(zhì)量 的圖形可視化與界面設(shè)計 便捷的與其他程序和語言接口的功能 MATLAB 語言 在各國高校與研究單位起著重大的作用 它是一種集數(shù)值運算 符號運算 可 視化建模 仿真和圖形處理等多種功能 在量子力學(xué)中 可以利用其幫助初學(xué) 者理解量子力學(xué)與經(jīng)典力學(xué)截然不同的思維方式和觀念 理解微觀粒子的波粒 二象性 下面將以一維勢阱問題 波函數(shù)和概率密度曲線問題為例講述 MATLAB 在量子力學(xué)中的應(yīng)用 二 用二 用 MATLABMATLAB 語言求解一維無限深勢阱問題語言求解一維無限深勢阱問題 如圖 1 所示 設(shè)想一粒子處在勢能為的力場中 并沿 x 軸作一維運動 p E 粒子的勢能滿足下述邊界條件 p E 1 當(dāng)粒子在范圍內(nèi)時 0 ax 0 p E 2 當(dāng)及時 0 xax p E 這就是說粒子只能在寬度為 a 的兩個無限高勢壘壁之間自由運動 就像一 小球被限制在無限深的平底深谷中運動那樣 我們理想化了得勢阱曲線叫無限 深方形勢阱 因為粒子只限于沿 x 軸方向運動 故這個勢阱為一維無限深的方 形勢阱 簡稱一維方勢阱 有上述邊界條件已知 粒子在勢阱中得勢能 x 與時間無關(guān) 且 0 因 p E p E 此 由一般的薛定諤方程 1 粒子在無限深方勢阱中得定態(tài)薛定諤方程為 1 0 8 2 2 2 2 h mE dx d 式中 m 為粒子的質(zhì)量 E 為粒子的總能量 如令 k 為 2 2 2 8 h mE k 2 則上式可寫成 3 0 2 2 2 k dx d 根據(jù)邊界條件 x 0 時 則可以利用 MATLAB 求解微分方程 0 0 利用 MATLAB 語言求解此方程程序如下 y dsolve D2y k 2 y y 0 0 x 求方程 16 p E 0 kxsin 0 0 a x 圖圖 1 1 一維無限深方勢阱中得粒子一維無限深方勢阱中得粒子 運行結(jié)果 y C1 sin k x 又根據(jù)邊界條件 x a 時 此時式 16 的解為kaCasin 1 4 0sin 1 kxCa 一般說來 A 可不為零 故 有0sin ka nka n 1 2 3 上式也可寫成 a n k 將上式與式 4 相比較 可得勢中粒子可能的能量值為 5 2 2 2 8ma h nE 式中 n 為量子數(shù) 表明粒子的能量只能取離散的值 當(dāng) n 1 時 勢阱中粒 3 子的能量為 n 2 3 4 時 4 9 16 這就是說 一維無 2 2 1 8ma h E 1 E 1 E 1 E 限深方勢阱中粒子的能量是量子化的 下面在來確定常數(shù) 由于粒子被限制在和 的勢阱中 因此 1 C 0 xax 按歸一化條件 粒子在此區(qū)間內(nèi)出現(xiàn)的概率總和為 1 即 6 aa dx 0 2 0 1 或 7 1sin 0 2 2 1 dxx a n C a 令 則上式左側(cè)積分為dx a d a x 8 aA aA d a na C a 2 2 0 2 2 1 2 1 2 sin 于是 可得 a C 2 1 這樣 式 7 所表現(xiàn)得波函數(shù)即為 9 x a n a x sin 2 ax 0 由此可得 能量為 E 所表示的粒子在勢阱中得概率密度為 10 sin 2 22 x a n a x 下面用 MATLAB 語言求解電子的各能級能量 波函數(shù)曲線和概率密度曲線 程序如下 function E shor m a 建立函數(shù)文件 n 1 10 量子數(shù) n h 6 63 1e 34 普朗克常量 E n 2 h 2 8 m a 2 n 能級的能量值 x 0 1 0 1e 12 a x 的值 subplot 4 2 1 分割繪圖區(qū)域 第一個子圖 y1 sqrt 2 a sin pi x a n 1 的波函數(shù) plot x y1 繪制 n 1 的波函數(shù)圖象 title n 1 給 n 1 的破函數(shù)曲線加標(biāo)題 subplot 4 2 3 分割繪圖區(qū)域 第三個子圖 y1 sqrt 2 a sin 2 pi x a n 2 的波函數(shù) plot x y1 繪制 n 2 的波函數(shù)圖象 4 title n 2 給 n 2 的破函數(shù)曲線加標(biāo)題 subplot 4 2 5 分割繪圖區(qū)域 第五個子圖 y1 sqrt 2 a sin 3 pi x a n 3 的波函數(shù) plot x y1 繪制 n 3 的波函數(shù)圖象 title n 3 給 n 3 的破函數(shù)曲線加標(biāo)題 subplot 4 2 7 分割繪圖區(qū)域 第七個子圖 y1 sqrt 2 a sin 4 pi x a n 4 的波函數(shù) plot x y1 繪制 n 4 的波函數(shù)圖象 title n 4 給 n 4 的波函數(shù)曲線加標(biāo)題 subplot 4 2 2 分割繪圖區(qū)域 第二個子圖 y2 2 a sin pi x a 2 n 1 的概率密度曲線 plot x y2 繪制 n 1 的概率密度曲線 title n 1 給 n 1 的概率函數(shù)曲線加標(biāo)題 subplot 4 2 4 分割繪圖區(qū)域 第四個子圖 y2 2 a sin 2 pi x a 2 n 2 的概率密度函數(shù) plot x y2 繪制 n 2 的概率密度曲線 title n 2 給 n 2 的概率函數(shù)曲線加標(biāo)題 subplot 4 2 6 分割繪圖區(qū)域 第六個子圖 y2 2 a sin 3 pi x a 2 n 3 的概率密度函數(shù) plot x y2 繪制 n 3 的概率密度曲線 title n 3 給 n 3 的概率函數(shù)曲線加標(biāo)題 subplot 4 2 8 分割繪圖區(qū)域 第八個子圖 y2 2 a sin 3 pi x a 2 n 4 的概率密度函數(shù) plot x y2 繪制 n 4 的概率密度曲線 title n 4 給 n 4 的概率函數(shù)曲線加標(biāo)題 假如有一電子在寬度為 0 02nm 的一維方勢阱中 則其 調(diào)kgm 31 101 9 用 shor 函數(shù)求解其各能級能量 波函數(shù)曲線和概率密度曲線 程序如下 a 0 2 1e 9 a 為勢阱的寬度 m 9 1 1e 31 粒子的質(zhì)量 E shor m a 調(diào)用 shor 函數(shù) 運行結(jié)果如下 E 1 0e 015 0 0015 0 0060 0 0136 0 0242 0 0377 0 0543 0 0740 0 0966 0 1223 0 1510 運行結(jié)果的圖像如圖 2 所示 5 圖圖 2 2 在一維無限深方勢阱中 粒子的能級 波函數(shù)和概率密度在一維無限深方勢阱中 粒子的能級 波函數(shù)和概率密度 通過以上程序?qū)αＷ釉谝痪S無限深方勢阱中運動的圖形描述 數(shù)據(jù)分析 我們發(fā)現(xiàn) 一維無限深方勢阱中粒子的能量是量子化的 為 1 0e 015 0 0015 1 0e 015 0 0060 1 0e 015 0 0136 1 0e 015 0 0242 一些列量子化的能量值 粒子在勢阱各處的概率密度并不是均勻分布的 隨量 子數(shù)而改變 當(dāng)量子數(shù) n 1 時 粒子在勢阱中部 x a 2 附近出現(xiàn)得概率最大 而在兩端出現(xiàn)的概率為零 隨著量子數(shù) n 的增大 概率密度分布曲線的峰值個 數(shù)也增多 例如 n 2 有兩個峰值 n 3 有三個峰值 而且兩相鄰峰值間的距 離將縮小得很小 彼此靠的很近 非常接近粒子在勢阱中各處概率處處相等的 情況 下面我們再來求解一個一維運動粒子的波函數(shù)和概率密度曲線問題 三 用三 用 MATLABMATLAB 語言求解一維運動粒子的波函數(shù)曲線 概率密度曲線語言求解一維運動粒子的波函數(shù)曲線 概率密度曲線 已知一維運動粒子的波函數(shù)為 11 00 0 x xAxe x x 有歸一化條件 aa dx 0 2 0 1 或 6 12 1 0 222 dxexA x 上式積分可得 13 3 2 0 222 4 A dxexA x 于是可得 14 3 2 A 這樣 式 12 表示的波函數(shù)即為 15 00 02 3 x xxe x x 由此可得 一維運動粒子的概念密度函數(shù)為 16 00 04 223 x xex x x 下面用 MATLAB 語言求解一維運動粒子的波函數(shù)曲線和概率密度曲線 程序如下 r input please input r 輸入 r 的值 a input please input a 請輸入 a 的值 x 0 0 0001 a 設(shè)置位置變量 x 的值 y1 4 r 3 x 2 exp 2 r x x 0 的概率密度函數(shù) y2 2 sqrt r 3 x exp r x x 0 的波函數(shù) y3 max y1 求概率密度最大值 x1 x find y1 y3 求概率密度最大值對應(yīng)的 x 值 subplot 1 2 1 分割繪圖區(qū)域 第一個子圖 plot x y1 x1 y3 p 繪制概率密度曲線 并標(biāo)記最大概率點 title 概率密度曲線 subplot 1 2 2 分割繪圖區(qū)域 第二個子圖 plot x y2 繪制波函數(shù)曲線 title 波函數(shù)曲線 運行結(jié)果如下 當(dāng) 1 時的運行結(jié)果如下 圖像如圖 3 所示 please input r2 please input a6 y3 1 0827 x1 0 5000 當(dāng) 0 5 時的運行結(jié)果如下 運行圖像如圖 4 所示 please input r0 5 please input a20 7 y3 0 2707 x1 2 由以上程序的數(shù)據(jù)分析和圖形顯示可知 當(dāng) 1 時 粒子在 x 1 處出現(xiàn)的 概率最大 概率為 0 5413 當(dāng)當(dāng) 0 5 時 粒子在 x 2 處出現(xiàn)的概率最大 概率為 0 2707 給定一個 就能求出一個概率最大值和概率最大 10 的 x 值 再次說明一維運動的粒子在空間各處出現(xiàn)的概率密度是不均勻的 隨 的改變而改變 越大 最大概率密度峰值也遠大 這與經(jīng)典力學(xué)很不相同 按照經(jīng)典力學(xué)粒子在空間各處運動是不受限制的 粒子在空間各處出現(xiàn)的概率 亦應(yīng)是相等的 但在量子力學(xué)中 這樣的結(jié)論明顯不成立 空間粒子在各處出 現(xiàn)的概率明顯不均勻 從圖 3 和圖 4 我們可以看出 空間粒子只有可能出現(xiàn)在 某一小區(qū)域內(nèi) 在其它區(qū)域出現(xiàn)的概率為零 圖圖 3 1 時一維運動粒子的波函數(shù)曲線和概率密度曲線時一維運動粒子的波函數(shù)曲線和概率密度曲線 8 圖圖 4 0 5 時一維運動粒子的波函數(shù)曲線和概率密度曲線時一維運動粒子的波函數(shù)曲線和概率密度曲線 四 結(jié)論四 結(jié)論 從以上利用 MATLAB 語言對兩個以為運動粒子波函數(shù)和概率密度曲線的分析 我們不難得出以下結(jié)論 1 MATLAB 語言強大的繪圖功能可以很有效地幫助我們理解量子力學(xué) 理解一 維運動粒子 的波粒二象性 在研究一維運動粒子問題時 使用的繪圖函數(shù)為 plot x y 這 種高層繪圖操作簡單便捷 方便實用 2 MATLAB 語言具有豐富的符號運算 對于求解一些復(fù)雜的微分方程帶來了很 大的便捷 本文采用了符號函數(shù) dslove 求解一維薛定諤方程的解析解 3 MATLAB 語言具有強大的矩陣運算功能 在作一些較繁雜的數(shù)學(xué)運算時 能 給我們帶來極大的方便 節(jié)約時間和精力 本文用之求粒子的能級快捷準(zhǔn)確 而且求出的能級數(shù)巨多 達到了手工無法達到的計算效果 五 課程體會五 課程體會 經(jīng)過一學(xué)期緊張而有序的課程學(xué)習(xí) 我收獲良多 我深刻體會到 MATLAB 對 于我們通信專業(yè)的重要性 MATLAB 語言相對與同類程序語言具有功能強 效率 高 簡單易學(xué)等特點 在許多領(lǐng)域得到廣泛應(yīng)用 目前最流行的幾種科學(xué)計算 機軟件各具特點 而且都在不斷地發(fā)展 版本不斷涌現(xiàn) 但其中影響最大 流 行最廣的當(dāng)屬 MATLAB 語言 因此學(xué)好 MATLAB 語言實屬必要之舉 但學(xué)校開的 MATLAB 的課程課時較少 我們課后花費了大量時間學(xué)習(xí) 經(jīng)過短時間的學(xué)習(xí) 我知道了 MATLAB 有很好的編輯環(huán)境 其有 MATLAB 桌 面和命令窗口 歷史命令窗口 編輯器和調(diào)試器 路徑搜索和用于用戶瀏覽幫 9 助 工作空間 文件的瀏覽器等工具可以方便的使用 MATLAB 的函數(shù)和文件 它 有簡單的易用的程序語言 它包含控制語句 函數(shù) 數(shù)據(jù)結(jié)構(gòu) 輸入和輸出和 面向?qū)ο缶幊烫攸c 它有強大的科學(xué)計算機數(shù)據(jù)處理能力 它是是一個包含大 量計算算法的集合 它有出色的圖形處理功能 MATLAB 自產(chǎn)生之日起就具有方 便的數(shù)據(jù)可視化功能 以將向量和距陣用圖形表現(xiàn)出來 并且可以對圖形進行 標(biāo)注和打印 高層次的作圖包括二維和三維的可視化 圖象處理 動畫和表達 式作圖可用于科學(xué)計算和工程繪圖等特點 下面我具體敘述我的學(xué)習(xí)體會 在第二章 我們學(xué)習(xí)了 MATLAB 數(shù)據(jù)及其運 算 MATLAB 數(shù)據(jù)類型比其他語言更為豐富 除數(shù)值型 字符型等基本數(shù)據(jù)類型 外 還有結(jié)構(gòu)體 單元體等更為復(fù)雜的對象 豐富的數(shù)據(jù)類型 增強了 MATLAB 的數(shù)據(jù)表達能力 給應(yīng)用帶來了很大的方便 并且 MATLAB 在運用數(shù)據(jù)時能自動 判別數(shù)據(jù)類型 不用先定義后使用 又減少了編程的難度和工作量 因此學(xué)習(xí) 起來倍感輕松 在第三章 我們學(xué)習(xí)了 MATLAB 矩陣分析與處理 MATLAB 的矩陣運算功能 非常豐富 應(yīng)用也非常廣泛 許多含有矩陣運算的復(fù)雜計算問題 在 MATLAB 中 很容易得到解決 在學(xué)習(xí)這一章時 深刻體會到 MATLAB 的計算功能只巨強 在 學(xué)習(xí)線性代數(shù)時 許許多多的計算特別繁瑣 并且手工計算容易出錯 有了 MATLAB 之后 矩陣的運算就簡單多了 在第四章我們學(xué)習(xí)了 MATLAB 一些簡單的程序設(shè)計 知道 MATLAB 有循環(huán)結(jié) 構(gòu) 選擇結(jié)構(gòu)和順序結(jié)構(gòu) 其中選擇結(jié)構(gòu)有 if else 結(jié)構(gòu) if 結(jié)構(gòu) 多 分支 if 語句 switch 語句 try 語句 循環(huán)結(jié)構(gòu)中有 for 語句和 while 語句 這些語言的語法規(guī)則雖然和其它語言有所區(qū)別 但其邏輯思維最基本相同 我 們還學(xué)習(xí)了 MATLAB 的向量運算 這種運算可以代替循