problem with Triangular distibution
Afficher commentaires plus anciens
i want to do Triangular distibution i found this online it does not seem to be working
in addtion doesmy code follow the attached 

clc
number_of_runs= 2;
x= 0.02;
syms n
syms k x
if d < 0.2
for j=(1:number_of_runs)
for k=(1:40)
for x = 0.02:.0018:0.2
pd1 = makedist('Triangular');
pdi = makedist('Triangular','a',0.05,'b',.10,'c',.15);
pdL = makedist('Triangular','a',19500,'b',20000,'c',20500);
pda = makedist('Triangular','a',5,'b',10,'c',15);
pdb = makedist('Triangular','a',1,'b',3,'c',5);
pdc = makedist('Triangular','a',1,'b',1.5,'c',2);
pdR = makedist('Triangular','a',1.69*10^-8,'b',1.7*10^-8,'c',1.7*10^-8);
pdD = makedist('Triangular','a',0.075,'b',0.08,'c',0.085);
pdP = makedist('Triangular','a',0.045,'b',0.05,'c',0.055);
pdDe = makedist('Triangular','a',8850,'b',8900,'c',8950);
pdV = makedist('Triangular','a',66000,'b',69000,'c',73000);
pdMC = makedist('Triangular','a',2.92,'b',3,'c',3.15);
pdFC = makedist('Triangular','a',500000,'b',1000000,'c',2000000);
pdT = makedist('Triangular','a',8550,'b',8600,'c',8750);
cdf1 = cdf(pd1,x);
cdfi = cdf(pdi,x);
cdfL = cdf(pdL,x);
cdfa = cdf(pda,x);
cdfb = cdf(pdb,x);
cdfc = cdf(pdc,x);
cdfR = cdf(pdR,x);
cdfD = cdf(pdD,x);
cdfP = cdf(pdP,x);
cdfDe = cdf(pdDe,x);
cdfV = cdf(pdV,x);
cdfMC = cdf(pdMC,x);
cdfFC = cdf(pdFC,x);
cdfT = cdf(pdT,x);
Io=36500*d - 160;
Rev= (3*1*Io*cdfV*cdfD*10^(-3)*(((1+cdfa)^2*(1+cdfb))^(k-1))*cdfT)/((1+cdfi)^k);
Fc=((1+cdfa)^2*(1+cdfc))^(k-1);
Fb=((1+cdfa)^2*(1+cdfb))^(k-1);
CL = (3*Io^2*cdfR*cdfL*((cdfT*cdfP*Fc)+(cdfD*Fb))/((1+cdfi)^k));
CI= cdfFC+ (cdfMC * cdfDe);
CT=CL + CI;
PV= (Rev -CT)/((1+cdfi)^k);
U= PV;
E = U*(1/number_of_runs);
end
disp (E)
disp (d)
d= d+0.0018;
end
end
Réponses (0)
Catégories
En savoir plus sur Propagation and Channel Models dans Centre d'aide et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!