Index exceeds the number of array elements (1).

close all

clear all

clc

%下列為參數

R = 0.308;

v = 2.29357 * (10^14);%輸出光的頻率

h = 6.626 * (10^-34);%單位是J*s 普郎克常數

q = 1.6 * (10 ^ -19);%單位為coulomb 電子的帶電量 老師講義42頁參數

d = 0.2 * (10 ^ -3);%單位為cm active region厚度 老師講義42頁參數

w = 2 * (10 ^ -3);%單位為cm active region寬度 老師講義42頁參數

L = 250 * (10 ^ -3);%單位為cm active region長度 老師講義42頁參數

alpha_int = 40; %單位為1/cm Internal absorption coefficient 老師講義42頁參數

V = d * w * L ;%單位為cm**3 active region的體積 老師講義42頁參數相乘算出來的

v_g = 8.57 *(10 ^ 9);%單位為cm/sec 群速度 老師講義42頁參數

tao_p = 1.6 *( 10 ^ -12);%光子壽命時間 老師講義42頁參數

tao_c = 2 * (10 ^ -9);%載子壽命時間 老師講義42頁參數

gama = 0.3; %光限制的因素 老師講義42頁參數

sigma_g = 2.5 * (10 ^ -16);%單位為cm**2 不同gain的係數 老師講義42頁參數

N_tr = 10 ^ 18 ;%單位為 1/cm**3 轉換載子密度 老師講義42頁參數

aslo_NL = 2 *(10^-17);%單位為 cm**3 gain的抑制係數 老師講義42頁參數

G = gama * v_g * sigma_g; %老師講義44頁參數

A = 10 ^ 8;%單位為1/s nonradiative recombination rate (老師講義28頁)

B = 10 ^ -10;%單位為(cm**3)/s radiative recombination rate (老師講義28頁)

C = 3 * (10 ^-29);%單位為(cm**6)/s Auger recombination rate (老師講義28頁)

beta_sp = [10^-3,10^-4,10^-5]; %自發輻射的因素

F1 = 0.9; %fractional power out of front facet MIT66頁

alpha_m = (1/L)*(log(1/R)); %mirror loss parameter MIT36頁

%下列計算順序N S P

power=[];%用來儲存P的矩陣

for beta = beta_sp%分別帶入三種beta值

current = [];%用來儲存電流的範圍的矩陣

NAns = [];%用來儲存N四次方程式解出的四個根矩陣

NAns1 = [];%用來儲存N四次方程式解出的四個根中的第一個根矩陣

NAns2 = [];%用來儲存N四次方程式解出的四個根中的第二個根矩陣

NAns3 = [];%用來儲存N四次方程式解出的四個根中的第三個根矩陣

NAns4 = [];%用來儲存N四次方程式解出的四個根中的第四個根矩陣

for I=0.01:0.01:30%電流的範圍

%下面解四次方程式先解出N值

coe0 = G*C;%N方程式四次係數

coe1=-1*C*(G*N_tr+1/tao_p)+B*G*(1-beta);%N方程式三次係數

coe2=-1*B*(G*N_tr+1/tao_p)+B*G*beta*N_tr+A*G;%N方程式兩次係數

coe3=-1*A(G*N_tr+1/tao_p);%N方程式一次係數

coe4=(I/(q*V))*(G*N_tr+1/tao_p);%N方程式常數項係數

eq = [coe0 coe1 coe2 coe3 coe4];%方程式矩陣

NAns = roots(eq);%求解分程式的根

NAns1 = [NAns1,NAns(1)];%儲存N四次方程式解出的四個根中的第一個根

NAns2 = [NAns2,NAns(2)];%儲存N四次方程式解出的四個根中的第一個根

NAns3 = [NAns3,NAns(3)];%儲存N四次方程式解出的四個根中的第一個根

NAns4 = [NAns4,NAns(4)];%儲存N四次方程式解出的四個根中的第一個根

current = [current,I];%儲存電流範圍的矩陣

end %下面解出S值

g = sigma_g*(NAns4-N_tr);%算出g的值

S = -1*(gama.*beta*B*(N.^2))./(gama*v_g.*g-(1/tao_p));%將算出的N跟g帶入算出S

V = V*10^-6;%轉換V的單位 從cm^3 變m^3

P = F1 * v_g *alpha_m*h*v.*S*V;%將算出的S帶入 算出P

power(i,:)=P;%儲存P的值

i = i+1;

end

power1 = power(1,:);

power2 = power(2,:);

power3 = power(3,:);

%畫圖一 Y軸為對數

figure(1);%圖一

semilogy(current,power1,'b',current,power2,'g',current,power3,'r');%Y以指數為底的畫圖

title('Lasers L-I Curve');%加整個圖的標題

xlabel('current(mA)');%加X軸標題

ylabel('Power(mW)');%加Y軸標題

legend('beta=1e-3','beta=1e-4','beta=1e-5','Location','northwest');%加圖標

grid on;%加格線

%畫圖二

figure(2);%圖二

plot(current,power1,'b',current,power2,'g',current,power3,'r')

title('Lasers L-I Curve');%加整個圖的標題

xlabel('current(mA)');%加X軸標題

ylabel('Power(mW)');%加Y軸標題

legend('beta=1e-3','beta=1e-4','beta=1e-5','Location','northwest');%加圖標

grid on;%加格線

Answer 1

[Walter Roberson]

coe3=-1*A(G*N_tr+1/tao_p);%N方程式一次係數

The A() is indexing A. You need A*()