| 查看: 1582 | 回复: 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里面设置计算的精度可以改进,但仍然有些密度矩阵的结果不正确.不知道问题出在哪里,怎么解决. 谢谢! |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有65人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关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
