MATLAB GUI code the Active Receptor should increased when at TstartLigand but im plot its decreased . Secondly how to get rid of this matlab ui.control even

Question

Ehtisham 2024-2-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2078781-matlab-gui-code-the-active-receptor-should-increased-when-at-tstartligand-but-im-plot-its-decreased

编辑： Cris LaPierre 2024-2-9

MATLAB GUI code the Active Receptor should increased when at TstartLigand but im plot its decreased . Secondly how to get rid of this matlab ui.control even(shared screenshot) involve in code.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Ehtisham 2024-2-7

移动：Cris LaPierre 2024-2-7

i did not edited the question i just added the comment beacuse i did a spelling mistake. As i said kinldy help me out to trace the fault in the Sim.m file script.

Cris LaPierre 2024-2-7

There was a screenshot and question about how to remove edit field labels from the app propertes. That is the only question answered so far, and it no longer appears in your question.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Cris LaPierre 2024-2-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2078781-matlab-gui-code-the-active-receptor-should-increased-when-at-tstartligand-but-im-plot-its-decreased#answer_1403836

The items you have circled are the Edit Field labels. An edit field component consists of an edit field and its label.

To remove the label, switch to the App Canvas, select the component, then position the cursor over the label and select it. You can then delete it using the delete or backspace key. Here I have deleted KBONEditFieldLabel.

As for the plotting, you'll need to tell us how to use your app in order to see what is going on.

6 个评论
显示 4更早的评论隐藏 4更早的评论

Cris LaPierre 2024-2-7

在 MATLAB Online 中打开

Sim.m

Thanks, that is helpful. The script you shared creates the exact same plot, with the same issue. Before digging into App Designer, I'd focus on getting the code correct in the script first.

% Define inputs

Kf1Max = 100;

L_Active = 0.5;

L_min = 0.05;

TauKFON = -1;

TauKFOFF = -1;

Kb1Max = 80;

Kb1Min = 10;

TauKBON = -0.01;

TauKBOFF = -0.01;

T_startLigand = 100;

T_endLigand = 1000;

timespan = [0 3000];

Non_Active_Receptor_concentration = 100;

Active_Receptor_concentration = 0;

% Call Sim function

[t, Non_Active_Receptor_concentration, Active_Receptor_concentration, peak_value, time_to_peak, ...

steady_state, time_to_50_percent, peak_to_steady_state_ratio, Delta, KF_LMaxA, KF_LMaxB, kf_L, kb_1] ...

= Sim (Kf1Max, L_Active, L_min, TauKFON, TauKFOFF, Kb1Max, Kb1Min, TauKBON, TauKBOFF, ...

T_startLigand, T_endLigand, timespan, Non_Active_Receptor_concentration, Active_Receptor_concentration);

% Create plot

plot(t, Non_Active_Receptor_concentration , 'r');

hold on

plot(t, Active_Receptor_concentration, 'b--');

hold off

legend('Non\_Active Receptor' , 'Active\_Receptor');

title('Activation of Receptor'); % Set a title for the first plot

Cris LaPierre 2024-2-9

在 MATLAB Online 中打开

After some reverse engineering, I believe the root cause is the difference between TauKFON and TauKFOFF (both -1) and TauKBON and TauKBOFF (both -0.01). The result is that kf_L is a step function while kb_1 is damped both turning on and off. Compare the blue and red curves below.

t = 0:3000;

T_startLigand = 100;

T_endLigand = 1000;

%% kf_L

Kf1Max = 100;

L_Active = 0.5;

L_min = 0.05;

TauKFON = -1;

TauKFOFF = -61;

KF_LMaxA = Kf1Max* (L_min / (L_min + 1));

KF_LMaxB = Kf1Max* (L_Active / (L_Active + 1));

kf_L = zeros(1,length(t));

kf_L(t < T_startLigand) = KF_LMaxA;

% Increasing phase: Exponential growth from KF_LMaxA to KF_LMaxB

kf_L(t >= T_startLigand & t < T_endLigand) = KF_LMaxB - (KF_LMaxB - KF_LMaxA) * exp(TauKFON * (0:((T_endLigand-1)-T_startLigand)));

