Tag Archives: matlab
How to speed up this code?
Dear All,
I have this code but when I run it it takes one day and half to run. How can I speed it up?
Best regards,
Ara
clear all;
clc;
% Set the folder path where the .nc files are located
folderPath = ‘E:datapodtc_apr’;
% Get a list of all NetCDF files in the folder
fileList = dir(fullfile(folderPath, ‘*_nc’));
numFiles = numel(fileList);
% Initialize the data structures
data = cell(numFiles, 1);
dateArray = [];
typeArray = cell(numFiles, 1);
% Loop through each file
for fileIndex = 1:numFiles
% Read the NetCDF file
filePath = fullfile(folderPath, fileList(fileIndex).name);
% Read the data from the NetCDF file
ncinfo_struct = ncinfo(filePath);
if isfield(ncinfo_struct, ‘Variables’)
variable_names = {ncinfo_struct.Variables.Name};
if all(ismember({‘time’, ‘TEC’, ‘S4’, ‘RFI’, ‘elevation’, ‘occheight’, ‘caL1_SNR’, ‘pL2_SNR’, ‘x_LEO’, ‘y_LEO’, ‘z_LEO’, ‘x_GPS’, ‘y_GPS’, ‘z_GPS’}, variable_names))
% Extract the date and type information from the file name
[~, filename, ~] = fileparts(fileList(fileIndex).name);
dateStr = regexp(filename, ‘d{4}.d{3}’, ‘match’, ‘once’);
typeStr = regexp(filename, ‘Gd{2}|Rd{2}’, ‘match’, ‘once’);
% Read the data from the NetCDF file
data{fileIndex}.time = ncread(filePath, ‘time’);
data{fileIndex}.TEC = ncread(filePath, ‘TEC’);
data{fileIndex}.S4 = ncread(filePath, ‘S4’);
data{fileIndex}.RFI = ncread(filePath, ‘RFI’);
data{fileIndex}.elevation = ncread(filePath, ‘elevation’);
data{fileIndex}.occheight = ncread(filePath, ‘occheight’);
data{fileIndex}.caL1_SNR = ncread(filePath, ‘caL1_SNR’);
data{fileIndex}.pL2_SNR = ncread(filePath, ‘pL2_SNR’);
data{fileIndex}.x_LEO = ncread(filePath, ‘x_LEO’);
data{fileIndex}.y_LEO = ncread(filePath, ‘y_LEO’);
data{fileIndex}.z_LEO = ncread(filePath, ‘z_LEO’);
data{fileIndex}.x_GPS = ncread(filePath, ‘x_GPS’);
data{fileIndex}.y_GPS = ncread(filePath, ‘y_GPS’);
data{fileIndex}.z_GPS = ncread(filePath, ‘z_GPS’);
% Store the date and type information
dateArray = [dateArray, str2double(dateStr)];
typeArray{fileIndex} = typeStr;
else
% Skip this file and move on to the next one
fprintf(‘File "%s" does not contain all the required variables. Skipping this file.n’, fileList(fileIndex).name);
continue;
end
else
% Skip this file and move on to the next one
fprintf(‘File "%s" does not contain the "Variables" field. Skipping this file.n’, fileList(fileIndex).name);
continue;
end
end
% Sort the data by date and type
[sortedDates, sortedIndices] = sort(dateArray);
sortedTypes = cellfun(@(x) x, typeArray(sortedIndices), ‘UniformOutput’, false);
% Create the sorted data structures
sortedData = cell(numFiles, 1);
for i = 1:numFiles
sortedData{i} = data{sortedIndices(i)};
end
% Save the data to a .mat file
save(‘podtc_apr_data.mat’, ‘sortedData’, ‘sortedDates’, ‘sortedTypes’);Dear All,
I have this code but when I run it it takes one day and half to run. How can I speed it up?
Best regards,
Ara
clear all;
clc;
% Set the folder path where the .nc files are located
folderPath = ‘E:datapodtc_apr’;
% Get a list of all NetCDF files in the folder
fileList = dir(fullfile(folderPath, ‘*_nc’));
numFiles = numel(fileList);
% Initialize the data structures
data = cell(numFiles, 1);
dateArray = [];
typeArray = cell(numFiles, 1);
% Loop through each file
for fileIndex = 1:numFiles
% Read the NetCDF file
filePath = fullfile(folderPath, fileList(fileIndex).name);
% Read the data from the NetCDF file
ncinfo_struct = ncinfo(filePath);
if isfield(ncinfo_struct, ‘Variables’)
variable_names = {ncinfo_struct.Variables.Name};
if all(ismember({‘time’, ‘TEC’, ‘S4’, ‘RFI’, ‘elevation’, ‘occheight’, ‘caL1_SNR’, ‘pL2_SNR’, ‘x_LEO’, ‘y_LEO’, ‘z_LEO’, ‘x_GPS’, ‘y_GPS’, ‘z_GPS’}, variable_names))
% Extract the date and type information from the file name
[~, filename, ~] = fileparts(fileList(fileIndex).name);
dateStr = regexp(filename, ‘d{4}.d{3}’, ‘match’, ‘once’);
typeStr = regexp(filename, ‘Gd{2}|Rd{2}’, ‘match’, ‘once’);
% Read the data from the NetCDF file
data{fileIndex}.time = ncread(filePath, ‘time’);
data{fileIndex}.TEC = ncread(filePath, ‘TEC’);
data{fileIndex}.S4 = ncread(filePath, ‘S4’);
data{fileIndex}.RFI = ncread(filePath, ‘RFI’);
data{fileIndex}.elevation = ncread(filePath, ‘elevation’);
data{fileIndex}.occheight = ncread(filePath, ‘occheight’);
data{fileIndex}.caL1_SNR = ncread(filePath, ‘caL1_SNR’);
data{fileIndex}.pL2_SNR = ncread(filePath, ‘pL2_SNR’);
data{fileIndex}.x_LEO = ncread(filePath, ‘x_LEO’);
data{fileIndex}.y_LEO = ncread(filePath, ‘y_LEO’);
data{fileIndex}.z_LEO = ncread(filePath, ‘z_LEO’);
data{fileIndex}.x_GPS = ncread(filePath, ‘x_GPS’);
data{fileIndex}.y_GPS = ncread(filePath, ‘y_GPS’);
data{fileIndex}.z_GPS = ncread(filePath, ‘z_GPS’);
% Store the date and type information
dateArray = [dateArray, str2double(dateStr)];
typeArray{fileIndex} = typeStr;
else
% Skip this file and move on to the next one
fprintf(‘File "%s" does not contain all the required variables. Skipping this file.n’, fileList(fileIndex).name);
continue;
end
else
% Skip this file and move on to the next one
fprintf(‘File "%s" does not contain the "Variables" field. Skipping this file.n’, fileList(fileIndex).name);
continue;
end
end
% Sort the data by date and type
[sortedDates, sortedIndices] = sort(dateArray);
sortedTypes = cellfun(@(x) x, typeArray(sortedIndices), ‘UniformOutput’, false);
% Create the sorted data structures
sortedData = cell(numFiles, 1);
for i = 1:numFiles
sortedData{i} = data{sortedIndices(i)};
end
% Save the data to a .mat file
save(‘podtc_apr_data.mat’, ‘sortedData’, ‘sortedDates’, ‘sortedTypes’); Dear All,
I have this code but when I run it it takes one day and half to run. How can I speed it up?
Best regards,
Ara
clear all;
clc;
% Set the folder path where the .nc files are located
folderPath = ‘E:datapodtc_apr’;
% Get a list of all NetCDF files in the folder
fileList = dir(fullfile(folderPath, ‘*_nc’));
numFiles = numel(fileList);
% Initialize the data structures
data = cell(numFiles, 1);
dateArray = [];
typeArray = cell(numFiles, 1);
% Loop through each file
for fileIndex = 1:numFiles
% Read the NetCDF file
filePath = fullfile(folderPath, fileList(fileIndex).name);
% Read the data from the NetCDF file
ncinfo_struct = ncinfo(filePath);
if isfield(ncinfo_struct, ‘Variables’)
variable_names = {ncinfo_struct.Variables.Name};
if all(ismember({‘time’, ‘TEC’, ‘S4’, ‘RFI’, ‘elevation’, ‘occheight’, ‘caL1_SNR’, ‘pL2_SNR’, ‘x_LEO’, ‘y_LEO’, ‘z_LEO’, ‘x_GPS’, ‘y_GPS’, ‘z_GPS’}, variable_names))
% Extract the date and type information from the file name
[~, filename, ~] = fileparts(fileList(fileIndex).name);
dateStr = regexp(filename, ‘d{4}.d{3}’, ‘match’, ‘once’);
typeStr = regexp(filename, ‘Gd{2}|Rd{2}’, ‘match’, ‘once’);
% Read the data from the NetCDF file
data{fileIndex}.time = ncread(filePath, ‘time’);
data{fileIndex}.TEC = ncread(filePath, ‘TEC’);
data{fileIndex}.S4 = ncread(filePath, ‘S4’);
data{fileIndex}.RFI = ncread(filePath, ‘RFI’);
data{fileIndex}.elevation = ncread(filePath, ‘elevation’);
data{fileIndex}.occheight = ncread(filePath, ‘occheight’);
data{fileIndex}.caL1_SNR = ncread(filePath, ‘caL1_SNR’);
data{fileIndex}.pL2_SNR = ncread(filePath, ‘pL2_SNR’);
data{fileIndex}.x_LEO = ncread(filePath, ‘x_LEO’);
data{fileIndex}.y_LEO = ncread(filePath, ‘y_LEO’);
data{fileIndex}.z_LEO = ncread(filePath, ‘z_LEO’);
data{fileIndex}.x_GPS = ncread(filePath, ‘x_GPS’);
data{fileIndex}.y_GPS = ncread(filePath, ‘y_GPS’);
data{fileIndex}.z_GPS = ncread(filePath, ‘z_GPS’);
% Store the date and type information
dateArray = [dateArray, str2double(dateStr)];
typeArray{fileIndex} = typeStr;
else
% Skip this file and move on to the next one
fprintf(‘File "%s" does not contain all the required variables. Skipping this file.n’, fileList(fileIndex).name);
continue;
end
else
% Skip this file and move on to the next one
fprintf(‘File "%s" does not contain the "Variables" field. Skipping this file.n’, fileList(fileIndex).name);
continue;
end
end
% Sort the data by date and type
[sortedDates, sortedIndices] = sort(dateArray);
sortedTypes = cellfun(@(x) x, typeArray(sortedIndices), ‘UniformOutput’, false);
% Create the sorted data structures
sortedData = cell(numFiles, 1);
for i = 1:numFiles
sortedData{i} = data{sortedIndices(i)};
end
% Save the data to a .mat file
save(‘podtc_apr_data.mat’, ‘sortedData’, ‘sortedDates’, ‘sortedTypes’); speed up run MATLAB Answers — New Questions
what is the meaning of error “Index exceeds the number of array elements. Index must not exceed 8.”
function Latest_code
%%
for Pr = [1 2 3 4]
% Define constants
M = 0.5;
B = 0.5;
C = 0.2;
A = 0.1;
R = 0.1;
E = 0.1;
K = 0.5;
W = 0.1;
N = 0.1;
T = 0.3;
Q = 0.1;
N1 = 0.1;
N2 = 0.1;
S = 0.1;
d = 1;
m = 1;
E1 = 0.1;
Pe = 0.2;
S1 = 0.1;
M1 = 0.1;
O = 0.1;
% Defining initial guess
solinit = bvpinit(linspace(0, 2, 50), [1 0 0 0 1 1 1 1]);
% Solving the boundary value problem
sol = bvp4c(@bvpexam2, @bcexam2, solinit);
x = sol.x;
y = sol.y;
% Plotting of velocity
figure(2)
plot(x, y(4, :), ‘linewidth’, 1)
ylim([0 1])
xlim([0 9])
hold on
xlabel(‘eta’, ‘FontWeight’, ‘bold’, ‘FontSize’, 16)
ylabel(‘theta(eta)’, ‘FontWeight’, ‘bold’, ‘FontSize’, 16)
end
% System of first order ODEs
function yvector = bvpexam2(~, y)
yy1 = (1 / (1 + (1 / B))) * (2 * y(2) * y(2) – y(1) * y(3) + M * y(2) + K * y(2) – W * (y(4) – N * y(7) – R));
yy2 = -(1 / (((4 / 3) * R) + 1)) * (Pr * (y(1) * y(5) – y(2) * y(4) – T * y(2) + Q * y(4) + N1 * y(8) * y(5) + N2 * y(5) * y(5) + E * (1 + 1 / B) * y(3) * y(3) + (E * y(2) * y(2)) * (M + K)));
yy3 = -S * (y(1) * y(8) – y(2) * y(7) – C * y(2) – A * y(7) * ((1 + d * y(4)) ^ m) * exp(-E1 / (1 + d * y(4)))) – (N2 / N1) * y(6);
yy4 = (Pe * (y(9) * (y(10) + O) + y(11) * y(8))) – (S1 * ((y(11) – y(10)) * y(1))) + (S1 * M1 * y(2));
yvector = [ y(2); y(3); yy1; y(5); yy2; y(8); yy3 ; y(11); yy4 ];
end
% Residual of the boundary conditions
function residual = bcexam2(y0, yinf)
residual = [y0(1)- 1; y0(2) – 1; y0(4) – 1 +T ; y0(7)-1+C ; y0(10)-1+M1 ; yinf(2); yinf(4); yinf(7); yinf(10)];
end
endfunction Latest_code
%%
for Pr = [1 2 3 4]
% Define constants
M = 0.5;
B = 0.5;
C = 0.2;
A = 0.1;
R = 0.1;
E = 0.1;
K = 0.5;
W = 0.1;
N = 0.1;
T = 0.3;
Q = 0.1;
N1 = 0.1;
N2 = 0.1;
S = 0.1;
d = 1;
m = 1;
E1 = 0.1;
Pe = 0.2;
S1 = 0.1;
M1 = 0.1;
O = 0.1;
% Defining initial guess
solinit = bvpinit(linspace(0, 2, 50), [1 0 0 0 1 1 1 1]);
% Solving the boundary value problem
sol = bvp4c(@bvpexam2, @bcexam2, solinit);
x = sol.x;
y = sol.y;
% Plotting of velocity
figure(2)
plot(x, y(4, :), ‘linewidth’, 1)
ylim([0 1])
xlim([0 9])
hold on
xlabel(‘eta’, ‘FontWeight’, ‘bold’, ‘FontSize’, 16)
ylabel(‘theta(eta)’, ‘FontWeight’, ‘bold’, ‘FontSize’, 16)
end
% System of first order ODEs
function yvector = bvpexam2(~, y)
yy1 = (1 / (1 + (1 / B))) * (2 * y(2) * y(2) – y(1) * y(3) + M * y(2) + K * y(2) – W * (y(4) – N * y(7) – R));
yy2 = -(1 / (((4 / 3) * R) + 1)) * (Pr * (y(1) * y(5) – y(2) * y(4) – T * y(2) + Q * y(4) + N1 * y(8) * y(5) + N2 * y(5) * y(5) + E * (1 + 1 / B) * y(3) * y(3) + (E * y(2) * y(2)) * (M + K)));
yy3 = -S * (y(1) * y(8) – y(2) * y(7) – C * y(2) – A * y(7) * ((1 + d * y(4)) ^ m) * exp(-E1 / (1 + d * y(4)))) – (N2 / N1) * y(6);
yy4 = (Pe * (y(9) * (y(10) + O) + y(11) * y(8))) – (S1 * ((y(11) – y(10)) * y(1))) + (S1 * M1 * y(2));
yvector = [ y(2); y(3); yy1; y(5); yy2; y(8); yy3 ; y(11); yy4 ];
end
% Residual of the boundary conditions
function residual = bcexam2(y0, yinf)
residual = [y0(1)- 1; y0(2) – 1; y0(4) – 1 +T ; y0(7)-1+C ; y0(10)-1+M1 ; yinf(2); yinf(4); yinf(7); yinf(10)];
end
end function Latest_code
%%
for Pr = [1 2 3 4]
% Define constants
M = 0.5;
B = 0.5;
C = 0.2;
A = 0.1;
R = 0.1;
E = 0.1;
K = 0.5;
W = 0.1;
N = 0.1;
T = 0.3;
Q = 0.1;
N1 = 0.1;
N2 = 0.1;
S = 0.1;
d = 1;
m = 1;
E1 = 0.1;
Pe = 0.2;
S1 = 0.1;
M1 = 0.1;
O = 0.1;
% Defining initial guess
solinit = bvpinit(linspace(0, 2, 50), [1 0 0 0 1 1 1 1]);
% Solving the boundary value problem
sol = bvp4c(@bvpexam2, @bcexam2, solinit);
x = sol.x;
y = sol.y;
% Plotting of velocity
figure(2)
plot(x, y(4, :), ‘linewidth’, 1)
ylim([0 1])
xlim([0 9])
hold on
xlabel(‘eta’, ‘FontWeight’, ‘bold’, ‘FontSize’, 16)
ylabel(‘theta(eta)’, ‘FontWeight’, ‘bold’, ‘FontSize’, 16)
end
% System of first order ODEs
function yvector = bvpexam2(~, y)
yy1 = (1 / (1 + (1 / B))) * (2 * y(2) * y(2) – y(1) * y(3) + M * y(2) + K * y(2) – W * (y(4) – N * y(7) – R));
yy2 = -(1 / (((4 / 3) * R) + 1)) * (Pr * (y(1) * y(5) – y(2) * y(4) – T * y(2) + Q * y(4) + N1 * y(8) * y(5) + N2 * y(5) * y(5) + E * (1 + 1 / B) * y(3) * y(3) + (E * y(2) * y(2)) * (M + K)));
yy3 = -S * (y(1) * y(8) – y(2) * y(7) – C * y(2) – A * y(7) * ((1 + d * y(4)) ^ m) * exp(-E1 / (1 + d * y(4)))) – (N2 / N1) * y(6);
yy4 = (Pe * (y(9) * (y(10) + O) + y(11) * y(8))) – (S1 * ((y(11) – y(10)) * y(1))) + (S1 * M1 * y(2));
yvector = [ y(2); y(3); yy1; y(5); yy2; y(8); yy3 ; y(11); yy4 ];
end
% Residual of the boundary conditions
function residual = bcexam2(y0, yinf)
residual = [y0(1)- 1; y0(2) – 1; y0(4) – 1 +T ; y0(7)-1+C ; y0(10)-1+M1 ; yinf(2); yinf(4); yinf(7); yinf(10)];
end
end fluidmechanics MATLAB Answers — New Questions
Why do I receive an error that the MATLAB Parallel Server service launched but failed to start?
When starting mjs on my Windows machine I receive the following error:
wrapper | Starting the MATLAB Parallel Server service…
wrapper | The MATLAB Parallel Server service was launched, but failed to start.
When checking the mjs-service.log I see the following error:
ERROR: STATUS | wrapper | 2014/07/22 09:47:17 | Starting the MATLAB Parallel Server service…
STATUS | wrapper | 2014/07/22 09:47:17 | –> Wrapper Started as Service
STATUS | wrapper | 2014/07/22 09:47:17 | Launching a JVM…
INFO | jvm 1 | 2014/07/22 09:47:17 | Wrapper (Version 3.1.2) https://wrapper.tanukisoftware.com/doc/english/home.html
INFO | jvm 1 | 2014/07/22 09:47:17 |
INFO | jvm 1 | 2014/07/22 09:47:17 |
INFO | jvm 1 | 2014/07/22 09:47:17 | WrapperSimpleApp: Encountered an error running main: java.lang.NoClassDefFoundError: com/mathworks/resource_core/BaseMsgID
INFO | jvm 1 | 2014/07/22 09:47:17 | java.lang.NoClassDefFoundError: com/mathworks/resource_core/BaseMsgID
INFO | jvm 1 | 2014/07/22 09:47:17 | at com.mathworks.toolbox.distcomp.control.Starter.run(Starter.java:241)
INFO | jvm 1 | 2014/07/22 09:47:17 | at com.mathworks.toolbox.distcomp.control.Starter.main(Starter.java:723)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.reflect.Method.invoke(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at org.tanukisoftware.wrapper.WrapperSimpleApp.run(WrapperSimpleApp.java:136)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.Thread.run(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | Caused by: java.lang.ClassNotFoundException: com.mathworks.resource_core.BaseMsgID
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.net.URLClassLoader$1.run(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.net.URLClassLoader$1.run(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.security.AccessController.doPrivileged(Native Method)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.net.URLClassLoader.findClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.ClassLoader.loadClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.misc.Launcher$AppClassLoader.loadClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.ClassLoader.loadClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | … 8 more
STATUS | wrapper | 2014/07/22 09:47:19 | <– Wrapper Stopped
ERROR | wrapper | 2014/07/22 09:47:21 | The MATLAB Parallel Server service was launched, but failed to start.When starting mjs on my Windows machine I receive the following error:
wrapper | Starting the MATLAB Parallel Server service…
wrapper | The MATLAB Parallel Server service was launched, but failed to start.
When checking the mjs-service.log I see the following error:
ERROR: STATUS | wrapper | 2014/07/22 09:47:17 | Starting the MATLAB Parallel Server service…
STATUS | wrapper | 2014/07/22 09:47:17 | –> Wrapper Started as Service
STATUS | wrapper | 2014/07/22 09:47:17 | Launching a JVM…
INFO | jvm 1 | 2014/07/22 09:47:17 | Wrapper (Version 3.1.2) https://wrapper.tanukisoftware.com/doc/english/home.html
INFO | jvm 1 | 2014/07/22 09:47:17 |
INFO | jvm 1 | 2014/07/22 09:47:17 |
INFO | jvm 1 | 2014/07/22 09:47:17 | WrapperSimpleApp: Encountered an error running main: java.lang.NoClassDefFoundError: com/mathworks/resource_core/BaseMsgID
INFO | jvm 1 | 2014/07/22 09:47:17 | java.lang.NoClassDefFoundError: com/mathworks/resource_core/BaseMsgID
INFO | jvm 1 | 2014/07/22 09:47:17 | at com.mathworks.toolbox.distcomp.control.Starter.run(Starter.java:241)
INFO | jvm 1 | 2014/07/22 09:47:17 | at com.mathworks.toolbox.distcomp.control.Starter.main(Starter.java:723)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.reflect.Method.invoke(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at org.tanukisoftware.wrapper.WrapperSimpleApp.run(WrapperSimpleApp.java:136)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.Thread.run(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | Caused by: java.lang.ClassNotFoundException: com.mathworks.resource_core.BaseMsgID
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.net.URLClassLoader$1.run(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.net.URLClassLoader$1.run(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.security.AccessController.doPrivileged(Native Method)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.net.URLClassLoader.findClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.ClassLoader.loadClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.misc.Launcher$AppClassLoader.loadClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.ClassLoader.loadClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | … 8 more
STATUS | wrapper | 2014/07/22 09:47:19 | <– Wrapper Stopped
ERROR | wrapper | 2014/07/22 09:47:21 | The MATLAB Parallel Server service was launched, but failed to start. When starting mjs on my Windows machine I receive the following error:
wrapper | Starting the MATLAB Parallel Server service…
wrapper | The MATLAB Parallel Server service was launched, but failed to start.
When checking the mjs-service.log I see the following error:
ERROR: STATUS | wrapper | 2014/07/22 09:47:17 | Starting the MATLAB Parallel Server service…
STATUS | wrapper | 2014/07/22 09:47:17 | –> Wrapper Started as Service
STATUS | wrapper | 2014/07/22 09:47:17 | Launching a JVM…
INFO | jvm 1 | 2014/07/22 09:47:17 | Wrapper (Version 3.1.2) https://wrapper.tanukisoftware.com/doc/english/home.html
INFO | jvm 1 | 2014/07/22 09:47:17 |
INFO | jvm 1 | 2014/07/22 09:47:17 |
INFO | jvm 1 | 2014/07/22 09:47:17 | WrapperSimpleApp: Encountered an error running main: java.lang.NoClassDefFoundError: com/mathworks/resource_core/BaseMsgID
INFO | jvm 1 | 2014/07/22 09:47:17 | java.lang.NoClassDefFoundError: com/mathworks/resource_core/BaseMsgID
INFO | jvm 1 | 2014/07/22 09:47:17 | at com.mathworks.toolbox.distcomp.control.Starter.run(Starter.java:241)
INFO | jvm 1 | 2014/07/22 09:47:17 | at com.mathworks.toolbox.distcomp.control.Starter.main(Starter.java:723)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.reflect.Method.invoke(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at org.tanukisoftware.wrapper.WrapperSimpleApp.run(WrapperSimpleApp.java:136)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.Thread.run(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | Caused by: java.lang.ClassNotFoundException: com.mathworks.resource_core.BaseMsgID
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.net.URLClassLoader$1.run(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.net.URLClassLoader$1.run(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.security.AccessController.doPrivileged(Native Method)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.net.URLClassLoader.findClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.ClassLoader.loadClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at sun.misc.Launcher$AppClassLoader.loadClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | at java.lang.ClassLoader.loadClass(Unknown Source)
INFO | jvm 1 | 2014/07/22 09:47:17 | … 8 more
STATUS | wrapper | 2014/07/22 09:47:19 | <– Wrapper Stopped
ERROR | wrapper | 2014/07/22 09:47:21 | The MATLAB Parallel Server service was launched, but failed to start. MATLAB Answers — New Questions
Cannot send PWM signals over telemetry with Connected IO to Pixhawk
I am currently able to send PWM signals over USB to my Holybro Pixhawk 6x to run my drone’s motors in Connected IO mode from Simulink. However, I need to run it while communicating over telemetry. I have a pair of telemetry modules (successfully paired to each other). One is on the Pixhawk /dev/ttyS1 port, and the other is on a Windows host computer COM5 port.
When I look in “Configuration Parameters -> Hardware Implementation -> Hardware board settings -> Target hardware resources” -> "Connected IO", I am unable to change the “Hardware board serial port” from “/dev/ttyACM0” (the Pixhawk-side USB port) to “/dev/ttyS1”, (the Pixhawk-side telemetry unit).I am currently able to send PWM signals over USB to my Holybro Pixhawk 6x to run my drone’s motors in Connected IO mode from Simulink. However, I need to run it while communicating over telemetry. I have a pair of telemetry modules (successfully paired to each other). One is on the Pixhawk /dev/ttyS1 port, and the other is on a Windows host computer COM5 port.
When I look in “Configuration Parameters -> Hardware Implementation -> Hardware board settings -> Target hardware resources” -> "Connected IO", I am unable to change the “Hardware board serial port” from “/dev/ttyACM0” (the Pixhawk-side USB port) to “/dev/ttyS1”, (the Pixhawk-side telemetry unit). I am currently able to send PWM signals over USB to my Holybro Pixhawk 6x to run my drone’s motors in Connected IO mode from Simulink. However, I need to run it while communicating over telemetry. I have a pair of telemetry modules (successfully paired to each other). One is on the Pixhawk /dev/ttyS1 port, and the other is on a Windows host computer COM5 port.
When I look in “Configuration Parameters -> Hardware Implementation -> Hardware board settings -> Target hardware resources” -> "Connected IO", I am unable to change the “Hardware board serial port” from “/dev/ttyACM0” (the Pixhawk-side USB port) to “/dev/ttyS1”, (the Pixhawk-side telemetry unit). px4, pixhawk, connected, io, mode, external MATLAB Answers — New Questions
partial differential equation solve
hi.i need help to write this equation and solve it in matlab.thanks.hi.i need help to write this equation and solve it in matlab.thanks. hi.i need help to write this equation and solve it in matlab.thanks. partial equation MATLAB Answers — New Questions
PV Array unable to connect with Blue N-Channel MOSFET
Object : I trying to simulate a MOSFET Based IV Tracer for PV Module.
Problem :
There have 2 type of MOSFET, first is the Black MOSFET which is under a same category with PV Array and only can adjust a very few parameter. While the second is Blue MOSFET which isn’t under a same category with PV Array but adjust much more parameter that from a real MOSFET datasheet, and the model converted from LTSpice lib also under this category.
Black MOSFET is able to connect with PV Array but Blue MOSFET cant. The reason i sent this email is because the Black MOSFET can’t simualte a normal MOSFET Characteristic IV Curve but Blue one can.
Question / Ask :
Is it possible to connect the N-Channel Mosfet(Blue) with PV-Array ? (I tried both PS-Simulink Converter)Object : I trying to simulate a MOSFET Based IV Tracer for PV Module.
Problem :
There have 2 type of MOSFET, first is the Black MOSFET which is under a same category with PV Array and only can adjust a very few parameter. While the second is Blue MOSFET which isn’t under a same category with PV Array but adjust much more parameter that from a real MOSFET datasheet, and the model converted from LTSpice lib also under this category.
Black MOSFET is able to connect with PV Array but Blue MOSFET cant. The reason i sent this email is because the Black MOSFET can’t simualte a normal MOSFET Characteristic IV Curve but Blue one can.
Question / Ask :
Is it possible to connect the N-Channel Mosfet(Blue) with PV-Array ? (I tried both PS-Simulink Converter) Object : I trying to simulate a MOSFET Based IV Tracer for PV Module.
Problem :
There have 2 type of MOSFET, first is the Black MOSFET which is under a same category with PV Array and only can adjust a very few parameter. While the second is Blue MOSFET which isn’t under a same category with PV Array but adjust much more parameter that from a real MOSFET datasheet, and the model converted from LTSpice lib also under this category.
Black MOSFET is able to connect with PV Array but Blue MOSFET cant. The reason i sent this email is because the Black MOSFET can’t simualte a normal MOSFET Characteristic IV Curve but Blue one can.
Question / Ask :
Is it possible to connect the N-Channel Mosfet(Blue) with PV-Array ? (I tried both PS-Simulink Converter) simulink, simscape, pv array, mosfet, connection, problem, matlab MATLAB Answers — New Questions
How to set specific colors in margin of errorbars line?
I use errorbar funtion in order to put errorbars in a plot. I would like to set a different color to the margin of the line.
I use
errorbar(x,y,z,’r-‘,’LineWidth’,3)
I would like for example to have red line with black margins..
Could you help me?I use errorbar funtion in order to put errorbars in a plot. I would like to set a different color to the margin of the line.
I use
errorbar(x,y,z,’r-‘,’LineWidth’,3)
I would like for example to have red line with black margins..
Could you help me? I use errorbar funtion in order to put errorbars in a plot. I would like to set a different color to the margin of the line.
I use
errorbar(x,y,z,’r-‘,’LineWidth’,3)
I would like for example to have red line with black margins..
Could you help me? errorbars, colors, set, matlab MATLAB Answers — New Questions
how can I detect and fill gapes in data
hello everyone.. im working with snow data in which present in one colon time and in other colon other parameter… and I have some missing value with pass of 7 minutes..now I would like to fill this gaps of 7 minutes with the interval of time so the last interval of time after this gaps… can someone help me please?hello everyone.. im working with snow data in which present in one colon time and in other colon other parameter… and I have some missing value with pass of 7 minutes..now I would like to fill this gaps of 7 minutes with the interval of time so the last interval of time after this gaps… can someone help me please? hello everyone.. im working with snow data in which present in one colon time and in other colon other parameter… and I have some missing value with pass of 7 minutes..now I would like to fill this gaps of 7 minutes with the interval of time so the last interval of time after this gaps… can someone help me please? data MATLAB Answers — New Questions
When using Matlab engine and calling from multiple ranks, does each rank use one Matlab license?
I’m trying to write a python code that will call a Matlab function compiled for python. I need to know if, when running this code using multiple MPI processes, will each need one instance of Matlab? Will this be a potential bottleneck in opening multiple instances of Matlab?I’m trying to write a python code that will call a Matlab function compiled for python. I need to know if, when running this code using multiple MPI processes, will each need one instance of Matlab? Will this be a potential bottleneck in opening multiple instances of Matlab? I’m trying to write a python code that will call a Matlab function compiled for python. I need to know if, when running this code using multiple MPI processes, will each need one instance of Matlab? Will this be a potential bottleneck in opening multiple instances of Matlab? python, engine, license MATLAB Answers — New Questions
puzzles on using mutiple “addOptional” in function definitions
How to use multiple |addOptional| in function definitions?
|addOptional| is to add positional arguments. So I want to use it
in the following way.
* First I define the function and use two
|addOptional| there.
* Then I thought I can call in the following
two ways:
|myfunc(‘Jim’, 12)| and |myfunc(‘Jim’, 12, ‘cm’)|.
However, the latter one throws an error:
_The argument ‘cm’ is a string and does not match any parameter names. It failed validation
for the argument ‘units’._
But I can use as |myfunc(‘Jim’, 12, ‘units’, ‘cm’)|, which seems the second |addOptional| behaves as |addParameter|, no longer like positional arguments. What’s the mistake I made? Did I misunderstand something about |addOptional|? Many thanks!
Here is the definition of |myfunc()|:
function myfunc(name, varargin)
par = inputParser;
defaultHeight = 1;
defaultUnits = ‘inches’;
defaultShape = ‘rectangle’;
expectedShape = {‘square’, ‘rectangle’, ‘parallelogram’};
addRequired(par, ‘name’, @ischar);
addOptional(par, ‘height’, defaultHeight, @isnumeric);
addOptional(par, ‘units’, defaultUnits);
% addParameter(par, ‘shape’, defaultShape, …
% @(x) any(validatestring(x,expectedShape)));
parse(par, name, varargin{:});
disp(par.Results)
———
*Updates*:
I found that if the positional arguments are NOT literal string (maybe called as character array in helpdoc, the type shown by whos is char, e.g. ‘any words’), there would be no problems. For example, use |myfunc(‘Jim’, 3, 4)| which is defined as
function myfunc(name, varargin)
p = inputParser;
addRequired(p, ‘name’);
addOptional(p, ‘i’, 1);
addOptional(p, ‘j’, 2);
parse(p, name, varargin{:});
disp(p.Results)
The MATLAB help doc doesn’t mention anything about the limitations in data types of additional positional arguments. Does it imply this is a bug?
*Info*:
* MATLAB: 2016a
* OS: OS X SierraHow to use multiple |addOptional| in function definitions?
|addOptional| is to add positional arguments. So I want to use it
in the following way.
* First I define the function and use two
|addOptional| there.
* Then I thought I can call in the following
two ways:
|myfunc(‘Jim’, 12)| and |myfunc(‘Jim’, 12, ‘cm’)|.
However, the latter one throws an error:
_The argument ‘cm’ is a string and does not match any parameter names. It failed validation
for the argument ‘units’._
But I can use as |myfunc(‘Jim’, 12, ‘units’, ‘cm’)|, which seems the second |addOptional| behaves as |addParameter|, no longer like positional arguments. What’s the mistake I made? Did I misunderstand something about |addOptional|? Many thanks!
Here is the definition of |myfunc()|:
function myfunc(name, varargin)
par = inputParser;
defaultHeight = 1;
defaultUnits = ‘inches’;
defaultShape = ‘rectangle’;
expectedShape = {‘square’, ‘rectangle’, ‘parallelogram’};
addRequired(par, ‘name’, @ischar);
addOptional(par, ‘height’, defaultHeight, @isnumeric);
addOptional(par, ‘units’, defaultUnits);
% addParameter(par, ‘shape’, defaultShape, …
% @(x) any(validatestring(x,expectedShape)));
parse(par, name, varargin{:});
disp(par.Results)
———
*Updates*:
I found that if the positional arguments are NOT literal string (maybe called as character array in helpdoc, the type shown by whos is char, e.g. ‘any words’), there would be no problems. For example, use |myfunc(‘Jim’, 3, 4)| which is defined as
function myfunc(name, varargin)
p = inputParser;
addRequired(p, ‘name’);
addOptional(p, ‘i’, 1);
addOptional(p, ‘j’, 2);
parse(p, name, varargin{:});
disp(p.Results)
The MATLAB help doc doesn’t mention anything about the limitations in data types of additional positional arguments. Does it imply this is a bug?
*Info*:
* MATLAB: 2016a
* OS: OS X Sierra How to use multiple |addOptional| in function definitions?
|addOptional| is to add positional arguments. So I want to use it
in the following way.
* First I define the function and use two
|addOptional| there.
* Then I thought I can call in the following
two ways:
|myfunc(‘Jim’, 12)| and |myfunc(‘Jim’, 12, ‘cm’)|.
However, the latter one throws an error:
_The argument ‘cm’ is a string and does not match any parameter names. It failed validation
for the argument ‘units’._
But I can use as |myfunc(‘Jim’, 12, ‘units’, ‘cm’)|, which seems the second |addOptional| behaves as |addParameter|, no longer like positional arguments. What’s the mistake I made? Did I misunderstand something about |addOptional|? Many thanks!
Here is the definition of |myfunc()|:
function myfunc(name, varargin)
par = inputParser;
defaultHeight = 1;
defaultUnits = ‘inches’;
defaultShape = ‘rectangle’;
expectedShape = {‘square’, ‘rectangle’, ‘parallelogram’};
addRequired(par, ‘name’, @ischar);
addOptional(par, ‘height’, defaultHeight, @isnumeric);
addOptional(par, ‘units’, defaultUnits);
% addParameter(par, ‘shape’, defaultShape, …
% @(x) any(validatestring(x,expectedShape)));
parse(par, name, varargin{:});
disp(par.Results)
———
*Updates*:
I found that if the positional arguments are NOT literal string (maybe called as character array in helpdoc, the type shown by whos is char, e.g. ‘any words’), there would be no problems. For example, use |myfunc(‘Jim’, 3, 4)| which is defined as
function myfunc(name, varargin)
p = inputParser;
addRequired(p, ‘name’);
addOptional(p, ‘i’, 1);
addOptional(p, ‘j’, 2);
parse(p, name, varargin{:});
disp(p.Results)
The MATLAB help doc doesn’t mention anything about the limitations in data types of additional positional arguments. Does it imply this is a bug?
*Info*:
* MATLAB: 2016a
* OS: OS X Sierra addoptional, inputparser MATLAB Answers — New Questions
Does MATLAB work with Qualcomm arm processors? For example, Microsoft surface pro 11 with X Elite processor, and Windows 11 on Arm.
Qualcomm has released its latest X Elite processor and I really like its long battery life and low carbon footprint. I would like to ask the community if matlab runs well on windows on arm platform through an emulator? Will it affect normal ploting and coding? OS version:the latest Windows 11 on Arm.
Thank you!Qualcomm has released its latest X Elite processor and I really like its long battery life and low carbon footprint. I would like to ask the community if matlab runs well on windows on arm platform through an emulator? Will it affect normal ploting and coding? OS version:the latest Windows 11 on Arm.
Thank you! Qualcomm has released its latest X Elite processor and I really like its long battery life and low carbon footprint. I would like to ask the community if matlab runs well on windows on arm platform through an emulator? Will it affect normal ploting and coding? OS version:the latest Windows 11 on Arm.
Thank you! windows 11 on arm MATLAB Answers — New Questions
Flexible Body Model Builder – Equations used inside the package
Hi,
I am experimenting with the Flexible body builder app and have some questions regarding the math behind it’s implementation.
StiffnessMatrix: I understand that the stiffness matrix in the output ROM data is calculated from the Young’s modulus and the poison’s ratio. But I want to change this young’s modulus with every run of the simulation without building the flexible body ROM from scratch. For this, I first export the ROM data with E = 1MPa, define the flexible body using the exported ROM data in simscape, and multiply the stiffness matrix with E_flex which is the variable I want to change with every simulation run.
This is because K = EA/L and since E is a scalar, it can be scaled up or down with just a simple scalar multiplication.
I did some experimentation of this usage by exporting stiffness matrices for E = 1MPa, 2MPa,…etc and I can verify that
det ([StiffnessMat.2MPa]-2*[StiffnessMat.1MPa]) = 0
It would be great if you could verify this usage.
Number of Fixed-Interface Normal Modes: What is the optimal number of frequency modes to retain for a given geometry? It would depend on the geometry and the loading conditions themselves. Is there a way to assess what’s a good estimate for this?
Algorithm for model Generation: Is there a quantitative difference between the included integro-differential modeling framework and the PDE toolbox? What is the difference between the two methods. What is the recommended method?
It’d be great if you could provide some insights into these areas.Hi,
I am experimenting with the Flexible body builder app and have some questions regarding the math behind it’s implementation.
StiffnessMatrix: I understand that the stiffness matrix in the output ROM data is calculated from the Young’s modulus and the poison’s ratio. But I want to change this young’s modulus with every run of the simulation without building the flexible body ROM from scratch. For this, I first export the ROM data with E = 1MPa, define the flexible body using the exported ROM data in simscape, and multiply the stiffness matrix with E_flex which is the variable I want to change with every simulation run.
This is because K = EA/L and since E is a scalar, it can be scaled up or down with just a simple scalar multiplication.
I did some experimentation of this usage by exporting stiffness matrices for E = 1MPa, 2MPa,…etc and I can verify that
det ([StiffnessMat.2MPa]-2*[StiffnessMat.1MPa]) = 0
It would be great if you could verify this usage.
Number of Fixed-Interface Normal Modes: What is the optimal number of frequency modes to retain for a given geometry? It would depend on the geometry and the loading conditions themselves. Is there a way to assess what’s a good estimate for this?
Algorithm for model Generation: Is there a quantitative difference between the included integro-differential modeling framework and the PDE toolbox? What is the difference between the two methods. What is the recommended method?
It’d be great if you could provide some insights into these areas. Hi,
I am experimenting with the Flexible body builder app and have some questions regarding the math behind it’s implementation.
StiffnessMatrix: I understand that the stiffness matrix in the output ROM data is calculated from the Young’s modulus and the poison’s ratio. But I want to change this young’s modulus with every run of the simulation without building the flexible body ROM from scratch. For this, I first export the ROM data with E = 1MPa, define the flexible body using the exported ROM data in simscape, and multiply the stiffness matrix with E_flex which is the variable I want to change with every simulation run.
This is because K = EA/L and since E is a scalar, it can be scaled up or down with just a simple scalar multiplication.
I did some experimentation of this usage by exporting stiffness matrices for E = 1MPa, 2MPa,…etc and I can verify that
det ([StiffnessMat.2MPa]-2*[StiffnessMat.1MPa]) = 0
It would be great if you could verify this usage.
Number of Fixed-Interface Normal Modes: What is the optimal number of frequency modes to retain for a given geometry? It would depend on the geometry and the loading conditions themselves. Is there a way to assess what’s a good estimate for this?
Algorithm for model Generation: Is there a quantitative difference between the included integro-differential modeling framework and the PDE toolbox? What is the difference between the two methods. What is the recommended method?
It’d be great if you could provide some insights into these areas. flexible body, flexible body model builder, craig bampton, reduced order modeling, rom MATLAB Answers — New Questions
“Unable to fetch some archives” when running sudo apt-get install matlab-rpi
I’m tyrying to install matlab-rpi on my pi, and I have several errors that begin with "Failed to fetch", and the error at the bottom says "E: Unable to fetch some archives, maybe run apt-get update or try with –fix-missing?". I am running Raspbian bullseye.I’m tyrying to install matlab-rpi on my pi, and I have several errors that begin with "Failed to fetch", and the error at the bottom says "E: Unable to fetch some archives, maybe run apt-get update or try with –fix-missing?". I am running Raspbian bullseye. I’m tyrying to install matlab-rpi on my pi, and I have several errors that begin with "Failed to fetch", and the error at the bottom says "E: Unable to fetch some archives, maybe run apt-get update or try with –fix-missing?". I am running Raspbian bullseye. raspberry pi MATLAB Answers — New Questions
run the system for 24 hours in short period
hello,
so i have model in simulink, it runs for Total time=2 seconds, and sampling time for microgrid=1e-6, sampling time for model predictive controller=0.1; the thing is when i run the simulation, it takes about 3 to 4 minutes to finish because of the sampling time, when in fact the Total time is 2 seconds, and i want the axe "X" to show results for 24 hours, so what can i do?
obiously i cannot run it for 24 hours, it will take few days to finish,
thanks in advancehello,
so i have model in simulink, it runs for Total time=2 seconds, and sampling time for microgrid=1e-6, sampling time for model predictive controller=0.1; the thing is when i run the simulation, it takes about 3 to 4 minutes to finish because of the sampling time, when in fact the Total time is 2 seconds, and i want the axe "X" to show results for 24 hours, so what can i do?
obiously i cannot run it for 24 hours, it will take few days to finish,
thanks in advance hello,
so i have model in simulink, it runs for Total time=2 seconds, and sampling time for microgrid=1e-6, sampling time for model predictive controller=0.1; the thing is when i run the simulation, it takes about 3 to 4 minutes to finish because of the sampling time, when in fact the Total time is 2 seconds, and i want the axe "X" to show results for 24 hours, so what can i do?
obiously i cannot run it for 24 hours, it will take few days to finish,
thanks in advance simulation time, sampling time, 24 hours MATLAB Answers — New Questions
Long simulation times with flexible simscape
Hello everybody,
I’ve been using flexible multibody simscape for a few days.
I have found the simulations to be very slow both using the default bodies and using Reduce Order. is this a common problem?
how can i fix it?Hello everybody,
I’ve been using flexible multibody simscape for a few days.
I have found the simulations to be very slow both using the default bodies and using Reduce Order. is this a common problem?
how can i fix it? Hello everybody,
I’ve been using flexible multibody simscape for a few days.
I have found the simulations to be very slow both using the default bodies and using Reduce Order. is this a common problem?
how can i fix it? simscape, simulation, simscape multibody, simscape flexible MATLAB Answers — New Questions
Error sending data to ThingSpeak – Data Intervals??
I am sending data from sensors on a XIAO esp32 S3 – there are 6 fields of data – I have timed to send every 5 minutes. When I initiate the code it sends 2 or 3 sets, but then I get an error message saying ‘Error sending data to ThingSpeak’ – and it doesnt send any more data unless I reset the esp , yet the data is being collected as shown on the serial monitor. I have played around witht the time intervals but no joy. Is it a function of thinkspeak or should I investigate the code?I am sending data from sensors on a XIAO esp32 S3 – there are 6 fields of data – I have timed to send every 5 minutes. When I initiate the code it sends 2 or 3 sets, but then I get an error message saying ‘Error sending data to ThingSpeak’ – and it doesnt send any more data unless I reset the esp , yet the data is being collected as shown on the serial monitor. I have played around witht the time intervals but no joy. Is it a function of thinkspeak or should I investigate the code? I am sending data from sensors on a XIAO esp32 S3 – there are 6 fields of data – I have timed to send every 5 minutes. When I initiate the code it sends 2 or 3 sets, but then I get an error message saying ‘Error sending data to ThingSpeak’ – and it doesnt send any more data unless I reset the esp , yet the data is being collected as shown on the serial monitor. I have played around witht the time intervals but no joy. Is it a function of thinkspeak or should I investigate the code? thingspeak, data, fields MATLAB Answers — New Questions
Missing gpu.lib in R2024a
I am trying to compile CUDA version of Matconvnet on R2024a.
Mex reporting an error of missing gpu.lib and indeed I can’t find it in R2024a. It however, exists in the previous versions, for example R2023b.
Is there a replacement for gpu.lib in R2024a?
Thank you,
IliaI am trying to compile CUDA version of Matconvnet on R2024a.
Mex reporting an error of missing gpu.lib and indeed I can’t find it in R2024a. It however, exists in the previous versions, for example R2023b.
Is there a replacement for gpu.lib in R2024a?
Thank you,
Ilia I am trying to compile CUDA version of Matconvnet on R2024a.
Mex reporting an error of missing gpu.lib and indeed I can’t find it in R2024a. It however, exists in the previous versions, for example R2023b.
Is there a replacement for gpu.lib in R2024a?
Thank you,
Ilia r2024a, gpu.lib MATLAB Answers — New Questions
Using a specific wavelet transform on new ECG signals
I’m working with ECG signals and am trying to use a wavelet technique to reduce some of the noise in various data sets. I was able to use a ‘Continuous Wavelet 1-D using FFT’ technique in the Wavelet Analyzer toolkit to achieve a good response (at least on one sample). I have exported the CWTFT struct to a variable called CWTS and want to apply the same wavelet function on a number of other pieces of data.
My overall goal would be to have a command or function that applies the same wavelet transform I generated to new input data.
I am somewhat of a novice in the wavelet / signal processing realm, so if there are smarter or better ways to approach this, please let me know. I tried several random wavelets and this option seemed to do the best for me.
How I generated the wavelet transformed data: 1-D/FFT -> ‘dog’, 6 parameter, linear analysis, linear synthesis, scale limited to first 16 options
Here are some of the details:
CWTS Contents: CWTS.cfs, CWTS.scales, CWTS.frequencies, CWTS.omega, CWTS.meanSIG, CWTS.dt, CWTS.wav
Here are the pre / post transform plots of the data:
Pre-Transform/Raw Data
<</matlabcentral/answers/uploaded_files/71656/PreWavelet.jpg>>
Post-Transform
<</matlabcentral/answers/uploaded_files/71657/post-wavelet.jpg>>
Thanks for any help you guys can offer!I’m working with ECG signals and am trying to use a wavelet technique to reduce some of the noise in various data sets. I was able to use a ‘Continuous Wavelet 1-D using FFT’ technique in the Wavelet Analyzer toolkit to achieve a good response (at least on one sample). I have exported the CWTFT struct to a variable called CWTS and want to apply the same wavelet function on a number of other pieces of data.
My overall goal would be to have a command or function that applies the same wavelet transform I generated to new input data.
I am somewhat of a novice in the wavelet / signal processing realm, so if there are smarter or better ways to approach this, please let me know. I tried several random wavelets and this option seemed to do the best for me.
How I generated the wavelet transformed data: 1-D/FFT -> ‘dog’, 6 parameter, linear analysis, linear synthesis, scale limited to first 16 options
Here are some of the details:
CWTS Contents: CWTS.cfs, CWTS.scales, CWTS.frequencies, CWTS.omega, CWTS.meanSIG, CWTS.dt, CWTS.wav
Here are the pre / post transform plots of the data:
Pre-Transform/Raw Data
<</matlabcentral/answers/uploaded_files/71656/PreWavelet.jpg>>
Post-Transform
<</matlabcentral/answers/uploaded_files/71657/post-wavelet.jpg>>
Thanks for any help you guys can offer! I’m working with ECG signals and am trying to use a wavelet technique to reduce some of the noise in various data sets. I was able to use a ‘Continuous Wavelet 1-D using FFT’ technique in the Wavelet Analyzer toolkit to achieve a good response (at least on one sample). I have exported the CWTFT struct to a variable called CWTS and want to apply the same wavelet function on a number of other pieces of data.
My overall goal would be to have a command or function that applies the same wavelet transform I generated to new input data.
I am somewhat of a novice in the wavelet / signal processing realm, so if there are smarter or better ways to approach this, please let me know. I tried several random wavelets and this option seemed to do the best for me.
How I generated the wavelet transformed data: 1-D/FFT -> ‘dog’, 6 parameter, linear analysis, linear synthesis, scale limited to first 16 options
Here are some of the details:
CWTS Contents: CWTS.cfs, CWTS.scales, CWTS.frequencies, CWTS.omega, CWTS.meanSIG, CWTS.dt, CWTS.wav
Here are the pre / post transform plots of the data:
Pre-Transform/Raw Data
<</matlabcentral/answers/uploaded_files/71656/PreWavelet.jpg>>
Post-Transform
<</matlabcentral/answers/uploaded_files/71657/post-wavelet.jpg>>
Thanks for any help you guys can offer! wavelet, ecg, transform, signal processing, fft MATLAB Answers — New Questions
Why is my solution incorrect here? How to implement the boundary conditions correctly
I am trying to solve for a 1D BVP with Neumann boundary conditions. I am using chebyshev derivative and so to implement these BCs I am initially replacing the first/last row of the second derivative to the values of first/last row of the first derivative and set them to zero so that u'(-1)=u'(1)=0 (i.e. first derivative at left and right boundries = 0). I am not sure if that’s what I am doing in the code though because I get completely wrong numerical solution.
clearvars; clc; close all;
N=10;%
[D,x] = cheb(N); D2 = D^2;
%D2(1,:) = D(1,:); %Right
%D2(N+1,:) = D(N+1,:); %Neumann BC at left endpoint (x=-1), esp w/ D
D2([1 N+1],:) = D([1 N+1],:);
%
BC=-D([1 N+1],[1 N+1])D([1 N+1],2:N);
a=-1; b=1;
u =(1/6)*x .*(6*a*b-3*a*x-3*b*x+2* x.^2);
%Exact Solution
figure
plot(x,u);
n =ones(size(u));
dndx = D *n;
dudx = D *u;
prod1 = dndx .* dudx;
du2dx = D * dudx;
prod2 = n .*du2dx;
invN = (1./n) ;
source = prod1 + prod2;
uold = ones(size(u));
max_iter =500;
err_max = 1e-8;
for iterations = 1:max_iter
phikmax_old = (max(abs(uold)));
duoldx =D *uold;
dnudx = dndx .* duoldx;
ffh = source;
RHS = ffh – dnudx;
Stilde =invN .* RHS;
unew = D2Stilde; %
%Enforce BCs
unew([1 N+1],:) = BC*unew(2:N,:);
phikmax = (max(abs(unew)));
if phikmax < err_max
it_error = err_max /2;
else
it_error = abs( phikmax – phikmax_old) / phikmax;
end
if it_error < err_max
break;
end
uold = unew;
end
figure
plot(x,unew,’–rs’,x,(u));
legend(‘Num solution’,’Exact solution’)
Function Cheb is
function [ D, x ] = cheb ( N )
if ( N == 0 )
D = 0.0;
x = 1.0;
return
end
x = cos ( pi * ( 0 : N ) / N )’;
c = [ 2.0; ones(N-1,1); 2.0 ] .* (-1.0).^(0:N)’;
X = repmat ( x, 1, N + 1 );
dX = X – X’;
% Set the off diagonal entries.
D =( c * (1.0 ./ c )’ ) ./ ( dX + ( eye ( N + 1 ) ) );
% Diagonal entries.
D = D – diag ( sum ( D’ ) );
return
endI am trying to solve for a 1D BVP with Neumann boundary conditions. I am using chebyshev derivative and so to implement these BCs I am initially replacing the first/last row of the second derivative to the values of first/last row of the first derivative and set them to zero so that u'(-1)=u'(1)=0 (i.e. first derivative at left and right boundries = 0). I am not sure if that’s what I am doing in the code though because I get completely wrong numerical solution.
clearvars; clc; close all;
N=10;%
[D,x] = cheb(N); D2 = D^2;
%D2(1,:) = D(1,:); %Right
%D2(N+1,:) = D(N+1,:); %Neumann BC at left endpoint (x=-1), esp w/ D
D2([1 N+1],:) = D([1 N+1],:);
%
BC=-D([1 N+1],[1 N+1])D([1 N+1],2:N);
a=-1; b=1;
u =(1/6)*x .*(6*a*b-3*a*x-3*b*x+2* x.^2);
%Exact Solution
figure
plot(x,u);
n =ones(size(u));
dndx = D *n;
dudx = D *u;
prod1 = dndx .* dudx;
du2dx = D * dudx;
prod2 = n .*du2dx;
invN = (1./n) ;
source = prod1 + prod2;
uold = ones(size(u));
max_iter =500;
err_max = 1e-8;
for iterations = 1:max_iter
phikmax_old = (max(abs(uold)));
duoldx =D *uold;
dnudx = dndx .* duoldx;
ffh = source;
RHS = ffh – dnudx;
Stilde =invN .* RHS;
unew = D2Stilde; %
%Enforce BCs
unew([1 N+1],:) = BC*unew(2:N,:);
phikmax = (max(abs(unew)));
if phikmax < err_max
it_error = err_max /2;
else
it_error = abs( phikmax – phikmax_old) / phikmax;
end
if it_error < err_max
break;
end
uold = unew;
end
figure
plot(x,unew,’–rs’,x,(u));
legend(‘Num solution’,’Exact solution’)
Function Cheb is
function [ D, x ] = cheb ( N )
if ( N == 0 )
D = 0.0;
x = 1.0;
return
end
x = cos ( pi * ( 0 : N ) / N )’;
c = [ 2.0; ones(N-1,1); 2.0 ] .* (-1.0).^(0:N)’;
X = repmat ( x, 1, N + 1 );
dX = X – X’;
% Set the off diagonal entries.
D =( c * (1.0 ./ c )’ ) ./ ( dX + ( eye ( N + 1 ) ) );
% Diagonal entries.
D = D – diag ( sum ( D’ ) );
return
end I am trying to solve for a 1D BVP with Neumann boundary conditions. I am using chebyshev derivative and so to implement these BCs I am initially replacing the first/last row of the second derivative to the values of first/last row of the first derivative and set them to zero so that u'(-1)=u'(1)=0 (i.e. first derivative at left and right boundries = 0). I am not sure if that’s what I am doing in the code though because I get completely wrong numerical solution.
clearvars; clc; close all;
N=10;%
[D,x] = cheb(N); D2 = D^2;
%D2(1,:) = D(1,:); %Right
%D2(N+1,:) = D(N+1,:); %Neumann BC at left endpoint (x=-1), esp w/ D
D2([1 N+1],:) = D([1 N+1],:);
%
BC=-D([1 N+1],[1 N+1])D([1 N+1],2:N);
a=-1; b=1;
u =(1/6)*x .*(6*a*b-3*a*x-3*b*x+2* x.^2);
%Exact Solution
figure
plot(x,u);
n =ones(size(u));
dndx = D *n;
dudx = D *u;
prod1 = dndx .* dudx;
du2dx = D * dudx;
prod2 = n .*du2dx;
invN = (1./n) ;
source = prod1 + prod2;
uold = ones(size(u));
max_iter =500;
err_max = 1e-8;
for iterations = 1:max_iter
phikmax_old = (max(abs(uold)));
duoldx =D *uold;
dnudx = dndx .* duoldx;
ffh = source;
RHS = ffh – dnudx;
Stilde =invN .* RHS;
unew = D2Stilde; %
%Enforce BCs
unew([1 N+1],:) = BC*unew(2:N,:);
phikmax = (max(abs(unew)));
if phikmax < err_max
it_error = err_max /2;
else
it_error = abs( phikmax – phikmax_old) / phikmax;
end
if it_error < err_max
break;
end
uold = unew;
end
figure
plot(x,unew,’–rs’,x,(u));
legend(‘Num solution’,’Exact solution’)
Function Cheb is
function [ D, x ] = cheb ( N )
if ( N == 0 )
D = 0.0;
x = 1.0;
return
end
x = cos ( pi * ( 0 : N ) / N )’;
c = [ 2.0; ones(N-1,1); 2.0 ] .* (-1.0).^(0:N)’;
X = repmat ( x, 1, N + 1 );
dX = X – X’;
% Set the off diagonal entries.
D =( c * (1.0 ./ c )’ ) ./ ( dX + ( eye ( N + 1 ) ) );
% Diagonal entries.
D = D – diag ( sum ( D’ ) );
return
end neumann, boundary, condition MATLAB Answers — New Questions
What is the most computationally efficient factorization of a matrix A?
I have an exam in 24 hours for my Linear Algebra and Geometry course and we have a section with MATLAB related questions. In a simulation test, I was asked to "solve the linear system Ax=b by the solution of two triangular systems, using the most computationally efficient factorization of a matrix A". What does this mean? What is "the solution of two triangular systems"? How would I choose the most computationally efficient factorization of a matrix?
The question I am talking about is the following, word for word:
Let Ax = b, a linear system of order 18, where A is a symmetric, tridiagonal matrix with all elements equal to 6 on the main diagonal and to 3 on the superior and inferior codiagonals. The elements of b are linearly spaced numbers in [8,11]. Compute the eigenvalues of A and, based on their properties, solve the linear system Ax = b by the solution of two triangular systems, using the most computationally efficient factorization of A. The 1-norm of the vector obtained by summing the vector that is the solution of the inferior triangular system and the one that is the solution of the superior triangular system is, approximately:
A) 8.4419e+01
B) 3.4638e+00
C) 7.5558e-01
D) 9.5155e-01
E) 6.5559e+01
Now, we went over all this with a friend and even asked chatgpt to give us an idea. In our curriculum, the main factorizations we learned about are the PA=LU, the Choleski, SVD and the QR factorization. We thought maybe after seeing all the eigenvalues are positive, and thus the matrix positive definite, the choleski should be used. But we aren’t really sure. Why would the others be less efficient? With the code chatgpt gave us, tweaked a little here and there, we got the answer right using the choleski. But we really aren’t sure we understand the wording of the question and the systems at play. Here is the code in question:
Thank you in advance for an answer to such a long question.
clear, clc
n = 18;
A = eye(n)*6 + diag(ones(n-1,1), 1)*3 + diag(ones(n-1,1), -1)*3;
b = linspace(8,11,18);
eigens = eig(A)
R = chol(A);
disp("R’*R*x = b so y = R*x = b")
Y = R’b’;
x = RY;
finalVector = x + Y;
norm(finalVector, 1)I have an exam in 24 hours for my Linear Algebra and Geometry course and we have a section with MATLAB related questions. In a simulation test, I was asked to "solve the linear system Ax=b by the solution of two triangular systems, using the most computationally efficient factorization of a matrix A". What does this mean? What is "the solution of two triangular systems"? How would I choose the most computationally efficient factorization of a matrix?
The question I am talking about is the following, word for word:
Let Ax = b, a linear system of order 18, where A is a symmetric, tridiagonal matrix with all elements equal to 6 on the main diagonal and to 3 on the superior and inferior codiagonals. The elements of b are linearly spaced numbers in [8,11]. Compute the eigenvalues of A and, based on their properties, solve the linear system Ax = b by the solution of two triangular systems, using the most computationally efficient factorization of A. The 1-norm of the vector obtained by summing the vector that is the solution of the inferior triangular system and the one that is the solution of the superior triangular system is, approximately:
A) 8.4419e+01
B) 3.4638e+00
C) 7.5558e-01
D) 9.5155e-01
E) 6.5559e+01
Now, we went over all this with a friend and even asked chatgpt to give us an idea. In our curriculum, the main factorizations we learned about are the PA=LU, the Choleski, SVD and the QR factorization. We thought maybe after seeing all the eigenvalues are positive, and thus the matrix positive definite, the choleski should be used. But we aren’t really sure. Why would the others be less efficient? With the code chatgpt gave us, tweaked a little here and there, we got the answer right using the choleski. But we really aren’t sure we understand the wording of the question and the systems at play. Here is the code in question:
Thank you in advance for an answer to such a long question.
clear, clc
n = 18;
A = eye(n)*6 + diag(ones(n-1,1), 1)*3 + diag(ones(n-1,1), -1)*3;
b = linspace(8,11,18);
eigens = eig(A)
R = chol(A);
disp("R’*R*x = b so y = R*x = b")
Y = R’b’;
x = RY;
finalVector = x + Y;
norm(finalVector, 1) I have an exam in 24 hours for my Linear Algebra and Geometry course and we have a section with MATLAB related questions. In a simulation test, I was asked to "solve the linear system Ax=b by the solution of two triangular systems, using the most computationally efficient factorization of a matrix A". What does this mean? What is "the solution of two triangular systems"? How would I choose the most computationally efficient factorization of a matrix?
The question I am talking about is the following, word for word:
Let Ax = b, a linear system of order 18, where A is a symmetric, tridiagonal matrix with all elements equal to 6 on the main diagonal and to 3 on the superior and inferior codiagonals. The elements of b are linearly spaced numbers in [8,11]. Compute the eigenvalues of A and, based on their properties, solve the linear system Ax = b by the solution of two triangular systems, using the most computationally efficient factorization of A. The 1-norm of the vector obtained by summing the vector that is the solution of the inferior triangular system and the one that is the solution of the superior triangular system is, approximately:
A) 8.4419e+01
B) 3.4638e+00
C) 7.5558e-01
D) 9.5155e-01
E) 6.5559e+01
Now, we went over all this with a friend and even asked chatgpt to give us an idea. In our curriculum, the main factorizations we learned about are the PA=LU, the Choleski, SVD and the QR factorization. We thought maybe after seeing all the eigenvalues are positive, and thus the matrix positive definite, the choleski should be used. But we aren’t really sure. Why would the others be less efficient? With the code chatgpt gave us, tweaked a little here and there, we got the answer right using the choleski. But we really aren’t sure we understand the wording of the question and the systems at play. Here is the code in question:
Thank you in advance for an answer to such a long question.
clear, clc
n = 18;
A = eye(n)*6 + diag(ones(n-1,1), 1)*3 + diag(ones(n-1,1), -1)*3;
b = linspace(8,11,18);
eigens = eig(A)
R = chol(A);
disp("R’*R*x = b so y = R*x = b")
Y = R’b’;
x = RY;
finalVector = x + Y;
norm(finalVector, 1) efficiency, factorization MATLAB Answers — New Questions