用Matlab实现随机区间的划分

[问题描述]:
随机区间的划分,有兴趣可以考虑一下问题:
将长度位len的区间划分成n段,每一段的长度是0~m的一个随机值.(len,n,m是给定的值,并n*m>len).

或者可以描述为:
产生一个数列使其满足一下要求(len,n,m是给定的值,并n*m>len)：
1、数列有n项；
2、数列的每一项使0～m的一个随机值；
3、此数列的各项之和是len；
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%matlab解决方案
%我用matlab便的程序，请高手点评
%将长度为len的区间随机分成n块,最大块不超过m（n×m>len)
%算法描述:
%1、先产生前n－1项a1、a2、a3、...、a(n-1)，每一项为0～m的随机数；
%2、下面确定最后一项：
% (1)令last = len - sum(n-1);sum(n-1)表示前n-1项的和
% (2)测试last的值;
% (3)如果last满足:last>0 且 last<m,则last就是最后一项,程序结束;
% (4)如果last<0,则从前n-1项中随机挑取一项，假设是第k项,
% 从中取出0～ak的随机数s，让last = last + s,
% (5)如果last>m,则从前n-1项中随机选取一项ak,从last取出0~(m-ak)的随机数s,
% last = last - s;
% (6)重复第(2)步;
function f = rand_seq(len,n,m)
f = m*rand(1,n) %产生1行n列的随机矩阵，每项不超过m
s = len - (sum(f) - f(n)); %最后剩余s
count = 0;
index = 0;
r = 0;
while s - m > 1e-18 || s < 1e-18
while s < 1e-18
index = int16(rand*(n-1) + 1);
t = f(index); %取第i项的值
r = rand;
f(index) = t*r; %从第i项中拿出t*(1-r)加到s里
s = s + t*(1-r);
count = count + 1;
end
while s - m > 1e-18
index = int16(rand*(n-1) + 1);
t = f(index); %取第i项的值
r = rand;
f(index) = t + (m-t)*r; %向第i项中加(m-t)*r
s = s - (m - t)*r;
count = count + 1;
end
end
f(n) = s; %最后一项确定
count %输出循环次数