So I have the funcion:
function y=function_for_Ex6a(x)
% This is a simple function
y=4*(1-cos(x)).^2+x.^3;
end
And I want to use a for loop to run through this function 100 times with an initial input of x=4. I need the output value of y to replace the input value for each interation. Any ideas? Please don’t just send links to other websites. Thanks.So I have the funcion:
function y=function_for_Ex6a(x)
% This is a simple function
y=4*(1-cos(x)).^2+x.^3;
end
And I want to use a for loop to run through this function 100 times with an initial input of x=4. I need the output value of y to replace the input value for each interation. Any ideas? Please don’t just send links to other websites. Thanks. So I have the funcion:
function y=function_for_Ex6a(x)
% This is a simple function
y=4*(1-cos(x)).^2+x.^3;
end
And I want to use a for loop to run through this function 100 times with an initial input of x=4. I need the output value of y to replace the input value for each interation. Any ideas? Please don’t just send links to other websites. Thanks. looping variables MATLAB Answers — New Questions

​

Hi,
i would like to forecast the Cl parameter for 10 years ahead. But keep getting error : Index in position 2 exceeds array bounds. Index must not exceed 1. by the way , when running the code, i have remove the heading, words from the excel. but still getting the error
Error in preparets (line 317)
xi = xx(:,FBS+((1-net.numInputDelays):0));
Error in timee (line 25)
[X,inputStates,layerStates,T] = preparets(net,inputSeries,{},targetSeries);

Please help to check which part that i did wrong.
data = readmatrix ("104CPr.csv");
Input = data(:,1:(end-1));
Target = data(:,end);
X = con2seq(Input)’;
T = con2seq(Target)’;
%% Data Preparation
N = 10; % Multi-step ahead prediction
% Input and target series are divided in two groups of data:
% 1st group: used to train the network
inputSeries = X(1:end-N);
targetSeries = T(1:end-N);
%% Network architecture
delay = 2;
inputDelays = 1:2;
feedbackDelays = 1:2;
hiddenLayerSize = 10;
trainFcn = ‘trainbr’;
% Network Creation
net = narxnet(inputDelays,feedbackDelays,hiddenLayerSize,’open’, trainFcn);
[X,inputStates,layerStates,T] = preparets(net,inputSeries,{},targetSeries);
net.divideParam.trainRatio = 70/100;
net.divideParam.valRatio = 15/100;
net.divideParam.testRatio = 15/100;
[net,Y] = train(net,X,T,inputStates,layerStates);
outputs = net(X,inputStates,layerStates);
errors = gsubtract(T,Y);
performance = perform(net,T,Y);
view(net)
delay = 2;
netc = closeloop(net);
netc.name = [net.name ‘ – Closed Loop’];
view(netc)
[xc,xic,aic,tc] = preparets(netc,inputSeries,{},targetSeries);
yc = netc(xc,xic,aic);
closedLoopPerformance = perform(netc,tc,yc);
% Early Prediction Network
nets = removedelay(net);
nets.name = [net.name ‘ – Predict One Step Ahead’];
view(nets)
[xs,xis,ais,ts] = preparets(nets,inputSeries,{},targetSeries);
ys = nets(xs,xis,ais);
earlyPredictPerformance = perform(nets,ts,ys);
%%5. Multi-step ahead prediction
inputSeriesPred = [inputSeries(end-delay+1:end),inputSeriesVal];
targetSeriesPred = [targetSeries(end-delay+1:end), con2seq(nan(1,N))];
[Xs,Xi,Ai,Ts] = preparets(netc,inputSeriesPred,{},targetSeriesPred);
yPred = netc(Xs,Xi,Ai);
perf = perform(net,yPred,targetSeriesVal);
figure;
plot([cell2mat(targetSeries),nan(1,N);
nan(1,length(targetSeries)),cell2mat(yPred);
nan(1,length(targetSeries)),cell2mat(targetSeriesVal)]’)
legend(‘Original Targets’,’Network Predictions’,’Expected Outputs’);Hi,
i would like to forecast the Cl parameter for 10 years ahead. But keep getting error : Index in position 2 exceeds array bounds. Index must not exceed 1. by the way , when running the code, i have remove the heading, words from the excel. but still getting the error
Error in preparets (line 317)
xi = xx(:,FBS+((1-net.numInputDelays):0));
Error in timee (line 25)
[X,inputStates,layerStates,T] = preparets(net,inputSeries,{},targetSeries);

