24小时热门版块排行榜

返回列表

long_3415

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 342.5
帖子: 39
在线: 90小时
虫号: 2087574
注册: 2012-10-26
专业: 光学

[交流] MATLAB模拟电磁诱导吸收图形

各位路过的大神，我在MATLAB软件中需要模拟电磁诱导产生的介质吸收图形，但是结果显示只是竖直方向y轴为0不知道为什么，请路过的各位大神给予指导啊，程序代码如下
k=1;q=1;f=1;gamma=1;
Delta_2=0,Delta=1,Omega_1=1.5,Omega_2=-2,C=500;
gamma_21=3,gamma_31=4,gamma_32=4,gamma_41=5,gamma_42=5;
gamma_2=gamma_21,gamma_3=gamma_31+gamma_32,gamma_4=gamma_41+gamma_42;
tau_41=gamma_4,tau_14=gamma_4,tau_42=gamma_4+gamma_2,tau_24=gamma_4+gamma_2,tau_43=gamma_4+gamma_3,tau_34=gamma_4+gamma_3;
tau_31=gamma_3,tau_13=gamma_3,tau_32=gamma_3+gamma_2,tau_23=gamma_3+gamma_2;
tau_21=gamma_2,tau_12=gamma_2;
x=-100*gamma:0.1*gamma:100*gamma;
y=zeros(1,length(x));
%.........relational expression..........
for n=1:length(x)
n
Delta_1=1.5*x(n)
rho_31=zeros(9,n);
rho_41=zeros(13,n);
A=[0,0,-Omega_1*i,-k*Omega_1*i,0,gamma_21,0,0,conj(Omega_1)*i,0,gamma_31,0,k*conj(Omega_1)*i,0,0,gamma_41;
   0,(Delta_2-Delta_1)*i-(1/2)*tau_12,-Omega_2*i,-q*Omega_2*i,0,0,0,0,0,conj(Omega_1)*i,0,0,0,k*conj(Omega_1)*i,0,0;
   -conj(Omega_1)*i,-conj(Omega_2)*i,-Delta_1*i-(1/2)*tau_13,0,0,0,0,0,0,0,conj(Omega_1)*i,0,0,0,k*conj(Omega_1)*i,0;
   -k*conj(Omega_1)*i,-q*conj(Omega_2)*i,0,(Delta-Delta_1)*i-(1/2)*tau_14,0,0,0,0,0,0,0,conj(Omega_1)*i,0,0,0,k*conj(Omega_1)*i;
   0,0,0,0,(Delta_1-Delta_2)*i-(1/2)*tau_21,0,-Omega_1*i,-k*Omega_1*i,conj(Omega_2)*i,0,0,0,q*conj(Omega_2)*i,0,0,0;
   0,0,0,0,0,-gamma_2,-Omega_2*i,-q*Omega_2*i,0,conj(Omega_2)*i,gamma_32,0,0,q*conj(Omega_2)*i,0,gamma_42;
   0,0,0,0,-conj(Omega_1)*i,-conj(Omega_2)*i,-Delta_2*i-(1/2)*tau_23,0,0,0,conj(Omega_2)*i,0,0,0,q*conj(Omega_2)*i,0;
   0,0,0,0,-k*conj(Omega_1)*i,-q*conj(Omega_2)*i,0,(Delta-Delta_2)*i-(1/2)*tau_24,0,0,0,conj(Omega_2)*i,0,0,0,q*conj(Omega_2)*i;
   Omega_1*i,0,0,0,Omega_2*i,0,0,0,Delta_1*i-(1/2)*tau_31,0,-Omega_1*i,-k*Omega_1*i,0,0,0,0;
   0,Omega_1*i,0,0,0,Omega_2*i,0,0,0,Delta_2*i-(1/2)*tau_32,-Omega_2*i,-q*Omega_2*i,0,0,0,0;
   0,0,Omega_1*i,0,0,0,Omega_2*i,0,-conj(Omega_1)*i,-conj(Omega_2)*i,-gamma_3,0,0,0,0,0;
   0,0,0,Omega_1*i,0,0,0,Omega_2*i,-k*conj(Omega_1)*i,-q*conj(Omega_2)*i,0,Delta*i-(1/2)*tau_34,0,0,0,0;
   k*Omega_1*i,0,0,0,q*Omega_2*i,0,0,0,0,0,0,0,(Delta_1-Delta)*i-(1/2)*tau_41,0,-Omega_1*i,-k*Omega_1*i;
   0,k*Omega_1*i,0,0,0,q*Omega_2*i,0,0,0,0,0,0,0,(Delta_2-Delta)*i-(1/2)*tau_42,-Omega_2*i,-q*Omega_2*i;
   0,0,k*Omega_1*i,0,0,0,q*Omega_2*i,0,0,0,0,0,-conj(Omega_1)*i,-conj(Omega_2)*i,-Delta*i-(1/2)*tau_43,0;
   1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1];
B=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1];
rho_31=inv(A)*B;
D(:,n)=rho_31;
rho_41=inv(A)*B;
F(:,n)=rho_41;
%...........E(:,n)=zeros(1,length(x));
E(:,n)=Delta_1-i*C*(rho_31+f*rho_41);
end
y=imag(D(9,