我想要matlab代码 并且omega = rand(130, 1) *pi-pi/2;

在MATLAB中，如果你想模拟一个包含130个振子的Kuramoto网络，并且每个振子的固有频率ω是均匀分布在`(-π/2, π/2)`之间的随机值，你可以使用以下代码：

% 初始化参数
num_oscillators = 130;
max_coupling = 1; % 考虑最大耦合强度
steps = 500; % 迭代次数
coupling_range = linspace(0, max_coupling, steps); % 耦合强度范围

% 生成随机固有频率
omega = (rand(num_oscillators, 1) * pi - pi/2);  % 每个振子的固有频率在 (-π/2, π/2)

% 初始化序参数数组
order_parameter = zeros(size(coupling_range));

for i = 1:length(coupling_range)
    % 实现Kuramoto方程的数值解
    % 通常使用欧拉法或其他数值积分方法
    % 这里简化了示例，假设有一个内置函数计算序参数
    order_parameter(i) = calculateOrderParameter(omega, coupling_range(i), num_oscillators);
end

% 可能的计算序参数函数（这里仅作示例）
function op = calculateOrderParameter(omega, coupling, num_oscillators)
    % 进行实际的Kuramoto模型模拟计算序参数...
    % 在真实应用中这部分可能更复杂，这里省略
    op = mean(abs(cos(omega - circshift(omega, 1))));
end

% 绘制序参数随耦合强度的变化
plot(coupling_range, order_parameter, 'LineWidth', 2);
xlabel('Coupling Strength');
ylabel('Order Parameter');
title(sprintf('Kuramoto Network with %d Oscillators, Random Omega', num_oscillators));
grid on;