Please help to check which part that i did wrong.
data = readmatrix ("104CPr.csv");
Input = data(:,1:(end-1));
Target = data(:,end);
X = con2seq(Input)’;
T = con2seq(Target)’;
%% Data Preparation
N = 10; % Multi-step ahead prediction
% Input and target series are divided in two groups of data:
% 1st group: used to train the network
inputSeries = X(1:end-N);
targetSeries = T(1:end-N);
%% Network architecture
delay = 2;
inputDelays = 1:2;
feedbackDelays = 1:2;
hiddenLayerSize = 10;
trainFcn = ‘trainbr’;
% Network Creation
net = narxnet(inputDelays,feedbackDelays,hiddenLayerSize,’open’, trainFcn);
[X,inputStates,layerStates,T] = preparets(net,inputSeries,{},targetSeries);
net.divideParam.trainRatio = 70/100;
net.divideParam.valRatio = 15/100;
net.divideParam.testRatio = 15/100;
[net,Y] = train(net,X,T,inputStates,layerStates);
outputs = net(X,inputStates,layerStates);
errors = gsubtract(T,Y);
performance = perform(net,T,Y);
view(net)
delay = 2;
netc = closeloop(net);
netc.name = [net.name ‘ – Closed Loop’];
view(netc)
[xc,xic,aic,tc] = preparets(netc,inputSeries,{},targetSeries);
yc = netc(xc,xic,aic);
closedLoopPerformance = perform(netc,tc,yc);
% Early Prediction Network
nets = removedelay(net);
nets.name = [net.name ‘ – Predict One Step Ahead’];
view(nets)
[xs,xis,ais,ts] = preparets(nets,inputSeries,{},targetSeries);
ys = nets(xs,xis,ais);
earlyPredictPerformance = perform(nets,ts,ys);
%%5. Multi-step ahead prediction
inputSeriesPred = [inputSeries(end-delay+1:end),inputSeriesVal];
targetSeriesPred = [targetSeries(end-delay+1:end), con2seq(nan(1,N))];
[Xs,Xi,Ai,Ts] = preparets(netc,inputSeriesPred,{},targetSeriesPred);
yPred = netc(Xs,Xi,Ai);
perf = perform(net,yPred,targetSeriesVal);
figure;
plot([cell2mat(targetSeries),nan(1,N);
nan(1,length(targetSeries)),cell2mat(yPred);
nan(1,length(targetSeries)),cell2mat(targetSeriesVal)]’)
legend(‘Original Targets’,’Network Predictions’,’Expected Outputs’); Hi,
i would like to forecast the Cl parameter for 10 years ahead. But keep getting error : Index in position 2 exceeds array bounds. Index must not exceed 1. by the way , when running the code, i have remove the heading, words from the excel. but still getting the error
Error in preparets (line 317)
xi = xx(:,FBS+((1-net.numInputDelays):0));
Error in timee (line 25)
[X,inputStates,layerStates,T] = preparets(net,inputSeries,{},targetSeries);

