Hi All!
I have a heatmap (attached below, made using imagesc) which represents the average result of 10 simulations. I would like to include the standard deviation for each parameter combination as well. What I would like to make is a 3D bar plot where the height of the bars corresponds to the standard deviation of a specific parameter combination and the color of the bars is the same as in the 2D heatmap. How could I achieve this?
Thank you for your help in advance!Hi All!
I have a heatmap (attached below, made using imagesc) which represents the average result of 10 simulations. I would like to include the standard deviation for each parameter combination as well. What I would like to make is a 3D bar plot where the height of the bars corresponds to the standard deviation of a specific parameter combination and the color of the bars is the same as in the 2D heatmap. How could I achieve this?
Thank you for your help in advance! Hi All!
I have a heatmap (attached below, made using imagesc) which represents the average result of 10 simulations. I would like to include the standard deviation for each parameter combination as well. What I would like to make is a 3D bar plot where the height of the bars corresponds to the standard deviation of a specific parameter combination and the color of the bars is the same as in the 2D heatmap. How could I achieve this?
Thank you for your help in advance! bar plot, 3d, heatmap, 2d, matlab, colormap MATLAB Answers — New Questions

​

Hi,
I have a couple of image files in ‘.webp’ format, how to read them in Matlab?Hi,
I have a couple of image files in ‘.webp’ format, how to read them in Matlab? Hi,
I have a couple of image files in ‘.webp’ format, how to read them in Matlab? webp, image file MATLAB Answers — New Questions

​

In C-Code there is a way to define a parameter using the #ifndef pre-processor directive. Is there a way to replicate this mechanism in Simulink, using data dictionaries?

For example, assume a vehicle system consisting of various subsystems. For this example assume a brake control system and an engine control system. To perform their calculations, both systems require to know the parameter "VehicleMass", so they define it as a parameter in their models. However, when they get integrated in a higher level vehicle model, the "VehicleMass" should best be described on the vehicle model level, to avoid incompatible / inconsistent definitions. But if it is define on vehicle model level, from my understanding it is currently required to first delete the lower level definitions, in order to avoid re-definition failures.
That causes another issue though, because if either of those lower level models is meant to be re-used, it can not be guaranteed that the top level integration model will always provide this value. What if the lower level model is meant to be standalone?
In C-Code, #ifndef checks if a value is already defined. If yes, that higher value is used. If no, a definition is done at this level. This mechanism prevents the issue described above.

In addition to parameters, I am experiencing the same issue with datatype definitions, so the lower level models use a given enum, but that enum might also be available in an integration environment.In C-Code there is a way to define a parameter using the #ifndef pre-processor directive. Is there a way to replicate this mechanism in Simulink, using data dictionaries?

For example, assume a vehicle system consisting of various subsystems. For this example assume a brake control system and an engine control system. To perform their calculations, both systems require to know the parameter "VehicleMass", so they define it as a parameter in their models. However, when they get integrated in a higher level vehicle model, the "VehicleMass" should best be described on the vehicle model level, to avoid incompatible / inconsistent definitions. But if it is define on vehicle model level, from my understanding it is currently required to first delete the lower level definitions, in order to avoid re-definition failures.
That causes another issue though, because if either of those lower level models is meant to be re-used, it can not be guaranteed that the top level integration model will always provide this value. What if the lower level model is meant to be standalone?
In C-Code, #ifndef checks if a value is already defined. If yes, that higher value is used. If no, a definition is done at this level. This mechanism prevents the issue described above.

In addition to parameters, I am experiencing the same issue with datatype definitions, so the lower level models use a given enum, but that enum might also be available in an integration environment. In C-Code there is a way to define a parameter using the #ifndef pre-processor directive. Is there a way to replicate this mechanism in Simulink, using data dictionaries?

For example, assume a vehicle system consisting of various subsystems. For this example assume a brake control system and an engine control system. To perform their calculations, both systems require to know the parameter "VehicleMass", so they define it as a parameter in their models. However, when they get integrated in a higher level vehicle model, the "VehicleMass" should best be described on the vehicle model level, to avoid incompatible / inconsistent definitions. But if it is define on vehicle model level, from my understanding it is currently required to first delete the lower level definitions, in order to avoid re-definition failures.
That causes another issue though, because if either of those lower level models is meant to be re-used, it can not be guaranteed that the top level integration model will always provide this value. What if the lower level model is meant to be standalone?
In C-Code, #ifndef checks if a value is already defined. If yes, that higher value is used. If no, a definition is done at this level. This mechanism prevents the issue described above.

In addition to parameters, I am experiencing the same issue with datatype definitions, so the lower level models use a given enum, but that enum might also be available in an integration environment. parameter, simulink, data dictionary MATLAB Answers — New Questions

​

