In a first set of measurements from my VNA (vector network analyser), I was able to collect .s2p files which I could plot the S-parameters from. In a later set of measurements, I somehow only had the option to save the files as either .s1p. Can I not gain the S21 and s12 plots from these? Matlab won’t let me import the data either, is this something that Matlab is not designed to do or am I doing something wrong?In a first set of measurements from my VNA (vector network analyser), I was able to collect .s2p files which I could plot the S-parameters from. In a later set of measurements, I somehow only had the option to save the files as either .s1p. Can I not gain the S21 and s12 plots from these? Matlab won’t let me import the data either, is this something that Matlab is not designed to do or am I doing something wrong? In a first set of measurements from my VNA (vector network analyser), I was able to collect .s2p files which I could plot the S-parameters from. In a later set of measurements, I somehow only had the option to save the files as either .s1p. Can I not gain the S21 and s12 plots from these? Matlab won’t let me import the data either, is this something that Matlab is not designed to do or am I doing something wrong? s-parameters, scattering parameters, vna data, vector network analyser data, .s1p MATLAB Answers — New Questions

​

Good morning,

I’m using MATLAB R2020a Update 2. I have a tiledlayout of five (three by two) area plots and would like to use the sixth, currently empty, tile to add a global legend. I’ve already found this question, and understand that there’s no official, built-in way of doing this, but perhaps it’s possible to get creative.

What I’ve tried is adding a new, empty area plot to the sixth tile, using NaNs as the shares, then adding a legend to that and setting the axis object for that tile to invisible, as follows:

