clear;
clc;
fun=‘24726*x_1+45960*x_2-70686*x_3-70686*x_4‘;%%%%%%%%%功能函數形式
beta=[];
mark=[1;1;2;3];%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%隨機變量類型
for?t=1:1:50%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%所需求解時間點
????mu_1=1.3988*exp（-0.0195*（exp（0.75*t）-1）^2）*40;
????sigma_1=（0.318*0.75*t+0.9881）*5;
????mu_2=（1-0.000015*（0.75*t）^3）*400;
????sigma_2=（1+0.289*0.75*t）*20;
????mu_3=15;
????sigma_3=10;
????mu_4=40;
????sigma_4=5;
????muX=[mu_1;mu_2;mu_3;mu_4];
????sigmaX=[sigma_1;sigma_2;sigma_3;sigma_4];
????beta=[beta;JC_method_ty（muXsigmaXmarkfun）];
end
plot（1:1:50beta‘--‘‘linewidth‘2）;
hold?on
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%以上為不加固工況下的時變可靠度
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%以下為加固工況下的時變可靠度
bbeta=[];
fun=‘24726*x_1+45960*x_2-70686*x_3-70686*x_4‘;%%%%%%%%%功能函數形式
t=15;
????mu_1=1.3988*exp（-0.0195*（exp（0.75*t）-1）^2）*40;
????sigma_1=（0.318*0.75*t+0.9881）*5;
????mu_2=（1-0.000015*（0.75*t）^3）*400;
????sigma_2=（1+0.289*0.75*t）*20;
????mu_3=15;
????sigma_3=10;
????mu_4=40;
????sigma_4=5;
????muX=[mu_1;mu_2;mu_3;mu_4];
????sigmaX=[sigma_1;sigma_2;sigma_3;sigma_4];
????bbeta=[bbeta;JC_method_ty（muXsigmaXmarkfun）];
fun=‘1.005*24726*x_1+1.005*45960*x_2-70686*x_3-70686*x_4‘;%%%%%%%%%功能函數形式
for?t=15:1:55%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%加固點之后
????mu_1=1.3988*exp（-0.0195*（exp（0.75*（t-14））-1）^2）*40;
????sigma_1=（0.318*0.75*（t-14）+0.9881）*5;
????mu_2=（1-0.000015*（0.75*（t-14））^3）*400;
????sigma_2=（1+0.289*0.75*（t-14））*20;
????mu_3=15;
????sigma_3=10;
????mu_4=40;
????sigma_4=5;
????muX=[mu_1;mu_2;mu_3;mu_4];
????sigmaX=[sigma_1;sigma_2;sigma_3;sigma_4];
????bbeta=[bbeta;JC_method_ty（muXsigmaXmarkfun）];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%繪圖%%%%%%%%%%%%%%%%%%%%%%%%%
t=[1515:1:55];
plot（tbbeta‘r‘‘linewidth‘2）
xlabel（‘時間t（年）‘）
ylabel（‘可靠度\beta‘）
hold?on
t=[0:1:60];
bbeta=2.5*ones（161）;
plot（tbbeta‘k‘）
legend（‘不加固短柱‘‘加固短柱‘‘允許可靠度指標‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2105??2020-09-22?05:50??JC_Calculate\beta_time.m
?????文件????????1274??2020-09-22?05:50??JC_Calculate\beta_time1.m
?????文件????????3516??2020-09-22?05:50??JC_Calculate\JC_method_ty.m
?????文件????????1382??2020-09-22?05:50??JC_Calculate\│╠╨≥╦╡├≈.txt
?????文件??????138240??2020-09-22?05:50??JC_Calculate\╩╣╙├╩┘├?╝╞╦π╩╡└2.doc