|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear7 n5 I, I, P) h# O1 f7 E
clc
" _- `- j- S) z2 K2 I8 ^& umag=0:0.1:30;
4 I% Q" _6 D' n8 [" sfor i=1:length(mag)
2 e* @, {7 P, P! ]7 s3 Mtau_v0=4.0d0;
$ A' Y9 Y! ?) V4 `! k# P& Dgamma=40e-6;! t7 t. n7 Q7 w( g: ?6 M; S
aRF=1e5; & f; E: V0 \% t) n/ [# `: e
T=40;
8 {% q0 O+ {0 g% P8 {, }- I3 yTnd=150;& |/ p% _1 Q, U! b' [
Tx=14;* Z+ l8 U1 b0 M/ b4 ~3 @. ]# n
Txx=33;& o1 ]3 o0 G! V; L0 e* B7 h
g=4.57e-4;
. |& h$ {6 Y- _: A( c: hbeta=100;
3 w0 N9 s' n# N+ ^4 {" ITbd=1;8 D+ l. j' d) M. A0 I7 M
kb=8.617e-5;
8 K# M; w/ x: x- ~6 `. q5 z; _+ Z& `4 eg1=4.57e-4;
. y% B" e1 h9 J A* a/ Zmub=0.058e-3;, u' h/ A* r! F% Y
delta_s=1.0*mub*mag;/ o" V" Q- u$ ~0 o
delta_v=2.0*mub*mag;/ G: b+ z3 ] H3 Y$ F
delta_a=2.0*mub*mag;
- e* h6 q, y' ]8 ~7 g) s6 W, nE_vs1=2*delta_a;%表示两个亮激子之间的能量差, U# a+ b. M) i9 w+ y, B% _& K
E_vs=4*delta_s;%表示两个双激子之间的能量差
- v: x. t) @8 _
* ~' r( e4 z: \+ S7 n# G, gs1=(gamma*gamma)./(gamma^2+E_vs1.^2).*(1.0/tau_v0);%激子3 H9 t9 i: w1 V- N! [
s=(gamma*gamma)./(gamma^2+E_vs.^2).*(1.0/tau_v0);%双激子
- Z3 R' e0 y1 T: qif mag=0! g- O9 }8 {/ ?& J. ^
* A! o6 l. |( s, |
我这里想做一个判断,如果mag=0我就可以这样,如果mag不等于0可以那样,那么这段代码可以怎么改呢 | . I. H4 [: O1 G; }$ W/ y; u: p
|
|