但是代码里的randi函数使用错误，你修改一下。注意输入参数的数目要准确

您指出的确实有误。`randi`函数用于生成整数随机数，它需要两个参数：一个表示生成数目的范围（包括最小值和最大值），另一个表示随机数的取值范围。在上述代码中，`randi([min_disturbance, max_disturbance], t_end_disturbance - t_start_disturbance + 1, 1)`这一行有问题，因为它试图在一个一维向量上生成一个元素，但应该生成一个与扰动时段长度相同的向量。

正确的做法应该是：

disturbance(t_start_disturbance:t_end_disturbance) = randi([min_disturbance, max_disturbance], t_end_disturbance - t_start_disturbance + 1, 1, 'uniform', 0, 1); % 使用'uniform'选项确保均匀分布

这里增加了第三个参数1，使其生成一个与扰动时间长度相等的一维数组，并使用`'uniform'`选项确保随机数均匀分布在[min_disturbance, max_disturbance]这个范围内。

现在，`disturbance`应该是一个包含了随机扰动的正确形状的向量了。