Please help to check which part that i did wrong.
data = readmatrix ("104CPr.csv");
Input = data(:,1:(end-1));
Target = data(:,end);
X = con2seq(Input)’;
T = con2seq(Target)’;
%% Data Preparation
N = 10; % Multi-step ahead prediction
% Input and target series are divided in two groups of data:
% 1st group: used to train the network
inputSeries = X(1:end-N);
targetSeries = T(1:end-N);
%% Network architecture
delay = 2;
inputDelays = 1:2;
feedbackDelays = 1:2;
hiddenLayerSize = 10;
trainFcn = ‘trainbr’;
% Network Creation
net = narxnet(inputDelays,feedbackDelays,hiddenLayerSize,’open’, trainFcn);
[X,inputStates,layerStates,T] = preparets(net,inputSeries,{},targetSeries);
net.divideParam.trainRatio = 70/100;
net.divideParam.valRatio = 15/100;
net.divideParam.testRatio = 15/100;
[net,Y] = train(net,X,T,inputStates,layerStates);
outputs = net(X,inputStates,layerStates);
errors = gsubtract(T,Y);
performance = perform(net,T,Y);
view(net)
delay = 2;
netc = closeloop(net);
netc.name = [net.name ‘ – Closed Loop’];
view(netc)
[xc,xic,aic,tc] = preparets(netc,inputSeries,{},targetSeries);
yc = netc(xc,xic,aic);
closedLoopPerformance = perform(netc,tc,yc);
% Early Prediction Network
nets = removedelay(net);
nets.name = [net.name ‘ – Predict One Step Ahead’];
view(nets)
[xs,xis,ais,ts] = preparets(nets,inputSeries,{},targetSeries);
ys = nets(xs,xis,ais);
earlyPredictPerformance = perform(nets,ts,ys);
%%5. Multi-step ahead prediction
inputSeriesPred = [inputSeries(end-delay+1:end),inputSeriesVal];
targetSeriesPred = [targetSeries(end-delay+1:end), con2seq(nan(1,N))];
[Xs,Xi,Ai,Ts] = preparets(netc,inputSeriesPred,{},targetSeriesPred);
yPred = netc(Xs,Xi,Ai);
perf = perform(net,yPred,targetSeriesVal);
figure;
plot([cell2mat(targetSeries),nan(1,N);
nan(1,length(targetSeries)),cell2mat(yPred);
nan(1,length(targetSeries)),cell2mat(targetSeriesVal)]’)
legend(‘Original Targets’,’Network Predictions’,’Expected Outputs’); ann, time series MATLAB Answers — New Questions

​

Hello all,

I need to create the subsystem which needs to be optionally Enable/Disabled before simulation starts but there si no need to be enabled/disable during the simulation …
What are the possiblities to create a such subsystem in Simulink?

Thank you very much for your suggestion.Hello all,

I need to create the subsystem which needs to be optionally Enable/Disabled before simulation starts but there si no need to be enabled/disable during the simulation …
What are the possiblities to create a such subsystem in Simulink?

Thank you very much for your suggestion. Hello all,

I need to create the subsystem which needs to be optionally Enable/Disabled before simulation starts but there si no need to be enabled/disable during the simulation …
What are the possiblities to create a such subsystem in Simulink?

Thank you very much for your suggestion. enabled subsystem MATLAB Answers — New Questions

​

Hello, I am using Simulink to generate C code and burn it to my DSP28388 to test Ethernet. The test program and code generation result are shown in the figure. It can be seen that the code generation result reported errors, the reason is that some structure is not defined, thank you very much for providing me with some ideas to solve the problem.Hello, I am using Simulink to generate C code and burn it to my DSP28388 to test Ethernet. The test program and code generation result are shown in the figure. It can be seen that the code generation result reported errors, the reason is that some structure is not defined, thank you very much for providing me with some ideas to solve the problem. Hello, I am using Simulink to generate C code and burn it to my DSP28388 to test Ethernet. The test program and code generation result are shown in the figure. It can be seen that the code generation result reported errors, the reason is that some structure is not defined, thank you very much for providing me with some ideas to solve the problem. code generation, 28388, simulink MATLAB Answers — New Questions

​

Hi,
I have updated my install from R2021a to R2022a. Using the RL toolbox when running the episode manager with the following code in R2022a, when I go to stop the training early, via "Stop Training" in episode manager, the training does not stop, it seems the only way to actual stop the current training early is via the "stop" button on the "run" tab in live editor. Has anyone encounered this in R2022a? and will the agent buffer still be saved via just stopping the script rather than stopping training via EM?
T = 0.066;
maxepisodes = 100;
maxsteps = ceil(T/Ts_agent)
trainingOpts = rlTrainingOptions(…
‘MaxEpisodes’,maxepisodes, …
‘MaxStepsPerEpisode’,maxsteps, …
‘StopTrainingCriteria’,"EpisodeCount",…
‘StopTrainingValue’,maxepisodes,…
‘ScoreAveragingWindowLength’,100,…
‘SaveAgentCriteria’,"EpisodeCount",…
‘SaveAgentValue’,maxepisodes,…
‘SaveAgentDirectory’,"savedAgents",…
‘Verbose’,true,…
‘Plots’,’training-progress’)

Cheers,

