蒙特卡洛求可靠度的一个问题。

看了论坛上的一个程序：

function [bb]=kkdmc

mm=1000;

nn=mm^2;

u=[155.7 30 40];%矩阵U第一列代表抗力R,第二列代表SG,第三列代表SQ

p=[155.7*0.15 30*0.11 40*0.15];%矩阵p的三列分别为R,SG,SQ的标准差

a1=log(lognrnd(u(1),p(1),mm,mm));%对数正态分布随机抽样

a3=normrnd(u(2),p(2),mm,mm);%正态分布随机抽样

aaa=pi/(p(3)*sqrt(6));%m为极值I型当量正态化过程中的一个常数

bbb=u(3)-0.5772/aaa;%n为极值I型当量正态化过程中的另一个常数

a=rand(mm);

a4=(aaa*bbb-log(log(1./a)))/aaa;

l=a1-a3-a4;%功能函数

o=find(l<0);%失效的次数的位置

s=length(o);%失效的次数

bb=norminv(1-s/mm^2)%可靠度指标

end

问题是nn=mm^2;表示什么意思？有什么作用？

初学者，请大家多多指点，谢谢！

Answer 1

[monej]

nn = mm^2 表示nn的值等于mm的平方。从你给出的这段程序来看，没有用到nn，看不出这一语句有何作用。