am working on a MATLAB project where I need to pass multiple arguments from two structs (n and d) to a function called mi_decode. Currently, I am passing each field individually, which makes the code quite verbose and hard to maintain. Here is a snippet of my current code:
function results = spectral_mi(labels, data, n, d)

d.pl_ u int8
d.prr_k cell
d.is_cnrm logical = true
d.is_chk logical = true
d.is_fbk logical = true
d.is_x logical = false
d.is_pl logical = false
d.is_prr logical = false

n.cl uint16 = numel(unique(labels))
n.t uint8 = size(data{1}, 2)
n.k uint8 = numel(data)
n.trials uint8 = numel(labels)
n.f uint8 = fix(size(data{1}, 2)/2)+1
n.pl uint16

% … other code …

pow_per_f = cell(1, n.f);
phase_per_f = cell(1, n.f);
for f = 1:n.f
pow_per_f{f} = squeeze(power(abs(st_data(:, f, :)), 2));
phase_per_f{f} = squeeze(angle(st_data(:, f, :)));
end
[mi_, zmi_, pmi_] = mi_decode( …
labels, pow_per_f{:}, …
"is_x", d.is_x, …
"t", n.t, …
"k", n.f, …
"pl", n.pl, …
"is_cnrm", d.is_cnrm, …
"is_chk", d.is_chk, …
"is_fbk", d.is_fbk, …
"is_pl", d.is_pl, …
"is_prr", d.is_prr);
[phase_mi_, phase_zmi_, phase_pmi_] = mi_decode( …
labels, phase_per_f{:}, …
"is_x", d.is_x, …
"t", n.t, …
"k", n.f, …
"pl", n.pl, …
"is_cnrm", d.is_cnrm, …
"is_chk", d.is_chk, …
"is_fbk", d.is_fbk, …
"is_pl", d.is_pl, …
"is_prr", d.is_prr);
end
Is there a more concise way to pass all the fields from the structs n and d to the mi_decode function without listing each field individually? Any suggestions or best practices for handling this in MATLAB would be greatly appreciated.am working on a MATLAB project where I need to pass multiple arguments from two structs (n and d) to a function called mi_decode. Currently, I am passing each field individually, which makes the code quite verbose and hard to maintain. Here is a snippet of my current code:
function results = spectral_mi(labels, data, n, d)

d.pl_ u int8
d.prr_k cell
d.is_cnrm logical = true
d.is_chk logical = true
d.is_fbk logical = true
d.is_x logical = false
d.is_pl logical = false
d.is_prr logical = false

n.cl uint16 = numel(unique(labels))
n.t uint8 = size(data{1}, 2)
n.k uint8 = numel(data)
n.trials uint8 = numel(labels)
n.f uint8 = fix(size(data{1}, 2)/2)+1
n.pl uint16

% … other code …

pow_per_f = cell(1, n.f);
phase_per_f = cell(1, n.f);
for f = 1:n.f
pow_per_f{f} = squeeze(power(abs(st_data(:, f, :)), 2));
phase_per_f{f} = squeeze(angle(st_data(:, f, :)));
end
[mi_, zmi_, pmi_] = mi_decode( …
labels, pow_per_f{:}, …
"is_x", d.is_x, …
"t", n.t, …
"k", n.f, …
"pl", n.pl, …
"is_cnrm", d.is_cnrm, …
"is_chk", d.is_chk, …
"is_fbk", d.is_fbk, …
"is_pl", d.is_pl, …
"is_prr", d.is_prr);
[phase_mi_, phase_zmi_, phase_pmi_] = mi_decode( …
labels, phase_per_f{:}, …
"is_x", d.is_x, …
"t", n.t, …
"k", n.f, …
"pl", n.pl, …
"is_cnrm", d.is_cnrm, …
"is_chk", d.is_chk, …
"is_fbk", d.is_fbk, …
"is_pl", d.is_pl, …
"is_prr", d.is_prr);
end
Is there a more concise way to pass all the fields from the structs n and d to the mi_decode function without listing each field individually? Any suggestions or best practices for handling this in MATLAB would be greatly appreciated. am working on a MATLAB project where I need to pass multiple arguments from two structs (n and d) to a function called mi_decode. Currently, I am passing each field individually, which makes the code quite verbose and hard to maintain. Here is a snippet of my current code:
function results = spectral_mi(labels, data, n, d)

d.pl_ u int8
d.prr_k cell
d.is_cnrm logical = true
d.is_chk logical = true
d.is_fbk logical = true
d.is_x logical = false
d.is_pl logical = false
d.is_prr logical = false

n.cl uint16 = numel(unique(labels))
n.t uint8 = size(data{1}, 2)
n.k uint8 = numel(data)
n.trials uint8 = numel(labels)
n.f uint8 = fix(size(data{1}, 2)/2)+1
n.pl uint16

