matlab怎么使用两个循环,急急急！matlab含有两个for循环，程序不知道哪儿错了，求指教！...

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

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);