程序分析和总结.docx

《程序分析和总结.docx》由会员分享，可在线阅读，更多相关《程序分析和总结.docx（3页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、模型1 clear all; clc;p=xlsread(T bia.xlsx *, T sheetl1, 1 m3:s31*);r=xlsread ( T bla . xlsx 1 , 1 sheetl 1 , 1 e3 : k31 * ) ; %00AQoAaex4-I aEAypAEyY x=r; %OOA06Ax4-I a5SIipAEyY %23.nEIOEa for m=l:29 px=0;f=0;for i=4:7 px=p (m, i) *x (m, i) +px;end for j=l:3x (m, j ) = (1-px)*r(m,j)/p (m,j); end1=(11(

2、k)*104*35)/sqrt (-log(0.5); for i=l:7u (m, i) =l-exp(-(x(m,i)/l)八2);j (m, i) =p (m, i) *u (m, i); f=f+p(m,i)*u(m,i);endF(m)=f; end 模型1 clear all clc;F=0;x=zeros(1,7); for n=29:60 for m=5:7 p (1)=l/nchoosek(n,m);p(2)=nchoosek(m,m-1)/nchoosek(n,m); p(3)=nchoosek(mz m-1)*nchoosek(n-m-1z1)/nchoosek(n,m)

3、; p (4)=nchoosek(m,m-2)nchoosek(n-m-1,1)/nchoosek(n,m); p(5)=nchoosek(m,m-2)nchoosek(n-m-1,2)/nchoosek(n,m); p(6)=nchoosek(m,m-3)*nchoosek(n-m-1,2)/nchoosek(n,m); p(7)=nchoosek(m,m-3)nchoosek(n-m-1,3)/nchoosek(n,m);for rl=0.5:0.01:0.8for r2=(l-rl)/2) :0.01: (0.98-rl) for x4=150:10:250r3=l-rl-r2 r= r

4、l, r2z r3;x(4)=x4;x(5)=30;x(6)=10;x(7)=5;pxl=0; for i=4:7 pxl=p (i)*x (i)+pxl;end for i=l:3 x(i) = (1-pxl)*r(i)/p (i);end f=0; 1= (3.2244*104*35)/sqrt(-log(0.5);for i=l:3u (i) =l-exp (- ( (x (i) /I) A2); f=f+p(i)*u (i);end if fFF=f;J(1)=F;J(2)=m;J(3)=n;J(4)=rl;J(5)=r2;J(6)=r3;J(7)=x(1); J(8)=x(2);J(

5、9)=x(3);J(10)=x(4);endendendendendend模型2clear all;clc;p=x_srad ( * bia . xlsx 1 r 1 sheetl , * m3 : s311 ;r=xlsread ( ! bta . xlsx *, 1 sheetl1 , 1 e3 : k31 ,) ; %00AQAaex4-l aEAy)jiAEy;Yll=xlsread(Bookl.xlsx1,sheet2,1b2:h21);pl=xlsread(1Bookl.xlsx1z 1 sheet2, 1b4:h41);x=r; %OOA06Ax4-Ia5SIipAEyY%23,

6、nEI0Eafor m=l:29F (m) =0;for k=l:7px=0;f=0;for i=4:7px=p(m,i)*x(m,i)+px;endfor j=l:3x(m, j) = (1-px)*r(m,j)/p (m,j);end1=(11(k)*35)/sqrt(-log (0.5);for i=l:7u (m, i) =l-exp (- ( (x (m, i) /I) A2);f=f+p(mAi)*u(m,i);endF(m)=f*pl(k)+F(m);endend clear all clc;ll=xlsread(1Bookl.xlsx *, f sheet21, 1b2:h21

7、); pl=xlsread(1Bookl.xlsx I 1sheet2, 1b4:h41); F=0;x=zeros(1,7); for n=29:60for m=5:7p (1)=l/nchoosek(n,m);p(2)=nchoosek(mzm-l)/nchoosek(n,m);p(3)=nchoosek(mf m-1) p(4)=nchoosek(m,m-2) p(5)=nchoosek(mr m-2) p(6)=nchoosek(m,m-3) p(7)=nchoosek(m,m-3)*nchoosek(n-m-1,1)/nchoosek(n,m);*nchoosek(n-m-1f1)/

8、nchoosek(n,m);*nchoosek(n-m-1A 2)/nchoosek(n,m);*nchoosek(n-m-1,2)/nchoosek(n,m);*nchoosek(n-m-1z 3)/nchoosek(n,m);for rl=0.5:0.01:0.8for r2=( (1-rl)/2) :0.01: (0.98-rl) for x4=150:20:250fl = 0;for k=l:7r3=l-rl-r2;r=r1,r2,r3;x (4)=x4;x(5)=30;x(6)=20;x(7)=2;pxl=0;for i=4:7 pxl=p(i)*x(i)+pxl; end for i=l:3x(i) = (1-pxl)*r(i)/p (i); endf=0;1=(11 (k)*35)/sqrt(-log(0.5); for i=l:7u(i) =l-exp(-(x (i)/I)A2);f=f+p(i)*u (i); endfl=f*pl(k)+f1; endif flFF=fl;J(1)=F; J(2)=m; J(3)=n; J(4)=rl; J(5)=r2; J(6)=r3; J(7)=x(l);J(8)=x(2);J(9)=x(3);J(10)=x(4);end end end end endend