版塊導航: 正在加載中...

調(diào)劑小程序

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (6591)
>考研 (2333)
>導師招生 (2066)
>文獻求助 (359)
>蟲友互識 (301)
>碩博家園 (169)
>考博 (155)
>博后之家 (146)
>休閑灌水 (117)
>招聘信息布告欄 (112)
>論文投稿 (106)
>基金申請 (87)
>找工作 (46)
>公派出國 (41)
>綠色求助(高懸賞) (40)
>高分子 (39)

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

phases

金蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 591.3
帖子: 3
在線: 30.9小時
蟲號: 443657
注冊: 2007-10-27
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 用matlab的ode45求解Liouville von-neumann方程的問題

最近為了求解Liouville von-neumann方程用matlab編了一個小程序.對于一個兩能級系統(tǒng),初態(tài)和系統(tǒng)的演化的Hamiltonian已知,用ODE45求解在任意時刻的兩能級系統(tǒng)的狀態(tài).代碼如下

% define the function
function dXdt = Liouvelle(t,X)
% Delta=5;
% epsilon=0.2; % driving strength
% Omega_d=10; % driving frequency
% Omega_r=10; % frequency of cavity

H=2*pi*[0.2 5; 5 -0.2];
%H=2*pi*[Delta -epsilon*cos(2*pi*Omega_d*t); -epsilon*cos(2*pi*Omega_d*t) -Delta]; %System Hamiltonian

X = reshape(X, size(H));

dXdt = -i*(H*X-X*H); %Liouville von-neumann equation
dXdt = dXdt(

;

% call function to solve the problem

clear

% initial state
X2=[0.5 0.5;0.5 0.5];
psi2=sqrt(2)/2*[1 1]';
input=psi2;
%options = odeset('RelTol',1e-8,'AbsTol',1e-10);
[T X]= ode45(@Liouvelle,[0:0.1:5],X2);

[m n] = size(X);
for j=1:m
XX(:,:,j)=reshape(X(j,

,size(X2));

F(j)=input'*XX(:,:,j)*input; % Compute the fidelity, Nielsen 'quantum computation and
% quantum information page 409 eq(9.60)

end

possibility=squeeze(XX(1,1,

);
figure(2); plot(T,F); hold on;

現(xiàn)在問題是發(fā)現(xiàn),對于有些時刻的系統(tǒng)態(tài)密度矩陣結(jié)果不正確.它們的非對角元的乘積竟然大于0.25,而我得到的應該是純態(tài)的密度矩陣,非對角元的乘積應該小于0.25.如果在ODE里面設置計算的精度可以改進,但仍然有些密度矩陣的結(jié)果不正確.不知道問題出在哪里,怎么解決.

謝謝!

回復此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有126人回復
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！已經(jīng)有1人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

有關(guān)MATLAB數(shù)據(jù)的問題已經(jīng)有6人回復
matlab 求指點動力學方程擬合過程中導數(shù)的獲取已經(jīng)有10人回復
matlab的函數(shù)調(diào)用問題已經(jīng)有3人回復
MATLAB中產(chǎn)生高斯白噪聲，涉及到awgn和wgn函數(shù) 已經(jīng)有3人回復
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復
積分方程和微分方程已經(jīng)有4人回復
求matlab求解一個方程組的問題已經(jīng)有6人回復
求高人指點用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復
用MATLAB ode45求解2階微分方程已經(jīng)有4人回復
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復
【求助完畢】用ode逆向求解微分方程問題已經(jīng)有4人回復
【求助】積分微分方程matlab求解已經(jīng)有6人回復
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復
【求助】幫我編個程啊~C語言，C++，Matlab，VB...都行~ 已經(jīng)有14人回復

1樓 2011-04-22 21:44:49

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wuhengbin

★
ben_ladeng(金幣-1): 求助帖，請不要純表 2011-04-23 14:16:35

2樓2011-04-23 13:47:03

已閱回復此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 phases 的主題更新