[討論] 新手初試算solar cell的IB detaile balance 請教高手迴圈已回收
小弟我剛觸碰MATLAB這好玩又實用的程式
老師這次咪聽交代我們跑論文的圖做比對
為此我已經兩天沒睡了 真的覺得會用程式很重要!
但是一直碰到問題 燈燈聲讓我的信心一再受到打擊XD
以下是我的code
想要把UCI在<=EG-EI的條件下 迴圈算200次
但是一直出現 ??? Undefined function or method 'mpower' for input arguments of type
'function_handle'.
請問一下高手哪裡錯了 或是哪裡需要改進 或是交流MSN 指導小弟 學費不會少@@
q=1.6e-19;
PS=73488034; %W/m^2
TA=300; %K
TS=6000; %K
K=8.6173324e-5; %eV/K
H=4.135667516e-15; %eV.s
C=3e8; %m/s
EI=0.5; %eV
EG=1.48; %eV
for i=1:200;
UCI=((EG-EI)/200)*i;
UIV=0.45;
A=(2*pi)/((H^3)*(C^2));
B=@(x)(x.^2)./(exp(x./(K*TS))-1);
C=@(x)(x.^2)./(exp((x-UCI-UIV)./(K*TA))-1);
D=@(x)(x.^2)./(exp((x-UCI)./(K*TA))-1);
E=@(x)(x.^2)./(exp((x-UIV)./(K*TA))-1);
F=@(x) (B(x)-C(x));
G=@(x) (B(x)-D(x));
H=@(x) (B(x)-E(x));
Q1=quadgk(F,EG,inf);
Q2=quadgk(G,(EG-EI),EG);
Q3=quadgk(H,EI,(EG-EI));
Differeence=abs(Q3-Q2);
Dlog(i,1)=Differeence;
end
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 115.43.213.57
推
10/01 18:33, , 1F
10/01 18:33, 1F
→
10/01 18:46, , 2F
10/01 18:46, 2F
→
10/01 18:48, , 3F
10/01 18:48, 3F
→
10/01 18:52, , 4F
10/01 18:52, 4F
推
10/01 19:01, , 5F
10/01 19:01, 5F