| 查看: 1336 | 回复: 1 | ||
phases金虫 (初入文坛)
|
[求助]
用matlab的ode45求解Liouville von-neumann方程的问题
|
|
最近为了求解Liouville von-neumann方程用matlab编了一个小程序.对于一个两能级系统,初态和系统的演化的Hamiltonian已知,用ODE45求解在任意时刻的两能级系统的状态.代码如下 % 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; 现在问题是发现,对于有些时刻的系统态密度矩阵结果不正确.它们的非对角元的乘积竟然大于0.25,而我得到的应该是纯态的密度矩阵,非对角元的乘积应该小于0.25.如果在ODE里面设置计算的精度可以改进,但仍然有些密度矩阵的结果不正确.不知道问题出在哪里,怎么解决. 谢谢! |
回复此楼» 猜你喜欢
第一性原理计算方向2026级博士申请 PRB*1，四级484
已经有1人回复
-
求助VISSIM破解版软件
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有92人回复
-
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有17人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
有关MATLAB数据的问题
已经有6人回复- matlab 求指点 动力学方程拟合过程中导数的获取
已经有10人回复
- matlab的函数调用问题
已经有3人回复
- MATLAB中产生高斯白噪声，涉及到awgn和wgn函数
已经有3人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 积分方程和微分方程
已经有4人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助完毕】用ode逆向求解微分方程问题
已经有4人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】帮我编个程啊~C语言，C++，Matlab，VB...都行~
已经有14人回复
★
ben_ladeng(金币-1): 求助帖，请不要纯表 2011-04-23 14:16:35
ben_ladeng(金币-1): 求助帖，请不要纯表 2011-04-23 14:16:35
 |
2楼2011-04-23 13:47:03