% … other code …

pow_per_f = cell(1, n.f);
phase_per_f = cell(1, n.f);
for f = 1:n.f
pow_per_f{f} = squeeze(power(abs(st_data(:, f, :)), 2));
phase_per_f{f} = squeeze(angle(st_data(:, f, :)));
end
[mi_, zmi_, pmi_] = mi_decode( …
labels, pow_per_f{:}, …
"is_x", d.is_x, …
"t", n.t, …
"k", n.f, …
"pl", n.pl, …
"is_cnrm", d.is_cnrm, …
"is_chk", d.is_chk, …
"is_fbk", d.is_fbk, …
"is_pl", d.is_pl, …
"is_prr", d.is_prr);
[phase_mi_, phase_zmi_, phase_pmi_] = mi_decode( …
labels, phase_per_f{:}, …
"is_x", d.is_x, …
"t", n.t, …
"k", n.f, …
"pl", n.pl, …
"is_cnrm", d.is_cnrm, …
"is_chk", d.is_chk, …
"is_fbk", d.is_fbk, …
"is_pl", d.is_pl, …
"is_prr", d.is_prr);
end
Is there a more concise way to pass all the fields from the structs n and d to the mi_decode function without listing each field individually? Any suggestions or best practices for handling this in MATLAB would be greatly appreciated. struct, functions, arguments MATLAB Answers — New Questions

​

Hello,

I’m thinking about grey-box modelling of heat diffusion in a rod with the 1d heat equation.
There is a interesting example at mathworls side:
https://de.mathworks.com/help/ident/ug/estimate-continuous-time-grey-box-model-for-heat-diffusion.html?status=SUCCESS
I would adapt the example (image: original) and change it a little bit to include a a convection to air at the left end of the rod (image: adapted).

Here is the code for the original version with the insulation at one end:
function [A,B,C,D,K,x0] = heatd(kappa,htf,T,Ngrid,L,temp)
% ODE file parameterizing the heat diffusion model

% kappa (first parameter) – heat diffusion coefficient
% htf (second parameter) – heat transfer coefficient
% at the far end of rod

% Auxiliary variables for computing state-space matrices:
% Ngrid: Number of points in the space-discretization
% L: Length of the rod
% temp: Initial room temperature (uniform)

% Compute space interval
deltaL = L/Ngrid;

% A matrix
A = zeros(Ngrid,Ngrid);
for kk = 2:Ngrid-1
A(kk,kk-1) = 1;
A(kk,kk) = -2;
A(kk,kk+1) = 1;
end

% Boundary condition on insulated end
A(1,1) = -1; A(1,2) = 1;
A(Ngrid,Ngrid-1) = 1;
A(Ngrid,Ngrid) = -1;
A = A*kappa/deltaL/deltaL;

% B matrix
B = zeros(Ngrid,1);
B(Ngrid,1) = htf/deltaL;

% C matrix
C = zeros(1,Ngrid);
C(1,1) = 1;

% D matrix (fixed to zero)
D = 0;

% K matrix: fixed to zero
K = zeros(Ngrid,1);

% Initial states: fixed to room temperature
x0 = temp*ones(Ngrid,1);
This section should be altered, but I’m not sure how to do it to include a known htc for air at the right end.
% Boundary condition on insulated end
A(1,1) = -1; A(1,2) = 1;
A(Ngrid,Ngrid-1) = 1;
A(Ngrid,Ngrid) = -1;
A = A*kappa/deltaL/deltaL;
Then the idgrey object is constructed:
m = idgrey(‘heatd’,{0.27 1},’c’,{10,1,22});
Another thing I don’t get is the structure of data. Unfortunately there is no explanation or example shown.
How can I implement my temperature data (for three points as two column table (t,T1) , (t,T2), (t,T3)) in the greyest?
me = greyest(data,m)
If I get it right, the parameters heat diffusion coefficient and heat transfer coefficient at the right end of rod are estimated, but I have to guess them as initial starting point? To do this the transient temperature data is required, unfortunately the example is incomplete.Hello,

I’m thinking about grey-box modelling of heat diffusion in a rod with the 1d heat equation.
There is a interesting example at mathworls side:
https://de.mathworks.com/help/ident/ug/estimate-continuous-time-grey-box-model-for-heat-diffusion.html?status=SUCCESS
I would adapt the example (image: original) and change it a little bit to include a a convection to air at the left end of the rod (image: adapted).

Here is the code for the original version with the insulation at one end:
function [A,B,C,D,K,x0] = heatd(kappa,htf,T,Ngrid,L,temp)
% ODE file parameterizing the heat diffusion model

% kappa (first parameter) – heat diffusion coefficient
% htf (second parameter) – heat transfer coefficient
% at the far end of rod

