Model Predictive Control for Subcarrier Power Allocation in a Simulated Multi-User System

Project info

• Field of studyTelecommunications
• LevelDoctoral dissertation (S3)
• Method typeControl Systems
• Estimated build time± 7 working days
• Toolsquadprog, optimoptions (Optimization Toolbox); c2d (Control System Toolbox) opsional
• PriceUS$31 - US$200

Research goal

Design and simulate MPC so that the total data rate tracks its target through the subcarrier power, accounting for the channel attenuation.

Variables

Example result

• 12Prediction horizon
• 0.219Tracking RMSE
• 1.114IAE
• 1.893Max |u|
• 2.857Saturated steps (%)
• 4.594Overshoot (%)

MPC with a 12-step horizon keeps Total data rate tracking the reference at RMSE 0.219 and 4.6% overshoot while honouring the Subcarrier power limits (3% of steps saturated); longer horizons reduce RMSE down to 0.219.

Computed from illustrative synthetic data, not from a real study.

About Basic Model Predictive Control

Basic Model Predictive Control solves an optimisation problem at each step to find the best control sequence over a prediction horizon while respecting constraints. It suits systems with actuator or state limits that cannot be ignored.

Steps

1. 1Define the discrete state-space model
2. 2Choose the prediction horizon, weights Q and R, and control limits
3. 3Build the prediction matrices (Phi and Gamma)
4. 4Solve the quadratic program with quadprog at each step
5. 5Apply only the first control move (receding horizon)
6. 6Plot the states and control signal

Assumptions to check

• The prediction model is sufficiently accurate
• The optimisation problem is feasible at each step
• Computation time per step is shorter than the sampling period

Reading and reporting the results

Report the horizon, weights, limits, and state and control responses that satisfy the constraints. Compare with unconstrained LQR to show the added value of MPC.

Example MATLAB code

A short starting script with this project's variable names. The delivered version is adapted to your data and fully commented.

Ts = 0.1;
A = [1 Ts; 0 1];  B = [0.5*Ts^2; Ts];             % double integrator
nx = 2; Np = 15;
Q = diag([10 1]);  R = 0.1;  umax = 1;

Phi = zeros(nx*Np, nx);  Gam = zeros(nx*Np, Np);
for i = 1:Np
    Phi((i-1)*nx+(1:nx), :) = A^i;
    for j = 1:i
        Gam((i-1)*nx+(1:nx), j) = A^(i-j) * B;
    end
end
Qb = kron(eye(Np), Q);  Rb = R * eye(Np);
H = 2 * (Gam' * Qb * Gam + Rb);  H = (H + H') / 2;

Nsim = 80;
subcarrier_power = [2; 0];                                  % initial state
X = zeros(nx, Nsim+1);  X(:,1) = subcarrier_power;  U = zeros(1, Nsim);
opts = optimoptions('quadprog', 'Display', 'off');
for k = 1:Nsim
    f = 2 * Gam' * Qb * Phi * X(:,k);
    u = quadprog(H, f, [], [], [], [], -umax*ones(Np,1), umax*ones(Np,1), [], opts);
    U(k) = u(1);                                  % receding horizon
    X(:,k+1) = A * X(:,k) + B * U(k);
end
total_data_rate = X(1,:);                                   % position trajectory
fprintf('Final position: %.4f, max |u|: %.3f\n', total_data_rate(end), max(abs(U)));

t = (0:Nsim) * Ts;
figure;
subplot(2,1,1); plot(t, total_data_rate); ylabel('Position'); grid on;
subplot(2,1,2); stairs(t(1:end-1), U); ylabel('Control u'); xlabel('Time (s)'); grid on;

Packages and prices

Choose how much help you need. Each higher package includes everything in the lower ones. Prices exclude 11% VAT.

• Payment: 50% down payment before work starts, balance after the deliverables are handed over
• Tax: Prices exclude 11% VAT
• Data and academic integrity: Results are computed from your own research data; you write and stand behind your manuscript, we help with the analysis and your understanding of it
• Revisions: Up to 2 revisions are included in the price (Package 3: 3 rounds)

Frequently asked questions