我的第一個for循環(huán)不循環(huán)

zeng119

第一個for循環(huán)不循環(huán)，我沒法把圖繪制在一起，諸君有什么辦法嘛


InE = 0.01*10^(-3);  %mA
DnB = 2;  %cm^2/s
WB = 50*10^(-7);  %cm
q = 1.6e-19;  %元電荷C
NB0 = 1.5e16;  % /cm^3                    參量定義
Wb=linspace(0,50*10^(-7),50);            %為了提高競速，將基區(qū)寬度劃分50等分
for n = 0: 2: 8                              %主循環(huán)，n分別循環(huán)0、2、4、6、8

end
for value = 1:50                         %循環(huán)求積分
    syms x;
    f1 = @ (x) NB0.*(1-x./WB);
    c1(value) = integral(f1,Wb(value),WB);
end

for value = 1:50                          %循環(huán)求積分
    syms x;
    f2 = @ (x) NB0.*exp(-n.*x./WB);
    c2(value) = integral(f2,Wb(value),WB);
end

z1=vpa(c1,5);                             %5位精度
z2=vpa(c2,5);                             %5位精度
x=linspace(0,50*10^(-7),50);              %X分割，函數定義
majority1=NB0.*(1-x./WB);                 %計算線性區(qū)雜質濃度
majority2=NB0.*exp(-n.*x./WB);            %計算指數區(qū)雜質濃度
nB1=InE.*z1./(q*DnB*majority1);           %計算線性區(qū)少子濃度
nB2=InE.*z2./(q*DnB*majority2);           %計算指數區(qū)少子濃度

%繪制圖形，4張子圖，上邊為第一種分布，下邊為第二種分布

%***********線性基區(qū)少子濃度分布圖********************
subplot(1,2,1);
plot(x,nB1);
title('線性基區(qū)少子濃度分布圖');
xlabel('基區(qū)寬度/cm');
ylabel('少子濃度/cm^3');
hold on;

%*************緩變基區(qū)少子濃度分布圖****************************
subplot(1,2,2);
plot(x,nB2);
title('緩變基區(qū)少子濃度分布圖');
text(WB./10,nB2(5),['n=' num2str(n) ]);
xlabel('基區(qū)寬度/cm');
ylabel('少子濃度/cm^3');
hold on;