TADA(8) TADA(16) TADA(1) TADA(2) TADA(1) TADA(2) TADA(3) TADA(4) TADA(5) TADA(16) TADA(0) TADA(0,8) TADA(3,8) TADA(334,8) TADA(33,8) TADA(47,8) %-- 3/29/07 12:07 PM --% fuzzy mojladownik d: cd d: cd Studenci\flak\ fuzzy mojladownik fuzzy fuzzy flakownik run sum1 run sym1 run sym1.m sym1 help fibnd help find %-- 4/03/07 10:25 AM --% fuzzy wahadlo fiswahadla wahadlo wahadlo(fiswahadlo) wahadlo(fiswahadla) wahadlo fiswahadla %-- 4/17/07 10:18 AM --% bayes PW = 0.002; PU = 0.5; PT = 0.06; PS_WU = 0.3; PS_nWU = 0.2; PS_WnU = 0.07; PS_nWnU = 0.02; PA_T = 0.1; PA_nT = 0.01; PZ_SA = 0.95; PZ_nSA = 0.90; PZ_Sn = 0.40; PZ_nSnA = 0.11; PW = 0.002; PW(0.002) PW = 0.002; function ps_anal(){ ps_anal = PS_WU*PW*PU + PS_nWU*(1-PW)*PU + PS_WnU*PW*(1-PU) + PS_nWnU*(1-PW)*(1-PU) for i=1 to 10000 for i=1,10000 zW = rand<pW end zW = rand < PW zW = rand(1,10000) < PW sum(ans) sum(zW) zW = rand(1,100000) < PW zU = rand(1,100000) < PU zT = rand(1,100000) < PT zS_WU = rand(1,100000) < PS_WU zS_WnU = rand(1,100000) < PS_WnU zS_nWU = rand(1,100000) < PS_nWU zS_nWnU = rand(1,100000) < PS_nWnU ps_bayes = sum(zS_WU) + sum(zS_WnU) + sum(zS_nWU) + sum(zS_nWnU) ps_bayes/100000 zW = rand(1,100000) < PW; zU = rand(1,100000) < PU; zT = rand(1,100000) < PT; zS_WU = rand(1,100000) < PS_WU; zS_WnU = rand(1,100000) < PS_WnU; zS_nWU = rand(1,100000) < PS_nWU; zS_nWnU = rand(1,100000) < PS_nWnU; sum(zW.*zS_WU + zU.*zWU).*zS_WU) + sum(zW.*zS_WnU + (1-zU).*zS_WnU) + sum((1-zW).*zS_nWU + zU.*zS_nWU) + sum((1-zW).*zS_nWnU + (1-zU).*zS_nWnU) zW = rand(1,100000) < PW; zU = rand(1,100000) < PU; zT = rand(1,100000) < PT; zS_WU = rand(1,100000) < PS_WU; zS_WnU = rand(1,100000) < PS_WnU; zS_nWU = rand(1,100000) < PS_nWU; zS_nWnU = rand(1,100000) < PS_nWnU; sum(zW.*zS_WU + zU.*zWU) + sum(zW.*zS_WnU + (1-zU).*zS_WnU) + sum((1-zW).*zS_nWU + zU.*zS_nWU) + sum((1-zW).*zS_nWnU + (1-zU).*zS_nWnU) zW = rand(1,100000) < PW; zU = rand(1,100000) < PU; zT = rand(1,100000) < PT; zS_WU = rand(1,100000) < PS_WU; zS_WnU = rand(1,100000) < PS_WnU; zS_nWU = rand(1,100000) < PS_nWU; zS_nWnU = rand(1,100000) < PS_nWnU; sum(zW.*zS_WU + zU.*zS?WU) + sum(zW.*zS_WnU + (1-zU).*zS_WnU) + sum((1-zW).*zS_nWU + zU.*zS_nWU) + sum((1-zW).*zS_nWnU + (1-zU).*zS_nWnU) zW = rand(1,100000) < PW; zU = rand(1,100000) < PU; zT = rand(1,100000) < PT; zS_WU = rand(1,100000) < PS_WU; zS_WnU = rand(1,100000) < PS_WnU; zS_nWU = rand(1,100000) < PS_nWU; zS_nWnU = rand(1,100000) < PS_nWnU; sum(zW.*zS_WU + zU.*zS_WU) + sum(zW.*zS_WnU + (1-zU).*zS_WnU) + sum((1-zW).*zS_nWU + zU.*zS_nWU) + sum((1-zW).*zS_nWnU + (1-zU).*zS_nWnU) sum(zW.*zS_WU + zU.*zS_WU) + sum(zW.*zS_WnU + (1-zU).*zS_WnU) + sum((1-zW).*zS_nWU + zU.*zS_nWU) + sum((1-zW).*zS_nWnU + (1-zU).*zS_nWnU)/100000 zW = rand(1,100000) < PW; zU = rand(1,100000) < PU; zT = rand(1,100000) < PT; zS_WU = rand(1,100000) < PS_WU; zS_WnU = rand(1,100000) < PS_WnU; zS_nWU = rand(1,100000) < PS_nWU; zS_nWnU = rand(1,100000) < PS_nWnU; sum(zW.*zS_WU.*zU) + sum(zW.*zS_WnU.*(1-zU)) + sum((1-zW).*zS_nWU.*zU) + sum((1-zW).*zS_nWnU.*(1-zU)) sum(zW.*zS_WU.*zU) + sum(zW.*zS_WnU.*(1-zU)) + sum((1-zW).*zS_nWU.*zU) + sum((1-zW).*zS_nWnU.*(1-zU))/100000 ans/100000 sum(zW.*zS_WU.*zU) + sum(zW.*zS_WnU.*(1-zU)) + sum((1-zW).*zS_nWU.*zU) + sum((1-zW).*zS_nWnU.*(1-zU))/100000 sum(zW.*zS_WU.*zU) + sum(zW.*zS_WnU.*(1-zU)) + sum((1-zW).*zS_nWU.*zU) + sum((1-zW).*zS_nWnU.*(1-zU)) ans/100000 PZ_SnA = 0.40; zW = rand(1,100000) < PW; zU = rand(1,100000) < PU; zT = rand(1,100000) < PT; zS_WU = rand(1,100000) < PS_WU; zS_WnU = rand(1,100000) < PS_WnU; zS_nWU = rand(1,100000) < PS_nWU; zS_nWnU = rand(1,100000) < PS_nWnU; zA_T = rand(1,100000) < PA_T; zA_nT = rand(1,100000) < PA_nT; zZ_SA = rand(1,100000) < PZ_SA; zZ_nSA = rand(1,100000) < PZ_nSA; zZ_SnA = rand(1,100000) < PZ_SnA; zZ_nSnA = rand(1,100000) < PZ_nSnA; zA = zT.*zA_T + zA_nT.*(1-zT) zS = zW.*zS_WU.*zU + zW.*zS_WnU.*(1-zU) + (1-zW).*zS_nWU.*zU + (1-zW).*zS_nWnU.*(1-zU) zZ = (zS.*zZ_SA.*zA + zW.*zZ_SnA.*(1-zA) + (1-zS).*zZ_nSA.*zA + (1-zA).*zZ_nSnA.*(1-zA) zZ_U = zZ.*zU zW = rand(1,100000) < PW; zU = rand(1,100000) < PU; zT = rand(1,100000) < PT; zS_WU = rand(1,100000) < PS_WU; zS_WnU = rand(1,100000) < PS_WnU; zS_nWU = rand(1,100000) < PS_nWU; zS_nWnU = rand(1,100000) < PS_nWnU; zA_T = rand(1,100000) < PA_T; zA_nT = rand(1,100000) < PA_nT; zZ_SA = rand(1,100000) < PZ_SA; zZ_nSA = rand(1,100000) < PZ_nSA; zZ_SnA = rand(1,100000) < PZ_SnA; zZ_nSnA = rand(1,100000) < PZ_nSnA; zA = zT.*zA_T + zA_nT.*(1-zT); zS = zW.*zS_WU.*zU + zW.*zS_WnU.*(1-zU) + (1-zW).*zS_nWU.*zU + (1-zW).*zS_nWnU.*(1-zU); zZ = (zS.*zZ_SA.*zA + zW.*zZ_SnA.*(1-zA) + (1-zS).*zZ_nSA.*zA + (1-zA).*zZ_nSnA.*(1-zA); zZ_U = zZ.*zU; zW = rand(1,100000) < PW; zU = rand(1,100000) < PU; zT = rand(1,100000) < PT; zS_WU = rand(1,100000) < PS_WU; zS_WnU = rand(1,100000) < PS_WnU; zS_nWU = rand(1,100000) < PS_nWU; zS_nWnU = rand(1,100000) < PS_nWnU; zA_T = rand(1,100000) < PA_T; zA_nT = rand(1,100000) < PA_nT; zZ_SA = rand(1,100000) < PZ_SA; zZ_nSA = rand(1,100000) < PZ_nSA; zZ_SnA = rand(1,100000) < PZ_SnA; zZ_nSnA = rand(1,100000) < PZ_nSnA; zA = zT.*zA_T + zA_nT.*(1-zT); zS = zW.*zS_WU.*zU + zW.*zS_WnU.*(1-zU) + (1-zW).*zS_nWU.*zU + (1-zW).*zS_nWnU.*(1-zU); zZ = zS.*zZ_SA.*zA + zW.*zZ_SnA.*(1-zA) + (1-zS).*zZ_nSA.*zA + (1-zA).*zZ_nSnA.*(1-zA); zZ_U = zZ.*zU; sum(zZ_U) ans/100000 zZ_UT = ZZ.*zU.*zT zZ_UT = zZ.*zU.*zT sum (zZ_UT) ans / 100000 zW = rand(1,100000) < PW; zU = rand(1,100000) < PU; zT = rand(1,100000) < PT; zS_WU = rand(1,100000) < PS_WU; zS_WnU = rand(1,100000) < PS_WnU; zS_nWU = rand(1,100000) < PS_nWU; zS_nWnU = rand(1,100000) < PS_nWnU; zA_T = rand(1,100000) < PA_T; zA_nT = rand(1,100000) < PA_nT; zZ_SA = rand(1,100000) < PZ_SA; zZ_nSA = rand(1,100000) < PZ_nSA; zZ_SnA = rand(1,100000) < PZ_SnA; zZ_nSnA = rand(1,100000) < PZ_nSnA; zA = zT.*zA_T + zA_nT.*(1-zT); zS = zW.*zS_WU.*zU + zW.*zS_WnU.*(1-zU) + (1-zW).*zS_nWU.*zU + (1-zW).*zS_nWnU.*(1-zU); zZ = zS.*zZ_SA.*zA + zW.*zZ_SnA.*(1-zA) + (1-zS).*zZ_nSA.*zA + (1-zA).*zZ_nSnA.*(1-zA); ilZ = sum(zZ) ilU = sum(zU) ilU/ilZ ilZ/ilU sum(zZ_U) pZU = ans/100000 pZ = ilZ/100000 pZU/pZ zW = rand(1,100000) < PW; zU = rand(1,100000) < PU; zT = rand(1,100000) < PT; zS_WU = rand(1,100000) < PS_WU; zS_WnU = rand(1,100000) < PS_WnU; zS_nWU = rand(1,100000) < PS_nWU; zS_nWnU = rand(1,100000) < PS_nWnU; zA_T = rand(1,100000) < PA_T; zA_nT = rand(1,100000) < PA_nT; zZ_SA = rand(1,100000) < PZ_SA; zZ_nSA = rand(1,100000) < PZ_nSA; zZ_SnA = rand(1,100000) < PZ_SnA; zZ_nSnA = rand(1,100000) < PZ_nSnA; zA = zT.*zA_T + zA_nT.*(1-zT); zS = zW.*zS_WU.*zU + zW.*zS_WnU.*(1-zU) + (1-zW).*zS_nWU.*zU + (1-zW).*zS_nWnU.*(1-zU); zZ = zS.*zZ_SA.*zA + zW.*zZ_SnA.*(1-zA) + (1-zS).*zZ_nSA.*zA + (1-zA).*zZ_nSnA.*(1-zA); zZ_U = zZ.*zU ilZU = sum(zZ_U) ilZ = sum(zZ) pZU = ilZU/100000 pzZ = ilZ/100000 pZU/pZZ pZU/pzZ %-- 4/21/07 11:40 AM --% cd D:\Studenci\Sidzi\ addpath 'D:\studenci\Sidzi\scheduling\' addpath 'D:\studenci\sidzi\grTheory\' addpath 'D:\studenci\Sidzi\scheduling\' t1=task('t1',5,2,10,2) plot(t1) t1=task('t1',5,2,INF,2) t1=task('t1',5,2,inf,2) t1=task('t1',2,2,inf,2) t2=task('t2',3,0,inf,3) t3=task('t3',4,3,inf,4) t4=task('t4',2,0,inf,5) T=taskset(t1,t2,y3,t4) T=taskset(t1,t2,t3,t4) T=taskset{[t1,t2,t3,t4]} T=taskset([t1,t2,t3,t4]) plot(T) T.procTime T.proctime T.proctime[1]=1 T.proctime(1)=1 T.procTime(1)=1 addpath 'D:\studenci\Sidzi\scheduling\' addpath 'D:\studenci\sidzi\grTheory\' T.procTime(1)=1 T.ProcTime(1)=1 p=problem('P|rj|Cmax') h=listsch(T,p,2) plot(T) T2=T st_t=[3 0 3 0] pr_t=[1 3 4 2] pr=[2 2 1 1] add_schedule(T2,'uszeregowanie niepodzielne',st_t,pr_t,pr);plot(T2) plot(T) plot(T2) Tprec=taskset([t1 t2 t3], [0 1 0; 0 0 1; 0 0 0]) plot(Tprec) Tprec=taskset([t1 t2 t3], [0 1 0; 0 0 1; 1 0 0]) plot(Tprec) Tprec=taskset([t1 t2 t3], [0 1 0; 0 0 1; 0 0 0]) Tprec=taskset([t1 t2 t3], [0 1 0; 1 0 1; 0 0 0]) plot(Tprec) Tprec=taskset([t1 t2 t3], [0 1 0; 0 0 1; 0 0 0]) plot(Tprec) st_t={[2] [0 3] [3] [0]} pr_t={[1] [2 1] [4] [2]} pr=[ 2 2 1 1] add_schedule(T2,'uszeregowanie podzielne',st_t,pr_t,pr); plot(T2) graphedit taskset ? help taskset T3 = taskset(g1) plot(T3) clear cls clean clrscr help screen pj=[2 2 3 1 4 5 2 1 2 3] rj=[0 1 3 4 5 0 2 3 4 5] dj=[5 4 6 5 3 6 3 4 5 8] p=problem('P2|rj|Cmax'); help taskset graphedit t1=task('t1'2,0,inf,5) t1=task('t1',2,0,inf,5) t2=task('t2',2,1,inf,4) t3=task('t3',3,3,inf,6) t4=task('t4',1,5,inf,3) t5=task('t5',4,0,inf,6) t6=task('t6',5,2,inf,3) pj=[2 2 3 1 4 5 2 1 2 3] rj=[0 1 3 4 5 0 2 3 4 5] dj=[5 4 6 5 3 6 3 4 5 8] pj=[2 2 3 1 4 5 2 1 2 3]; rj=[0 1 3 4 5 0 2 3 4 5]; dj=[5 4 6 5 3 6 3 4 5 8] m={pj,rj,dj} m m(,1) m(1,1) m=[pj,rj,dj] pj={2,2,3,1,4,5,2,1,2,3} rj={0,1,3,4,5,0,2,3,4,5} dj={5,4,6,5,3,6,3,4,5,8} pj={2,2,3,1,4,5,2,1,2,3}; rj={0,1,3,4,5,0,2,3,4,5}; dj={5,4,6,5,3,6,3,4,5,8} m={pj,rj,dj} m m(0,0) m(1,1) m={pj rj dj} t1=task(2,0,inf,5) t2=task(2,1,inf,4) t3=task(3,3,inf,6) t3=task(1,4,inf,5) t3=task(4,5,inf,3) t3=task(5,0,inf,6) t3=task(2,2,inf,3) t3=task(1,3,inf,4) t4=task(1,4,inf,5); t5=task(4,5,inf,3); t6=task(5,0,inf,6); t7=task(2,2,inf,3); t8=task(1,3,inf,4); t9=task(2,4,inf,5); t10=task(3,5,inf,8); T=taskset(t1,t2,t3,t4,t5,t6,t7,t8,t9,t10);h=listsch(T,p,