for i=1:n_load p_load(i,:)=normrnd(mu_load_p(i),sigma_load_p(i),1,m); q_load(i,:)=normrnd(mu_load_q(i),sigma_load_q(i),1,m); end

这段代码使用了for循环，对每个负荷节点进行了m次抽样，得到了每次抽样的有功功率和无功功率值，存储在p_load和q_load中。其中，对第i个负荷节点，使用normrnd函数基于上一段代码中计算得到的mu_load_p(i)和sigma_load_p(i)来生成有功功率抽样值，同时也基于mu_load_q(i)和sigma_load_q(i)来生成无功功率抽样值。这里采用了Matlab中的normrnd函数，用于生成符合正态分布的随机数。

这段代码是电力系统概率分析中的抽样部分，用于对每个负荷节点进行多次抽样，得到相应的有功功率和无功功率抽样值，为后续的概率分布计算和分析提供了基础。