There are four cases that I am trying to incorporate as multiple IF statements. Here are the four cases below. The first case is the one I used in my code to generate the output.
dx(1,1) > 0 && dx(3,1) > 0
dx(1,1) > 0 || dx(3,1) > 0
dx(1,1) > 0 && dx(3,1) <= 0
dx(1,1) <= 0 && dx(3,1) > 0
Below is the code that I used to generate the first case:
close all
clear all
clc

tspan = [0 100];
y0 = 0.7; % initial value of state variable x1
r0 = 0.7; % initial value of state variable x2
w0 = 0.8; % initial value of state variable x3
x0 = [y0; r0; w0];
[t, x] = ode45(@odefcn, tspan, x0);

%% Plot results
figure
subplot(3,1,1);
plot(t, x(:,1)); grid on
xlabel(‘Time’), ylabel(‘Output’);

subplot(3,1,2);
plot(t, x(:,2)); grid on
xlabel(‘Time’), ylabel(‘Policy Rate’);

subplot(3,1,3);
plot(t, x(:,3)); grid on
xlabel(‘Time’), ylabel(‘Wage Share’);
y0 = [0.7; 0.7; 0.8];

%% System of three differential equations
function dx = odefcn(t, x)
% definitions
y = x(1);
r = x(2);
w = x(3);

% parameters
alpha = 1.0;
beta = 1.0;
gamma = 0.6;
delta = 0.6;
mu = 0.1;
lambda = 0.6;
theta = 0.4;
omega = 0.4;
sigma = 0.1;
tau = 0.1;

% ODEs
dx(1,1) = alpha * y – beta * r * y;
dx(3,1) = – theta * w + lambda * y * w – mu * w * w;

% Asymmetrical Reaction Function
if dx(1,1) > 0 && dx(3,1) > 0
dx(2,1) = – omega * r + gamma * w * r + sigma * w * r + delta * y * r + tau * y * r;
else
dx(2,1) = – omega * r + gamma * w * r + delta * y * r;
end
endThere are four cases that I am trying to incorporate as multiple IF statements. Here are the four cases below. The first case is the one I used in my code to generate the output.
dx(1,1) > 0 && dx(3,1) > 0
dx(1,1) > 0 || dx(3,1) > 0
dx(1,1) > 0 && dx(3,1) <= 0
dx(1,1) <= 0 && dx(3,1) > 0
Below is the code that I used to generate the first case:
close all
clear all
clc

tspan = [0 100];
y0 = 0.7; % initial value of state variable x1
r0 = 0.7; % initial value of state variable x2
w0 = 0.8; % initial value of state variable x3
x0 = [y0; r0; w0];
[t, x] = ode45(@odefcn, tspan, x0);

%% Plot results
figure
subplot(3,1,1);
plot(t, x(:,1)); grid on
xlabel(‘Time’), ylabel(‘Output’);

subplot(3,1,2);
plot(t, x(:,2)); grid on
xlabel(‘Time’), ylabel(‘Policy Rate’);

subplot(3,1,3);
plot(t, x(:,3)); grid on
xlabel(‘Time’), ylabel(‘Wage Share’);
y0 = [0.7; 0.7; 0.8];

%% System of three differential equations
function dx = odefcn(t, x)
% definitions
y = x(1);
r = x(2);
w = x(3);

% parameters
alpha = 1.0;
beta = 1.0;
gamma = 0.6;
delta = 0.6;
mu = 0.1;
lambda = 0.6;
theta = 0.4;
omega = 0.4;
sigma = 0.1;
tau = 0.1;

% ODEs
dx(1,1) = alpha * y – beta * r * y;
dx(3,1) = – theta * w + lambda * y * w – mu * w * w;

% Asymmetrical Reaction Function
if dx(1,1) > 0 && dx(3,1) > 0
dx(2,1) = – omega * r + gamma * w * r + sigma * w * r + delta * y * r + tau * y * r;
else
dx(2,1) = – omega * r + gamma * w * r + delta * y * r;
end
end There are four cases that I am trying to incorporate as multiple IF statements. Here are the four cases below. The first case is the one I used in my code to generate the output.
dx(1,1) > 0 && dx(3,1) > 0
dx(1,1) > 0 || dx(3,1) > 0
dx(1,1) > 0 && dx(3,1) <= 0
dx(1,1) <= 0 && dx(3,1) > 0
Below is the code that I used to generate the first case:
close all
clear all
clc

