5耦合率方程如何建模?

2视图(30天)

显示旧的评论

威廉•罗伊尔 2017年9月4日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/355140-how-to-model-5-coupled-rate-equations

评论道: 卡兰吉尔 2017年9月5日

你好,我是试图利用dsolve解决5耦合速率方程组。检查我的阴谋粒子的总数来看,以确保它保持不变。如果phonon_dist设置为0这是发生了什么代码的工作原理。但这删除一个水平的系统,因为它只是填充和de-populated条款,包括phonon_dist。

当phonon_dist 0.5所有我看到的是N_ex腐烂,我看到没有其他瞬态行为。我不知道这是为什么。

我有多次推导速度方程和双检查矩阵是正确的。我可以看到所有的平衡方面,负面的平衡是积极的方面。所以我不确定粒子的人口变化的原因,我不能失去他们的我。

解决方案万博尤文图斯由Matlab计算非常长。是Matlab解决这个大问题都处理不了呢?万博尤文图斯还是我盲目,不能看到一个明显的错误代码吗?

提前谢谢。

                         如果真正的
                        
                         % %定义生命周期变量
                        
                         T_th = 10 e-11;
                        
                         T_cap = 10 e-9;
                        
                         T_ph = 10 e-20;
                        
                         T_r2p = 10 e-8;
                        
                         T_nr = 10 e;
                        
                         T_rcb = 10 e-8;

%定义的兴奋声子温度T phonon_dist = 5 e 1;

%设置问题,表示与符号函数u和v sym N_ex (t) N_cb (t) N_2p (t) N_1s (t) N_gr (t)

%建立微分方程。与N_ex N_cb、N_2p N_1s N_gr

一个= [(1 / T_th) 0 0 0 0;(1 / T_th) (1 / T_cap 2 * phonon_dist / T_ph 1 / T_rcb) (phonon_dist * 1 / T_ph) (phonon_dist * 1 / T_ph) 0;0 (1 / T_cap + phonon_dist * 1 / T_ph) (2 * phonon_dist / T_ph 1 / T_r2p 1 / T_nr) (phonon_dist * 1 / T_ph) 0;0 (phonon_dist * 1 / T_ph) (phonon_dist * 1 / T_ph) (2 * phonon_dist / T_ph 1 / T_nr) 0;0 (1 / T_rcb) (1 / T_nr + 1 / T_r2p) (1 / T_nr) 0];

% ()

B = [N_ex;N_cb;N_2p;N_1s;N_gr];

%设置微分eqn ordinary_diff_eqns = diff (B) = = A * B;

%创建初始条件Initial_conds = B (0) = = (1000;0;0;0;0);

%解决微分方程使用上面定义的开始条件[Soln_Struct] = dsolve (ordinary_diff_eqns Initial_conds);

%提取个人水平的解决方案。万博尤文图斯dsolve解决命令然后%将按字母顺序排列的解决方案需要提取的解决方案万博尤文图斯。N_ex_Soln = Soln_Struct (1) .N_ex;N_cb_Soln = Soln_Struct (1) .N_cb;N_2p_Soln = Soln_Struct (1) .N_2p;N_1s_Soln = Soln_Struct (1) .N_1s;N_gr_Soln = Soln_Struct (1) .N_gr;

Total_electrons = N_ex_Soln + N_cb_Soln + N_2p_Soln + N_1s_Soln + N_gr_Soln;

%画出解决方案看时间依万博尤文图斯赖性。% clf删除从当前所有图形对象clf抓住网格图标题(V6——时间表= 5平台以及)fplot (N_ex_Soln) xlim([0、5平台以及])ylim([0, 1050])包含(时间)ylabel(人口)fplot (N_cb_Soln) fplot (N_2p_Soln) fplot (N_1s_Soln) fplot (N_gr_Soln) fplot (Total_electrons)传说(“N_{前任}’,‘N_ {cb}’,‘N_ {2 p}’,‘N_ {1}’,‘N_ {gr}’,‘总电子’,“位置”,“最好”)