example_usage.m

% Example usage of the functions "nPSO_model" and "create_mixture_gaussian_gamma_pdf".
% The script generates and plots four examples:
% - nPSO model with 4 communities and default Gaussian mixture distribution
% - nPSO model with 4 communities and custom Gaussian mixture distribution
% - nPSO model with 8 communities and custom mixture distribution (Gaussian and Gamma mixture distribution)
% - PSO model

N = 100;
m = 4;
T = 0.1;
gamma = 3;
plot_flag = 1;

% nPSO model with 4 communities and default Gaussian mixture distribution
distr = 4;
[x_nPSO1, coords_nPSO1, comm_nPSO1, d_nPSO1] = nPSO_model(N, m, T, gamma, distr, plot_flag);

% nPSO model with 4 communities and custom Gaussian mixture distribution
C = 4; mu = (0:C-1).*(2*pi/C); sigma = (2*pi/C/6)^2; p = rand(1,C);
distr = gmdistribution(mu', sigma, p);
[x_nPSO2, coords_nPSO2, comm_nPSO2, d_nPSO2] = nPSO_model(N, m, T, gamma, distr, plot_flag);

% nPSO model with 8 communities and custom mixture distribution (Gaussian and Gamma mixture distribution)
distr = create_mixture_gaussian_gamma_pdf(8);
[x_nPSO3, coords_nPSO3, comm_nPSO3, d_nPSO3] = nPSO_model(N, m, T, gamma, distr, plot_flag);

% PSO model
distr = 0;
[x_PSO, coords_PSO, ~, d_PSO] = nPSO_model(N, m, T, gamma, distr, plot_flag);