lifetime.m

来自「MATSNL is a package of MATLAB M-files fo」· M 代码 · 共 26 行

function [u,p,Mt,Me,Pe]=lifetime(d)

for k=1:1000
a=1;
Tc=uniform_gen(1, 60);%2sec
y=uniform_gen(0.01, 100);%2sec
z=uniform_gen(1*10^-3, 10);%100msec
P0=uniform_gen(0,100); %0.01mW
P3=uniform_gen(300,1000); %315mW
P4=uniform_gen(100, 300); %132mW
Cr=6.63*10^-3; %6.63uJ
Cp=uniform_gen(1,100);%43mJ
e=0; %100msec
L=uniform_gen(Tc*d, 10000*Tc*d);
lambda=1./L;

Mt=y+a*z-(Tc*d.*(Tc*d+2*e)./(2*Tc*d+2*e));
Me=y*P4+a*(z*P3+Cr)-(Tc*d.*(Tc*d+2*e)./(2*Tc*d+2*e))*P4;
Pe=P0+d*(P4-P0)+Cp./Tc;
u=min(1,(d+e./Tc));
p=(Pe+lambda.*u.*Me)./(1+lambda.*u.*Mt);
if min(Pe)<Me/Mt
plot(d,p,'r');
end
hold on;
end