% …
ax = nexttile;
area([NaN NaN], NaN(2, 4));
leg = legend({‘tar’, ‘sta’, ‘ext’, ‘dmp’);
leg.Location = ‘none’;
leg.Interpreter = ‘latex’;
leg.FontSize = 16;
ax.Visible = false;

This works in principle, but leaves a horizontal line where the X axis would be:

This only happens when I use area(), not e.g. plot(). Perhaps it’s a bug, perhaps it’s intentional and due to the way area() works. In any case it’s not what I want, but I can’t get rid of it (without breaking other things in the process).

I had the idea of making the children of this axis, i.e. the areas of the area plot, invisible as well, like so:

for i = 1:length(ax.Children)
ax.Children(i).Visible = false;
end

But while this gets rid of the horizontal line it also grays out the legend entries:

This behavior in turn is known and expected and apparently cannot be changed. The workaround suggested by a Mathworks staffer in the linked question is to plot NaNs; but that’s what I’m already doing and what’s leaving the horizontal line, due to area()’s quirks.

Can anyone help?

I’m not hung up on specifically creating an invisible area(), this is merely the best (first, only) idea I had for fudging a global legend. If anyone can make this approach work, that’d be wonderful. If anyone has another idea of how to achieve a similar effect, that’d be wonderful as well.

Thank you very much!Good morning,

I’m using MATLAB R2020a Update 2. I have a tiledlayout of five (three by two) area plots and would like to use the sixth, currently empty, tile to add a global legend. I’ve already found this question, and understand that there’s no official, built-in way of doing this, but perhaps it’s possible to get creative.

What I’ve tried is adding a new, empty area plot to the sixth tile, using NaNs as the shares, then adding a legend to that and setting the axis object for that tile to invisible, as follows:

% …
ax = nexttile;
area([NaN NaN], NaN(2, 4));
leg = legend({‘tar’, ‘sta’, ‘ext’, ‘dmp’);
leg.Location = ‘none’;
leg.Interpreter = ‘latex’;
leg.FontSize = 16;
ax.Visible = false;

This works in principle, but leaves a horizontal line where the X axis would be:

This only happens when I use area(), not e.g. plot(). Perhaps it’s a bug, perhaps it’s intentional and due to the way area() works. In any case it’s not what I want, but I can’t get rid of it (without breaking other things in the process).

I had the idea of making the children of this axis, i.e. the areas of the area plot, invisible as well, like so:

for i = 1:length(ax.Children)
ax.Children(i).Visible = false;
end

But while this gets rid of the horizontal line it also grays out the legend entries:

This behavior in turn is known and expected and apparently cannot be changed. The workaround suggested by a Mathworks staffer in the linked question is to plot NaNs; but that’s what I’m already doing and what’s leaving the horizontal line, due to area()’s quirks.

Can anyone help?

I’m not hung up on specifically creating an invisible area(), this is merely the best (first, only) idea I had for fudging a global legend. If anyone can make this approach work, that’d be wonderful. If anyone has another idea of how to achieve a similar effect, that’d be wonderful as well.

Thank you very much! Good morning,

I’m using MATLAB R2020a Update 2. I have a tiledlayout of five (three by two) area plots and would like to use the sixth, currently empty, tile to add a global legend. I’ve already found this question, and understand that there’s no official, built-in way of doing this, but perhaps it’s possible to get creative.

What I’ve tried is adding a new, empty area plot to the sixth tile, using NaNs as the shares, then adding a legend to that and setting the axis object for that tile to invisible, as follows:

% …
ax = nexttile;
area([NaN NaN], NaN(2, 4));
leg = legend({‘tar’, ‘sta’, ‘ext’, ‘dmp’);
leg.Location = ‘none’;
leg.Interpreter = ‘latex’;
leg.FontSize = 16;
ax.Visible = false;

This works in principle, but leaves a horizontal line where the X axis would be:

This only happens when I use area(), not e.g. plot(). Perhaps it’s a bug, perhaps it’s intentional and due to the way area() works. In any case it’s not what I want, but I can’t get rid of it (without breaking other things in the process).

I had the idea of making the children of this axis, i.e. the areas of the area plot, invisible as well, like so:

for i = 1:length(ax.Children)
ax.Children(i).Visible = false;
end

But while this gets rid of the horizontal line it also grays out the legend entries:

This behavior in turn is known and expected and apparently cannot be changed. The workaround suggested by a Mathworks staffer in the linked question is to plot NaNs; but that’s what I’m already doing and what’s leaving the horizontal line, due to area()’s quirks.

Can anyone help?

I’m not hung up on specifically creating an invisible area(), this is merely the best (first, only) idea I had for fudging a global legend. If anyone can make this approach work, that’d be wonderful. If anyone has another idea of how to achieve a similar effect, that’d be wonderful as well.

Thank you very much! tiledlayout, legend MATLAB Answers — New Questions

​

I am using File Log blocks in my Simulink Real-Time (SLRT) model (available since R2020b) to log data on the Speedgoat target computer file system. After the real-time application stops, if the development computer running MATLAB is connected to the Speedgoat, File Log data is automatically transferred and loaded on the Simulation Data Inspector (SDI). 
What are the different ways to turn off this auto-import behavior for File Logs?I am using File Log blocks in my Simulink Real-Time (SLRT) model (available since R2020b) to log data on the Speedgoat target computer file system. After the real-time application stops, if the development computer running MATLAB is connected to the Speedgoat, File Log data is automatically transferred and loaded on the Simulation Data Inspector (SDI). 
What are the different ways to turn off this auto-import behavior for File Logs? I am using File Log blocks in my Simulink Real-Time (SLRT) model (available since R2020b) to log data on the Speedgoat target computer file system. After the real-time application stops, if the development computer running MATLAB is connected to the Speedgoat, File Log data is automatically transferred and loaded on the Simulation Data Inspector (SDI). 
What are the different ways to turn off this auto-import behavior for File Logs? autoimportfilelog, filelog MATLAB Answers — New Questions

​

% [P] * [TPM]^n* [R]= CR
disp(‘Create an array with seven elements in a single row:’)
disp(‘>> P = [1 0 0 0 0 0 0]’)
p = [1 0 0 0 0 0 0]

disp(‘Create an array with seven elements in a single column:’)
disp(‘>> R = [9; 8; 7; 6; 5; 4; 3]’)
R = [9; 8; 7; 6; 5; 4; 3]

CR = [0.45]

Solve (TPM)^20=CR/ (p*R)

I would like to find the matrics [TPM] with size 7*7? How to do it ?% [P] * [TPM]^n* [R]= CR
disp(‘Create an array with seven elements in a single row:’)
disp(‘>> P = [1 0 0 0 0 0 0]’)
p = [1 0 0 0 0 0 0]

disp(‘Create an array with seven elements in a single column:’)
disp(‘>> R = [9; 8; 7; 6; 5; 4; 3]’)
R = [9; 8; 7; 6; 5; 4; 3]

CR = [0.45]

Solve (TPM)^20=CR/ (p*R)

I would like to find the matrics [TPM] with size 7*7? How to do it ? % [P] * [TPM]^n* [R]= CR
disp(‘Create an array with seven elements in a single row:’)
disp(‘>> P = [1 0 0 0 0 0 0]’)
p = [1 0 0 0 0 0 0]

disp(‘Create an array with seven elements in a single column:’)
disp(‘>> R = [9; 8; 7; 6; 5; 4; 3]’)
R = [9; 8; 7; 6; 5; 4; 3]

CR = [0.45]

Solve (TPM)^20=CR/ (p*R)

I would like to find the matrics [TPM] with size 7*7? How to do it ? matrix, solve, equation, matrix manipulation, matrix array, matrices, matlab, mathematics MATLAB Answers — New Questions

​

Hi, I try to call the eigenvectors of the respective eigenvalue, but that part in the code gives an error:
"Unrecognize function or variable ‘eigvecs_theta’

I tried eig and eigvec instead of eigvecs_theta, but nothing works. What did I miss here? Thanks

% Discretization parameters
theta_min = 0;
theta_max = pi;
M = 100; % Number of grid points
theta = linspace(theta_min, theta_max, M);
dtheta = theta(2) – theta(1);

% Initialize the Theta(theta) vector
Theta = zeros(M, 1);

% Set up the finite difference matrix
B = zeros(M, M);
for j = 3:M-2
B(j, j-2) = -1 / (2 * dtheta^3);
B(j, j-1) = 2 / (dtheta^3);
B(j, j) = -2 / (dtheta^3) + l * (l + 1);
B(j, j+1) = 2 / (dtheta^3);
B(j, j+2) = -1 / (2 * dtheta^3);
end

% Apply boundary conditions (example: Theta(0) = 0 and Theta(pi) = 0)
B(1,1) = 1;
B(M,M) = 1;

% Solve the eigenvalue problem
[~, D_theta] = eig(B);

% The eigenvalues are the diagonal elements of D_theta
eigenvalues_theta = diag(D_theta);

% The solution Theta(theta) corresponds to the eigenvector with the desired eigenvalue
Theta = eig(:, idx);

% Plot the solution
plot(theta, Theta);
xlabel(‘theta’);
ylabel(‘Theta(theta)’);
title(‘Angular Wavefunction’);
endHi, I try to call the eigenvectors of the respective eigenvalue, but that part in the code gives an error:
"Unrecognize function or variable ‘eigvecs_theta’

I tried eig and eigvec instead of eigvecs_theta, but nothing works. What did I miss here? Thanks

% Discretization parameters
theta_min = 0;
theta_max = pi;
M = 100; % Number of grid points
theta = linspace(theta_min, theta_max, M);
dtheta = theta(2) – theta(1);

% Initialize the Theta(theta) vector
Theta = zeros(M, 1);

% Set up the finite difference matrix
B = zeros(M, M);
for j = 3:M-2
B(j, j-2) = -1 / (2 * dtheta^3);
B(j, j-1) = 2 / (dtheta^3);
B(j, j) = -2 / (dtheta^3) + l * (l + 1);
B(j, j+1) = 2 / (dtheta^3);
B(j, j+2) = -1 / (2 * dtheta^3);
end

% Apply boundary conditions (example: Theta(0) = 0 and Theta(pi) = 0)
B(1,1) = 1;
B(M,M) = 1;

% Solve the eigenvalue problem
[~, D_theta] = eig(B);

% The eigenvalues are the diagonal elements of D_theta
eigenvalues_theta = diag(D_theta);

% The solution Theta(theta) corresponds to the eigenvector with the desired eigenvalue
Theta = eig(:, idx);

% Plot the solution
plot(theta, Theta);
xlabel(‘theta’);
ylabel(‘Theta(theta)’);
title(‘Angular Wavefunction’);
end Hi, I try to call the eigenvectors of the respective eigenvalue, but that part in the code gives an error:
"Unrecognize function or variable ‘eigvecs_theta’

I tried eig and eigvec instead of eigvecs_theta, but nothing works. What did I miss here? Thanks

% Discretization parameters
theta_min = 0;
theta_max = pi;
M = 100; % Number of grid points
theta = linspace(theta_min, theta_max, M);
dtheta = theta(2) – theta(1);

% Initialize the Theta(theta) vector
Theta = zeros(M, 1);

% Set up the finite difference matrix
B = zeros(M, M);
for j = 3:M-2
B(j, j-2) = -1 / (2 * dtheta^3);
B(j, j-1) = 2 / (dtheta^3);
B(j, j) = -2 / (dtheta^3) + l * (l + 1);
B(j, j+1) = 2 / (dtheta^3);
B(j, j+2) = -1 / (2 * dtheta^3);
end

% Apply boundary conditions (example: Theta(0) = 0 and Theta(pi) = 0)
B(1,1) = 1;
B(M,M) = 1;

% Solve the eigenvalue problem
[~, D_theta] = eig(B);

% The eigenvalues are the diagonal elements of D_theta
eigenvalues_theta = diag(D_theta);

% The solution Theta(theta) corresponds to the eigenvector with the desired eigenvalue
Theta = eig(:, idx);

% Plot the solution
plot(theta, Theta);
xlabel(‘theta’);
ylabel(‘Theta(theta)’);
title(‘Angular Wavefunction’);
end eigenvalues, eigenvectors MATLAB Answers — New Questions

​

I get this error trying to solve a system of PDEs, and I do not know if such system is solvable with ‘pdepe’. The equations are:

energy balance equation, and:

mass balance equations. Initial conditions:

Boundary conditions:

The variables in which I am solving are however T,, knowing that:
, I can write initial condition for , and I have also a value for .
I tried writing the equations in the form required by ‘pdepe’. Also, note that many parameters depends on the vector , IN particular, the term ‘v’ includes a partial spatial derivative w.r.t. P, which i transformed in a spatial derivative w.r.t. u(1) and u(2) with the above formula. Note that many constants are imported with the file.m including the initial conditions that are: This is the code of the 3 functions used inside the pde call (sorry is a bit long and messy):
function [c,f,s] = pdefun(x,t,u,dudx)
% load data
run(‘HM_data.m’);

% compute epsilon , F and dfdt
F = (hm.rho_s_in-u(3))/(u(3)*hm.tau_p – hm.rho_s_in*hm.w_max);
epsilon = 1 – (1 – hm.epsilon_0)*((1 + hm.tau_p*F)/(1 + hm.tau_a*F));
P = u(2)*u(1)*hm.R_gas/hm.MH2;
P_eq_a = (hm.C0_a*hm.C1_a*(F)^(hm.C2_a)/(1+hm.C1_a*(F)^(hm.C2_a)) + hm.C3_a*F + exp(hm.C4_a*(F-hm.C5_a)))*exp(-hm.K_a*((1/u(1))-(1/303)));
k_a = (hm.kappa_a/(1 – epsilon)) * exp(-hm.E_a/(hm.R_gas*u(1))) * log(P/P_eq_a);
dfdt = k_a*(1-F);

% compute diagonal matrix c
Cps = 6000*(3.1*hm.R + 10.04*hm.x_max*F)/(hm.Ms + 6*hm.x_max*F);
rho_C_eff = epsilon*u(2)*hm.Cpg + (1- epsilon)*u(3)*Cps;
c = [rho_C_eff; epsilon; 1-epsilon];

% COMPUTE VECTOR f
% velocity
Dp = hm.D_in*(1+ F*hm.tau_p)^(1/3);
Kp = (Dp^2)*(epsilon^3)/(150*(1-epsilon)^2);
v_cost = -Kp*hm.R_gas/(hm.mu_g*hm.MH2);
v = v_cost*(dudx(1)*u(2) + dudx(2)*u(1));

% effective thermal conductivity
N = 3.08/epsilon – 1.13;
Fn = 4/3 * hm.E_prime * hm.R^(1/2) * hm.dv^(1.5);
Hv = hm.c1*hm.dv^(hm.c2);
Rs = 0.565*Hv*hm.dv/(hm.ks*Fn);
a_H = (0.75*Fn*hm.R/hm.E_prime)^(1/3);
a_LH = [];
if epsilon <= 0.47 && epsilon >= 0.01
a_LH = 1.605/sqrt(epsilon);
elseif epsilon > 0.47 && epsilon <= 1
a_LH = 3.51 – 2.51*epsilon;
else
error(‘error epsilon’);
end
a_L = a_LH*a_H;
R_L = 1/(2*hm.ks*a_L);

Pmax = 2/pi*hm.E_prime*(hm.dv/hm.R)^(0.5);
Hc = hm.c1*(1.62*hm.dv)^(hm.c2);
a1 = erfcinv(2*Pmax/Hc);
a2 = erfcinv(0.03*Pmax/Hc) – a1;
coef_a = 2*hm.b/Dp;
coef_c = -(6*(hm.gamma-1)/(9*hm.gamma-5))* (hm.kg_ref*hm.MH2/(u(1)*u(2)*hm.R_gas))*(hm.MH2*u(1)/(2*hm.kb))^(0.5);

l_m = (-1 + sqrt(1 – 4*coef_a*coef_c))/(2*coef_a);
kg = hm.kg_ref/(1+2*hm.b*l_m/Dp);
M = (((2-hm.alpha_T1)/hm.alpha_T1)+((2-hm.alpha_T2)/hm.alpha_T2))*(2*hm.gamma/(1+hm.gamma))*(1/hm.Pr)*l_m;
R_g = sqrt(2)*hm.sigma*a2/(pi * kg * a_L^2 * log(1+ a2/(a1+M/(sqrt(2)*hm.sigma))));
L = (hm.gamma+1)*3*Dp/((9*hm.gamma-5)*4*l_m*sqrt(pi));
R_G = 1/(2*pi*kg*Dp*(0.5*log(1+L) + log(1+sqrt(L)) + 1/(1+sqrt(L)) – 1));
R_mic = Rs*R_g/(Rs+R_g);
R_c_inv = 1/(R_mic+R_L) + 1/R_G;

k_eff = N*(1-epsilon)*R_c_inv/(pi*Dp);

f = [k_eff*dudx(1); -u(2)*v; 0];

% compute vector s
m_dot_a = (1-epsilon)*(hm.rho_sat-u(3))*dfdt;
s = [u(2)*hm.Cpg*v*dudx(1) + m_dot_a*hm.delta_H ; -m_dot_a+hm.phi_abs; m_dot_a];

end

function u0 = icfun(x)
run(‘HM_data.m’);
u0 = hm.ic;
end

function [pl,ql,pr,qr] = bcfun(xl, ul, xr, ur, t)
run(‘HM_data.m’);
Nu_abs = 0.3+((0.62*(hm.Re_a^(0.5))*(hm.Pr_a^(1/3)))/(1+(0.4/hm.Pr_a)^(2/3))^(1/4)*((1+(hm.Re_a/282000)^(5/8))^(4/5)));
h_f = Nu_abs*hm.kf/hm.D_tank;
pl = [0 ; 0; 0];
ql = [1; 0; 0];
pr = [h_f*(ur(1)-hm.Tf_a); 0; 0];
qr = [-1; 0; 0];
endI get this error trying to solve a system of PDEs, and I do not know if such system is solvable with ‘pdepe’. The equations are:

energy balance equation, and:

mass balance equations. Initial conditions:

Boundary conditions:

The variables in which I am solving are however T,, knowing that:
, I can write initial condition for , and I have also a value for .
I tried writing the equations in the form required by ‘pdepe’. Also, note that many parameters depends on the vector , IN particular, the term ‘v’ includes a partial spatial derivative w.r.t. P, which i transformed in a spatial derivative w.r.t. u(1) and u(2) with the above formula. Note that many constants are imported with the file.m including the initial conditions that are: This is the code of the 3 functions used inside the pde call (sorry is a bit long and messy):
function [c,f,s] = pdefun(x,t,u,dudx)
% load data
run(‘HM_data.m’);

% compute epsilon , F and dfdt
F = (hm.rho_s_in-u(3))/(u(3)*hm.tau_p – hm.rho_s_in*hm.w_max);
epsilon = 1 – (1 – hm.epsilon_0)*((1 + hm.tau_p*F)/(1 + hm.tau_a*F));
P = u(2)*u(1)*hm.R_gas/hm.MH2;
P_eq_a = (hm.C0_a*hm.C1_a*(F)^(hm.C2_a)/(1+hm.C1_a*(F)^(hm.C2_a)) + hm.C3_a*F + exp(hm.C4_a*(F-hm.C5_a)))*exp(-hm.K_a*((1/u(1))-(1/303)));
k_a = (hm.kappa_a/(1 – epsilon)) * exp(-hm.E_a/(hm.R_gas*u(1))) * log(P/P_eq_a);
dfdt = k_a*(1-F);

% compute diagonal matrix c
Cps = 6000*(3.1*hm.R + 10.04*hm.x_max*F)/(hm.Ms + 6*hm.x_max*F);
rho_C_eff = epsilon*u(2)*hm.Cpg + (1- epsilon)*u(3)*Cps;
c = [rho_C_eff; epsilon; 1-epsilon];

% COMPUTE VECTOR f
% velocity
Dp = hm.D_in*(1+ F*hm.tau_p)^(1/3);
Kp = (Dp^2)*(epsilon^3)/(150*(1-epsilon)^2);
v_cost = -Kp*hm.R_gas/(hm.mu_g*hm.MH2);
v = v_cost*(dudx(1)*u(2) + dudx(2)*u(1));

% effective thermal conductivity
N = 3.08/epsilon – 1.13;
Fn = 4/3 * hm.E_prime * hm.R^(1/2) * hm.dv^(1.5);
Hv = hm.c1*hm.dv^(hm.c2);
Rs = 0.565*Hv*hm.dv/(hm.ks*Fn);
a_H = (0.75*Fn*hm.R/hm.E_prime)^(1/3);
a_LH = [];
if epsilon <= 0.47 && epsilon >= 0.01
a_LH = 1.605/sqrt(epsilon);
elseif epsilon > 0.47 && epsilon <= 1
a_LH = 3.51 – 2.51*epsilon;
else
error(‘error epsilon’);
end
a_L = a_LH*a_H;
R_L = 1/(2*hm.ks*a_L);

Pmax = 2/pi*hm.E_prime*(hm.dv/hm.R)^(0.5);
Hc = hm.c1*(1.62*hm.dv)^(hm.c2);
a1 = erfcinv(2*Pmax/Hc);
a2 = erfcinv(0.03*Pmax/Hc) – a1;
coef_a = 2*hm.b/Dp;
coef_c = -(6*(hm.gamma-1)/(9*hm.gamma-5))* (hm.kg_ref*hm.MH2/(u(1)*u(2)*hm.R_gas))*(hm.MH2*u(1)/(2*hm.kb))^(0.5);

l_m = (-1 + sqrt(1 – 4*coef_a*coef_c))/(2*coef_a);
kg = hm.kg_ref/(1+2*hm.b*l_m/Dp);
M = (((2-hm.alpha_T1)/hm.alpha_T1)+((2-hm.alpha_T2)/hm.alpha_T2))*(2*hm.gamma/(1+hm.gamma))*(1/hm.Pr)*l_m;
R_g = sqrt(2)*hm.sigma*a2/(pi * kg * a_L^2 * log(1+ a2/(a1+M/(sqrt(2)*hm.sigma))));
L = (hm.gamma+1)*3*Dp/((9*hm.gamma-5)*4*l_m*sqrt(pi));
R_G = 1/(2*pi*kg*Dp*(0.5*log(1+L) + log(1+sqrt(L)) + 1/(1+sqrt(L)) – 1));
R_mic = Rs*R_g/(Rs+R_g);
R_c_inv = 1/(R_mic+R_L) + 1/R_G;

k_eff = N*(1-epsilon)*R_c_inv/(pi*Dp);

f = [k_eff*dudx(1); -u(2)*v; 0];

% compute vector s
m_dot_a = (1-epsilon)*(hm.rho_sat-u(3))*dfdt;
s = [u(2)*hm.Cpg*v*dudx(1) + m_dot_a*hm.delta_H ; -m_dot_a+hm.phi_abs; m_dot_a];

end

function u0 = icfun(x)
run(‘HM_data.m’);
u0 = hm.ic;
end

function [pl,ql,pr,qr] = bcfun(xl, ul, xr, ur, t)
run(‘HM_data.m’);
Nu_abs = 0.3+((0.62*(hm.Re_a^(0.5))*(hm.Pr_a^(1/3)))/(1+(0.4/hm.Pr_a)^(2/3))^(1/4)*((1+(hm.Re_a/282000)^(5/8))^(4/5)));
h_f = Nu_abs*hm.kf/hm.D_tank;
pl = [0 ; 0; 0];
ql = [1; 0; 0];
pr = [h_f*(ur(1)-hm.Tf_a); 0; 0];
qr = [-1; 0; 0];
end I get this error trying to solve a system of PDEs, and I do not know if such system is solvable with ‘pdepe’. The equations are:

energy balance equation, and:

mass balance equations. Initial conditions:

Boundary conditions:

The variables in which I am solving are however T,, knowing that:
, I can write initial condition for , and I have also a value for .
I tried writing the equations in the form required by ‘pdepe’. Also, note that many parameters depends on the vector , IN particular, the term ‘v’ includes a partial spatial derivative w.r.t. P, which i transformed in a spatial derivative w.r.t. u(1) and u(2) with the above formula. Note that many constants are imported with the file.m including the initial conditions that are: This is the code of the 3 functions used inside the pde call (sorry is a bit long and messy):
function [c,f,s] = pdefun(x,t,u,dudx)
% load data
run(‘HM_data.m’);

% compute epsilon , F and dfdt
F = (hm.rho_s_in-u(3))/(u(3)*hm.tau_p – hm.rho_s_in*hm.w_max);
epsilon = 1 – (1 – hm.epsilon_0)*((1 + hm.tau_p*F)/(1 + hm.tau_a*F));
P = u(2)*u(1)*hm.R_gas/hm.MH2;
P_eq_a = (hm.C0_a*hm.C1_a*(F)^(hm.C2_a)/(1+hm.C1_a*(F)^(hm.C2_a)) + hm.C3_a*F + exp(hm.C4_a*(F-hm.C5_a)))*exp(-hm.K_a*((1/u(1))-(1/303)));
k_a = (hm.kappa_a/(1 – epsilon)) * exp(-hm.E_a/(hm.R_gas*u(1))) * log(P/P_eq_a);
dfdt = k_a*(1-F);

% compute diagonal matrix c
Cps = 6000*(3.1*hm.R + 10.04*hm.x_max*F)/(hm.Ms + 6*hm.x_max*F);
rho_C_eff = epsilon*u(2)*hm.Cpg + (1- epsilon)*u(3)*Cps;
c = [rho_C_eff; epsilon; 1-epsilon];

% COMPUTE VECTOR f
% velocity
Dp = hm.D_in*(1+ F*hm.tau_p)^(1/3);
Kp = (Dp^2)*(epsilon^3)/(150*(1-epsilon)^2);
v_cost = -Kp*hm.R_gas/(hm.mu_g*hm.MH2);
v = v_cost*(dudx(1)*u(2) + dudx(2)*u(1));

% effective thermal conductivity
N = 3.08/epsilon – 1.13;
Fn = 4/3 * hm.E_prime * hm.R^(1/2) * hm.dv^(1.5);
Hv = hm.c1*hm.dv^(hm.c2);
Rs = 0.565*Hv*hm.dv/(hm.ks*Fn);
a_H = (0.75*Fn*hm.R/hm.E_prime)^(1/3);
a_LH = [];
if epsilon <= 0.47 && epsilon >= 0.01
a_LH = 1.605/sqrt(epsilon);
elseif epsilon > 0.47 && epsilon <= 1
a_LH = 3.51 – 2.51*epsilon;
else
error(‘error epsilon’);
end
a_L = a_LH*a_H;
R_L = 1/(2*hm.ks*a_L);

Pmax = 2/pi*hm.E_prime*(hm.dv/hm.R)^(0.5);
Hc = hm.c1*(1.62*hm.dv)^(hm.c2);
a1 = erfcinv(2*Pmax/Hc);
a2 = erfcinv(0.03*Pmax/Hc) – a1;
coef_a = 2*hm.b/Dp;
coef_c = -(6*(hm.gamma-1)/(9*hm.gamma-5))* (hm.kg_ref*hm.MH2/(u(1)*u(2)*hm.R_gas))*(hm.MH2*u(1)/(2*hm.kb))^(0.5);

l_m = (-1 + sqrt(1 – 4*coef_a*coef_c))/(2*coef_a);
kg = hm.kg_ref/(1+2*hm.b*l_m/Dp);
M = (((2-hm.alpha_T1)/hm.alpha_T1)+((2-hm.alpha_T2)/hm.alpha_T2))*(2*hm.gamma/(1+hm.gamma))*(1/hm.Pr)*l_m;
R_g = sqrt(2)*hm.sigma*a2/(pi * kg * a_L^2 * log(1+ a2/(a1+M/(sqrt(2)*hm.sigma))));
L = (hm.gamma+1)*3*Dp/((9*hm.gamma-5)*4*l_m*sqrt(pi));
R_G = 1/(2*pi*kg*Dp*(0.5*log(1+L) + log(1+sqrt(L)) + 1/(1+sqrt(L)) – 1));
R_mic = Rs*R_g/(Rs+R_g);
R_c_inv = 1/(R_mic+R_L) + 1/R_G;

k_eff = N*(1-epsilon)*R_c_inv/(pi*Dp);

f = [k_eff*dudx(1); -u(2)*v; 0];

% compute vector s
m_dot_a = (1-epsilon)*(hm.rho_sat-u(3))*dfdt;
s = [u(2)*hm.Cpg*v*dudx(1) + m_dot_a*hm.delta_H ; -m_dot_a+hm.phi_abs; m_dot_a];

end

function u0 = icfun(x)
run(‘HM_data.m’);
u0 = hm.ic;
end

function [pl,ql,pr,qr] = bcfun(xl, ul, xr, ur, t)
run(‘HM_data.m’);
Nu_abs = 0.3+((0.62*(hm.Re_a^(0.5))*(hm.Pr_a^(1/3)))/(1+(0.4/hm.Pr_a)^(2/3))^(1/4)*((1+(hm.Re_a/282000)^(5/8))^(4/5)));
h_f = Nu_abs*hm.kf/hm.D_tank;
pl = [0 ; 0; 0];
ql = [1; 0; 0];
pr = [h_f*(ur(1)-hm.Tf_a); 0; 0];
qr = [-1; 0; 0];
end pde, pdepe MATLAB Answers — New Questions

​

Hi. I am trying to use manual variant sink programatically using a m-script and there is an error for the parameter. Can anyone suggest me what parameter is to be used for manual variant sink as in case of manual switch it is ‘sw’?
Thanks in advance.Hi. I am trying to use manual variant sink programatically using a m-script and there is an error for the parameter. Can anyone suggest me what parameter is to be used for manual variant sink as in case of manual switch it is ‘sw’?
Thanks in advance. Hi. I am trying to use manual variant sink programatically using a m-script and there is an error for the parameter. Can anyone suggest me what parameter is to be used for manual variant sink as in case of manual switch it is ‘sw’?
Thanks in advance. manual variant sink MATLAB Answers — New Questions

​

I would like to import google maps into matlab by using mapping toolbox in order to edit a specific city region of the map with dynamic view.Then everytime i click on a building, matlab should be able to display some informations, which i would have provided.I would like to import google maps into matlab by using mapping toolbox in order to edit a specific city region of the map with dynamic view.Then everytime i click on a building, matlab should be able to display some informations, which i would have provided. I would like to import google maps into matlab by using mapping toolbox in order to edit a specific city region of the map with dynamic view.Then everytime i click on a building, matlab should be able to display some informations, which i would have provided. google map, mapping toolbox MATLAB Answers — New Questions

​

Hi.
I need some help in predefining a test procedure for a remote Oven with ON/OFF Control. A step response is unsuitable because the oven is overpowered. Will a simple ON till the temperature reaches a typical working temperature followed by on OFF be sufficient? Can the collected data then be used for analysis to determine controller constants?
If that ON/OFF test is inadequate, will multiple ON/OFF’s be better as this would fall into the category of variable input?
Thanks.Hi.
I need some help in predefining a test procedure for a remote Oven with ON/OFF Control. A step response is unsuitable because the oven is overpowered. Will a simple ON till the temperature reaches a typical working temperature followed by on OFF be sufficient? Can the collected data then be used for analysis to determine controller constants?
If that ON/OFF test is inadequate, will multiple ON/OFF’s be better as this would fall into the category of variable input?
Thanks. Hi.
I need some help in predefining a test procedure for a remote Oven with ON/OFF Control. A step response is unsuitable because the oven is overpowered. Will a simple ON till the temperature reaches a typical working temperature followed by on OFF be sufficient? Can the collected data then be used for analysis to determine controller constants?
If that ON/OFF test is inadequate, will multiple ON/OFF’s be better as this would fall into the category of variable input?
Thanks. pid, pid constants, on/off control MATLAB Answers — New Questions

​

Hi. I’m trying to use matlabFunction to run solve() statements, so the app will be deployed fine. I have this example code:
if ~isdeployed()
syms x
eq = solve(x + 5 == 8, x);
fun = matlabFunction(eq, ‘File’, ‘SolveTest.m’);
clear x eq fun
else
fun = @SolveTest;
end

app.EditField.Value = SolveTest;
it works fine inside appdesigner, but when I try to deploy it by ‘deploytool’, I get this error:
Warning: In "test.mlapp", "syms" are excluded from packaging for the MATLAB Runtime environment according to the MATLAB Compiler license. Either remove the file or function from your code, or use the MATLAB function "isdeployed" to ensure the function is not invoked in the deployed component.
how can it be fixed?
I’m using R2022b.Hi. I’m trying to use matlabFunction to run solve() statements, so the app will be deployed fine. I have this example code:
if ~isdeployed()
syms x
eq = solve(x + 5 == 8, x);
fun = matlabFunction(eq, ‘File’, ‘SolveTest.m’);
clear x eq fun
else
fun = @SolveTest;
end

app.EditField.Value = SolveTest;
it works fine inside appdesigner, but when I try to deploy it by ‘deploytool’, I get this error:
Warning: In "test.mlapp", "syms" are excluded from packaging for the MATLAB Runtime environment according to the MATLAB Compiler license. Either remove the file or function from your code, or use the MATLAB function "isdeployed" to ensure the function is not invoked in the deployed component.
how can it be fixed?
I’m using R2022b. Hi. I’m trying to use matlabFunction to run solve() statements, so the app will be deployed fine. I have this example code:
if ~isdeployed()
syms x
eq = solve(x + 5 == 8, x);
fun = matlabFunction(eq, ‘File’, ‘SolveTest.m’);
clear x eq fun
else
fun = @SolveTest;
end

app.EditField.Value = SolveTest;
it works fine inside appdesigner, but when I try to deploy it by ‘deploytool’, I get this error:
Warning: In "test.mlapp", "syms" are excluded from packaging for the MATLAB Runtime environment according to the MATLAB Compiler license. Either remove the file or function from your code, or use the MATLAB function "isdeployed" to ensure the function is not invoked in the deployed component.
how can it be fixed?
I’m using R2022b. matlabfunction, appdesigner, solve MATLAB Answers — New Questions

​