% Decreasing phase: Exponential decay from KF_LMaxA to KF_LMaxB

kf_end = KF_LMaxB - (KF_LMaxB - KF_LMaxA) * exp(TauKFON * (T_endLigand - T_startLigand));

kf_L(t >= T_endLigand) = KF_LMaxA + (kf_end - KF_LMaxA) * exp(TauKFOFF * ((T_endLigand:t(end))-T_endLigand));

%% kb_1

Kb1Max = 80;

Kb1Min = 10;

TauKBON = -0.01;

TauKBOFF = -0.01;

kb_1 = zeros(1,length(t));

kb_1(t < T_startLigand) = Kb1Min;

% Increasing phase: Exponential growth from Kb1Min to Kb1Max

kb_1(t >= T_startLigand & t < T_endLigand) = Kb1Max - (Kb1Max - Kb1Min) * exp(TauKBON * (0:((T_endLigand-1)-T_startLigand)));

% Decreasing phase: Exponential decay from Kb1Max to Kb1Min

kb_end = Kb1Max - (Kb1Max - Kb1Min) * exp(TauKBON * (T_endLigand - T_startLigand));

kb_1(t >= T_endLigand) = Kb1Min + (kb_end - Kb1Min) * exp(TauKBOFF * ((T_endLigand:t(end))-T_endLigand));

yyaxis left

plot(t,kf_L)

yyaxis right

plot(t,kb_1)

legend("kf_L","kb_1",'interpreter','none')

title('Reaction Rates')

That difference in response time means that the forward rate is going 100% as soon as it is turned on while the backward rate takes some time to get there. Based on how you defined the reactions in ode_LR, kf_L is the rate non-active receptors turn into active receptors, and kb_1 is the reverse rate.

Non_Active_Receptor_concentration = 100;

Active_Receptor_concentration = 0;

y0 = [Non_Active_Receptor_concentration, Active_Receptor_concentration];

[t, y] = ode45(@(t,y) ode_LR(t,y,kf_L,kb_1), t, y0);

figure

plot(t,y)

legend('[Non-Active Receptors]','[Active Receptors]')

title('Rerceptor Concentration')

% I've modififed your code for demo purposes only. This is not 'better'

function dydt = ode_LR(t,y,kf_L,kb_1)

t0 = 0:length(kf_L)-1;

kf = interp1(t0,kf_L,t);

kb = interp1(t0,kb_1,t);

dydt = [-kf*y(1) + kb * y(2);

kf * y(1) - kb * y(2)];

end

The sharp change at t=100 is because you have receptors becoming active while none are becoming non-active. This causes a rapid depletion of your non-active concentration, and a rapid increase in your active concentration. After some time, kb_1 increases to a point where the two are more or less in equilibrium until they turn off.

I can't say if this is correct or not. The code is performing as written, and makes intuitive sense to me. I'd start by making sure the rates and other constants are correct, paying close attention to units.

Ehtisham 2024-2-9

编辑：Ehtisham 2024-2-9

I want that this sharp change should maintain till time TendLigand beacuse the kb_1 is slower reaction than the Kf_L .

Cris LaPierre 2024-2-9

编辑：Cris LaPierre 2024-2-9

There is no MATLAB problem in your code. I suspect you have either not implemented your model correctly, or you are using the wrong values. However, I have no way of knowing what to change without more information from you.

You can guess-and-check different values if that is all that remains to be done.

请先登录，再进行评论。

MATLAB GUI code the Active Receptor should increased when at TstartLigand but im plot its decreased . Secondly how to get rid of this matlab ui.control even

3 个评论
显示 1更早的评论隐藏 1更早的评论

采纳的回答

6 个评论
显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

MATLAB GUI code the Active Receptor should increased when at TstartLigand but im plot its decreased . Secondly how to get rid of this matlab ui.control even

3 个评论 显示 1更早的评论隐藏 1更早的评论

采纳的回答

6 个评论 显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

3 个评论
显示 1更早的评论隐藏 1更早的评论

6 个评论
显示 4更早的评论隐藏 4更早的评论