% Auxiliary variables for computing state-space matrices:
% Ngrid: Number of points in the space-discretization
% L: Length of the rod
% temp: Initial room temperature (uniform)

% Compute space interval
deltaL = L/Ngrid;

% A matrix
A = zeros(Ngrid,Ngrid);
for kk = 2:Ngrid-1
A(kk,kk-1) = 1;
A(kk,kk) = -2;
A(kk,kk+1) = 1;
end

% Boundary condition on insulated end
A(1,1) = -1; A(1,2) = 1;
A(Ngrid,Ngrid-1) = 1;
A(Ngrid,Ngrid) = -1;
A = A*kappa/deltaL/deltaL;

% B matrix
B = zeros(Ngrid,1);
B(Ngrid,1) = htf/deltaL;

% C matrix
C = zeros(1,Ngrid);
C(1,1) = 1;

% D matrix (fixed to zero)
D = 0;

% K matrix: fixed to zero
K = zeros(Ngrid,1);

% Initial states: fixed to room temperature
x0 = temp*ones(Ngrid,1);
This section should be altered, but I’m not sure how to do it to include a known htc for air at the right end.
% Boundary condition on insulated end
A(1,1) = -1; A(1,2) = 1;
A(Ngrid,Ngrid-1) = 1;
A(Ngrid,Ngrid) = -1;
A = A*kappa/deltaL/deltaL;
Then the idgrey object is constructed:
m = idgrey(‘heatd’,{0.27 1},’c’,{10,1,22});
Another thing I don’t get is the structure of data. Unfortunately there is no explanation or example shown.
How can I implement my temperature data (for three points as two column table (t,T1) , (t,T2), (t,T3)) in the greyest?
me = greyest(data,m)
If I get it right, the parameters heat diffusion coefficient and heat transfer coefficient at the right end of rod are estimated, but I have to guess them as initial starting point? To do this the transient temperature data is required, unfortunately the example is incomplete. Hello,

I’m thinking about grey-box modelling of heat diffusion in a rod with the 1d heat equation.
There is a interesting example at mathworls side:
https://de.mathworks.com/help/ident/ug/estimate-continuous-time-grey-box-model-for-heat-diffusion.html?status=SUCCESS
I would adapt the example (image: original) and change it a little bit to include a a convection to air at the left end of the rod (image: adapted).

Here is the code for the original version with the insulation at one end:
function [A,B,C,D,K,x0] = heatd(kappa,htf,T,Ngrid,L,temp)
% ODE file parameterizing the heat diffusion model

% kappa (first parameter) – heat diffusion coefficient
% htf (second parameter) – heat transfer coefficient
% at the far end of rod

% Auxiliary variables for computing state-space matrices:
% Ngrid: Number of points in the space-discretization
% L: Length of the rod
% temp: Initial room temperature (uniform)

% Compute space interval
deltaL = L/Ngrid;

% A matrix
A = zeros(Ngrid,Ngrid);
for kk = 2:Ngrid-1
A(kk,kk-1) = 1;
A(kk,kk) = -2;
A(kk,kk+1) = 1;
end

% Boundary condition on insulated end
A(1,1) = -1; A(1,2) = 1;
A(Ngrid,Ngrid-1) = 1;
A(Ngrid,Ngrid) = -1;
A = A*kappa/deltaL/deltaL;

% B matrix
B = zeros(Ngrid,1);
B(Ngrid,1) = htf/deltaL;

% C matrix
C = zeros(1,Ngrid);
C(1,1) = 1;

% D matrix (fixed to zero)
D = 0;

% K matrix: fixed to zero
K = zeros(Ngrid,1);

% Initial states: fixed to room temperature
x0 = temp*ones(Ngrid,1);
This section should be altered, but I’m not sure how to do it to include a known htc for air at the right end.
% Boundary condition on insulated end
A(1,1) = -1; A(1,2) = 1;
A(Ngrid,Ngrid-1) = 1;
A(Ngrid,Ngrid) = -1;
A = A*kappa/deltaL/deltaL;
Then the idgrey object is constructed:
m = idgrey(‘heatd’,{0.27 1},’c’,{10,1,22});
Another thing I don’t get is the structure of data. Unfortunately there is no explanation or example shown.
How can I implement my temperature data (for three points as two column table (t,T1) , (t,T2), (t,T3)) in the greyest?
me = greyest(data,m)
If I get it right, the parameters heat diffusion coefficient and heat transfer coefficient at the right end of rod are estimated, but I have to guess them as initial starting point? To do this the transient temperature data is required, unfortunately the example is incomplete. grey box, system identification, database, data import MATLAB Answers — New Questions

​