PatrickHi,
I have updated my install from R2021a to R2022a. Using the RL toolbox when running the episode manager with the following code in R2022a, when I go to stop the training early, via "Stop Training" in episode manager, the training does not stop, it seems the only way to actual stop the current training early is via the "stop" button on the "run" tab in live editor. Has anyone encounered this in R2022a? and will the agent buffer still be saved via just stopping the script rather than stopping training via EM?
T = 0.066;
maxepisodes = 100;
maxsteps = ceil(T/Ts_agent)
trainingOpts = rlTrainingOptions(…
‘MaxEpisodes’,maxepisodes, …
‘MaxStepsPerEpisode’,maxsteps, …
‘StopTrainingCriteria’,"EpisodeCount",…
‘StopTrainingValue’,maxepisodes,…
‘ScoreAveragingWindowLength’,100,…
‘SaveAgentCriteria’,"EpisodeCount",…
‘SaveAgentValue’,maxepisodes,…
‘SaveAgentDirectory’,"savedAgents",…
‘Verbose’,true,…
‘Plots’,’training-progress’)

Cheers,

Patrick Hi,
I have updated my install from R2021a to R2022a. Using the RL toolbox when running the episode manager with the following code in R2022a, when I go to stop the training early, via "Stop Training" in episode manager, the training does not stop, it seems the only way to actual stop the current training early is via the "stop" button on the "run" tab in live editor. Has anyone encounered this in R2022a? and will the agent buffer still be saved via just stopping the script rather than stopping training via EM?
T = 0.066;
maxepisodes = 100;
maxsteps = ceil(T/Ts_agent)
trainingOpts = rlTrainingOptions(…
‘MaxEpisodes’,maxepisodes, …
‘MaxStepsPerEpisode’,maxsteps, …
‘StopTrainingCriteria’,"EpisodeCount",…
‘StopTrainingValue’,maxepisodes,…
‘ScoreAveragingWindowLength’,100,…
‘SaveAgentCriteria’,"EpisodeCount",…
‘SaveAgentValue’,maxepisodes,…
‘SaveAgentDirectory’,"savedAgents",…
‘Verbose’,true,…
‘Plots’,’training-progress’)

Cheers,

Patrick reinforcement learning MATLAB Answers — New Questions

​

i was trying to make a predictive code algorithm however i recieve the same error every time .

code below.
function [Sa,Sb,Sc] = control(T_ref,sflux_ref,wm,i_meas)
% Variables defined in the parameters file
global Ts Rs Lr Lm Ls p tr kr r_sigma t_sigma lambda v states
persistent x_opt Fs %Fsp1 Tp1 g
if isempty(x_opt), x_opt = 1; end
if isempty(Fs), Fs = 0 + 0i*1; end

T_ref = T_ref(1) + 1j*T_ref(2);
sflux_ref = sflux_ref(1) + 1j*sflux(2);
i_meas = i_meas(1) + 1j*i_meas(2);
wm = wm(1) + 1j*wm(2);
% Stator flux estimate
Fs = Fs + Ts*(v(x_opt) – Rs*i_meas);
% Rotor flux estimate
Fr = Lr/Lm*Fs+i_meas*(Lm-Lr*Ls/Lm);
g_opt = 1e10;
for i = 1:8
% i-th voltage vector for current prediction
v_o1 = v(i);
% Stator flux prediction at instant k+1
Fsp1 = Fs + Ts*v_o1 – Rs*Ts*i_meas;
% Stator current prediction at instant k+1
Isp1 = (1+Ts/t_sigma)*i_meas+Ts/(t_sigma+Ts)*(1/r_sigma*((kr/tr-kr*1i*wm)*Fr+v_o1));
% Torque prediction at instant k+1
Tp1 = 3/2*p*imag(conj(Fsp1)*Isp1);
% Cost function
g = abs(T_ref – Tp1)+ lambda*abs(sflux_ref-abs(Fsp1));
if (g<g_opt)
g_opt = g;
x_opt = i;
end
endi was trying to make a predictive code algorithm however i recieve the same error every time .

code below.
function [Sa,Sb,Sc] = control(T_ref,sflux_ref,wm,i_meas)
% Variables defined in the parameters file
global Ts Rs Lr Lm Ls p tr kr r_sigma t_sigma lambda v states
persistent x_opt Fs %Fsp1 Tp1 g
if isempty(x_opt), x_opt = 1; end
if isempty(Fs), Fs = 0 + 0i*1; end