tspan = [0 100];
y0 = 0.7; % initial value of state variable x1
r0 = 0.7; % initial value of state variable x2
w0 = 0.8; % initial value of state variable x3
x0 = [y0; r0; w0];
[t, x] = ode45(@odefcn, tspan, x0);

%% Plot results
figure
subplot(3,1,1);
plot(t, x(:,1)); grid on
xlabel(‘Time’), ylabel(‘Output’);

subplot(3,1,2);
plot(t, x(:,2)); grid on
xlabel(‘Time’), ylabel(‘Policy Rate’);

subplot(3,1,3);
plot(t, x(:,3)); grid on
xlabel(‘Time’), ylabel(‘Wage Share’);
y0 = [0.7; 0.7; 0.8];

%% System of three differential equations
function dx = odefcn(t, x)
% definitions
y = x(1);
r = x(2);
w = x(3);

% parameters
alpha = 1.0;
beta = 1.0;
gamma = 0.6;
delta = 0.6;
mu = 0.1;
lambda = 0.6;
theta = 0.4;
omega = 0.4;
sigma = 0.1;
tau = 0.1;

% ODEs
dx(1,1) = alpha * y – beta * r * y;
dx(3,1) = – theta * w + lambda * y * w – mu * w * w;

% Asymmetrical Reaction Function
if dx(1,1) > 0 && dx(3,1) > 0
dx(2,1) = – omega * r + gamma * w * r + sigma * w * r + delta * y * r + tau * y * r;
else
dx(2,1) = – omega * r + gamma * w * r + delta * y * r;
end
end ode, differential equations, if statement MATLAB Answers — New Questions

​

