结构可靠度理论代码共享
本文将记录《结构可靠度理论》的一些代码及课堂笔记，采用Matlab编程 1 概率分布密度（pdf、cdf） 自己编写函数 clc clear close all % unif a=1; b=3; x=1:0.1:3; for i=1:length(x) fx(i)=1/(b-a); Fx(i)=(x(i)-a)/(b-a); end h1=figure(1); ax1=axes('Parent',h1,'Position',[0.1,0.4,0.35,0.35],'FontSize',10); plot(x,fx,'r-') xlabel('x','fontsize',10) ylabel('pdf','FontSize',10) ax2=axes('Parent',h1,'Position',[0.55,0.4,0.35,0.35],'FontSize',10); plot(x,Fx,'r-') xlabel('x','fontsize',10) ylabel('cdf','FontSize',10) %------- mux=1;sigx=1; x=-2:0.1:4; syms t for i=1:length(x) tx=x(i); a1=1/(sigx*(2*pi)^0.5); b1=exp(-0.5*((tx-mux)/sigx)^2); fx(i)=a1*b1; c1=int(exp(-0.5*((t-mux)/sigx)^2),t,-inf,tx); Fx(i)=a1*c1; end h1=figure(2); ax1=axes('Parent',h1,'Position',[0.1,0.4,0.35,0.35],'FontSize',10); plot(x,fx,'r-') xlabel('x','fontsize',10) ylabel('pdf','FontSize',10) ax2=axes("Parent",h1,"Position",[0.55,0.4,0.35,0.35],"FontSize",10); plot(x,Fx,'r-') xlabel('x','fontsize',10) ylabel('cdf','FontSize',10) 调用Matlab自带函数 clc clear close all %----- a=1;b=3; x=-2:0.1:5; fx=unifpdf(x,a,b); Fx=unifcdf(x,a,b); figure(1) plot(x,fx,'-r') hold on plot(x,Fx,'--b') %--------- mux=0;sigx=1; x=-3:0.1:3; fx=normpdf(x,mux,sigx); Fx=normcdf(x,mux,sigx); figure(2) plot(x,fx,'r-') hold on plot(x,Fx,'r-') ...