T_ref = T_ref(1) + 1j*T_ref(2);
sflux_ref = sflux_ref(1) + 1j*sflux(2);
i_meas = i_meas(1) + 1j*i_meas(2);
wm = wm(1) + 1j*wm(2);
% Stator flux estimate
Fs = Fs + Ts*(v(x_opt) – Rs*i_meas);
% Rotor flux estimate
Fr = Lr/Lm*Fs+i_meas*(Lm-Lr*Ls/Lm);
g_opt = 1e10;
for i = 1:8
% i-th voltage vector for current prediction
v_o1 = v(i);
% Stator flux prediction at instant k+1
Fsp1 = Fs + Ts*v_o1 – Rs*Ts*i_meas;
% Stator current prediction at instant k+1
Isp1 = (1+Ts/t_sigma)*i_meas+Ts/(t_sigma+Ts)*(1/r_sigma*((kr/tr-kr*1i*wm)*Fr+v_o1));
% Torque prediction at instant k+1
Tp1 = 3/2*p*imag(conj(Fsp1)*Isp1);
% Cost function
g = abs(T_ref – Tp1)+ lambda*abs(sflux_ref-abs(Fsp1));
if (g<g_opt)
g_opt = g;
x_opt = i;
end
end i was trying to make a predictive code algorithm however i recieve the same error every time .

code below.
function [Sa,Sb,Sc] = control(T_ref,sflux_ref,wm,i_meas)
% Variables defined in the parameters file
global Ts Rs Lr Lm Ls p tr kr r_sigma t_sigma lambda v states
persistent x_opt Fs %Fsp1 Tp1 g
if isempty(x_opt), x_opt = 1; end
if isempty(Fs), Fs = 0 + 0i*1; end

T_ref = T_ref(1) + 1j*T_ref(2);
sflux_ref = sflux_ref(1) + 1j*sflux(2);
i_meas = i_meas(1) + 1j*i_meas(2);
wm = wm(1) + 1j*wm(2);
% Stator flux estimate
Fs = Fs + Ts*(v(x_opt) – Rs*i_meas);
% Rotor flux estimate
Fr = Lr/Lm*Fs+i_meas*(Lm-Lr*Ls/Lm);
g_opt = 1e10;
for i = 1:8
% i-th voltage vector for current prediction
v_o1 = v(i);
% Stator flux prediction at instant k+1
Fsp1 = Fs + Ts*v_o1 – Rs*Ts*i_meas;
% Stator current prediction at instant k+1
Isp1 = (1+Ts/t_sigma)*i_meas+Ts/(t_sigma+Ts)*(1/r_sigma*((kr/tr-kr*1i*wm)*Fr+v_o1));
% Torque prediction at instant k+1
Tp1 = 3/2*p*imag(conj(Fsp1)*Isp1);
% Cost function
g = abs(T_ref – Tp1)+ lambda*abs(sflux_ref-abs(Fsp1));
if (g<g_opt)
g_opt = g;
x_opt = i;
end
end ports, data manager, store memory MATLAB Answers — New Questions

​

Motor not responding after dumping and run the code in TI F28379d LaunchPad using BOOSTXL-3PhGaNInv booster pack also chacked the PMW pin no 37-40 but did not get any signal. Is there any way to observe the PWM signal from Lunchpad without using the motor to validate the code generation is working.Motor not responding after dumping and run the code in TI F28379d LaunchPad using BOOSTXL-3PhGaNInv booster pack also chacked the PMW pin no 37-40 but did not get any signal. Is there any way to observe the PWM signal from Lunchpad without using the motor to validate the code generation is working. Motor not responding after dumping and run the code in TI F28379d LaunchPad using BOOSTXL-3PhGaNInv booster pack also chacked the PMW pin no 37-40 but did not get any signal. Is there any way to observe the PWM signal from Lunchpad without using the motor to validate the code generation is working. pmsm, motor control, ti f28379d launchpad, boostxl-3phganinv MATLAB Answers — New Questions

​