Hi team. I have a simple script which takes in stock ticker from the my list and downloads histrorical data and then plots it on line chart (X axis is date, and Y axis is close price). I want to add scatter points on the dates when I bought securities. Example – If I bought stock A on 1/2/25 and 5/2/25 then I want my codes to show two dots on the line chart on these two dates. I know scatter needs second dimension so the x value will be the dates and y value should be what is there already on the line chart. I have attached a figure of what I want to achieve.
I have a list of securities (as shown in the image column "Name"), I would like to know best way to create a list for the dates that I want as scatters (as shown in image column "Purchased (scatter)"). Please advise how to create a list of dates which can have single or multiple dates for each securities (as shown in the image below). So when loop runs for security A it should then check for all dates for A and plot them as scatter before moving to security B in second iteraton. So total number of main loop iterations would be equal to number of securities I have in Name column.
My current code takes name from the list, downloads historical prices and then plots them on line chart in Loop.
Any help would be appreciated.
kind regards
AmitHi team. I have a simple script which takes in stock ticker from the my list and downloads histrorical data and then plots it on line chart (X axis is date, and Y axis is close price). I want to add scatter points on the dates when I bought securities. Example – If I bought stock A on 1/2/25 and 5/2/25 then I want my codes to show two dots on the line chart on these two dates. I know scatter needs second dimension so the x value will be the dates and y value should be what is there already on the line chart. I have attached a figure of what I want to achieve.
I have a list of securities (as shown in the image column "Name"), I would like to know best way to create a list for the dates that I want as scatters (as shown in image column "Purchased (scatter)"). Please advise how to create a list of dates which can have single or multiple dates for each securities (as shown in the image below). So when loop runs for security A it should then check for all dates for A and plot them as scatter before moving to security B in second iteraton. So total number of main loop iterations would be equal to number of securities I have in Name column.
My current code takes name from the list, downloads historical prices and then plots them on line chart in Loop.
Any help would be appreciated.
kind regards
Amit Hi team. I have a simple script which takes in stock ticker from the my list and downloads histrorical data and then plots it on line chart (X axis is date, and Y axis is close price). I want to add scatter points on the dates when I bought securities. Example – If I bought stock A on 1/2/25 and 5/2/25 then I want my codes to show two dots on the line chart on these two dates. I know scatter needs second dimension so the x value will be the dates and y value should be what is there already on the line chart. I have attached a figure of what I want to achieve.
I have a list of securities (as shown in the image column "Name"), I would like to know best way to create a list for the dates that I want as scatters (as shown in image column "Purchased (scatter)"). Please advise how to create a list of dates which can have single or multiple dates for each securities (as shown in the image below). So when loop runs for security A it should then check for all dates for A and plot them as scatter before moving to security B in second iteraton. So total number of main loop iterations would be equal to number of securities I have in Name column.
My current code takes name from the list, downloads historical prices and then plots them on line chart in Loop.
Any help would be appreciated.
kind regards
Amit scatter and line chart superimposed in one loop MATLAB Answers — New Questions

​

