for k = 1:Nc % Nc马氏链的条数 msg = fprintf('***Sampling from seed no. %g/%g', k, Nc); %作用是什么? 得到的是一个 n/100 的现象 theta_l(k,:) = F_seeds{:,j-1}(k,:); % 每个子集种子储存 for l = 1:Ns-1 %Ns从每个马尔可夫链模拟的样本数 为甚要 ‘-1’是不是应为下面有个加一了 theta_l(l+1,:) = MMA(theta_l,S_rnd,pi_pdf,g,b(j-1)); %【这个地方就用到MMA函数了 end theta_jk{k} = theta_l; fprintf(repmat('\b',1,msg)); end theta_j = cell2mat(theta_jk');,,需要花括号或点索引表达式中的一个输出,但结果有 10 个。 出错 SS (第 55 行) theta_l(k,:) = F_seeds{:,j-1}(k,:); % 每个子集种子储存
时间: 2024-04-23 12:23:40 浏览: 174
c文件读写函数
5星 · 资源好评率100%
这段代码是一个MATLAB语言中的for循环,其中k表示循环的计数器,从1到Nc,Nc为马氏链的条数。msg是一个输出的提示信息,作用是告诉用户当前正在采样第几个种子,输出格式为***Sampling from seed no. k/Nc。theta_l是一个矩阵,用来储存每个子集的种子,其中第一行为第k个种子,后面的行是通过MMA函数生成的样本。Ns是从每个马尔可夫链模拟的样本数,这里为什么要减1,应该是因为第一行已经是种子了,后面只
阅读全文