I have Navisys GR-401 gps receiver which works just fine with SiRF-Live. I know it uses COM3, but I am unable to read data with matlab.
with:
s = serialport(‘COM3’,9600);
gps = gpsdev(s,’OutputFormat’,"timetable");
and
read(gps)
gives error:
Obtaining the required NMEA sentences (RMC, GGA and GSA) is taking longer than expected. Verify the physical connection of the GPS device to the specified serial port and the configuration of the GPS device.
also
gps = gpsReceiver(‘COM3’);
and
read(gps)
gives error
Warning: Unsuccessful read: The input buffer was filled before the Terminator was reached.
‘serial’ unable to read all requested data. For more information on possible reasons, see Serial Read Warnings.

any ideas what is the culprit here?I have Navisys GR-401 gps receiver which works just fine with SiRF-Live. I know it uses COM3, but I am unable to read data with matlab.
with:
s = serialport(‘COM3’,9600);
gps = gpsdev(s,’OutputFormat’,"timetable");
and
read(gps)
gives error:
Obtaining the required NMEA sentences (RMC, GGA and GSA) is taking longer than expected. Verify the physical connection of the GPS device to the specified serial port and the configuration of the GPS device.
also
gps = gpsReceiver(‘COM3’);
and
read(gps)
gives error
Warning: Unsuccessful read: The input buffer was filled before the Terminator was reached.
‘serial’ unable to read all requested data. For more information on possible reasons, see Serial Read Warnings.

any ideas what is the culprit here? I have Navisys GR-401 gps receiver which works just fine with SiRF-Live. I know it uses COM3, but I am unable to read data with matlab.
with:
s = serialport(‘COM3’,9600);
gps = gpsdev(s,’OutputFormat’,"timetable");
and
read(gps)
gives error:
Obtaining the required NMEA sentences (RMC, GGA and GSA) is taking longer than expected. Verify the physical connection of the GPS device to the specified serial port and the configuration of the GPS device.
also
gps = gpsReceiver(‘COM3’);
and
read(gps)
gives error
Warning: Unsuccessful read: The input buffer was filled before the Terminator was reached.
‘serial’ unable to read all requested data. For more information on possible reasons, see Serial Read Warnings.

any ideas what is the culprit here? gps, gr-401 MATLAB Answers — New Questions

​

I have a logical array that is 100×1 logical. sometimes it contains a mix of numbers and sometimes it contains all 0s.

How would i go about making an if statement that checks to see if all 100×1 is 0s and if all 0s = do something
I thought i was correct in using any command but im not sure that works correctly for me..

Thanks!I have a logical array that is 100×1 logical. sometimes it contains a mix of numbers and sometimes it contains all 0s.

How would i go about making an if statement that checks to see if all 100×1 is 0s and if all 0s = do something
I thought i was correct in using any command but im not sure that works correctly for me..

Thanks! I have a logical array that is 100×1 logical. sometimes it contains a mix of numbers and sometimes it contains all 0s.

How would i go about making an if statement that checks to see if all 100×1 is 0s and if all 0s = do something
I thought i was correct in using any command but im not sure that works correctly for me..

Thanks! array, logic, if statement MATLAB Answers — New Questions

​

Hi all, what i want is a code that can actually place random carbon nanotubes in a cube domain so that i can use it in comsol livelink MATLAB bellow is a code that puts spheres which i found in this forum i need it based on that.. :
function P = GetRandomSpheres(nWant, Width, Radius)
% INPUT:
% nWant: Number of spheres
% Width: Dimension of 3d box as [1 x 3] double vector
% Radius: Radius of spheres
% OUTPUT:
% P: [nWant x 3] matrix, centers

P = zeros(nWant, 3);

