Tag Archives: matlab
Click disp() in the Icon Drawing Commands palette to add the command to the editable region.
Open the Mask Editor dialog and navigate to the Icon tab.
Click disp() in the Icon Drawing Commands palette to add the command to the editable region.
Replace the sample text with the variable from the previous task, label.
In the Property Editor at right, Change Run initialization to On.
Save the mask.Open the Mask Editor dialog and navigate to the Icon tab.
Click disp() in the Icon Drawing Commands palette to add the command to the editable region.
Replace the sample text with the variable from the previous task, label.
In the Property Editor at right, Change Run initialization to On.
Save the mask. Open the Mask Editor dialog and navigate to the Icon tab.
Click disp() in the Icon Drawing Commands palette to add the command to the editable region.
Replace the sample text with the variable from the previous task, label.
In the Property Editor at right, Change Run initialization to On.
Save the mask. display MATLAB Answers — New Questions
Rejected stereo image pairs in Stereo Camera Calibrator App
Hello all,
I’m trying to use the “Stereo Camera Calibrator App”.
I’m trying to import four pairs of images, but it only accepts one pair and rejects the other three. I followed the advice from the camera preparation page, but it didn’t do any good. I have no clue why it rejects most of the image pairs.
Please see the attached four images of each camera (left and right). The checkerboard size is 3 mm.
Thank you.Hello all,
I’m trying to use the “Stereo Camera Calibrator App”.
I’m trying to import four pairs of images, but it only accepts one pair and rejects the other three. I followed the advice from the camera preparation page, but it didn’t do any good. I have no clue why it rejects most of the image pairs.
Please see the attached four images of each camera (left and right). The checkerboard size is 3 mm.
Thank you. Hello all,
I’m trying to use the “Stereo Camera Calibrator App”.
I’m trying to import four pairs of images, but it only accepts one pair and rejects the other three. I followed the advice from the camera preparation page, but it didn’t do any good. I have no clue why it rejects most of the image pairs.
Please see the attached four images of each camera (left and right). The checkerboard size is 3 mm.
Thank you. image processing, stereo camera calibrator app MATLAB Answers — New Questions
Calculate the normal vector between two nodes in the space
Hi! I have two nodes A and B. These are positioned in space with the same value of X and Y but different Z. The normal in this case would be n=[0 0 1] (red plane parallel to the XY plane).
I would like to calculate the normal n ‘green’ in the case of nodes A and C. Is this possible?
A = [-33.24 -10.70 7.41];
B = [A(1,1), A(1,2), A(1,3)+5];
N = [A;B];
C = [A(1,1), A(1,2)+7, A(1,3)+5];
NN = [A;C];
figure
plot3(A(1,1),A(1,2),A(1,3),’k.’,’Markersize’,30);
hold on
plot3(B(1,1),B(1,2),B(1,3),’r.’,’Markersize’,30);
plot3(C(1,1),C(1,2),C(1,3),’g.’,’Markersize’,30);
plot3(N(:,1),N(:,2),N(:,3),’-k’,’LineWidth’,1);
plot3(NN(:,1),NN(:,2),NN(:,3),’-k’,’LineWidth’,1);
hold off
axis equal
(YZ view)Hi! I have two nodes A and B. These are positioned in space with the same value of X and Y but different Z. The normal in this case would be n=[0 0 1] (red plane parallel to the XY plane).
I would like to calculate the normal n ‘green’ in the case of nodes A and C. Is this possible?
A = [-33.24 -10.70 7.41];
B = [A(1,1), A(1,2), A(1,3)+5];
N = [A;B];
C = [A(1,1), A(1,2)+7, A(1,3)+5];
NN = [A;C];
figure
plot3(A(1,1),A(1,2),A(1,3),’k.’,’Markersize’,30);
hold on
plot3(B(1,1),B(1,2),B(1,3),’r.’,’Markersize’,30);
plot3(C(1,1),C(1,2),C(1,3),’g.’,’Markersize’,30);
plot3(N(:,1),N(:,2),N(:,3),’-k’,’LineWidth’,1);
plot3(NN(:,1),NN(:,2),NN(:,3),’-k’,’LineWidth’,1);
hold off
axis equal
(YZ view) Hi! I have two nodes A and B. These are positioned in space with the same value of X and Y but different Z. The normal in this case would be n=[0 0 1] (red plane parallel to the XY plane).
I would like to calculate the normal n ‘green’ in the case of nodes A and C. Is this possible?
A = [-33.24 -10.70 7.41];
B = [A(1,1), A(1,2), A(1,3)+5];
N = [A;B];
C = [A(1,1), A(1,2)+7, A(1,3)+5];
NN = [A;C];
figure
plot3(A(1,1),A(1,2),A(1,3),’k.’,’Markersize’,30);
hold on
plot3(B(1,1),B(1,2),B(1,3),’r.’,’Markersize’,30);
plot3(C(1,1),C(1,2),C(1,3),’g.’,’Markersize’,30);
plot3(N(:,1),N(:,2),N(:,3),’-k’,’LineWidth’,1);
plot3(NN(:,1),NN(:,2),NN(:,3),’-k’,’LineWidth’,1);
hold off
axis equal
(YZ view) normal, vector MATLAB Answers — New Questions
my Simcape Electrical model does not want to wire the reference (earth ) block
As you see below the reference block desnot want to connect to the volatge source or any other block. I am using MATLAB R2023b academic version. Can someone give me a clue what I have to look atAs you see below the reference block desnot want to connect to the volatge source or any other block. I am using MATLAB R2023b academic version. Can someone give me a clue what I have to look at As you see below the reference block desnot want to connect to the volatge source or any other block. I am using MATLAB R2023b academic version. Can someone give me a clue what I have to look at r2023b MATLAB Answers — New Questions
libstdc++.so.6: version `GLIBCXX_3.4.29′ not found
When I run the MATLAB on my arch linux computer, I met the following issue. I am using R2021b, it works fine on this arch linux computer before. This problems happens when I change the font of matlab and reboot.
MATLAB is selecting SOFTWARE OPENGL rendering.
MESA-LOADER: failed to open radeonsi: /home/deck/matlab/bin/glnxa64/../../sys/os/g
lnxa64/libstdc++.so.6: version `GLIBCXX_3.4.29′ not found (required by /usr/lib/dr
i/radeonsi_dri.so) (search paths /usr/lib/dri, suffix _dri)
failed to load driver: radeonsi
MESA-LOADER: failed to open kms_swrast: /home/deck/matlab/bin/glnxa64/../../sys/os
/glnxa64/libstdc++.so.6: version `GLIBCXX_3.4.29′ not found (required by /usr/lib/
dri/kms_swrast_dri.so) (search paths /usr/lib/dri, suffix _dri)
failed to load driver: kms_swrast
MESA-LOADER: failed to open swrast: /home/deck/matlab/bin/glnxa64/../../sys/os/gln
xa64/libstdc++.so.6: version `GLIBCXX_3.4.29′ not found (required by /usr/lib/dri/
swrast_dri.so) (search paths /usr/lib/dri, suffix _dri)
failed to load swrast driver
Gtk-Message: 00:00:13.007: Failed to load module "canberra-gtk-module"When I run the MATLAB on my arch linux computer, I met the following issue. I am using R2021b, it works fine on this arch linux computer before. This problems happens when I change the font of matlab and reboot.
MATLAB is selecting SOFTWARE OPENGL rendering.
MESA-LOADER: failed to open radeonsi: /home/deck/matlab/bin/glnxa64/../../sys/os/g
lnxa64/libstdc++.so.6: version `GLIBCXX_3.4.29′ not found (required by /usr/lib/dr
i/radeonsi_dri.so) (search paths /usr/lib/dri, suffix _dri)
failed to load driver: radeonsi
MESA-LOADER: failed to open kms_swrast: /home/deck/matlab/bin/glnxa64/../../sys/os
/glnxa64/libstdc++.so.6: version `GLIBCXX_3.4.29′ not found (required by /usr/lib/
dri/kms_swrast_dri.so) (search paths /usr/lib/dri, suffix _dri)
failed to load driver: kms_swrast
MESA-LOADER: failed to open swrast: /home/deck/matlab/bin/glnxa64/../../sys/os/gln
xa64/libstdc++.so.6: version `GLIBCXX_3.4.29′ not found (required by /usr/lib/dri/
swrast_dri.so) (search paths /usr/lib/dri, suffix _dri)
failed to load swrast driver
Gtk-Message: 00:00:13.007: Failed to load module "canberra-gtk-module" When I run the MATLAB on my arch linux computer, I met the following issue. I am using R2021b, it works fine on this arch linux computer before. This problems happens when I change the font of matlab and reboot.
MATLAB is selecting SOFTWARE OPENGL rendering.
MESA-LOADER: failed to open radeonsi: /home/deck/matlab/bin/glnxa64/../../sys/os/g
lnxa64/libstdc++.so.6: version `GLIBCXX_3.4.29′ not found (required by /usr/lib/dr
i/radeonsi_dri.so) (search paths /usr/lib/dri, suffix _dri)
failed to load driver: radeonsi
MESA-LOADER: failed to open kms_swrast: /home/deck/matlab/bin/glnxa64/../../sys/os
/glnxa64/libstdc++.so.6: version `GLIBCXX_3.4.29′ not found (required by /usr/lib/
dri/kms_swrast_dri.so) (search paths /usr/lib/dri, suffix _dri)
failed to load driver: kms_swrast
MESA-LOADER: failed to open swrast: /home/deck/matlab/bin/glnxa64/../../sys/os/gln
xa64/libstdc++.so.6: version `GLIBCXX_3.4.29′ not found (required by /usr/lib/dri/
swrast_dri.so) (search paths /usr/lib/dri, suffix _dri)
failed to load swrast driver
Gtk-Message: 00:00:13.007: Failed to load module "canberra-gtk-module" installation, linux MATLAB Answers — New Questions
How to enable “Snap to Tick” for a slider in App Designer?
Hi, I want to know how to make a slider move discretely rather than continuously. As an example, I want the slider to jump from 1 to 2 when I move it. In other words, I need to enable the "Snap to Tick" feature.
Thanks in advance.Hi, I want to know how to make a slider move discretely rather than continuously. As an example, I want the slider to jump from 1 to 2 when I move it. In other words, I need to enable the "Snap to Tick" feature.
Thanks in advance. Hi, I want to know how to make a slider move discretely rather than continuously. As an example, I want the slider to jump from 1 to 2 when I move it. In other words, I need to enable the "Snap to Tick" feature.
Thanks in advance. app designer, snap to tick, slider MATLAB Answers — New Questions
bifurcation for Fractional order
I need an example matlab code to draw a bifurcation diagram for a fractional orde system. (for parametre varying or FO q varying, as seen in figures)
I’ve done some research on the subject and found The code by @Marius-F. Danca (2023) can be found at this File Exchange link:
Matlab code for FO logistic map, MATLAB Central File Exchange. Retrieved October 24, 2023.
However, I can’t plot the bifurcation diagram for the FO case.
code for the graphical representation of the bifurcation proposed by Sam Chak
interval = [2.8 4.0];
accuracy = 0.001;
reps = 100; % number of repetitions
numtoplot = 50;
initidx = reps – (numtoplot – 1);
for r = interval(1):accuracy:interval(2)
x = 0.4;
xo = x;
for n = 2:reps
xn = r*xo*(1 – xo);
x = [x xn];
xo = xn;
end
plot(r*ones(numtoplot), x(initidx:reps), ‘.’, ‘MarkerSize’, 1, ‘Color’, ‘#be4d25’),
hold on
end
hold off
xlim([2.8 4])
title(‘Logistic Map’)
xlabel(‘itr’, ‘FontSize’, 14)
ylabel(‘itx’, ‘FontSize’, 14)
However, the fractional order need a workaround, i fond the code by @Marius-F. Danca
but I don’t know how it works, Can anyone kindly share *.m file
function [u]=FO_discrete(q,p,u0,n_max)
% Author: Marius-F. Danca 2021
u=zeros(n_max,1); % memory alocation
u(1)=u0+f(u0); % u(1) outside of the loops
for n=2:n_max; % $u(2),u(3),…,n(n_max)$
s=exp(gammaln(n-1+q)-gammaln(n))*f(u0);
for j=2:n;
s=s+exp(gammaln(n-j+q)-gammaln(n-j+1))*f(u(j-1));
end
u(n)=u0+s/gamma(q);
end
function ff=f(x); % logistic map
ff=p*x*(1-x);
end
endI need an example matlab code to draw a bifurcation diagram for a fractional orde system. (for parametre varying or FO q varying, as seen in figures)
I’ve done some research on the subject and found The code by @Marius-F. Danca (2023) can be found at this File Exchange link:
Matlab code for FO logistic map, MATLAB Central File Exchange. Retrieved October 24, 2023.
However, I can’t plot the bifurcation diagram for the FO case.
code for the graphical representation of the bifurcation proposed by Sam Chak
interval = [2.8 4.0];
accuracy = 0.001;
reps = 100; % number of repetitions
numtoplot = 50;
initidx = reps – (numtoplot – 1);
for r = interval(1):accuracy:interval(2)
x = 0.4;
xo = x;
for n = 2:reps
xn = r*xo*(1 – xo);
x = [x xn];
xo = xn;
end
plot(r*ones(numtoplot), x(initidx:reps), ‘.’, ‘MarkerSize’, 1, ‘Color’, ‘#be4d25’),
hold on
end
hold off
xlim([2.8 4])
title(‘Logistic Map’)
xlabel(‘itr’, ‘FontSize’, 14)
ylabel(‘itx’, ‘FontSize’, 14)
However, the fractional order need a workaround, i fond the code by @Marius-F. Danca
but I don’t know how it works, Can anyone kindly share *.m file
function [u]=FO_discrete(q,p,u0,n_max)
% Author: Marius-F. Danca 2021
u=zeros(n_max,1); % memory alocation
u(1)=u0+f(u0); % u(1) outside of the loops
for n=2:n_max; % $u(2),u(3),…,n(n_max)$
s=exp(gammaln(n-1+q)-gammaln(n))*f(u0);
for j=2:n;
s=s+exp(gammaln(n-j+q)-gammaln(n-j+1))*f(u(j-1));
end
u(n)=u0+s/gamma(q);
end
function ff=f(x); % logistic map
ff=p*x*(1-x);
end
end I need an example matlab code to draw a bifurcation diagram for a fractional orde system. (for parametre varying or FO q varying, as seen in figures)
I’ve done some research on the subject and found The code by @Marius-F. Danca (2023) can be found at this File Exchange link:
Matlab code for FO logistic map, MATLAB Central File Exchange. Retrieved October 24, 2023.
However, I can’t plot the bifurcation diagram for the FO case.
code for the graphical representation of the bifurcation proposed by Sam Chak
interval = [2.8 4.0];
accuracy = 0.001;
reps = 100; % number of repetitions
numtoplot = 50;
initidx = reps – (numtoplot – 1);
for r = interval(1):accuracy:interval(2)
x = 0.4;
xo = x;
for n = 2:reps
xn = r*xo*(1 – xo);
x = [x xn];
xo = xn;
end
plot(r*ones(numtoplot), x(initidx:reps), ‘.’, ‘MarkerSize’, 1, ‘Color’, ‘#be4d25’),
hold on
end
hold off
xlim([2.8 4])
title(‘Logistic Map’)
xlabel(‘itr’, ‘FontSize’, 14)
ylabel(‘itx’, ‘FontSize’, 14)
However, the fractional order need a workaround, i fond the code by @Marius-F. Danca
but I don’t know how it works, Can anyone kindly share *.m file
function [u]=FO_discrete(q,p,u0,n_max)
% Author: Marius-F. Danca 2021
u=zeros(n_max,1); % memory alocation
u(1)=u0+f(u0); % u(1) outside of the loops
for n=2:n_max; % $u(2),u(3),…,n(n_max)$
s=exp(gammaln(n-1+q)-gammaln(n))*f(u0);
for j=2:n;
s=s+exp(gammaln(n-j+q)-gammaln(n-j+1))*f(u(j-1));
end
u(n)=u0+s/gamma(q);
end
function ff=f(x); % logistic map
ff=p*x*(1-x);
end
end mar MATLAB Answers — New Questions
How can i print at app designer on a Text Area box, text with hyperlinks?
Hi,
How can i print in Text area box (with app designer), hyperlinks for files that located in local PC?
when im using the sprintf function rom conslue window i can see the print with hypelink, for example:
sprintf(‘<a href="matlab: winopen(”%s”)">%s</a>nn’, filepath, CurrentFileName),
but when im trying to print the same inside the text area its printing me the whole string:
How can i print hyperlink inside the text area?Hi,
How can i print in Text area box (with app designer), hyperlinks for files that located in local PC?
when im using the sprintf function rom conslue window i can see the print with hypelink, for example:
sprintf(‘<a href="matlab: winopen(”%s”)">%s</a>nn’, filepath, CurrentFileName),
but when im trying to print the same inside the text area its printing me the whole string:
How can i print hyperlink inside the text area? Hi,
How can i print in Text area box (with app designer), hyperlinks for files that located in local PC?
when im using the sprintf function rom conslue window i can see the print with hypelink, for example:
sprintf(‘<a href="matlab: winopen(”%s”)">%s</a>nn’, filepath, CurrentFileName),
but when im trying to print the same inside the text area its printing me the whole string:
How can i print hyperlink inside the text area? appdesigner, hyperlink MATLAB Answers — New Questions
Live Script opening in new window please help
When I open a new live script it opens it in a new window instead of splitting the screen in half vertically! Its drivng me crazy but I cant figure out what option in the toolbar editer it is.When I open a new live script it opens it in a new window instead of splitting the screen in half vertically! Its drivng me crazy but I cant figure out what option in the toolbar editer it is. When I open a new live script it opens it in a new window instead of splitting the screen in half vertically! Its drivng me crazy but I cant figure out what option in the toolbar editer it is. workspace layout MATLAB Answers — New Questions
how to azimuthally average a N*N data ?
I have a 291*291 velocity data from LES . I want to azimuthally average the data so that for any angle from the origin at a particular radial distance i get the same value of velocity.
please suggest the way to do it.
Following figure shows the pictorial reperesentation of velocity data in yz plane
here is the attached file for the dataI have a 291*291 velocity data from LES . I want to azimuthally average the data so that for any angle from the origin at a particular radial distance i get the same value of velocity.
please suggest the way to do it.
Following figure shows the pictorial reperesentation of velocity data in yz plane
here is the attached file for the data I have a 291*291 velocity data from LES . I want to azimuthally average the data so that for any angle from the origin at a particular radial distance i get the same value of velocity.
please suggest the way to do it.
Following figure shows the pictorial reperesentation of velocity data in yz plane
here is the attached file for the data azimuthally-average data MATLAB Answers — New Questions
How matlab reads ansys files
I used ansys_aas toolbox to link between ansys workbench and matlab but i’m confused of how can the matlab reads the chart in ansys and make optimization to get minimum error between desired and actual curvesI used ansys_aas toolbox to link between ansys workbench and matlab but i’m confused of how can the matlab reads the chart in ansys and make optimization to get minimum error between desired and actual curves I used ansys_aas toolbox to link between ansys workbench and matlab but i’m confused of how can the matlab reads the chart in ansys and make optimization to get minimum error between desired and actual curves ansys workbench, matlab MATLAB Answers — New Questions
problem with loops and conditional statements
I am writing a simple program with if loop inside a for loop but the program gets stopped in middle with error:
Index exceeds the number of array elements. Index must not exceed 1.
Error in newtrial (line 19)
Y(i) = Y(i) + 1;
and the plot is not generated as well even for values generated.
the code is:
clc
clear
x1 = 5
y1 = 3
x0 = 0
y0 = 0
dx = x1-x0
dy = y1-y0
e = dy-dx
X = x0
Y = y0
for i = 1:1:5;
X(i+1) = X(i) + 1
if e(i) >= 0
Y(i+1) = Y(i) + 1;
e(i+1) = e(i) – dx;
else
e(i+1) = e(i) + dy;
end
end
Y
plot(X,Y,’bo’)
axis([0 6 0 6])
I am new to this, so any help or correction in code is appreciated, thanks in advance.I am writing a simple program with if loop inside a for loop but the program gets stopped in middle with error:
Index exceeds the number of array elements. Index must not exceed 1.
Error in newtrial (line 19)
Y(i) = Y(i) + 1;
and the plot is not generated as well even for values generated.
the code is:
clc
clear
x1 = 5
y1 = 3
x0 = 0
y0 = 0
dx = x1-x0
dy = y1-y0
e = dy-dx
X = x0
Y = y0
for i = 1:1:5;
X(i+1) = X(i) + 1
if e(i) >= 0
Y(i+1) = Y(i) + 1;
e(i+1) = e(i) – dx;
else
e(i+1) = e(i) + dy;
end
end
Y
plot(X,Y,’bo’)
axis([0 6 0 6])
I am new to this, so any help or correction in code is appreciated, thanks in advance. I am writing a simple program with if loop inside a for loop but the program gets stopped in middle with error:
Index exceeds the number of array elements. Index must not exceed 1.
Error in newtrial (line 19)
Y(i) = Y(i) + 1;
and the plot is not generated as well even for values generated.
the code is:
clc
clear
x1 = 5
y1 = 3
x0 = 0
y0 = 0
dx = x1-x0
dy = y1-y0
e = dy-dx
X = x0
Y = y0
for i = 1:1:5;
X(i+1) = X(i) + 1
if e(i) >= 0
Y(i+1) = Y(i) + 1;
e(i+1) = e(i) – dx;
else
e(i+1) = e(i) + dy;
end
end
Y
plot(X,Y,’bo’)
axis([0 6 0 6])
I am new to this, so any help or correction in code is appreciated, thanks in advance. index error, loops MATLAB Answers — New Questions
Processing EMG data in MATLAB?
Hi all,
I am VERY new to MATLAB, so please excuse my ignorance if this is a very novice problem. I am analysing EMG signals, I have imported from my data from xlsx as time (X) and voltage (y).
After extracting data;
set(0, ‘defaultFigureToolbar’, ‘none’)
% RAW_DATA_EXTRACTION
%% BACK_SQUAT
%%% VISIT_ONE
T1_10rm_BS_1 = xlsread(‘T1_70%_10rm_BS_1.xlsx’,’Sheet3′,’A5:E77494′);
time = T1_10rm_BS_1(:,1);
BS_RF_1 = T1_10rm_BS_1(:,4);
BS_BF_1 = T1_10rm_BS_1(:,5);
T1_10rm_BS_2 = xlsread(‘T1_70%_10rm_BS_2.xlsx’, ‘Sheet4′,’A5:E78394’);
time = T1_10rm_BS_2(:,1);
BS_RF_2 = T1_10rm_BS_2(:,4);
BS_BF_2 = T1_10rm_BS_2(:,5);
running filters;
%FILTER
% Lowpass filter 500Hz
Freq = 1500;
NyqFreq = Freq/2;
Lowpass = 500;
Wn = Lowpass/NyqFreq;
[B,A] = butter (4,Wn,’low’);
% run filter
BS_RF_1_low_filtered = filtfilt(B,A,BS_RF_1);
BS_BF_1_low_filtered = filtfilt(B,A,BS_BF_1);
BS_RF_2_low_filtered = filtfilt(B,A,BS_RF_2);
BS_BF_2_low_filtered = filtfilt(B,A,BS_BF_2);
% Highpass filter 10Hz
Highpass = 10;
Wo = Highpass/NyqFreq;
[D,C] = butter (4,Wo,’high’);
% run filter
high_pass_BS_RF_1_filtdata = filtfilt(D,C,BS_RF_1_low_filtered);
high_pass_BS_BF_1_filtdata = filtfilt(D,C,BS_BF_1_low_filtered);
high_pass_BS_RF_2_filtdata = filtfilt(D,C,BS_RF_2_low_filtered);
high_pass_BS_BF_2_filtdata = filtfilt(D,C,BS_BF_2_low_filtered);
% Full wave rectification of EMG
filtered_BS_RF_1 = abs(high_pass_BS_RF_1_filtdata);
filtered_BS_BF_1 = abs(high_pass_BS_BF_1_filtdata);
filtered_BS_RF_2 = abs(high_pass_BS_RF_2_filtdata);
filtered_BS_BF_2 = abs(high_pass_BS_BF_2_filtdata);
% Linear envelope (2Hz lowpass)
LP = 2;
Wp = LP/NyqFreq;
[F,E] = butter (4,Wp,’low’);
linear_BS_RF_1 = filtfilt(F,E,filtered_BS_RF_1);
linear_BS_BF_1 = filtfilt(F,E,filtered_BS_BF_1);
linear_BS_RF_2 = filtfilt(F,E,filtered_BS_RF_2);
linear_BS_BF_2 = filtfilt(F,E,filtered_BS_BF_2);
and generating my plot (for BS_RF_1);
I am now struggling to try and find the area under the curve for my set of 10 repetitions.
I am looking to define mean activation and peak activation.
All constructive help welcome,
Cheers.Hi all,
I am VERY new to MATLAB, so please excuse my ignorance if this is a very novice problem. I am analysing EMG signals, I have imported from my data from xlsx as time (X) and voltage (y).
After extracting data;
set(0, ‘defaultFigureToolbar’, ‘none’)
% RAW_DATA_EXTRACTION
%% BACK_SQUAT
%%% VISIT_ONE
T1_10rm_BS_1 = xlsread(‘T1_70%_10rm_BS_1.xlsx’,’Sheet3′,’A5:E77494′);
time = T1_10rm_BS_1(:,1);
BS_RF_1 = T1_10rm_BS_1(:,4);
BS_BF_1 = T1_10rm_BS_1(:,5);
T1_10rm_BS_2 = xlsread(‘T1_70%_10rm_BS_2.xlsx’, ‘Sheet4′,’A5:E78394’);
time = T1_10rm_BS_2(:,1);
BS_RF_2 = T1_10rm_BS_2(:,4);
BS_BF_2 = T1_10rm_BS_2(:,5);
running filters;
%FILTER
% Lowpass filter 500Hz
Freq = 1500;
NyqFreq = Freq/2;
Lowpass = 500;
Wn = Lowpass/NyqFreq;
[B,A] = butter (4,Wn,’low’);
% run filter
BS_RF_1_low_filtered = filtfilt(B,A,BS_RF_1);
BS_BF_1_low_filtered = filtfilt(B,A,BS_BF_1);
BS_RF_2_low_filtered = filtfilt(B,A,BS_RF_2);
BS_BF_2_low_filtered = filtfilt(B,A,BS_BF_2);
% Highpass filter 10Hz
Highpass = 10;
Wo = Highpass/NyqFreq;
[D,C] = butter (4,Wo,’high’);
% run filter
high_pass_BS_RF_1_filtdata = filtfilt(D,C,BS_RF_1_low_filtered);
high_pass_BS_BF_1_filtdata = filtfilt(D,C,BS_BF_1_low_filtered);
high_pass_BS_RF_2_filtdata = filtfilt(D,C,BS_RF_2_low_filtered);
high_pass_BS_BF_2_filtdata = filtfilt(D,C,BS_BF_2_low_filtered);
% Full wave rectification of EMG
filtered_BS_RF_1 = abs(high_pass_BS_RF_1_filtdata);
filtered_BS_BF_1 = abs(high_pass_BS_BF_1_filtdata);
filtered_BS_RF_2 = abs(high_pass_BS_RF_2_filtdata);
filtered_BS_BF_2 = abs(high_pass_BS_BF_2_filtdata);
% Linear envelope (2Hz lowpass)
LP = 2;
Wp = LP/NyqFreq;
[F,E] = butter (4,Wp,’low’);
linear_BS_RF_1 = filtfilt(F,E,filtered_BS_RF_1);
linear_BS_BF_1 = filtfilt(F,E,filtered_BS_BF_1);
linear_BS_RF_2 = filtfilt(F,E,filtered_BS_RF_2);
linear_BS_BF_2 = filtfilt(F,E,filtered_BS_BF_2);
and generating my plot (for BS_RF_1);
I am now struggling to try and find the area under the curve for my set of 10 repetitions.
I am looking to define mean activation and peak activation.
All constructive help welcome,
Cheers. Hi all,
I am VERY new to MATLAB, so please excuse my ignorance if this is a very novice problem. I am analysing EMG signals, I have imported from my data from xlsx as time (X) and voltage (y).
After extracting data;
set(0, ‘defaultFigureToolbar’, ‘none’)
% RAW_DATA_EXTRACTION
%% BACK_SQUAT
%%% VISIT_ONE
T1_10rm_BS_1 = xlsread(‘T1_70%_10rm_BS_1.xlsx’,’Sheet3′,’A5:E77494′);
time = T1_10rm_BS_1(:,1);
BS_RF_1 = T1_10rm_BS_1(:,4);
BS_BF_1 = T1_10rm_BS_1(:,5);
T1_10rm_BS_2 = xlsread(‘T1_70%_10rm_BS_2.xlsx’, ‘Sheet4′,’A5:E78394’);
time = T1_10rm_BS_2(:,1);
BS_RF_2 = T1_10rm_BS_2(:,4);
BS_BF_2 = T1_10rm_BS_2(:,5);
running filters;
%FILTER
% Lowpass filter 500Hz
Freq = 1500;
NyqFreq = Freq/2;
Lowpass = 500;
Wn = Lowpass/NyqFreq;
[B,A] = butter (4,Wn,’low’);
% run filter
BS_RF_1_low_filtered = filtfilt(B,A,BS_RF_1);
BS_BF_1_low_filtered = filtfilt(B,A,BS_BF_1);
BS_RF_2_low_filtered = filtfilt(B,A,BS_RF_2);
BS_BF_2_low_filtered = filtfilt(B,A,BS_BF_2);
% Highpass filter 10Hz
Highpass = 10;
Wo = Highpass/NyqFreq;
[D,C] = butter (4,Wo,’high’);
% run filter
high_pass_BS_RF_1_filtdata = filtfilt(D,C,BS_RF_1_low_filtered);
high_pass_BS_BF_1_filtdata = filtfilt(D,C,BS_BF_1_low_filtered);
high_pass_BS_RF_2_filtdata = filtfilt(D,C,BS_RF_2_low_filtered);
high_pass_BS_BF_2_filtdata = filtfilt(D,C,BS_BF_2_low_filtered);
% Full wave rectification of EMG
filtered_BS_RF_1 = abs(high_pass_BS_RF_1_filtdata);
filtered_BS_BF_1 = abs(high_pass_BS_BF_1_filtdata);
filtered_BS_RF_2 = abs(high_pass_BS_RF_2_filtdata);
filtered_BS_BF_2 = abs(high_pass_BS_BF_2_filtdata);
% Linear envelope (2Hz lowpass)
LP = 2;
Wp = LP/NyqFreq;
[F,E] = butter (4,Wp,’low’);
linear_BS_RF_1 = filtfilt(F,E,filtered_BS_RF_1);
linear_BS_BF_1 = filtfilt(F,E,filtered_BS_BF_1);
linear_BS_RF_2 = filtfilt(F,E,filtered_BS_RF_2);
linear_BS_BF_2 = filtfilt(F,E,filtered_BS_BF_2);
and generating my plot (for BS_RF_1);
I am now struggling to try and find the area under the curve for my set of 10 repetitions.
I am looking to define mean activation and peak activation.
All constructive help welcome,
Cheers. onset, offset, emg, novice MATLAB Answers — New Questions
How to draw a semi-circular heat map?(怎么画半圆形的热力图)
How to draw such a heatmap, should I use ‘polaroplot’ or ‘colormap’?
怎么画出这样的热力图,应该使用polarplot还是colormap?How to draw such a heatmap, should I use ‘polaroplot’ or ‘colormap’?
怎么画出这样的热力图,应该使用polarplot还是colormap? How to draw such a heatmap, should I use ‘polaroplot’ or ‘colormap’?
怎么画出这样的热力图,应该使用polarplot还是colormap? polarplot, colormap MATLAB Answers — New Questions
Is it possible to have a string constant (output datatype: string) generated in UTF-8?
"Simulink strings support 256 characters of the ISO/IEC 8859-1 character set" according to Simulink Strings in the Help Center.
I’ve created a "string constant" block called "String ãºç" in a simulink model where I’ve inserted the value "exampleàñçº".
The block’s name is exported correctly in UTF-8 ("String ãºç"), but what I get in the exported file for the value is "example340361347272" (octal values for extended characters in ISO/IEC 8859-1).
I was wondering if there was a way to export the string value in UTF-8."Simulink strings support 256 characters of the ISO/IEC 8859-1 character set" according to Simulink Strings in the Help Center.
I’ve created a "string constant" block called "String ãºç" in a simulink model where I’ve inserted the value "exampleàñçº".
The block’s name is exported correctly in UTF-8 ("String ãºç"), but what I get in the exported file for the value is "example340361347272" (octal values for extended characters in ISO/IEC 8859-1).
I was wondering if there was a way to export the string value in UTF-8. "Simulink strings support 256 characters of the ISO/IEC 8859-1 character set" according to Simulink Strings in the Help Center.
I’ve created a "string constant" block called "String ãºç" in a simulink model where I’ve inserted the value "exampleàñçº".
The block’s name is exported correctly in UTF-8 ("String ãºç"), but what I get in the exported file for the value is "example340361347272" (octal values for extended characters in ISO/IEC 8859-1).
I was wondering if there was a way to export the string value in UTF-8. block parameters, simulink, simulink coder, matlab coder, embedded coder, string, datatype, string constant, model configuration, encoding, latin-1, iso/iec 8859-1, utf-8, unicode MATLAB Answers — New Questions
How to create the 3D mesh plot in Simulink
I am having simple m-code ..
[X,Y] = meshgrid(-8:.5:8);
R = sqrt(X.^2 + Y.^2) + eps;
Z = sin(R)./R;
figure(1)
mesh(Z)
and it displays
I want to do it in Simulink – Is it possible to do 3D mesh in SImulink – how?I am having simple m-code ..
[X,Y] = meshgrid(-8:.5:8);
R = sqrt(X.^2 + Y.^2) + eps;
Z = sin(R)./R;
figure(1)
mesh(Z)
and it displays
I want to do it in Simulink – Is it possible to do 3D mesh in SImulink – how? I am having simple m-code ..
[X,Y] = meshgrid(-8:.5:8);
R = sqrt(X.^2 + Y.^2) + eps;
Z = sin(R)./R;
figure(1)
mesh(Z)
and it displays
I want to do it in Simulink – Is it possible to do 3D mesh in SImulink – how? 3d mesh in simulink MATLAB Answers — New Questions
weboptions w/ Windows encoding
Hi
What’s the CharacterEncoding option for Windows?
Example only lists ‘US-ASCII’, ‘UTF-8’, ‘latin1’, ‘Shift_JIS’, and’ISO-8859-1′.Hi
What’s the CharacterEncoding option for Windows?
Example only lists ‘US-ASCII’, ‘UTF-8’, ‘latin1’, ‘Shift_JIS’, and’ISO-8859-1′. Hi
What’s the CharacterEncoding option for Windows?
Example only lists ‘US-ASCII’, ‘UTF-8’, ‘latin1’, ‘Shift_JIS’, and’ISO-8859-1′. matlab MATLAB Answers — New Questions
I’m doing an project in matlab software and my data set having more than 1 lakh images so can you prefer any best laptop for this processing in i7
I’m doing an project on Alzheimer’s disease with an osias dataset which is having nearly 1 lakh images so to process this i5 system is going too slow so prefer me which system is best to perform this operation and for future use. I’m planning to take i7 is this ok for working these dataset.I’m doing an project on Alzheimer’s disease with an osias dataset which is having nearly 1 lakh images so to process this i5 system is going too slow so prefer me which system is best to perform this operation and for future use. I’m planning to take i7 is this ok for working these dataset. I’m doing an project on Alzheimer’s disease with an osias dataset which is having nearly 1 lakh images so to process this i5 system is going too slow so prefer me which system is best to perform this operation and for future use. I’m planning to take i7 is this ok for working these dataset. #maltab, #new laptops, #alzhimer disease MATLAB Answers — New Questions
Error using contourf (line 57) The size of X must match the size of Z or the number of columns of Z.
clear all; clc;
% Parameters (units at m, C, W..)
L = 0.3; % Length of the cross section of square
H = 0.15; % Height of the cross section of square
T_flueinterior = 350; % Interior flue temperature
h_flueinterior = 100; % Interior flue convection coefficient
T_airexterior = 25; % Exterior air temperature
h_airexterior = 5; % Exterior air convection coefficient
k = 0.85; % Thermal conductivity of refractory brick
alpha = 5.5e-7; % Thermal diffusivity of refractory brick
Ti = 25; % Initial temperature of the wall
dx = 0.05; % Grid spacing at x direction
dy = 0.05; % Grid spacing at y direction
dt = 3600; % Time step in one hour in seconds
Nt = [1, 3, 10, 30]; % Time steps to output
% Converting Nt to seconds
Nt = Nt * 3600;
% Setting up the number of grid points in x and y direction
Nx = round(L/dx) + 1;
Ny = round(H/dy) + 1;
% Initializing temperature field
T = Ti * ones(Nx,Ny);
% Defining a coefficient to simplify the discretized FD equation
a = alpha * dt / dx^2;
% Time loop
for n = 1:max(Nt)/dt
T_new = T;
for j = 2:Ny-1
for i = 2:Nx-1
T_new(i,j) = T(i,j) + a * (T(i,j+1) + T(i,j-1) + T(i+1,j) + T(i-1,j) – 4 * T(i,j));
end
end
% Applying boundary conditions
% Left side flue
T_new(:,1) = T_flueinterior + (T(:,2) – T_flueinterior) / (1 + dx * h_flueinterior / k);
% Right side air
T_new(:,end) = T_airexterior + (T(:,end-1) – T_airexterior) / (1 + dx * h_airexterior / k);
% Top and bottom boundaries (adiabatic and insulated)
T_new(1,:) = T_new(2,:);
T_new(end, 🙂 = T_new(end-1, :);
% Re-update temperature field
T = T_new;
% Save results at the specified times
if any(n*dt == Nt)
figure;
contourf(linspace(0,L,Nx), linspace(0,H,Ny), T, 20, ‘LineColor’, ‘none’);
colorbar;
title([‘Temperature Distribution at t = ‘, num2str(n*dt/3600), ‘hours’]);
xlabel(‘x (in m)’);
ylabel(‘y (in m)’);
end
end
% The final temperature distribution
figure;
contourf(linspace(0,L,Nx), linspace(0,H,Ny), T, 20, ‘LineColor’, ‘none’);
title(‘Final Temperature Distribution’);
xlabel(‘x (in m)’);
ylabel(‘y (in m)’);
I am having the error in the title and am not sure how to fix it? I could use some helpclear all; clc;
% Parameters (units at m, C, W..)
L = 0.3; % Length of the cross section of square
H = 0.15; % Height of the cross section of square
T_flueinterior = 350; % Interior flue temperature
h_flueinterior = 100; % Interior flue convection coefficient
T_airexterior = 25; % Exterior air temperature
h_airexterior = 5; % Exterior air convection coefficient
k = 0.85; % Thermal conductivity of refractory brick
alpha = 5.5e-7; % Thermal diffusivity of refractory brick
Ti = 25; % Initial temperature of the wall
dx = 0.05; % Grid spacing at x direction
dy = 0.05; % Grid spacing at y direction
dt = 3600; % Time step in one hour in seconds
Nt = [1, 3, 10, 30]; % Time steps to output
% Converting Nt to seconds
Nt = Nt * 3600;
% Setting up the number of grid points in x and y direction
Nx = round(L/dx) + 1;
Ny = round(H/dy) + 1;
% Initializing temperature field
T = Ti * ones(Nx,Ny);
% Defining a coefficient to simplify the discretized FD equation
a = alpha * dt / dx^2;
% Time loop
for n = 1:max(Nt)/dt
T_new = T;
for j = 2:Ny-1
for i = 2:Nx-1
T_new(i,j) = T(i,j) + a * (T(i,j+1) + T(i,j-1) + T(i+1,j) + T(i-1,j) – 4 * T(i,j));
end
end
% Applying boundary conditions
% Left side flue
T_new(:,1) = T_flueinterior + (T(:,2) – T_flueinterior) / (1 + dx * h_flueinterior / k);
% Right side air
T_new(:,end) = T_airexterior + (T(:,end-1) – T_airexterior) / (1 + dx * h_airexterior / k);
% Top and bottom boundaries (adiabatic and insulated)
T_new(1,:) = T_new(2,:);
T_new(end, 🙂 = T_new(end-1, :);
% Re-update temperature field
T = T_new;
% Save results at the specified times
if any(n*dt == Nt)
figure;
contourf(linspace(0,L,Nx), linspace(0,H,Ny), T, 20, ‘LineColor’, ‘none’);
colorbar;
title([‘Temperature Distribution at t = ‘, num2str(n*dt/3600), ‘hours’]);
xlabel(‘x (in m)’);
ylabel(‘y (in m)’);
end
end
% The final temperature distribution
figure;
contourf(linspace(0,L,Nx), linspace(0,H,Ny), T, 20, ‘LineColor’, ‘none’);
title(‘Final Temperature Distribution’);
xlabel(‘x (in m)’);
ylabel(‘y (in m)’);
I am having the error in the title and am not sure how to fix it? I could use some help clear all; clc;
% Parameters (units at m, C, W..)
L = 0.3; % Length of the cross section of square
H = 0.15; % Height of the cross section of square
T_flueinterior = 350; % Interior flue temperature
h_flueinterior = 100; % Interior flue convection coefficient
T_airexterior = 25; % Exterior air temperature
h_airexterior = 5; % Exterior air convection coefficient
k = 0.85; % Thermal conductivity of refractory brick
alpha = 5.5e-7; % Thermal diffusivity of refractory brick
Ti = 25; % Initial temperature of the wall
dx = 0.05; % Grid spacing at x direction
dy = 0.05; % Grid spacing at y direction
dt = 3600; % Time step in one hour in seconds
Nt = [1, 3, 10, 30]; % Time steps to output
% Converting Nt to seconds
Nt = Nt * 3600;
% Setting up the number of grid points in x and y direction
Nx = round(L/dx) + 1;
Ny = round(H/dy) + 1;
% Initializing temperature field
T = Ti * ones(Nx,Ny);
% Defining a coefficient to simplify the discretized FD equation
a = alpha * dt / dx^2;
% Time loop
for n = 1:max(Nt)/dt
T_new = T;
for j = 2:Ny-1
for i = 2:Nx-1
T_new(i,j) = T(i,j) + a * (T(i,j+1) + T(i,j-1) + T(i+1,j) + T(i-1,j) – 4 * T(i,j));
end
end
% Applying boundary conditions
% Left side flue
T_new(:,1) = T_flueinterior + (T(:,2) – T_flueinterior) / (1 + dx * h_flueinterior / k);
% Right side air
T_new(:,end) = T_airexterior + (T(:,end-1) – T_airexterior) / (1 + dx * h_airexterior / k);
% Top and bottom boundaries (adiabatic and insulated)
T_new(1,:) = T_new(2,:);
T_new(end, 🙂 = T_new(end-1, :);
% Re-update temperature field
T = T_new;
% Save results at the specified times
if any(n*dt == Nt)
figure;
contourf(linspace(0,L,Nx), linspace(0,H,Ny), T, 20, ‘LineColor’, ‘none’);
colorbar;
title([‘Temperature Distribution at t = ‘, num2str(n*dt/3600), ‘hours’]);
xlabel(‘x (in m)’);
ylabel(‘y (in m)’);
end
end
% The final temperature distribution
figure;
contourf(linspace(0,L,Nx), linspace(0,H,Ny), T, 20, ‘LineColor’, ‘none’);
title(‘Final Temperature Distribution’);
xlabel(‘x (in m)’);
ylabel(‘y (in m)’);
I am having the error in the title and am not sure how to fix it? I could use some help optimization, error, matlab MATLAB Answers — New Questions
How to perform Uplink synchronization between two SDR’s?
Dear Everyone,
I would like to generate an Uplink signal using 5G toolbox from a UE (SDR-1) and send it to the gNB (SDR-2). How should I synchronize in between the two SDR’s so that the signal is recevied properly on the gNB side. Do I have to add preamble from the UE? As I did not find anything related to it. Thank you.Dear Everyone,
I would like to generate an Uplink signal using 5G toolbox from a UE (SDR-1) and send it to the gNB (SDR-2). How should I synchronize in between the two SDR’s so that the signal is recevied properly on the gNB side. Do I have to add preamble from the UE? As I did not find anything related to it. Thank you. Dear Everyone,
I would like to generate an Uplink signal using 5G toolbox from a UE (SDR-1) and send it to the gNB (SDR-2). How should I synchronize in between the two SDR’s so that the signal is recevied properly on the gNB side. Do I have to add preamble from the UE? As I did not find anything related to it. Thank you. 5g toolbox, 5g nr uplink MATLAB Answers — New Questions