该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
clear;
t=4;t_L=1*t;e=1;h=1;v_F=2*t;
N=50;
B=0.3*t;
V4=2*t;
epsilon_D=-4:0.4:4;
tau=0:0.1:5;
C4=zeros(length(epsilon_D),length(tau));
T4=zeros(length(epsilon_D),length(epsilon_D));
D4=zeros(1,length(epsilon_D));
Q4=zeros(50,50);
t_m5=zeros(length(epsilon_D),50);
z4=zeros(length(epsilon_D),length(tau));
f4=zeros(length(epsilon_D),length(tau));
g4=zeros(length(epsilon_D),length(tau));
for i=1:length(epsilon_D);
for j=1:length(tau);
for m=1:N;
t_m4=(2*i*(t_L)^2*sin(m*(e*V4)/(N*h*v_F)))/(t*exp(-2*i*m*(e*V4)/(N*h*v_F))*(-2*t*cos(m*(e*V4)/(N*h*v_F))-epsilon_D(i)+B)+2*(t_L)^2*exp(-i*m*(e*V4)/(N*h*v_F)));
t_m5(1,m)=t_m5(1,m)+(2*i*(t_L)^2*sin(m*(e*V4)/(N*h*v_F)))/(t*exp(-2*i*m*(e*V4)/(N*h*v_F))*(-2*t*cos(m*(e*V4)/(N*h*v_F))-epsilon_D(i)+B)+2*(t_L)^2*exp(-i*m*(e*V4)/(N*h*v_F)));
for n=1:N
t_n4=(2*i*(t_L)^2*sin(n*(e*V4)/(N*h*v_F)))/(t*exp(-2*i*n*(e*V4)/(N*h*v_F))*(-2*t*cos(n*(e*V4)/(N*h*v_F))-epsilon_D(i)+B)+2*(t_L)^2*exp(-i*n*(e*V4)/(N*h*v_F)));
if n==m
Q4(m,n)=1-tau(j)*((e*V4)/(2*N*h))*(t_m4*conj(t_m4)/pi);
else if n~=m
Q4(m,n)=-(conj(t_m4)*t_n4/pi)*(exp(-i.*tau*(n-m)*(e*V4)/(2*N*h))*sin(tau(j)*(n-m)*(e*V4)/(2*N*h)))./(n-m);
end
end
end
end
T4(1,i)=T4(1,i)+(t_m5)*((t_m5)');
D4(1,i)= D4(1,i)+1/(T4(1,i)*(e*V4)/(2*pi*N*h));
C4(i,j)=C4(i,j)+real(det(Q4));
end
end
tau=0:0.1:5;
for i=1:length(epsilon_D);
z4(i,:)= z4(i,:)+gradient(C4(i,:))./gradient(tau);
f4(i,:)=f4(i,:)+gradient(z4(i,:))./gradient(tau);
g4(i,:)=D4(1,i)*f4(i,:);
end
mesh(epsilon_D,tau,g);
最后
以上就是炙热大侠最近收集整理的关于matlab怎么使用两个循环,急急急!matlab含有两个for循环,程序不知道哪儿错了,求指教!...的全部内容,更多相关matlab怎么使用两个循环,急急急!matlab含有两个for循环,程序不知道哪儿错了,求指教!内容请搜索靠谱客的其他文章。
发表评论 取消回复