R2 = (2 * Radius) ^ 2; % Squared once instead of SQRT each time
W = Width – 2 * Radius; % Avoid interesction with borders
iLoop = 1; % Security break to avoid infinite loop
nValid = 0;
while nValid < nWant && iLoop < 1e12
newP = rand(1, 3) .* W + Radius;
% Auto-expanding, need Matlab >= R2016b. For earlier versions:
% Dist2 = sum(bsxfun(@minus, P(1:nValid, :), newP) .^ 2, 2);
Dist2 = sum((P(1:nValid, 🙂 – newP) .^ 2, 2);
if all(Dist2 > R2)
% Success: The new point does not touch existing sheres:
nValid = nValid + 1; % Append this point
P(nValid, 🙂 = newP;
end
iLoop = iLoop + 1;
end
% Stop if too few values have been found:
if nValid < nWant
error(‘Cannot find wanted number of points in %d iterations.’, iLoop)
end
end

n = 100;
R = 5;
P = GetRandomSpheres(n, [100, 100, 100], R);
figure
axes(‘NextPlot’, ‘add’, …
‘XLim’, [0, 100], ‘YLim’, [0, 100], ‘ZLim’, [0, 100]);
view(3);
[X, Y, Z] = sphere();
for k = 1:n
surf(X * R + P(k, 1), Y * R + P(k, 2), Z * R + P(k, 3));
endHi all, what i want is a code that can actually place random carbon nanotubes in a cube domain so that i can use it in comsol livelink MATLAB bellow is a code that puts spheres which i found in this forum i need it based on that.. :
function P = GetRandomSpheres(nWant, Width, Radius)
% INPUT:
% nWant: Number of spheres
% Width: Dimension of 3d box as [1 x 3] double vector
% Radius: Radius of spheres
% OUTPUT:
% P: [nWant x 3] matrix, centers

P = zeros(nWant, 3);

R2 = (2 * Radius) ^ 2; % Squared once instead of SQRT each time
W = Width – 2 * Radius; % Avoid interesction with borders
iLoop = 1; % Security break to avoid infinite loop
nValid = 0;
while nValid < nWant && iLoop < 1e12
newP = rand(1, 3) .* W + Radius;
% Auto-expanding, need Matlab >= R2016b. For earlier versions:
% Dist2 = sum(bsxfun(@minus, P(1:nValid, :), newP) .^ 2, 2);
Dist2 = sum((P(1:nValid, 🙂 – newP) .^ 2, 2);
if all(Dist2 > R2)
% Success: The new point does not touch existing sheres:
nValid = nValid + 1; % Append this point
P(nValid, 🙂 = newP;
end
iLoop = iLoop + 1;
end
% Stop if too few values have been found:
if nValid < nWant
error(‘Cannot find wanted number of points in %d iterations.’, iLoop)
end
end

n = 100;
R = 5;
P = GetRandomSpheres(n, [100, 100, 100], R);
figure
axes(‘NextPlot’, ‘add’, …
‘XLim’, [0, 100], ‘YLim’, [0, 100], ‘ZLim’, [0, 100]);
view(3);
[X, Y, Z] = sphere();
for k = 1:n
surf(X * R + P(k, 1), Y * R + P(k, 2), Z * R + P(k, 3));
end Hi all, what i want is a code that can actually place random carbon nanotubes in a cube domain so that i can use it in comsol livelink MATLAB bellow is a code that puts spheres which i found in this forum i need it based on that.. :
function P = GetRandomSpheres(nWant, Width, Radius)
% INPUT:
% nWant: Number of spheres
% Width: Dimension of 3d box as [1 x 3] double vector
% Radius: Radius of spheres
% OUTPUT:
% P: [nWant x 3] matrix, centers

P = zeros(nWant, 3);

R2 = (2 * Radius) ^ 2; % Squared once instead of SQRT each time
W = Width – 2 * Radius; % Avoid interesction with borders
iLoop = 1; % Security break to avoid infinite loop
nValid = 0;
while nValid < nWant && iLoop < 1e12
newP = rand(1, 3) .* W + Radius;
% Auto-expanding, need Matlab >= R2016b. For earlier versions:
% Dist2 = sum(bsxfun(@minus, P(1:nValid, :), newP) .^ 2, 2);
Dist2 = sum((P(1:nValid, 🙂 – newP) .^ 2, 2);
if all(Dist2 > R2)
% Success: The new point does not touch existing sheres:
nValid = nValid + 1; % Append this point
P(nValid, 🙂 = newP;
end
iLoop = iLoop + 1;
end
% Stop if too few values have been found:
if nValid < nWant
error(‘Cannot find wanted number of points in %d iterations.’, iLoop)
end
end

n = 100;
R = 5;
P = GetRandomSpheres(n, [100, 100, 100], R);
figure
axes(‘NextPlot’, ‘add’, …
‘XLim’, [0, 100], ‘YLim’, [0, 100], ‘ZLim’, [0, 100]);
view(3);
[X, Y, Z] = sphere();
for k = 1:n
surf(X * R + P(k, 1), Y * R + P(k, 2), Z * R + P(k, 3));
end matlab, random, java, comsol, geometry, matlab code, programming MATLAB Answers — New Questions

​