shuju=data_ieee30; index_load=shuju.bus(:,3)>0; n_load=sum(index_load); %%负荷节点数 load=shuju.bus(index_load,[1,3,4]); mu_load_p=shuju.bus(index_load,3); sigma_load_p=0.3*mu_load_p; %%负荷标准差为期望值30% mu_load_q=shuju.bus(index_load,4); sigma_load_q=0.3*mu_load_q; m=500; %%抽样数 p_load=zeros(n_load,m); p_loss=zeros(m,1); v_mc=zeros(30,m);

这段代码中，首先从IEEE 30节点系统数据中读取负荷节点的信息，即节点编号、有功功率和无功功率，并统计负荷节点数。然后，分别计算每个负荷节点的有功功率和无功功率的期望值和标准差，并设定抽样数m=500。

接下来，使用正态分布随机数生成器normrnd()生成m个负荷样本，其中有功功率和无功功率分别服从以相应期望值为均值、以相应标准差为标准差的正态分布。生成的负荷样本用于后续的潮流计算和线路损耗计算。同时，将每个负荷样本代入潮流计算中，计算出每个样本下的系统损耗和节点电压，以及每条线路的有功和无功功率流。最后统计出m个样本下的系统损耗和节点电压，以及每条线路的有功和无功功率流的期望值和方差。