| 查看: 610 | 回复: 3 | |||
[交流]
Matlab程序问题 已有2人参与
|
|
各位大神，我写了个程序怎么修改能把图画成如下的图啊？请各位大神给予指导啊。 %.........syms x z k t gamma Delta rho alpha Omega A B C D E gamma=1; Delta_2=0,gamma_1=5*gamma,gamma_2=5*gamma,rho=0,C=600, Omega_1=10*gamma*(1-rho^2)^(1/2);Omega_2=50*gamma*(1-rho^2)^(1/2); x=-100*gamma:0.1*gamma:100*gamma; %.........relational expression.......... z1=zeros(1,length(x)); for k=1:length(x) k Delta_1=0.1*gamma*k A=[0,-conj(Omega_1),0,Omega_1,2*gamma_1,0,0,0,0; -Omega_1,-Delta_1,-conj(Omega_2),0,Omega_1,0,0,0,0; 0,-Omega_2,-Delta_1-Delta_2,0,0,Omega_1,0,0,0; conj(Omega_1)*i,0,0,Delta_1*i-gamma_1,-conj(Omega_1)*i,0,Omega_2*i,2*rho*(gamma_1*gamma_2)^(1/2),0; 0,conj(Omega_1),0,-Omega_1,0,-conj(Omega_2),0,Omega_2,0; 0,0,conj(Omega_1),0,-Omega_2,-Delta_2,0,0,Omega_2; 0,0,0,conj(Omega_2)*i,0,0,Delta_1*i+Delta_2*i-gamma_2,conj(Omega_1)*i,0; 0,0,0,0,conj(Omega_2)*i,0,-Omega_1*i,Delta_2*i-gamma_1-gamma_2,-conj(Omega_2)*i; 1,0,0,0,1,0,0,0,1]; B=[0;0;0;0;0;0;0;0;1]; t=inv(A)*B; D(:,k)=t; %........E(:,k)=zeros(1,length(x)); E(:,k)=Omega_1+C*gamma_1*imag(t)-i*C*gamma_1*real(t); end z1=imag(D(4,  );%......................................... z2=zeros(1,length(x)); for k=1:length(x) k Delta_1=-0.1*gamma*k A=[0,-conj(Omega_1),0,Omega_1,2*gamma_1,0,0,0,0; -Omega_1,-Delta_1,-conj(Omega_2),0,Omega_1,0,0,0,0; 0,-Omega_2,-Delta_1-Delta_2,0,0,Omega_1,0,0,0; conj(Omega_1)*i,0,0,Delta_1*i-gamma_1,-conj(Omega_1)*i,0,Omega_2*i,2*rho*(gamma_1*gamma_2)^(1/2),0; 0,conj(Omega_1),0,-Omega_1,0,-conj(Omega_2),0,Omega_2,0; 0,0,conj(Omega_1),0,-Omega_2,-Delta_2,0,0,Omega_2; 0,0,0,conj(Omega_2)*i,0,0,Delta_1*i+Delta_2*i-gamma_2,conj(Omega_1)*i,0; 0,0,0,0,conj(Omega_2)*i,0,-Omega_1*i,Delta_2*i-gamma_1-gamma_2,-conj(Omega_2)*i; 1,0,0,0,1,0,0,0,1]; B=[0;0;0;0;0;0;0;0;1]; t=inv(A)*B; D(:,k)=t; %........E(:,k)=zeros(1,length(x)); E(:,k)=Omega_1+C*gamma_1*imag(t)-i*C*gamma_1*real(t); end z2=real(D(4, );plot(x,z1,x,z2,'LineWidth',1.5); xlabel('Delta_1/gamma'); ylabel('Im(D(4, )');axis([-100,100,-1,1]); legend('z1','z2'); grid off 其中上面图的Y轴中的字母表示是程序当中D(4，1)的虚部，程序中rho=0表示图中的P=0.  11.png |
回复此楼» 猜你喜欢
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有160人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab 程序问题
已经有15人回复
- Matlab程序问题
已经有2人回复
- matlab程序问题
已经有2人回复
- matlab程序问题求解答
已经有10人回复
- matlab程序问题求助
已经有1人回复
- 【求助】matlab程序问题
已经有0人回复
- TSP matlab程序问题 请各位大侠帮忙啊
已经有0人回复
- matlab程序问题，总是出错，麻烦大家看看问题在那里？
已经有6人回复
- matlab 程序问题
已经有2人回复
- matlab 程序问题
已经有4人回复
2楼2013-07-18 12:54:43
3楼2013-07-18 14:52:32
knightman
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 110.5
- 散金: 1080
- 帖子: 65
- 在线: 57.4小时
- 虫号: 496079
- 注册: 2008-01-16
- 专业: 电力电子学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
把你的程序发给我吧，我给你检查一下。 另，你想画的那个图，也给我。 在此网站上的文本，你自己看看，没法看啊，都转为图片了，很多的笑脸哦（ ）邮箱： knightman@163.com . |
4楼2013-07-19 21:58:05