Hello everyone ;
Please I need a help, i wrote a general code for determining natural frequencies and shape modes,but during the execution, i got a problem, no figure!!!
My pleasure if someone help me, please!! here is the code which i wrote it!!
%%%%—————-General Code——————-%%
clear
clc
close all
n=input(‘degrees of freedom’); %—-le nombre des masses—
M=zeros(n,n);
for i=1:n
M(i,i)=input(‘input la masse d elemnt’);
end
K=zeros(n,n);
disp(‘0 = pas deconnection’);
for i=1:n
K(i,i)=input(‘element de la raideur connecté à la terre’);
end
for i=1:n-1
for j=i+1:n
k=input(‘connection entre les elements de raideur’);
K(i,i)=K(i,i)+k;
K(i,j)=K(i,j)-k;
K(j,i)=K(j,i)-k;
K(j,j)=K(j,j)+k;
end
end
%%———-Calcul de vitesse angulaire naturelle———%%
[V,D]=eig(K,M);
omega=diag(D).^0.5;
freq_nat=omega/2*pi;
F = cell(numel(n),1);%%– forces’ vectors%%
%%%%%———Determination of eigenvalues——–%%%
C=zeros(n,n);
disp(‘0 = pas deconnection’);
for i=1:n
C(i,i)=input(‘element d amortisseur connecté à la terre’);
end
for i=1:n-1
for j=i+1:n
k=input(‘connection entre les elements damortisseur’);
C(i,i)=C(i,i)+k;
C(i,j)=C(i,j)-k;
C(j,i)=C(j,i)-k;
C(j,j)=C(j,j)+k;
end
end
%%%%————Fréquences propres amorties———%%%
A=[zeros(n,n),M;M,C];
B=[-M,zeros(n,n);zeros(n,n),K];
[V,D]=eig(full(B),full(A));
[omega,tmp]=sort(sqrt(diag(imag(D))),’ascend’);
freq_nat=omega/2/pi;
[mode,tmp]=sort((V),’ascend’);
%%%for i=1:n
%%V(:,i)=V(:,i)/V(n,i);
%end
for i=1:length(D)
zeta(i)=-(real(D(i)))/(imag(D(i)));
end
for i=1:length(D)
plot(freq_nat(i),zeta(i));
xlabel(‘frequences naturelles’);
ylabel(‘coefficient damortissement’);
end
for i=1:length(V)
plot(V(:,i));
figure()
hold on
title(‘mode (i)’)
hold off
xlabel(‘overturning moment’);
ylabel(”);
endHello everyone ;
Please I need a help, i wrote a general code for determining natural frequencies and shape modes,but during the execution, i got a problem, no figure!!!
My pleasure if someone help me, please!! here is the code which i wrote it!!
%%%%—————-General Code——————-%%
clear
clc
close all
n=input(‘degrees of freedom’); %—-le nombre des masses—
M=zeros(n,n);
for i=1:n
M(i,i)=input(‘input la masse d elemnt’);
end
K=zeros(n,n);
disp(‘0 = pas deconnection’);
for i=1:n
K(i,i)=input(‘element de la raideur connecté à la terre’);
end
for i=1:n-1
for j=i+1:n
k=input(‘connection entre les elements de raideur’);
K(i,i)=K(i,i)+k;
K(i,j)=K(i,j)-k;
K(j,i)=K(j,i)-k;
K(j,j)=K(j,j)+k;
end
end
%%———-Calcul de vitesse angulaire naturelle———%%
[V,D]=eig(K,M);
omega=diag(D).^0.5;
freq_nat=omega/2*pi;
F = cell(numel(n),1);%%– forces’ vectors%%
%%%%%———Determination of eigenvalues——–%%%
C=zeros(n,n);
disp(‘0 = pas deconnection’);
for i=1:n
C(i,i)=input(‘element d amortisseur connecté à la terre’);
end
for i=1:n-1
for j=i+1:n
k=input(‘connection entre les elements damortisseur’);
C(i,i)=C(i,i)+k;
C(i,j)=C(i,j)-k;
C(j,i)=C(j,i)-k;
C(j,j)=C(j,j)+k;
end
end
%%%%————Fréquences propres amorties———%%%
A=[zeros(n,n),M;M,C];
B=[-M,zeros(n,n);zeros(n,n),K];
[V,D]=eig(full(B),full(A));
[omega,tmp]=sort(sqrt(diag(imag(D))),’ascend’);
freq_nat=omega/2/pi;
[mode,tmp]=sort((V),’ascend’);
%%%for i=1:n
%%V(:,i)=V(:,i)/V(n,i);
%end
for i=1:length(D)
zeta(i)=-(real(D(i)))/(imag(D(i)));
end
for i=1:length(D)
plot(freq_nat(i),zeta(i));
xlabel(‘frequences naturelles’);
ylabel(‘coefficient damortissement’);
end
for i=1:length(V)
plot(V(:,i));
figure()
hold on
title(‘mode (i)’)
hold off
xlabel(‘overturning moment’);
ylabel(”);
end Hello everyone ;
Please I need a help, i wrote a general code for determining natural frequencies and shape modes,but during the execution, i got a problem, no figure!!!
My pleasure if someone help me, please!! here is the code which i wrote it!!
%%%%—————-General Code——————-%%
clear
clc
close all
n=input(‘degrees of freedom’); %—-le nombre des masses—
M=zeros(n,n);
for i=1:n
M(i,i)=input(‘input la masse d elemnt’);
end
K=zeros(n,n);
disp(‘0 = pas deconnection’);
for i=1:n
K(i,i)=input(‘element de la raideur connecté à la terre’);
end
for i=1:n-1
for j=i+1:n
k=input(‘connection entre les elements de raideur’);
K(i,i)=K(i,i)+k;
K(i,j)=K(i,j)-k;
K(j,i)=K(j,i)-k;
K(j,j)=K(j,j)+k;
end
end
%%———-Calcul de vitesse angulaire naturelle———%%
[V,D]=eig(K,M);
omega=diag(D).^0.5;
freq_nat=omega/2*pi;
F = cell(numel(n),1);%%– forces’ vectors%%
%%%%%———Determination of eigenvalues——–%%%
C=zeros(n,n);
disp(‘0 = pas deconnection’);
for i=1:n
C(i,i)=input(‘element d amortisseur connecté à la terre’);
end
for i=1:n-1
for j=i+1:n
k=input(‘connection entre les elements damortisseur’);
C(i,i)=C(i,i)+k;
C(i,j)=C(i,j)-k;
C(j,i)=C(j,i)-k;
C(j,j)=C(j,j)+k;
end
end
%%%%————Fréquences propres amorties———%%%
A=[zeros(n,n),M;M,C];
B=[-M,zeros(n,n);zeros(n,n),K];
[V,D]=eig(full(B),full(A));
[omega,tmp]=sort(sqrt(diag(imag(D))),’ascend’);
freq_nat=omega/2/pi;
[mode,tmp]=sort((V),’ascend’);
%%%for i=1:n
%%V(:,i)=V(:,i)/V(n,i);
%end
for i=1:length(D)
zeta(i)=-(real(D(i)))/(imag(D(i)));
end
for i=1:length(D)
plot(freq_nat(i),zeta(i));
xlabel(‘frequences naturelles’);
ylabel(‘coefficient damortissement’);
end
for i=1:length(V)
plot(V(:,i));
figure()
hold on
title(‘mode (i)’)
hold off
xlabel(‘overturning moment’);
ylabel(”);
end general code, natural frequencies, mode shapes, critical speeds MATLAB Answers — New Questions

