| 查看: 574 | 回复: 0 | |||
[交流]
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,  +F(13,);plot(x,y,'r-','LineWidth',1.5) 上面的循环程序中是一个四阶的矩阵循环，但是感觉程序就是其中的 rho_31=inv(A)*B; D(:,n)=rho_31; rho_41=inv(A)*B; F(:,n)=rho_41; 有问题，请路过的大神们给予指导啊 |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有94人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
各位好，我想问怎么把地震波数据输入MATLAB中显示图形呢
已经有5人回复
- 能否利用matlab画出如下图形，主要是红色框里的细节？
已经有4人回复
- Matlab中3维图形的绘制
已经有17人回复
- matlab图形用户界面的设计-gui命令
已经有7人回复
- 【已完结】matlab仿真图形（交流问题即有奖，呵呵）
已经有6人回复
- matlab画图问题: 怎样用小图标或自定义图形代替+ x . o等做图
已经有10人回复
- 求matlab图形保存方法
已经有11人回复
- matlab图形的横坐标其实坐标点不显示？
已经有3人回复
- matlab里有没有类似于origin合并图形的功能
已经有8人回复
- 求助 matlab正负并存的宽范围数据的图形如何画
已经有7人回复
- 请教Matlab图形保存问题
已经有5人回复
- Matlab中这种图形怎么画？
已经有10人回复
- [求助]matlab画一个非参数方程的三维图形
已经有3人回复
- 如何在matlab中绘制含有三个参数方程的三维图形
已经有6人回复
- 【求助】利用MATLAB求解的外围轮廓包络线为什么不与图形外围轮廓重合？怎样才可以精确
已经有8人回复
- 【求助】vb与matlab混合编程时图形显示窗口中，怎样使图形比窗口小一些？
已经有3人回复
- 【求助】这个MATLAB四维图形如何绘制，请教！！
已经有5人回复
- 【原创】MATLAB模拟李萨茹图形
已经有4人回复
- 【求助】Matlab 四维图形的绘制。
已经有12人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
