1.控制硬部件、集线器部件、接口部件等的状态机描述函数.docVIP

1.控制硬部件、集线器部件、接口部件等的状态机描述函数： function [state,md_time] = Series(n,lamda) %参数n表示该部件串联的个数，lamda表示错误概率的描述参数 %state表示当天该部件的状态，1表示正确，0表示错误;md_time表示需要修理元件的个数 P=1-exp(-24*lamda);%错误的概率 S=rand(1,n); S=(SP); md_time=sum(S); if md_time1 state=1; else state=0; end; 2. 控制软件部件的状态机描述函数： function [state,n] = software(lamda,P_auto,action1) %控制软部件的状态机 %lamda是描述错误概率的参数，P_auto是表示自己恢复的概率,action1表示是否采用措施1，0表示不采用，1表示采用 %n表示需要修理的天数天数，可能取值为0，1，3 P=exp(-24*lamda);%正确的概率 if action1==0%不采用措施1 if rand()P state=1;n=0; else state=0;n=3; end; else%采用措施1 if rand()P state=1;n=0; elseif rand()P_auto state=1;n=1; else state=0;n=3; end; end; 3. 包含控制硬部件，控制软部件和接口部件的串联线的描述函数： function state = kofn(k,n,lamda) %光隔部件的状态机 %参数k,n描述该k-of-n部件的参数，lamda表示错误概率的描述参数 %state表示当天该部件的状态，1表示正确，0表示错误 P_te=exp(-24*lamda);%正确的概率 te1=rand(1,n); te1=(te1P_te); te2=sum(te1); if te2=k state=1; else state=0; end 4. 包含控制硬部件，控制软部件，接口部件和光隔部件的串联线的描述函数： function [n,subfatal,k] = singleline1(action1,action2) %n是天数，k是手动修复次数，sub_fatal 子系统严重错误增量 lamda1=1/(5e4); lamda2=1/(5e3); lamda3=1/(5e3); k=0; subfatal=0; n=0; [a1,b1]=Series(20,lamda1); if a1==0 n=3; k=k+b1; if action2==0 subfatal=(rand()0.1);end;end; [a2,b2]=Series(20,lamda1); if a2==0 n=3; k=k+b2;end; [te1,te2]=software(lamda2,0.95,action1); if te2==3 n=3;k=k+1; elseif te1==0%状态是错误的 if n1 n=1; end end 5. 控制光隔部件等的状态机描述函数： function [n,subfatal,k] = singleline2(action1,action2) %n是天数，k是手动修复次数，sub_fatal 子系统严重错误增量 lamda3=1/(5e3); [n,subfatal,k] = singleline1(action1,action2); if kofn(3,5,lamda3)==0 if n1 n=1; end; end; 6．无双机热备份的代码： N=100;%样本个数 action1=0;%是否采用措施1；0表示不采用，1表示采用 action2=1;%是否采用措施2；0表示不采用，1表示采用 H=2920;%仿真天数 lamda1=1/(5e4); lamda2=1/(5e3); lamda3=1/(5e3); %最后结果变量 mendtime_mean=0;%平均人工修理次数 mean_error=0;%系统故障平均发生次数 fatal_error_mean=0;%重大故障平均发生次数 no_error_day_mean=0;%平均无故障运行总天数 continue_noerrorday_mean=0;%平均连续无故障运行的天数 continue_nofatalday_mean=0;%平均连续无重大故障运行的天数 %临时变量 fatal_error=zeros(1,N);%重大错误次数 e