​

Hi,

I’m currently stuck on performing a two-way repeated-measures Manova in Matlab using the RepeatedMeasuresModel. My data set consists of six units, from each of which I measured three properties, which might be highly correlated. Each unit was tested using two within-subject factors each of which has two levels (4 conditions in total). Thus, my data table consists of six rows and twelve columns (3 properties * 4 conditions).

I created the within-table as follows:

within = table({‘A’;’A’;’A’;’A’;’B’;’B’;’B’;’B’;’C’;’C’;’C’;’C’},…
{‘1′;’2′;’1′;’2′;’1′;’2′;’1′;’2′;’1′;’2′;’1′;’2’},…
{‘a’;’a’;’b’;’b’;’a’;’a’;’b’;’b’;’a’;’a’;’b’;’b’},…
‘VariableNames’,{‘Property’,’Condition1′,’Condition2′});

The table looks like:

Property Condition1 Condition2
________ __________ __________

‘A’ ‘1’ ‘a’
‘A’ ‘2’ ‘a’
‘A’ ‘1’ ‘b’
‘A’ ‘2’ ‘b’
‘B’ ‘1’ ‘a’
‘B’ ‘2’ ‘a’
‘B’ ‘1’ ‘b’
‘B’ ‘2’ ‘b’
‘C’ ‘1’ ‘a’
‘C’ ‘2’ ‘a’
‘C’ ‘1’ ‘b’
‘C’ ‘2’ ‘b’

The repeated-measures Model I created using

rm = fitrm(t,’M1-M12~1′,’withinDesign’,within);

where t is the data table and M1-M12 are the responses.

If I now use

manova(rm)

the output looks like this:

Statistic Value F RSquare df1 df2 pValue
_________ _______ ___ _______ ___ ___ ______

Pillai 0.54237 0 0.54237 6 0 NaN
Wilks NaN NaN NaN 6 0 NaN
Hotelling NaN NaN NaN 6 NaN NaN
Roy Inf NaN NaN 6 0 NaN

I don’t know, how I can specify, which of my within-subject factors are the different properties and which ones are the different conditions.

Help of any kind would be highly appreciated.

Best,
ChristianHi,

I’m currently stuck on performing a two-way repeated-measures Manova in Matlab using the RepeatedMeasuresModel. My data set consists of six units, from each of which I measured three properties, which might be highly correlated. Each unit was tested using two within-subject factors each of which has two levels (4 conditions in total). Thus, my data table consists of six rows and twelve columns (3 properties * 4 conditions).

I created the within-table as follows:

within = table({‘A’;’A’;’A’;’A’;’B’;’B’;’B’;’B’;’C’;’C’;’C’;’C’},…
{‘1′;’2′;’1′;’2′;’1′;’2′;’1′;’2′;’1′;’2′;’1′;’2’},…
{‘a’;’a’;’b’;’b’;’a’;’a’;’b’;’b’;’a’;’a’;’b’;’b’},…
‘VariableNames’,{‘Property’,’Condition1′,’Condition2′});

The table looks like:

Property Condition1 Condition2
________ __________ __________

‘A’ ‘1’ ‘a’
‘A’ ‘2’ ‘a’
‘A’ ‘1’ ‘b’
‘A’ ‘2’ ‘b’
‘B’ ‘1’ ‘a’
‘B’ ‘2’ ‘a’
‘B’ ‘1’ ‘b’
‘B’ ‘2’ ‘b’
‘C’ ‘1’ ‘a’
‘C’ ‘2’ ‘a’
‘C’ ‘1’ ‘b’
‘C’ ‘2’ ‘b’

The repeated-measures Model I created using

rm = fitrm(t,’M1-M12~1′,’withinDesign’,within);

where t is the data table and M1-M12 are the responses.

If I now use

manova(rm)