I am plotting orbital data, using plot3 to progressively lengthen the orbital track of a spacecraft.
plot3(earthMoon.UIAxes, x,y,z,"w’);
Each time I do this, another "Children’ is added to the UIAxes instance. Eventually, when the number of children approaches 200 or so, the graphical object management becomes unwieldy and my GUI stops responding to button pushes; even the Matlab STOP and CONTINUE buttons fail to respond.
I found a similar question and answer here:
https://www.mathworks.com/matlabcentral/answers/433845-maximum-number-of-uiaxes-in-app-designer-2017b?s_tid=prof_contriblnk
Is there a better way to manage the Children? I have tried thinning out my track data to reducte the number of points, but that results in jerky updates that don’t align with the orbital state. I have also looked at using splines, but I will still need to limit the number of current track points to something under 200.
I have another app for displaying radar data where I progressively delete the oldest children in the UIApp, so I never have more than about 50. That works well, but for orbital data I would like to display the whole trajectory from launch to current time. It appears I can’t have it both ways.I am plotting orbital data, using plot3 to progressively lengthen the orbital track of a spacecraft.
plot3(earthMoon.UIAxes, x,y,z,"w’);
Each time I do this, another "Children’ is added to the UIAxes instance. Eventually, when the number of children approaches 200 or so, the graphical object management becomes unwieldy and my GUI stops responding to button pushes; even the Matlab STOP and CONTINUE buttons fail to respond.
I found a similar question and answer here:
https://www.mathworks.com/matlabcentral/answers/433845-maximum-number-of-uiaxes-in-app-designer-2017b?s_tid=prof_contriblnk
Is there a better way to manage the Children? I have tried thinning out my track data to reducte the number of points, but that results in jerky updates that don’t align with the orbital state. I have also looked at using splines, but I will still need to limit the number of current track points to something under 200.
I have another app for displaying radar data where I progressively delete the oldest children in the UIApp, so I never have more than about 50. That works well, but for orbital data I would like to display the whole trajectory from launch to current time. It appears I can’t have it both ways. I am plotting orbital data, using plot3 to progressively lengthen the orbital track of a spacecraft.
plot3(earthMoon.UIAxes, x,y,z,"w’);
Each time I do this, another "Children’ is added to the UIAxes instance. Eventually, when the number of children approaches 200 or so, the graphical object management becomes unwieldy and my GUI stops responding to button pushes; even the Matlab STOP and CONTINUE buttons fail to respond.
I found a similar question and answer here:
https://www.mathworks.com/matlabcentral/answers/433845-maximum-number-of-uiaxes-in-app-designer-2017b?s_tid=prof_contriblnk
Is there a better way to manage the Children? I have tried thinning out my track data to reducte the number of points, but that results in jerky updates that don’t align with the orbital state. I have also looked at using splines, but I will still need to limit the number of current track points to something under 200.
I have another app for displaying radar data where I progressively delete the oldest children in the UIApp, so I never have more than about 50. That works well, but for orbital data I would like to display the whole trajectory from launch to current time. It appears I can’t have it both ways. plot3, uiaxes, uifigure MATLAB Answers — New Questions

​

I want to add semi-colons to multiple lines of MATLAB code automatically.I want to add semi-colons to multiple lines of MATLAB code automatically. I want to add semi-colons to multiple lines of MATLAB code automatically. semicolon, multilinesemicolon MATLAB Answers — New Questions

​

Hello, I want to specify intervals of increasing & decreasing of a function for a specific range; say -20-20. This is my code here:
clc, close, clear
syms x y real

dy=1-(4/x^2); x=~0; % This is the derivative of the function
% dy=(sin(x)-1)*(2*cos(x)+1); % (0<=x<=2*pi) % Code doesn’t work for this.

fplot(dy); hold on; axis([-5, 5, -5, 5]);
j=xline(0); set(j,’color’,’black’,’Linewidth’, 1.5)
k=yline(0); set(k,’color’,’black’,’Linewidth’, 1.5)
cp=solve(dy); % Finding critical points
y_cp=subs(dy,cp);
[~,idx] = unique(cp);
cp = cp(idx);
y_cp = y_cp(idx);
[num, den]=numden(dy); % Not a conventional way to find other critical points
cp2=solve(den==0);
disp(‘The critical points are:’)
for k =1:length(cp)
fprintf("(%.2f,%.2f)n",cp(k), y_cp(k))
plot(cp,y_cp,’*’, ‘color’, ‘black’);
end
for k =1:length(cp2)
fprintf("x = %.2f is a critical point, y is undefinedn",cp2(k))
end
% y is increasing if dy > 0, decreasing if dy < 0, if dy=0 it’s a critical point. Code below doesn’t work.
for i = -20:20
dy_sub_i=subs(dy,i);
end
if dy_sub_i>0
fprintf(‘f is increasing in %s-%s’, i(0), i);
elseif dy_sub_i<0
fprintf(‘f is decreasing in ‘);
endHello, I want to specify intervals of increasing & decreasing of a function for a specific range; say -20-20. This is my code here:
clc, close, clear
syms x y real

dy=1-(4/x^2); x=~0; % This is the derivative of the function
% dy=(sin(x)-1)*(2*cos(x)+1); % (0<=x<=2*pi) % Code doesn’t work for this.

fplot(dy); hold on; axis([-5, 5, -5, 5]);
j=xline(0); set(j,’color’,’black’,’Linewidth’, 1.5)
k=yline(0); set(k,’color’,’black’,’Linewidth’, 1.5)
cp=solve(dy); % Finding critical points
y_cp=subs(dy,cp);
[~,idx] = unique(cp);
cp = cp(idx);
y_cp = y_cp(idx);
[num, den]=numden(dy); % Not a conventional way to find other critical points
cp2=solve(den==0);
disp(‘The critical points are:’)
for k =1:length(cp)
fprintf("(%.2f,%.2f)n",cp(k), y_cp(k))
plot(cp,y_cp,’*’, ‘color’, ‘black’);
end
for k =1:length(cp2)
fprintf("x = %.2f is a critical point, y is undefinedn",cp2(k))
end
% y is increasing if dy > 0, decreasing if dy < 0, if dy=0 it’s a critical point. Code below doesn’t work.
for i = -20:20
dy_sub_i=subs(dy,i);
end
if dy_sub_i>0
fprintf(‘f is increasing in %s-%s’, i(0), i);
elseif dy_sub_i<0
fprintf(‘f is decreasing in ‘);
end Hello, I want to specify intervals of increasing & decreasing of a function for a specific range; say -20-20. This is my code here:
clc, close, clear
syms x y real

dy=1-(4/x^2); x=~0; % This is the derivative of the function
% dy=(sin(x)-1)*(2*cos(x)+1); % (0<=x<=2*pi) % Code doesn’t work for this.

fplot(dy); hold on; axis([-5, 5, -5, 5]);
j=xline(0); set(j,’color’,’black’,’Linewidth’, 1.5)
k=yline(0); set(k,’color’,’black’,’Linewidth’, 1.5)
cp=solve(dy); % Finding critical points
y_cp=subs(dy,cp);
[~,idx] = unique(cp);
cp = cp(idx);
y_cp = y_cp(idx);
[num, den]=numden(dy); % Not a conventional way to find other critical points
cp2=solve(den==0);
disp(‘The critical points are:’)
for k =1:length(cp)
fprintf("(%.2f,%.2f)n",cp(k), y_cp(k))
plot(cp,y_cp,’*’, ‘color’, ‘black’);
end
for k =1:length(cp2)
fprintf("x = %.2f is a critical point, y is undefinedn",cp2(k))
end
% y is increasing if dy > 0, decreasing if dy < 0, if dy=0 it’s a critical point. Code below doesn’t work.
for i = -20:20
dy_sub_i=subs(dy,i);
end
if dy_sub_i>0
fprintf(‘f is increasing in %s-%s’, i(0), i);
elseif dy_sub_i<0
fprintf(‘f is decreasing in ‘);
end for loop, fplot, fprintf MATLAB Answers — New Questions

​

I’m experiencing a problem with channel interference when using the Analog Input block of Simulink Desktop Real-Time (SDRT). Although it initially appears to be a case of ‘crosstalk’, the problem also occurs when measuring DC signals from very low impedance sources compared to the input impedance of the data acquisition boards used (National Instruments).
Test description
– I set up several simple voltage dividers (two 10 kΩ resistors in series) powered by an external 5 V supply.
– The voltage at the center of the dividers was measured using a single Analog Input block configured to access multiple channels in vector form [1:n].
– When I short-circuited the resistor connected to GND on the channel (1), its voltage dropped to zero as expected, but I observed a change of around 30 mV in the voltage measured on the channel (2).
– This behavior occurs consistently between any channel (n) and the subsequent channel (n+1).
Additional Tests
– I repeated the test using independent Analog Input blocks for each channel, and, in this case, the problem did not occur for any of the channels, demonstrating that the behavior is related to the vector configuration of the block.
– However, according to the SDRT user guide and the recommendation of MathWorks experts in the past, using a single Analog Input block for multiple signals, using vector notation for channel selection. This approach is necessary to maintain real-time operation, as using independent blocks introduces significant delays, even at low sampling frequencies.
Additional Information
– The same tests were carried out on two different computes (one slower and one faster and robut), with different cards of acquisition (PCI-6221 and PCIe-6323) and in two versions of MATLAB (2021b and 2024a).
– I also randomly changed the order of the input channels to rule out possible problems related to the sequence of the channels.
– The behavior persisted in all the configurations tested.
Requesting Assistance
Are there any known limitations or solutions to this behavior when using the Analog Input block in vector configuration? Considering that SDRT recommends this approach for multiple channels, would a technical review or guidance be possible to solve or mitigate the problem?
.I’m experiencing a problem with channel interference when using the Analog Input block of Simulink Desktop Real-Time (SDRT). Although it initially appears to be a case of ‘crosstalk’, the problem also occurs when measuring DC signals from very low impedance sources compared to the input impedance of the data acquisition boards used (National Instruments).
Test description
– I set up several simple voltage dividers (two 10 kΩ resistors in series) powered by an external 5 V supply.
– The voltage at the center of the dividers was measured using a single Analog Input block configured to access multiple channels in vector form [1:n].
– When I short-circuited the resistor connected to GND on the channel (1), its voltage dropped to zero as expected, but I observed a change of around 30 mV in the voltage measured on the channel (2).
– This behavior occurs consistently between any channel (n) and the subsequent channel (n+1).
Additional Tests
– I repeated the test using independent Analog Input blocks for each channel, and, in this case, the problem did not occur for any of the channels, demonstrating that the behavior is related to the vector configuration of the block.
– However, according to the SDRT user guide and the recommendation of MathWorks experts in the past, using a single Analog Input block for multiple signals, using vector notation for channel selection. This approach is necessary to maintain real-time operation, as using independent blocks introduces significant delays, even at low sampling frequencies.
Additional Information
– The same tests were carried out on two different computes (one slower and one faster and robut), with different cards of acquisition (PCI-6221 and PCIe-6323) and in two versions of MATLAB (2021b and 2024a).
– I also randomly changed the order of the input channels to rule out possible problems related to the sequence of the channels.
– The behavior persisted in all the configurations tested.
Requesting Assistance
Are there any known limitations or solutions to this behavior when using the Analog Input block in vector configuration? Considering that SDRT recommends this approach for multiple channels, would a technical review or guidance be possible to solve or mitigate the problem?
. I’m experiencing a problem with channel interference when using the Analog Input block of Simulink Desktop Real-Time (SDRT). Although it initially appears to be a case of ‘crosstalk’, the problem also occurs when measuring DC signals from very low impedance sources compared to the input impedance of the data acquisition boards used (National Instruments).
Test description
– I set up several simple voltage dividers (two 10 kΩ resistors in series) powered by an external 5 V supply.
– The voltage at the center of the dividers was measured using a single Analog Input block configured to access multiple channels in vector form [1:n].
– When I short-circuited the resistor connected to GND on the channel (1), its voltage dropped to zero as expected, but I observed a change of around 30 mV in the voltage measured on the channel (2).
– This behavior occurs consistently between any channel (n) and the subsequent channel (n+1).
Additional Tests
– I repeated the test using independent Analog Input blocks for each channel, and, in this case, the problem did not occur for any of the channels, demonstrating that the behavior is related to the vector configuration of the block.
– However, according to the SDRT user guide and the recommendation of MathWorks experts in the past, using a single Analog Input block for multiple signals, using vector notation for channel selection. This approach is necessary to maintain real-time operation, as using independent blocks introduces significant delays, even at low sampling frequencies.
Additional Information
– The same tests were carried out on two different computes (one slower and one faster and robut), with different cards of acquisition (PCI-6221 and PCIe-6323) and in two versions of MATLAB (2021b and 2024a).
– I also randomly changed the order of the input channels to rule out possible problems related to the sequence of the channels.
– The behavior persisted in all the configurations tested.
Requesting Assistance
Are there any known limitations or solutions to this behavior when using the Analog Input block in vector configuration? Considering that SDRT recommends this approach for multiple channels, would a technical review or guidance be possible to solve or mitigate the problem?
. #input interference, #simulink, #sdrt MATLAB Answers — New Questions

​

Hi All,
Forgive me if this topic of signal resetting in Simscape has been ask and solved countless times before. However, I am struggling to what seems a basic bit of logic in my Simscape model, and the solutons in the MATHWORKS examples I can find, seem to not reset the angular position signal by way of my trigger signal (I also have a rising edge signal).
Simply need to reset the angular position output from the Ideal Rotational Motion Sensor (converted to degrees) when my piston reaches its TDC position.
Using triggered subsystems to pass the signal through, or an integrator with a reset and manual swiches, are not permit the signal value to increase at each time-step once reset to 0.
Any pointer would be great on what should be a 5 minute task.
Thanks again.

Patrick

TDC Switch

Rising Edge SignalHi All,
Forgive me if this topic of signal resetting in Simscape has been ask and solved countless times before. However, I am struggling to what seems a basic bit of logic in my Simscape model, and the solutons in the MATHWORKS examples I can find, seem to not reset the angular position signal by way of my trigger signal (I also have a rising edge signal).
Simply need to reset the angular position output from the Ideal Rotational Motion Sensor (converted to degrees) when my piston reaches its TDC position.
Using triggered subsystems to pass the signal through, or an integrator with a reset and manual swiches, are not permit the signal value to increase at each time-step once reset to 0.
Any pointer would be great on what should be a 5 minute task.
Thanks again.

Patrick

TDC Switch

Rising Edge Signal Hi All,
Forgive me if this topic of signal resetting in Simscape has been ask and solved countless times before. However, I am struggling to what seems a basic bit of logic in my Simscape model, and the solutons in the MATHWORKS examples I can find, seem to not reset the angular position signal by way of my trigger signal (I also have a rising edge signal).
Simply need to reset the angular position output from the Ideal Rotational Motion Sensor (converted to degrees) when my piston reaches its TDC position.
Using triggered subsystems to pass the signal through, or an integrator with a reset and manual swiches, are not permit the signal value to increase at each time-step once reset to 0.
Any pointer would be great on what should be a 5 minute task.
Thanks again.

Patrick

TDC Switch

Rising Edge Signal simulink, simscape, stateflow MATLAB Answers — New Questions

​

Good evening,
I am writing code in Simulink for the TI f28p65 Launchpad of Texas instruments. Over all I use the c2000 library.
Following the tutorial https://de.mathworks.com/help/ti-c2000/ug/mat-file-sd-card-logging.html I am using a Scope block to log the data via the SPI-interface. Additionally I use SCIA-UART transmit to output the data in form of a string.
Sending the data via UART works fine in case I dont enable the sd card logging.
I have checked the correct setting of the pins as well as the funcitonality of the cables.

A similar issue occured already for the microcontroller, only it was programmed in c, such that no help was given for simulink. (https://e2e.ti.com/support/microcontrollers/c2000-microcontrollers-group/c2000/f/c2000-microcontrollers-forum/1080828/tms320f280025-uart-and-spi)

I am looking forward to hearing from you!

Cheers,
IsabellaGood evening,
I am writing code in Simulink for the TI f28p65 Launchpad of Texas instruments. Over all I use the c2000 library.
Following the tutorial https://de.mathworks.com/help/ti-c2000/ug/mat-file-sd-card-logging.html I am using a Scope block to log the data via the SPI-interface. Additionally I use SCIA-UART transmit to output the data in form of a string.
Sending the data via UART works fine in case I dont enable the sd card logging.
I have checked the correct setting of the pins as well as the funcitonality of the cables.

A similar issue occured already for the microcontroller, only it was programmed in c, such that no help was given for simulink. (https://e2e.ti.com/support/microcontrollers/c2000-microcontrollers-group/c2000/f/c2000-microcontrollers-forum/1080828/tms320f280025-uart-and-spi)

I am looking forward to hearing from you!

Cheers,
Isabella Good evening,
I am writing code in Simulink for the TI f28p65 Launchpad of Texas instruments. Over all I use the c2000 library.
Following the tutorial https://de.mathworks.com/help/ti-c2000/ug/mat-file-sd-card-logging.html I am using a Scope block to log the data via the SPI-interface. Additionally I use SCIA-UART transmit to output the data in form of a string.
Sending the data via UART works fine in case I dont enable the sd card logging.
I have checked the correct setting of the pins as well as the funcitonality of the cables.

A similar issue occured already for the microcontroller, only it was programmed in c, such that no help was given for simulink. (https://e2e.ti.com/support/microcontrollers/c2000-microcontrollers-group/c2000/f/c2000-microcontrollers-forum/1080828/tms320f280025-uart-and-spi)

I am looking forward to hearing from you!

Cheers,
Isabella c2000, simulink, sdcardlogging, sci MATLAB Answers — New Questions

​

I’ve tried running the "Infineon Buck Simscape Example" both online and locally and get the following error
View_system(‘BuckSystem’) unrecognized function or variable ‘view system"
I’m running R2023bI’ve tried running the "Infineon Buck Simscape Example" both online and locally and get the following error
View_system(‘BuckSystem’) unrecognized function or variable ‘view system"
I’m running R2023b I’ve tried running the "Infineon Buck Simscape Example" both online and locally and get the following error
View_system(‘BuckSystem’) unrecognized function or variable ‘view system"
I’m running R2023b infineon buck, view_system MATLAB Answers — New Questions

​

I am trying to control a DC to DC boost converter using FSC MPC to generate a 20 kW PV power but the MPC controller generated a maximum of 1135 W PV power and a maximum of 7184 W DC load power instead of 20 kW power. The PV voltage has a maximum of 59 V against 847 V boosted voltage with 21 A maximum PV current against a maximum of 27 A inductor current. Below is the MPC script used in the Matlab function block:
function S = MPC_Boost(iL, V_PV, Vout, iL_ref, wf_min, wf_max, Dref)
% MPC controller for boost converter with iL_ref from MPPT logic

%% System Parameters
Ts = 10e-6;
RL = 0.02;
L = 0.0199;
C = 3.1222e-04;
R_load = 18;

%% Limits
iL_max = 20;
iL_min = 0.2;
Vout_max = 420;
Vout_min = 360;
P_PV_max = 20000;
Vout_ref = 400;

%% Persistent memory
persistent S_old wf;
if isempty(S_old)
S_old = Dref;
wf = wf_min;
end

%% Sanitize Inputs
iL = abs(iL);
Vout = abs(Vout);
V_PV = max(V_PV, 50);

%% Duty Cycle Candidates
states = [Dref – 0.05; Dref; Dref + 0.05];
states = min(max(states, 0.3), 0.9); % Clamp

%% Cost Evaluation
g = zeros(length(states),1);
g_min = inf;
n_best = 1;

X_k = [iL; Vout];

for n = 1:length(states)
D = states(n);

% Discrete model
A_d = [1 – (RL * Ts / L), -(1 – D) * (Ts / L);
(1 – D) * (Ts / C), 1 – (Ts / (R_load * C))];
B_d = [Ts / L; 0];

% Predict state
X_k1 = A_d * X_k + B_d * V_PV;
iboost_k1 = X_k1(1);
Vout_k1 = X_k1(2);
P_PV_pred = V_PV * iboost_k1;

% Constraint check
if iboost_k1 < iL_min || iboost_k1 > iL_max || …
Vout_k1 < Vout_min || Vout_k1 > Vout_max || …
P_PV_pred <= 0 || P_PV_pred > P_PV_max || …
isnan(iboost_k1) || isnan(Vout_k1)
g(n) = inf;
continue;
end

% Cost terms
g_iboost = (iL_ref – iboost_k1)^2;
g_vout = (Vout_ref – Vout_k1)^2;
g_sw = (D ~= S_old) * wf;
g_dref = (D – Dref)^2;
g_pv = (P_PV_pred < 100) * 1000;

% Total cost
g(n) = g_iboost + 0.1 * g_vout + g_sw + g_dref + g_pv;

if g(n) < g_min
g_min = g(n);
n_best = n;
end
end

%% Apply Selected Duty Cycle
S = states(n_best);
S = min(max(S, 0.3), 0.9); % Clamp again

% Adaptive penalty update
if S ~= S_old
wf = min(wf * 1.10, wf_max);
else
wf = max(wf * 0.90, wf_min);
end

S_old = S;
endI am trying to control a DC to DC boost converter using FSC MPC to generate a 20 kW PV power but the MPC controller generated a maximum of 1135 W PV power and a maximum of 7184 W DC load power instead of 20 kW power. The PV voltage has a maximum of 59 V against 847 V boosted voltage with 21 A maximum PV current against a maximum of 27 A inductor current. Below is the MPC script used in the Matlab function block:
function S = MPC_Boost(iL, V_PV, Vout, iL_ref, wf_min, wf_max, Dref)
% MPC controller for boost converter with iL_ref from MPPT logic

%% System Parameters
Ts = 10e-6;
RL = 0.02;
L = 0.0199;
C = 3.1222e-04;
R_load = 18;

%% Limits
iL_max = 20;
iL_min = 0.2;
Vout_max = 420;
Vout_min = 360;
P_PV_max = 20000;
Vout_ref = 400;

%% Persistent memory
persistent S_old wf;
if isempty(S_old)
S_old = Dref;
wf = wf_min;
end

%% Sanitize Inputs
iL = abs(iL);
Vout = abs(Vout);
V_PV = max(V_PV, 50);

%% Duty Cycle Candidates
states = [Dref – 0.05; Dref; Dref + 0.05];
states = min(max(states, 0.3), 0.9); % Clamp

%% Cost Evaluation
g = zeros(length(states),1);
g_min = inf;
n_best = 1;

X_k = [iL; Vout];

for n = 1:length(states)
D = states(n);

% Discrete model
A_d = [1 – (RL * Ts / L), -(1 – D) * (Ts / L);
(1 – D) * (Ts / C), 1 – (Ts / (R_load * C))];
B_d = [Ts / L; 0];

% Predict state
X_k1 = A_d * X_k + B_d * V_PV;
iboost_k1 = X_k1(1);
Vout_k1 = X_k1(2);
P_PV_pred = V_PV * iboost_k1;

% Constraint check
if iboost_k1 < iL_min || iboost_k1 > iL_max || …
Vout_k1 < Vout_min || Vout_k1 > Vout_max || …
P_PV_pred <= 0 || P_PV_pred > P_PV_max || …
isnan(iboost_k1) || isnan(Vout_k1)
g(n) = inf;
continue;
end

% Cost terms
g_iboost = (iL_ref – iboost_k1)^2;
g_vout = (Vout_ref – Vout_k1)^2;
g_sw = (D ~= S_old) * wf;
g_dref = (D – Dref)^2;
g_pv = (P_PV_pred < 100) * 1000;

% Total cost
g(n) = g_iboost + 0.1 * g_vout + g_sw + g_dref + g_pv;

if g(n) < g_min
g_min = g(n);
n_best = n;
end
end

%% Apply Selected Duty Cycle
S = states(n_best);
S = min(max(S, 0.3), 0.9); % Clamp again

% Adaptive penalty update
if S ~= S_old
wf = min(wf * 1.10, wf_max);
else
wf = max(wf * 0.90, wf_min);
end

S_old = S;
end I am trying to control a DC to DC boost converter using FSC MPC to generate a 20 kW PV power but the MPC controller generated a maximum of 1135 W PV power and a maximum of 7184 W DC load power instead of 20 kW power. The PV voltage has a maximum of 59 V against 847 V boosted voltage with 21 A maximum PV current against a maximum of 27 A inductor current. Below is the MPC script used in the Matlab function block:
function S = MPC_Boost(iL, V_PV, Vout, iL_ref, wf_min, wf_max, Dref)
% MPC controller for boost converter with iL_ref from MPPT logic

%% System Parameters
Ts = 10e-6;
RL = 0.02;
L = 0.0199;
C = 3.1222e-04;
R_load = 18;

%% Limits
iL_max = 20;
iL_min = 0.2;
Vout_max = 420;
Vout_min = 360;
P_PV_max = 20000;
Vout_ref = 400;

%% Persistent memory
persistent S_old wf;
if isempty(S_old)
S_old = Dref;
wf = wf_min;
end

%% Sanitize Inputs
iL = abs(iL);
Vout = abs(Vout);
V_PV = max(V_PV, 50);

%% Duty Cycle Candidates
states = [Dref – 0.05; Dref; Dref + 0.05];
states = min(max(states, 0.3), 0.9); % Clamp

%% Cost Evaluation
g = zeros(length(states),1);
g_min = inf;
n_best = 1;

X_k = [iL; Vout];

for n = 1:length(states)
D = states(n);

% Discrete model
A_d = [1 – (RL * Ts / L), -(1 – D) * (Ts / L);
(1 – D) * (Ts / C), 1 – (Ts / (R_load * C))];
B_d = [Ts / L; 0];

% Predict state
X_k1 = A_d * X_k + B_d * V_PV;
iboost_k1 = X_k1(1);
Vout_k1 = X_k1(2);
P_PV_pred = V_PV * iboost_k1;

% Constraint check
if iboost_k1 < iL_min || iboost_k1 > iL_max || …
Vout_k1 < Vout_min || Vout_k1 > Vout_max || …
P_PV_pred <= 0 || P_PV_pred > P_PV_max || …
isnan(iboost_k1) || isnan(Vout_k1)
g(n) = inf;
continue;
end

% Cost terms
g_iboost = (iL_ref – iboost_k1)^2;
g_vout = (Vout_ref – Vout_k1)^2;
g_sw = (D ~= S_old) * wf;
g_dref = (D – Dref)^2;
g_pv = (P_PV_pred < 100) * 1000;

% Total cost
g(n) = g_iboost + 0.1 * g_vout + g_sw + g_dref + g_pv;

if g(n) < g_min
g_min = g(n);
n_best = n;
end
end

%% Apply Selected Duty Cycle
S = states(n_best);
S = min(max(S, 0.3), 0.9); % Clamp again

% Adaptive penalty update
if S ~= S_old
wf = min(wf * 1.10, wf_max);
else
wf = max(wf * 0.90, wf_min);
end

S_old = S;
end mpc control of solar pv MATLAB Answers — New Questions

​