the output looks like this:

Statistic Value F RSquare df1 df2 pValue
_________ _______ ___ _______ ___ ___ ______

Pillai 0.54237 0 0.54237 6 0 NaN
Wilks NaN NaN NaN 6 0 NaN
Hotelling NaN NaN NaN 6 NaN NaN
Roy Inf NaN NaN 6 0 NaN

I don’t know, how I can specify, which of my within-subject factors are the different properties and which ones are the different conditions.

Help of any kind would be highly appreciated.

Best,
Christian Hi,

I’m currently stuck on performing a two-way repeated-measures Manova in Matlab using the RepeatedMeasuresModel. My data set consists of six units, from each of which I measured three properties, which might be highly correlated. Each unit was tested using two within-subject factors each of which has two levels (4 conditions in total). Thus, my data table consists of six rows and twelve columns (3 properties * 4 conditions).

I created the within-table as follows:

within = table({‘A’;’A’;’A’;’A’;’B’;’B’;’B’;’B’;’C’;’C’;’C’;’C’},…
{‘1′;’2′;’1′;’2′;’1′;’2′;’1′;’2′;’1′;’2′;’1′;’2’},…
{‘a’;’a’;’b’;’b’;’a’;’a’;’b’;’b’;’a’;’a’;’b’;’b’},…
‘VariableNames’,{‘Property’,’Condition1′,’Condition2′});

The table looks like:

Property Condition1 Condition2
________ __________ __________

‘A’ ‘1’ ‘a’
‘A’ ‘2’ ‘a’
‘A’ ‘1’ ‘b’
‘A’ ‘2’ ‘b’
‘B’ ‘1’ ‘a’
‘B’ ‘2’ ‘a’
‘B’ ‘1’ ‘b’
‘B’ ‘2’ ‘b’
‘C’ ‘1’ ‘a’
‘C’ ‘2’ ‘a’
‘C’ ‘1’ ‘b’
‘C’ ‘2’ ‘b’

The repeated-measures Model I created using

rm = fitrm(t,’M1-M12~1′,’withinDesign’,within);

where t is the data table and M1-M12 are the responses.

If I now use

manova(rm)

the output looks like this:

Statistic Value F RSquare df1 df2 pValue
_________ _______ ___ _______ ___ ___ ______

Pillai 0.54237 0 0.54237 6 0 NaN
Wilks NaN NaN NaN 6 0 NaN
Hotelling NaN NaN NaN 6 NaN NaN
Roy Inf NaN NaN 6 0 NaN

I don’t know, how I can specify, which of my within-subject factors are the different properties and which ones are the different conditions.

Help of any kind would be highly appreciated.

Best,
Christian manova, repeatedmeasures MATLAB Answers — New Questions

​

few day before, i was unable to open desktop window on PC. I have tried any answer form Math Work, but still failed. than i decided to reinstall the computer and installed again. After installed desktop window, i logged my friend’s account (cause i can’t recive the verify mail). Everything is down.

here is the problem
if my desktop window can open is system promblem, which is easy to slove.
What if the problem is the account verify fail cause desktop window don’t show up. How to log out my account before open Matlab desktop window ?few day before, i was unable to open desktop window on PC. I have tried any answer form Math Work, but still failed. than i decided to reinstall the computer and installed again. After installed desktop window, i logged my friend’s account (cause i can’t recive the verify mail). Everything is down.

here is the problem
if my desktop window can open is system promblem, which is easy to slove.
What if the problem is the account verify fail cause desktop window don’t show up. How to log out my account before open Matlab desktop window ? few day before, i was unable to open desktop window on PC. I have tried any answer form Math Work, but still failed. than i decided to reinstall the computer and installed again. After installed desktop window, i logged my friend’s account (cause i can’t recive the verify mail). Everything is down.

here is the problem
if my desktop window can open is system promblem, which is easy to slove.
What if the problem is the account verify fail cause desktop window don’t show up. How to log out my account before open Matlab desktop window ? log in/out MATLAB Answers — New Questions

​

Hello,
how can I plot the first k terms of a Fourier series on the same graph (without showing the Gibb’s effect). I would also like to plot the amplitude of those terms over frequency. I want the graphs on the attached picture. How can I make them with Matlab?
Thanks in advanceHello,
how can I plot the first k terms of a Fourier series on the same graph (without showing the Gibb’s effect). I would also like to plot the amplitude of those terms over frequency. I want the graphs on the attached picture. How can I make them with Matlab?
Thanks in advance Hello,
how can I plot the first k terms of a Fourier series on the same graph (without showing the Gibb’s effect). I would also like to plot the amplitude of those terms over frequency. I want the graphs on the attached picture. How can I make them with Matlab?
Thanks in advance fourier, plotting MATLAB Answers — New Questions

​