Tag Archives: matlab
Matlab 2017b No Compilers Detected
I am trying to compile a simulink model I made to check the I/O on my Speedgoat but when I go to do mex -setup Matlab replies with "Error using mex, No supported compiler or SDK was found." This is confusing to me since I just installed visual studios per https://www.mathworks.com/matlabcentral/answers/348269-how-do-i-set-up-microsoft-visual-studio-2017-for-slrt. I have tried adding the Visual Studio to the Matlab file path but it has not fixed the problem. Has anyone encountered and solved this problem?I am trying to compile a simulink model I made to check the I/O on my Speedgoat but when I go to do mex -setup Matlab replies with "Error using mex, No supported compiler or SDK was found." This is confusing to me since I just installed visual studios per https://www.mathworks.com/matlabcentral/answers/348269-how-do-i-set-up-microsoft-visual-studio-2017-for-slrt. I have tried adding the Visual Studio to the Matlab file path but it has not fixed the problem. Has anyone encountered and solved this problem? I am trying to compile a simulink model I made to check the I/O on my Speedgoat but when I go to do mex -setup Matlab replies with "Error using mex, No supported compiler or SDK was found." This is confusing to me since I just installed visual studios per https://www.mathworks.com/matlabcentral/answers/348269-how-do-i-set-up-microsoft-visual-studio-2017-for-slrt. I have tried adding the Visual Studio to the Matlab file path but it has not fixed the problem. Has anyone encountered and solved this problem? mex compiler MATLAB Answers — New Questions
Matlab ODE Solvers – Paramater Question
Greetings,
I working with a system of differential equations that model amount of components as a reaction proceeds. I want to pass a variable parameter through one of matlab’s ode solvers (probably 45/15). This is in addition to time bounds, so it would be like ode45(function, times, [X1(0),X2(0),…X5(0)],… ). I want to find λ given something like X1(final) = 0.35. How to include this information?
Thank you (:Greetings,
I working with a system of differential equations that model amount of components as a reaction proceeds. I want to pass a variable parameter through one of matlab’s ode solvers (probably 45/15). This is in addition to time bounds, so it would be like ode45(function, times, [X1(0),X2(0),…X5(0)],… ). I want to find λ given something like X1(final) = 0.35. How to include this information?
Thank you (: Greetings,
I working with a system of differential equations that model amount of components as a reaction proceeds. I want to pass a variable parameter through one of matlab’s ode solvers (probably 45/15). This is in addition to time bounds, so it would be like ode45(function, times, [X1(0),X2(0),…X5(0)],… ). I want to find λ given something like X1(final) = 0.35. How to include this information?
Thank you (: ode45, differential equations, system, ode15 MATLAB Answers — New Questions
API를 사용하여 젠킨스에서 CLI를 구성하여 아래 f_PSGen.m을 실행하였습니다. 오류가 발생하여, 오류분석 부탁드립니다.
■■■■■■■■■■■■■■■■■■■■■ 소스 코드 ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
% f_PSGen
% Abstract: function for Polyspace for Simulink Model Only
%
%%% Syntax
% * [chrBFResult,chrCPResult]=f_SLPSGen(chrModelPath, chrProjectPath)
%
%%% Example
% [chrBFResult,chrCPResult]=f_SLPSGen(chrModelPath, chrProjectPath)
%
%%% Input
% * chrModelPath : Simulink Model Path
% * chrProjectPath : Simulink Project Path
%
%%% Output
% * chrBFResult : Bug Finder Result Path
% * chrCPResult : Code Prover Result Path
%
%%% Copyright
% ————————————————————————-
% Copyright 2023 Consulting Services, The MathWorks, Inc.
% ————————————————————————-
function [chrBFResult,chrCPResult]=f_PSGen(chrModelPath, chrProjectPath,chrLegacyCPath,chrTarget,chrCompiler)
arguments
chrModelPath (1,:) char {mustBeText}
chrProjectPath (1,:) char {mustBeText}
chrLegacyCPath (1,:) char {mustBeText}
chrTarget (1,:) char {mustBeText}
chrCompiler (1,:) char {mustBeText}
end
%Get Current Path for Code Verification
strCurrentPath = fileparts(mfilename(‘fullpath’));
fprintf(‘[PJK] Current path: %sn’, strCurrentPath);
%Incase Simulink model , then set up the code generation folder
if( exist(char(chrModelPath),’file’) && exist(char(chrProjectPath),’file’))
%% Load Simulink Project to get Cach and CodeGen Folder
hdlSimulinkProject = openProject(char(chrProjectPath));
chrCodeGenFolder = hdlSimulinkProject.SimulinkCodeGenFolder;
chrCacheFolder = hdlSimulinkProject.SimulinkCacheFolder;
fprintf(‘[PJK] Simulink project loaded. CodeGenFolder: %s, CacheFolder: %sn’, chrCodeGenFolder, chrCacheFolder);
%%Set Reference Project for Simulink Model
hdlProject = openProject(strCurrentPath);
addReference(hdlProject,fileparts(char(chrProjectPath)),"relative");
hdlProject.SimulinkCacheFolder = chrCacheFolder;
hdlProject.SimulinkCodeGenFolder = chrCodeGenFolder;
reload(hdlProject);
fprintf(‘[PJK] Reference project set. SimulinkCacheFolder: %s, SimulinkCodeGenFolder: %sn’, hdlProject.SimulinkCacheFolder, hdlProject.SimulinkCodeGenFolder);
%Check CodeGen Folder Existing
if( ~exist(chrCodeGenFolder,’dir’) )
chrModelPath = [];
end
else
chrModelPath = [];
end
%% Checking Legacy C Code Path Existing
if( ~exist(char(chrLegacyCPath),’dir’) )
chrLegacyCPath = [];
end
%%Create Polyspace Project & Run
hPSProject = Polyspace.PSProject(chrModelPath,chrLegacyCPath,chrTarget,chrCompiler);
fprintf(‘[PJK] Polyspace project created. ModelPath: %s, LegacyCPath: %s, Target: %s, Compiler: %sn’, chrModelPath, chrLegacyCPath, chrTarget, chrCompiler);
%%Run Polyspace for Bug Finder
[chrBFResult,chrCPResult] = hPSProject.run(Polyspace.PSMode.AllMode);
fprintf(‘[PJK] Polyspace project run. BFResult: %s, CPResult: %sn’, chrBFResult, chrCPResult);
%delete handle
hPSProject.delete();
fprintf(‘[PJK] Polyspace project deleted.n’);
if( exist(char(chrModelPath),’file’) && exist(char(chrProjectPath),’file’))
removeReference(hdlProject,fileparts(char(chrProjectPath)));
hdlProject.SimulinkCacheFolder = ”;
hdlProject.SimulinkCodeGenFolder = ”;
fprintf(‘[PJK] Reference removed. SimulinkCacheFolder: %s, SimulinkCodeGenFolder: %sn’, hdlProject.SimulinkCacheFolder, hdlProject.SimulinkCodeGenFolder);
end
end
■■■■■■■■■■■■■■■■■■■■■ Error Log ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
Generating MATLAB script with content:
cd ‘D:t2-master-work-ModelworkModel’;runAll;
#################### Starting command output ####################
[PolySpace_ICU_SU2idL_FL-ICU-HMID-23-07] $ cmd.exe /C .matlabYJKDTKesrun_matlab_command.bat "cd(‘D:JenkinsworkspacePolySpace_ICU_SU2idL_FL-ICU-HMID-23-07/.matlab/YJKDTKes’);command_lKfQMwWL"
———This MATLAB R2021a Version Supported ————-
[PJK] Current path: D:t2-master-work-ModelworkModelUtilitiesM08
[PJK] Polyspace project created. ModelPath: , LegacyCPath: D:t2-master-work-ModelworkModelLegacyCProjectSU2idL_FL-ICU-HMID-23-07, Target: powerpc, Compiler: greenhills
Picked up JAVA_TOOL_OPTIONS: -Dfile.encoding=UTF-8
{���� ��� �� ������ ����: rmdir
������ ���ŵ��� �ʾҽ��ϴ�.
���� ��: Polyspace.PSProject/saveLegacy (45�� ����)
rmdir(chrFolderPath,’s’);
���� ��: Polyspace.PSProject/run (63�� ����)
obj.saveLegacy(true);
���� ��: f_PSGen (68�� ����)
[chrBFResult,chrCPResult] = hPSProject.run(Polyspace.PSMode.AllMode);
���� ��: PSGen (53�� ����)
f_PSGen(tblPSList.ModelPath{index},
tblPSList.ProjectPath{index},tblPSList.LegacyCPath{index},tblPSList.Target{index},tblPSList.Compiler{index}
);
���� ��: runAll (8�� ����)
PSGen;
���� ��: command_lKfQMwWL (2�� ����)
cd ‘D:t2-master-work-ModelworkModel’;runAll;
}
Picked up JAVA_TOOL_OPTIONS: -Dfile.encoding=UTF-8
ERROR: MATLAB error Exit Status: 0x00000001
Build step ‘Run MATLAB Command’ changed build result to FAILURE
[PolySpace_ICU_SU2idL_FL-ICU-HMID-23-07] $ cmd /c call C:UsersADMINI~1AppDataLocalTempjenkins10177263456194605002.bat■■■■■■■■■■■■■■■■■■■■■ 소스 코드 ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
% f_PSGen
% Abstract: function for Polyspace for Simulink Model Only
%
%%% Syntax
% * [chrBFResult,chrCPResult]=f_SLPSGen(chrModelPath, chrProjectPath)
%
%%% Example
% [chrBFResult,chrCPResult]=f_SLPSGen(chrModelPath, chrProjectPath)
%
%%% Input
% * chrModelPath : Simulink Model Path
% * chrProjectPath : Simulink Project Path
%
%%% Output
% * chrBFResult : Bug Finder Result Path
% * chrCPResult : Code Prover Result Path
%
%%% Copyright
% ————————————————————————-
% Copyright 2023 Consulting Services, The MathWorks, Inc.
% ————————————————————————-
function [chrBFResult,chrCPResult]=f_PSGen(chrModelPath, chrProjectPath,chrLegacyCPath,chrTarget,chrCompiler)
arguments
chrModelPath (1,:) char {mustBeText}
chrProjectPath (1,:) char {mustBeText}
chrLegacyCPath (1,:) char {mustBeText}
chrTarget (1,:) char {mustBeText}
chrCompiler (1,:) char {mustBeText}
end
%Get Current Path for Code Verification
strCurrentPath = fileparts(mfilename(‘fullpath’));
fprintf(‘[PJK] Current path: %sn’, strCurrentPath);
%Incase Simulink model , then set up the code generation folder
if( exist(char(chrModelPath),’file’) && exist(char(chrProjectPath),’file’))
%% Load Simulink Project to get Cach and CodeGen Folder
hdlSimulinkProject = openProject(char(chrProjectPath));
chrCodeGenFolder = hdlSimulinkProject.SimulinkCodeGenFolder;
chrCacheFolder = hdlSimulinkProject.SimulinkCacheFolder;
fprintf(‘[PJK] Simulink project loaded. CodeGenFolder: %s, CacheFolder: %sn’, chrCodeGenFolder, chrCacheFolder);
%%Set Reference Project for Simulink Model
hdlProject = openProject(strCurrentPath);
addReference(hdlProject,fileparts(char(chrProjectPath)),"relative");
hdlProject.SimulinkCacheFolder = chrCacheFolder;
hdlProject.SimulinkCodeGenFolder = chrCodeGenFolder;
reload(hdlProject);
fprintf(‘[PJK] Reference project set. SimulinkCacheFolder: %s, SimulinkCodeGenFolder: %sn’, hdlProject.SimulinkCacheFolder, hdlProject.SimulinkCodeGenFolder);
%Check CodeGen Folder Existing
if( ~exist(chrCodeGenFolder,’dir’) )
chrModelPath = [];
end
else
chrModelPath = [];
end
%% Checking Legacy C Code Path Existing
if( ~exist(char(chrLegacyCPath),’dir’) )
chrLegacyCPath = [];
end
%%Create Polyspace Project & Run
hPSProject = Polyspace.PSProject(chrModelPath,chrLegacyCPath,chrTarget,chrCompiler);
fprintf(‘[PJK] Polyspace project created. ModelPath: %s, LegacyCPath: %s, Target: %s, Compiler: %sn’, chrModelPath, chrLegacyCPath, chrTarget, chrCompiler);
%%Run Polyspace for Bug Finder
[chrBFResult,chrCPResult] = hPSProject.run(Polyspace.PSMode.AllMode);
fprintf(‘[PJK] Polyspace project run. BFResult: %s, CPResult: %sn’, chrBFResult, chrCPResult);
%delete handle
hPSProject.delete();
fprintf(‘[PJK] Polyspace project deleted.n’);
if( exist(char(chrModelPath),’file’) && exist(char(chrProjectPath),’file’))
removeReference(hdlProject,fileparts(char(chrProjectPath)));
hdlProject.SimulinkCacheFolder = ”;
hdlProject.SimulinkCodeGenFolder = ”;
fprintf(‘[PJK] Reference removed. SimulinkCacheFolder: %s, SimulinkCodeGenFolder: %sn’, hdlProject.SimulinkCacheFolder, hdlProject.SimulinkCodeGenFolder);
end
end
■■■■■■■■■■■■■■■■■■■■■ Error Log ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
Generating MATLAB script with content:
cd ‘D:t2-master-work-ModelworkModel’;runAll;
#################### Starting command output ####################
[PolySpace_ICU_SU2idL_FL-ICU-HMID-23-07] $ cmd.exe /C .matlabYJKDTKesrun_matlab_command.bat "cd(‘D:JenkinsworkspacePolySpace_ICU_SU2idL_FL-ICU-HMID-23-07/.matlab/YJKDTKes’);command_lKfQMwWL"
———This MATLAB R2021a Version Supported ————-
[PJK] Current path: D:t2-master-work-ModelworkModelUtilitiesM08
[PJK] Polyspace project created. ModelPath: , LegacyCPath: D:t2-master-work-ModelworkModelLegacyCProjectSU2idL_FL-ICU-HMID-23-07, Target: powerpc, Compiler: greenhills
Picked up JAVA_TOOL_OPTIONS: -Dfile.encoding=UTF-8
{���� ��� �� ������ ����: rmdir
������ ���ŵ��� �ʾҽ��ϴ�.
���� ��: Polyspace.PSProject/saveLegacy (45�� ����)
rmdir(chrFolderPath,’s’);
���� ��: Polyspace.PSProject/run (63�� ����)
obj.saveLegacy(true);
���� ��: f_PSGen (68�� ����)
[chrBFResult,chrCPResult] = hPSProject.run(Polyspace.PSMode.AllMode);
���� ��: PSGen (53�� ����)
f_PSGen(tblPSList.ModelPath{index},
tblPSList.ProjectPath{index},tblPSList.LegacyCPath{index},tblPSList.Target{index},tblPSList.Compiler{index}
);
���� ��: runAll (8�� ����)
PSGen;
���� ��: command_lKfQMwWL (2�� ����)
cd ‘D:t2-master-work-ModelworkModel’;runAll;
}
Picked up JAVA_TOOL_OPTIONS: -Dfile.encoding=UTF-8
ERROR: MATLAB error Exit Status: 0x00000001
Build step ‘Run MATLAB Command’ changed build result to FAILURE
[PolySpace_ICU_SU2idL_FL-ICU-HMID-23-07] $ cmd /c call C:UsersADMINI~1AppDataLocalTempjenkins10177263456194605002.bat ■■■■■■■■■■■■■■■■■■■■■ 소스 코드 ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
% f_PSGen
% Abstract: function for Polyspace for Simulink Model Only
%
%%% Syntax
% * [chrBFResult,chrCPResult]=f_SLPSGen(chrModelPath, chrProjectPath)
%
%%% Example
% [chrBFResult,chrCPResult]=f_SLPSGen(chrModelPath, chrProjectPath)
%
%%% Input
% * chrModelPath : Simulink Model Path
% * chrProjectPath : Simulink Project Path
%
%%% Output
% * chrBFResult : Bug Finder Result Path
% * chrCPResult : Code Prover Result Path
%
%%% Copyright
% ————————————————————————-
% Copyright 2023 Consulting Services, The MathWorks, Inc.
% ————————————————————————-
function [chrBFResult,chrCPResult]=f_PSGen(chrModelPath, chrProjectPath,chrLegacyCPath,chrTarget,chrCompiler)
arguments
chrModelPath (1,:) char {mustBeText}
chrProjectPath (1,:) char {mustBeText}
chrLegacyCPath (1,:) char {mustBeText}
chrTarget (1,:) char {mustBeText}
chrCompiler (1,:) char {mustBeText}
end
%Get Current Path for Code Verification
strCurrentPath = fileparts(mfilename(‘fullpath’));
fprintf(‘[PJK] Current path: %sn’, strCurrentPath);
%Incase Simulink model , then set up the code generation folder
if( exist(char(chrModelPath),’file’) && exist(char(chrProjectPath),’file’))
%% Load Simulink Project to get Cach and CodeGen Folder
hdlSimulinkProject = openProject(char(chrProjectPath));
chrCodeGenFolder = hdlSimulinkProject.SimulinkCodeGenFolder;
chrCacheFolder = hdlSimulinkProject.SimulinkCacheFolder;
fprintf(‘[PJK] Simulink project loaded. CodeGenFolder: %s, CacheFolder: %sn’, chrCodeGenFolder, chrCacheFolder);
%%Set Reference Project for Simulink Model
hdlProject = openProject(strCurrentPath);
addReference(hdlProject,fileparts(char(chrProjectPath)),"relative");
hdlProject.SimulinkCacheFolder = chrCacheFolder;
hdlProject.SimulinkCodeGenFolder = chrCodeGenFolder;
reload(hdlProject);
fprintf(‘[PJK] Reference project set. SimulinkCacheFolder: %s, SimulinkCodeGenFolder: %sn’, hdlProject.SimulinkCacheFolder, hdlProject.SimulinkCodeGenFolder);
%Check CodeGen Folder Existing
if( ~exist(chrCodeGenFolder,’dir’) )
chrModelPath = [];
end
else
chrModelPath = [];
end
%% Checking Legacy C Code Path Existing
if( ~exist(char(chrLegacyCPath),’dir’) )
chrLegacyCPath = [];
end
%%Create Polyspace Project & Run
hPSProject = Polyspace.PSProject(chrModelPath,chrLegacyCPath,chrTarget,chrCompiler);
fprintf(‘[PJK] Polyspace project created. ModelPath: %s, LegacyCPath: %s, Target: %s, Compiler: %sn’, chrModelPath, chrLegacyCPath, chrTarget, chrCompiler);
%%Run Polyspace for Bug Finder
[chrBFResult,chrCPResult] = hPSProject.run(Polyspace.PSMode.AllMode);
fprintf(‘[PJK] Polyspace project run. BFResult: %s, CPResult: %sn’, chrBFResult, chrCPResult);
%delete handle
hPSProject.delete();
fprintf(‘[PJK] Polyspace project deleted.n’);
if( exist(char(chrModelPath),’file’) && exist(char(chrProjectPath),’file’))
removeReference(hdlProject,fileparts(char(chrProjectPath)));
hdlProject.SimulinkCacheFolder = ”;
hdlProject.SimulinkCodeGenFolder = ”;
fprintf(‘[PJK] Reference removed. SimulinkCacheFolder: %s, SimulinkCodeGenFolder: %sn’, hdlProject.SimulinkCacheFolder, hdlProject.SimulinkCodeGenFolder);
end
end
■■■■■■■■■■■■■■■■■■■■■ Error Log ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■
Generating MATLAB script with content:
cd ‘D:t2-master-work-ModelworkModel’;runAll;
#################### Starting command output ####################
[PolySpace_ICU_SU2idL_FL-ICU-HMID-23-07] $ cmd.exe /C .matlabYJKDTKesrun_matlab_command.bat "cd(‘D:JenkinsworkspacePolySpace_ICU_SU2idL_FL-ICU-HMID-23-07/.matlab/YJKDTKes’);command_lKfQMwWL"
———This MATLAB R2021a Version Supported ————-
[PJK] Current path: D:t2-master-work-ModelworkModelUtilitiesM08
[PJK] Polyspace project created. ModelPath: , LegacyCPath: D:t2-master-work-ModelworkModelLegacyCProjectSU2idL_FL-ICU-HMID-23-07, Target: powerpc, Compiler: greenhills
Picked up JAVA_TOOL_OPTIONS: -Dfile.encoding=UTF-8
{���� ��� �� ������ ����: rmdir
������ ���ŵ��� �ʾҽ��ϴ�.
���� ��: Polyspace.PSProject/saveLegacy (45�� ����)
rmdir(chrFolderPath,’s’);
���� ��: Polyspace.PSProject/run (63�� ����)
obj.saveLegacy(true);
���� ��: f_PSGen (68�� ����)
[chrBFResult,chrCPResult] = hPSProject.run(Polyspace.PSMode.AllMode);
���� ��: PSGen (53�� ����)
f_PSGen(tblPSList.ModelPath{index},
tblPSList.ProjectPath{index},tblPSList.LegacyCPath{index},tblPSList.Target{index},tblPSList.Compiler{index}
);
���� ��: runAll (8�� ����)
PSGen;
���� ��: command_lKfQMwWL (2�� ����)
cd ‘D:t2-master-work-ModelworkModel’;runAll;
}
Picked up JAVA_TOOL_OPTIONS: -Dfile.encoding=UTF-8
ERROR: MATLAB error Exit Status: 0x00000001
Build step ‘Run MATLAB Command’ changed build result to FAILURE
[PolySpace_ICU_SU2idL_FL-ICU-HMID-23-07] $ cmd /c call C:UsersADMINI~1AppDataLocalTempjenkins10177263456194605002.bat 젠킨스, cli, 오류 MATLAB Answers — New Questions
How to draw Domino tiles in a figure?
Hi everyone,
I’m here to ask you if anyone knows, how to draw a bunch of domino tiles in a figure. I already generated a figure with the pixels black, but know I want to draw domino tiles into this figure.
And basically I want a figure something like this: (the color os the domino tiles can be just one).
Thanks,
NunoHi everyone,
I’m here to ask you if anyone knows, how to draw a bunch of domino tiles in a figure. I already generated a figure with the pixels black, but know I want to draw domino tiles into this figure.
And basically I want a figure something like this: (the color os the domino tiles can be just one).
Thanks,
Nuno Hi everyone,
I’m here to ask you if anyone knows, how to draw a bunch of domino tiles in a figure. I already generated a figure with the pixels black, but know I want to draw domino tiles into this figure.
And basically I want a figure something like this: (the color os the domino tiles can be just one).
Thanks,
Nuno image processing MATLAB Answers — New Questions
Can anyone help me in minimizing the following integral equation and finding the roots ?
Here i am posting an integral equation for which i want the minimization code and methods.
Ω =λ^ 4 (σ 2 + ~π2 − v 2 ) 2 − cσ + νqT∫ (d ^3k)/ (2π) ^3 {ln(1 − nq(T, µ)) + ln(1 − nq¯(T, µ))} equation(1)
we have to minimize the above equation with respect to sigma(σ) and in the scond term,there is integration on the momentum vector k, where integration term also contain the sigma variable in the following ways
nq(T, µ) = 1/ (1 + exp((sqrt(k^2+g^2(σ^2+π^2)) − µ)/T ))
nq¯(T, µ))= 1/ (1 + exp((sqrt(k^2+g^2(σ^2+π^2)) +µ)/T ))
where all the other notation like pi,lambda,vq,T are the some parametrs which we can fix like
lambda=87.6 Mev
T = 500(temperature)
vq=12
g=3.8
mu=100
c=13.6
pi=0
so, we have to find ∂Ω/∂σ and corresponding value of sigma at which above equation (1) is minimized .Here i am posting an integral equation for which i want the minimization code and methods.
Ω =λ^ 4 (σ 2 + ~π2 − v 2 ) 2 − cσ + νqT∫ (d ^3k)/ (2π) ^3 {ln(1 − nq(T, µ)) + ln(1 − nq¯(T, µ))} equation(1)
we have to minimize the above equation with respect to sigma(σ) and in the scond term,there is integration on the momentum vector k, where integration term also contain the sigma variable in the following ways
nq(T, µ) = 1/ (1 + exp((sqrt(k^2+g^2(σ^2+π^2)) − µ)/T ))
nq¯(T, µ))= 1/ (1 + exp((sqrt(k^2+g^2(σ^2+π^2)) +µ)/T ))
where all the other notation like pi,lambda,vq,T are the some parametrs which we can fix like
lambda=87.6 Mev
T = 500(temperature)
vq=12
g=3.8
mu=100
c=13.6
pi=0
so, we have to find ∂Ω/∂σ and corresponding value of sigma at which above equation (1) is minimized . Here i am posting an integral equation for which i want the minimization code and methods.
Ω =λ^ 4 (σ 2 + ~π2 − v 2 ) 2 − cσ + νqT∫ (d ^3k)/ (2π) ^3 {ln(1 − nq(T, µ)) + ln(1 − nq¯(T, µ))} equation(1)
we have to minimize the above equation with respect to sigma(σ) and in the scond term,there is integration on the momentum vector k, where integration term also contain the sigma variable in the following ways
nq(T, µ) = 1/ (1 + exp((sqrt(k^2+g^2(σ^2+π^2)) − µ)/T ))
nq¯(T, µ))= 1/ (1 + exp((sqrt(k^2+g^2(σ^2+π^2)) +µ)/T ))
where all the other notation like pi,lambda,vq,T are the some parametrs which we can fix like
lambda=87.6 Mev
T = 500(temperature)
vq=12
g=3.8
mu=100
c=13.6
pi=0
so, we have to find ∂Ω/∂σ and corresponding value of sigma at which above equation (1) is minimized . optimization of equation containing integration MATLAB Answers — New Questions
Is FORMAT DEBUG still any useful?
Can someone has an explanation for this (R2020B)
clear
clc
x=1;
y=x;
x=x+1;
format debug
y
x
that produces
y =
Structure address = 23ebc936f80
m = 1
n = 1
pr = 23edbdd6660
1
x =
Structure address = 23ebc936f80
m = 1
n = 1
pr = 23edbdd6660
2
So x and y have the same structure address and data address, yet data are different (1 for y and 2 for x).
Finish the time where FORMAT DEBUG returns useful information for us.Can someone has an explanation for this (R2020B)
clear
clc
x=1;
y=x;
x=x+1;
format debug
y
x
that produces
y =
Structure address = 23ebc936f80
m = 1
n = 1
pr = 23edbdd6660
1
x =
Structure address = 23ebc936f80
m = 1
n = 1
pr = 23edbdd6660
2
So x and y have the same structure address and data address, yet data are different (1 for y and 2 for x).
Finish the time where FORMAT DEBUG returns useful information for us. Can someone has an explanation for this (R2020B)
clear
clc
x=1;
y=x;
x=x+1;
format debug
y
x
that produces
y =
Structure address = 23ebc936f80
m = 1
n = 1
pr = 23edbdd6660
1
x =
Structure address = 23ebc936f80
m = 1
n = 1
pr = 23edbdd6660
2
So x and y have the same structure address and data address, yet data are different (1 for y and 2 for x).
Finish the time where FORMAT DEBUG returns useful information for us. internal structure, data sharing, format debug MATLAB Answers — New Questions
Optimal allocation of DG in IEEE 33 bus system
I need the code for optimal placement and sizing of DGs for ieee 33 bus system using optimization algorithmI need the code for optimal placement and sizing of DGs for ieee 33 bus system using optimization algorithm I need the code for optimal placement and sizing of DGs for ieee 33 bus system using optimization algorithm pso, gwo, optimal dg allocation MATLAB Answers — New Questions
HOW TO CREATE TSUNAMI MODEL
be discovered:
A = 5
X = 0 – 600
t = 0 – 60 minute
long = 0,5 – 30 kilometers
count the lamda, sigma and plot psibe discovered:
A = 5
X = 0 – 600
t = 0 – 60 minute
long = 0,5 – 30 kilometers
count the lamda, sigma and plot psi be discovered:
A = 5
X = 0 – 600
t = 0 – 60 minute
long = 0,5 – 30 kilometers
count the lamda, sigma and plot psi #model #code #tsunami MATLAB Answers — New Questions
I have some questions about the minibatchsize attribute of PPO+LSTM
Hello everyone
I found this sentence when I searched the content of PPO+LSTM in the official MATLAB documentation:“For a PPO agent, the trajectory length is the MiniBatchSize property of its options object” ,But I have some doubts about this statement。When using PPO+LSTM, does the agent no longer need to sample minbatch from the current experience sequence? How to understand "For a PPO agent, the trajectory length is the MiniBatchSize property of its options object".
Hope to get your answerHello everyone
I found this sentence when I searched the content of PPO+LSTM in the official MATLAB documentation:“For a PPO agent, the trajectory length is the MiniBatchSize property of its options object” ,But I have some doubts about this statement。When using PPO+LSTM, does the agent no longer need to sample minbatch from the current experience sequence? How to understand "For a PPO agent, the trajectory length is the MiniBatchSize property of its options object".
Hope to get your answer Hello everyone
I found this sentence when I searched the content of PPO+LSTM in the official MATLAB documentation:“For a PPO agent, the trajectory length is the MiniBatchSize property of its options object” ,But I have some doubts about this statement。When using PPO+LSTM, does the agent no longer need to sample minbatch from the current experience sequence? How to understand "For a PPO agent, the trajectory length is the MiniBatchSize property of its options object".
Hope to get your answer drl MATLAB Answers — New Questions
Is my PPO+LSTM agent reasonable?
Hello everyone, I created PPO+LSTM agent, here is my code
cnet = [
sequenceInputLayer(numObs,"Name","name1")
fullyConnectedLayer(256,"Name","fc1")
tanhLayer("Name","tanh1")
fullyConnectedLayer(128,"Name","fc2")
tanhLayer("Name","tanh2")
fullyConnectedLayer(64,"Name","fc3")
tanhLayer("Name","tanh3")
lstmLayer(16,"Name","LSTM2")
reluLayer
fullyConnectedLayer(1,"Name","CriticOutput")];
criticdlnet = dlnetwork(cnet);
critic= rlValueFunction(criticdlnet,obsInfo);
anet = [
sequenceInputLayer(numObs,"Name","name1")
fullyConnectedLayer(256,"Name","fc1")
tanhLayer("Name","tanh1")
fullyConnectedLayer(128,"Name","fc2")
tanhLayer("Name","tanh2")
lstmLayer(16,"Name","LSTM1")
reluLayer("Name","relu2")
];
meanPath = [
fullyConnectedLayer(32,Name="meanPathIn")
tanhLayer
fullyConnectedLayer(numAct,"Name","mean")
tanhLayer("Name","meanPathOut")];
stdPath = [
fullyConnectedLayer(32,Name="stdPathIn")
reluLayer
fullyConnectedLayer(numAct,"Name","std")
softplusLayer("Name","stdPathOut")];
actordlnet = layerGraph(anet);
actordlnet = addLayers(actordlnet,meanPath);
actordlnet = addLayers(actordlnet,stdPath);
actordlnet = connectLayers(actordlnet,"relu2","meanPathIn/in");
actordlnet = connectLayers(actordlnet,"relu2","stdPathIn/in");
actordlnet = dlnetwork(actordlnet);
actor = rlContinuousGaussianActor(actordlnet,obsInfo,actInfo,…
"ActionMeanOutputNames","meanPathOut", …
"ActionStandardDeviationOutputNames","stdPathOut","ObservationInputNames","name1");
agentOptions=rlPPOAgentOptions("SampleTime",Ts,"DiscountFactor",0.992,"ExperienceHorizon",1024,"MiniBatchSize",64,"ClipFactor",0.2, …
"EntropyLossWeight",0.01,"NumEpoch",1,"AdvantageEstimateMethod","gae","GAEFactor",0.95, …
"NormalizedAdvantageMethod","current");
agent=rlPPOAgent(actor,critic,agentOptions);
Is my deep network structure and hyperparameters reasonable?
I hope I can get your help?Hello everyone, I created PPO+LSTM agent, here is my code
cnet = [
sequenceInputLayer(numObs,"Name","name1")
fullyConnectedLayer(256,"Name","fc1")
tanhLayer("Name","tanh1")
fullyConnectedLayer(128,"Name","fc2")
tanhLayer("Name","tanh2")
fullyConnectedLayer(64,"Name","fc3")
tanhLayer("Name","tanh3")
lstmLayer(16,"Name","LSTM2")
reluLayer
fullyConnectedLayer(1,"Name","CriticOutput")];
criticdlnet = dlnetwork(cnet);
critic= rlValueFunction(criticdlnet,obsInfo);
anet = [
sequenceInputLayer(numObs,"Name","name1")
fullyConnectedLayer(256,"Name","fc1")
tanhLayer("Name","tanh1")
fullyConnectedLayer(128,"Name","fc2")
tanhLayer("Name","tanh2")
lstmLayer(16,"Name","LSTM1")
reluLayer("Name","relu2")
];
meanPath = [
fullyConnectedLayer(32,Name="meanPathIn")
tanhLayer
fullyConnectedLayer(numAct,"Name","mean")
tanhLayer("Name","meanPathOut")];
stdPath = [
fullyConnectedLayer(32,Name="stdPathIn")
reluLayer
fullyConnectedLayer(numAct,"Name","std")
softplusLayer("Name","stdPathOut")];
actordlnet = layerGraph(anet);
actordlnet = addLayers(actordlnet,meanPath);
actordlnet = addLayers(actordlnet,stdPath);
actordlnet = connectLayers(actordlnet,"relu2","meanPathIn/in");
actordlnet = connectLayers(actordlnet,"relu2","stdPathIn/in");
actordlnet = dlnetwork(actordlnet);
actor = rlContinuousGaussianActor(actordlnet,obsInfo,actInfo,…
"ActionMeanOutputNames","meanPathOut", …
"ActionStandardDeviationOutputNames","stdPathOut","ObservationInputNames","name1");
agentOptions=rlPPOAgentOptions("SampleTime",Ts,"DiscountFactor",0.992,"ExperienceHorizon",1024,"MiniBatchSize",64,"ClipFactor",0.2, …
"EntropyLossWeight",0.01,"NumEpoch",1,"AdvantageEstimateMethod","gae","GAEFactor",0.95, …
"NormalizedAdvantageMethod","current");
agent=rlPPOAgent(actor,critic,agentOptions);
Is my deep network structure and hyperparameters reasonable?
I hope I can get your help? Hello everyone, I created PPO+LSTM agent, here is my code
cnet = [
sequenceInputLayer(numObs,"Name","name1")
fullyConnectedLayer(256,"Name","fc1")
tanhLayer("Name","tanh1")
fullyConnectedLayer(128,"Name","fc2")
tanhLayer("Name","tanh2")
fullyConnectedLayer(64,"Name","fc3")
tanhLayer("Name","tanh3")
lstmLayer(16,"Name","LSTM2")
reluLayer
fullyConnectedLayer(1,"Name","CriticOutput")];
criticdlnet = dlnetwork(cnet);
critic= rlValueFunction(criticdlnet,obsInfo);
anet = [
sequenceInputLayer(numObs,"Name","name1")
fullyConnectedLayer(256,"Name","fc1")
tanhLayer("Name","tanh1")
fullyConnectedLayer(128,"Name","fc2")
tanhLayer("Name","tanh2")
lstmLayer(16,"Name","LSTM1")
reluLayer("Name","relu2")
];
meanPath = [
fullyConnectedLayer(32,Name="meanPathIn")
tanhLayer
fullyConnectedLayer(numAct,"Name","mean")
tanhLayer("Name","meanPathOut")];
stdPath = [
fullyConnectedLayer(32,Name="stdPathIn")
reluLayer
fullyConnectedLayer(numAct,"Name","std")
softplusLayer("Name","stdPathOut")];
actordlnet = layerGraph(anet);
actordlnet = addLayers(actordlnet,meanPath);
actordlnet = addLayers(actordlnet,stdPath);
actordlnet = connectLayers(actordlnet,"relu2","meanPathIn/in");
actordlnet = connectLayers(actordlnet,"relu2","stdPathIn/in");
actordlnet = dlnetwork(actordlnet);
actor = rlContinuousGaussianActor(actordlnet,obsInfo,actInfo,…
"ActionMeanOutputNames","meanPathOut", …
"ActionStandardDeviationOutputNames","stdPathOut","ObservationInputNames","name1");
agentOptions=rlPPOAgentOptions("SampleTime",Ts,"DiscountFactor",0.992,"ExperienceHorizon",1024,"MiniBatchSize",64,"ClipFactor",0.2, …
"EntropyLossWeight",0.01,"NumEpoch",1,"AdvantageEstimateMethod","gae","GAEFactor",0.95, …
"NormalizedAdvantageMethod","current");
agent=rlPPOAgent(actor,critic,agentOptions);
Is my deep network structure and hyperparameters reasonable?
I hope I can get your help? drl ppo MATLAB Answers — New Questions
Real time Neural Network issue with “preparets”.
Hi everyone,
I am trying to execute a real time neural network and, as long as I know, I need the preparets() function to “prepare” the variables for the neural network.
net= load(‘NARX_CMON’);
X = tonndata(up,true,false);
T = tonndata(yk’,true,false);
%Execution of the net
[x,xi,ai] = preparets(net,X,{},T);
%where:
%net is the neural network
%X input data
%T target data
yp(k) = cell2mat(neural_function(x, xi, ai));
The problem here is the following.
As I said before, I need to execute preparets () function, and for that, I need the “net” network class variable. To load that variable, I use, as you can see in the following code lines, the function load() .
But as you can see in the error picture that I attach below, Matlab says me that this variable is unsupported and thus, preparets() doesn’t work.
What am I doing wrong? What is going on? How can I solve this problem?
Thanks a lot,Hi everyone,
I am trying to execute a real time neural network and, as long as I know, I need the preparets() function to “prepare” the variables for the neural network.
net= load(‘NARX_CMON’);
X = tonndata(up,true,false);
T = tonndata(yk’,true,false);
%Execution of the net
[x,xi,ai] = preparets(net,X,{},T);
%where:
%net is the neural network
%X input data
%T target data
yp(k) = cell2mat(neural_function(x, xi, ai));
The problem here is the following.
As I said before, I need to execute preparets () function, and for that, I need the “net” network class variable. To load that variable, I use, as you can see in the following code lines, the function load() .
But as you can see in the error picture that I attach below, Matlab says me that this variable is unsupported and thus, preparets() doesn’t work.
What am I doing wrong? What is going on? How can I solve this problem?
Thanks a lot, Hi everyone,
I am trying to execute a real time neural network and, as long as I know, I need the preparets() function to “prepare” the variables for the neural network.
net= load(‘NARX_CMON’);
X = tonndata(up,true,false);
T = tonndata(yk’,true,false);
%Execution of the net
[x,xi,ai] = preparets(net,X,{},T);
%where:
%net is the neural network
%X input data
%T target data
yp(k) = cell2mat(neural_function(x, xi, ai));
The problem here is the following.
As I said before, I need to execute preparets () function, and for that, I need the “net” network class variable. To load that variable, I use, as you can see in the following code lines, the function load() .
But as you can see in the error picture that I attach below, Matlab says me that this variable is unsupported and thus, preparets() doesn’t work.
What am I doing wrong? What is going on? How can I solve this problem?
Thanks a lot, net, neuralnetwork, matlab, preparets, neural, nnarx, mpc MATLAB Answers — New Questions
Why does readmatrix return empty matrix if text file contains quotes?
Function readmatrix fails to read data from text file when it contains quotes. Why? How to make readmatrix working?
Extended example:
let first text.txt be:
header1
header2
1
2
3
footer
then
readmatrix(‘text.txt’, ‘Range’, [3 1 5 1])
gives usual vector [1;2;3].
But whe text.txt is like this:
"
header2
1
2
3
"
readmatrix returns empty vector while dlmread works fine.
How to force readmatrix ignore quotes?Function readmatrix fails to read data from text file when it contains quotes. Why? How to make readmatrix working?
Extended example:
let first text.txt be:
header1
header2
1
2
3
footer
then
readmatrix(‘text.txt’, ‘Range’, [3 1 5 1])
gives usual vector [1;2;3].
But whe text.txt is like this:
"
header2
1
2
3
"
readmatrix returns empty vector while dlmread works fine.
How to force readmatrix ignore quotes? Function readmatrix fails to read data from text file when it contains quotes. Why? How to make readmatrix working?
Extended example:
let first text.txt be:
header1
header2
1
2
3
footer
then
readmatrix(‘text.txt’, ‘Range’, [3 1 5 1])
gives usual vector [1;2;3].
But whe text.txt is like this:
"
header2
1
2
3
"
readmatrix returns empty vector while dlmread works fine.
How to force readmatrix ignore quotes? readmatrix MATLAB Answers — New Questions
Error in Trying to Build a Custom Variable Resistor In MATLAB Simulink
Hey everyone, I’m in the process of trying to build a custom variable resistor using simulink. I created a code for a simscape component in MATLAB but it’s giving me an error saying "Invalid use of a MATLAB function in This Context" when I try inserting the simscape component in simulink. Shown below is my code:
The error happens on line 26 (the equation with variable R). Are there any tips or suggestions to fix my error?
component VoltageDependentResistor
% Voltage-dependent resistor
% This resistor’s value changes based on the voltage across it.
nodes
p = foundation.electrical.electrical; % positive terminal
n = foundation.electrical.electrical; % negative terminal
end
parameters
alpha = {1, ‘1/V’}; % Coefficient for voltage dependence
epsilon = 0.6;
sigma = 5.67E-8
end
variables
v = { 0, ‘V’ }; % Voltage across the resistor
i = { 0, ‘A’ }; % Current through the resistor
end
branches
i : p.i -> n.i; % Current from p to n
end
equations
v == p.v – n.v; % Voltage across the resistor
R == 1/((epsilon*sigma*(p.v^2+n.v^2)*(p.v+n.v))/(2-epsilon)); % Resistance as a function of radiative heat transfer
v == R * i; % Ohm’s Law
end
end
Thanks for your time and help,
Elias LynchHey everyone, I’m in the process of trying to build a custom variable resistor using simulink. I created a code for a simscape component in MATLAB but it’s giving me an error saying "Invalid use of a MATLAB function in This Context" when I try inserting the simscape component in simulink. Shown below is my code:
The error happens on line 26 (the equation with variable R). Are there any tips or suggestions to fix my error?
component VoltageDependentResistor
% Voltage-dependent resistor
% This resistor’s value changes based on the voltage across it.
nodes
p = foundation.electrical.electrical; % positive terminal
n = foundation.electrical.electrical; % negative terminal
end
parameters
alpha = {1, ‘1/V’}; % Coefficient for voltage dependence
epsilon = 0.6;
sigma = 5.67E-8
end
variables
v = { 0, ‘V’ }; % Voltage across the resistor
i = { 0, ‘A’ }; % Current through the resistor
end
branches
i : p.i -> n.i; % Current from p to n
end
equations
v == p.v – n.v; % Voltage across the resistor
R == 1/((epsilon*sigma*(p.v^2+n.v^2)*(p.v+n.v))/(2-epsilon)); % Resistance as a function of radiative heat transfer
v == R * i; % Ohm’s Law
end
end
Thanks for your time and help,
Elias Lynch Hey everyone, I’m in the process of trying to build a custom variable resistor using simulink. I created a code for a simscape component in MATLAB but it’s giving me an error saying "Invalid use of a MATLAB function in This Context" when I try inserting the simscape component in simulink. Shown below is my code:
The error happens on line 26 (the equation with variable R). Are there any tips or suggestions to fix my error?
component VoltageDependentResistor
% Voltage-dependent resistor
% This resistor’s value changes based on the voltage across it.
nodes
p = foundation.electrical.electrical; % positive terminal
n = foundation.electrical.electrical; % negative terminal
end
parameters
alpha = {1, ‘1/V’}; % Coefficient for voltage dependence
epsilon = 0.6;
sigma = 5.67E-8
end
variables
v = { 0, ‘V’ }; % Voltage across the resistor
i = { 0, ‘A’ }; % Current through the resistor
end
branches
i : p.i -> n.i; % Current from p to n
end
equations
v == p.v – n.v; % Voltage across the resistor
R == 1/((epsilon*sigma*(p.v^2+n.v^2)*(p.v+n.v))/(2-epsilon)); % Resistance as a function of radiative heat transfer
v == R * i; % Ohm’s Law
end
end
Thanks for your time and help,
Elias Lynch simulink, matlab, simscape MATLAB Answers — New Questions
plotting the bifurcation diagram of Duffing oscillations in the (x1, omega) plane?
function dx = duffing(t,x);
global alpha delta beta F omega
dx=[x(2);
-delta*x(2)-alpha*x(1)-beta*x(1)^3+F*cos(x(3));
omega];
% Define parameters
global alpha delta beta F omega
alpha=0.1;
delta=0.2;
beta=0.3;
F=0.5;
omega_range = 0.8:0.1:1; % Range of omega values
x0 = [0.1, 0.1, 0.1]; % Initial conditions
% Set simulation step
dt = 0.001;
% Presave a matrix of NaNs to save points of intersection
M = NaN * zeros(1000, length(omega_range));
pos = 0; % Indexing dummy variable
% Set ode options
options = odeset(‘RelTol’, 1e-5, ‘AbsTol’, 1e-5);
for omega = omega_range
omega % Print omega, just to track progress
pos = pos + 1; % Increase index
% Simulate the system
[t, x] = ode45(@duffing, 0:dt:1000, x0, options);
% Discard transient
index = t > 400;
X = x(index, :); % Save states without transient in a new vector
l = length(X);
% Save the final x(1) value
M(1:l, pos) = X(:, 1);
end
% Plot the bifurcation diagram
figure;
plot(omega_range, M, ‘.b’, ‘MarkerSize’, 2);
xlabel(‘omega’);
ylabel(‘x_1’);
title(‘Bifurcation Diagram of Duffing System (x_1 vs. omega)’);
grid on;
Can anyone help troubleshoot my code for plotting the bifurcation diagram of Duffing oscillations in the (x1, omega) plane? It’s not functioning as expected. Thanks!function dx = duffing(t,x);
global alpha delta beta F omega
dx=[x(2);
-delta*x(2)-alpha*x(1)-beta*x(1)^3+F*cos(x(3));
omega];
% Define parameters
global alpha delta beta F omega
alpha=0.1;
delta=0.2;
beta=0.3;
F=0.5;
omega_range = 0.8:0.1:1; % Range of omega values
x0 = [0.1, 0.1, 0.1]; % Initial conditions
% Set simulation step
dt = 0.001;
% Presave a matrix of NaNs to save points of intersection
M = NaN * zeros(1000, length(omega_range));
pos = 0; % Indexing dummy variable
% Set ode options
options = odeset(‘RelTol’, 1e-5, ‘AbsTol’, 1e-5);
for omega = omega_range
omega % Print omega, just to track progress
pos = pos + 1; % Increase index
% Simulate the system
[t, x] = ode45(@duffing, 0:dt:1000, x0, options);
% Discard transient
index = t > 400;
X = x(index, :); % Save states without transient in a new vector
l = length(X);
% Save the final x(1) value
M(1:l, pos) = X(:, 1);
end
% Plot the bifurcation diagram
figure;
plot(omega_range, M, ‘.b’, ‘MarkerSize’, 2);
xlabel(‘omega’);
ylabel(‘x_1’);
title(‘Bifurcation Diagram of Duffing System (x_1 vs. omega)’);
grid on;
Can anyone help troubleshoot my code for plotting the bifurcation diagram of Duffing oscillations in the (x1, omega) plane? It’s not functioning as expected. Thanks! function dx = duffing(t,x);
global alpha delta beta F omega
dx=[x(2);
-delta*x(2)-alpha*x(1)-beta*x(1)^3+F*cos(x(3));
omega];
% Define parameters
global alpha delta beta F omega
alpha=0.1;
delta=0.2;
beta=0.3;
F=0.5;
omega_range = 0.8:0.1:1; % Range of omega values
x0 = [0.1, 0.1, 0.1]; % Initial conditions
% Set simulation step
dt = 0.001;
% Presave a matrix of NaNs to save points of intersection
M = NaN * zeros(1000, length(omega_range));
pos = 0; % Indexing dummy variable
% Set ode options
options = odeset(‘RelTol’, 1e-5, ‘AbsTol’, 1e-5);
for omega = omega_range
omega % Print omega, just to track progress
pos = pos + 1; % Increase index
% Simulate the system
[t, x] = ode45(@duffing, 0:dt:1000, x0, options);
% Discard transient
index = t > 400;
X = x(index, :); % Save states without transient in a new vector
l = length(X);
% Save the final x(1) value
M(1:l, pos) = X(:, 1);
end
% Plot the bifurcation diagram
figure;
plot(omega_range, M, ‘.b’, ‘MarkerSize’, 2);
xlabel(‘omega’);
ylabel(‘x_1’);
title(‘Bifurcation Diagram of Duffing System (x_1 vs. omega)’);
grid on;
Can anyone help troubleshoot my code for plotting the bifurcation diagram of Duffing oscillations in the (x1, omega) plane? It’s not functioning as expected. Thanks! differential equations, nonlinear dynamics MATLAB Answers — New Questions
Find centroid of binary image
I want to know find centroid area in picture output show X,Y axis
Sample
centroid x 60 centroid y 40I want to know find centroid area in picture output show X,Y axis
Sample
centroid x 60 centroid y 40 I want to know find centroid area in picture output show X,Y axis
Sample
centroid x 60 centroid y 40 image processing, image analysis, image segmentation, regionprops MATLAB Answers — New Questions
Summing results at the end of a loop
I have a piece of code I want to execute 3 times, changing a set of two variables to different values each time, and then add the results of the three at the end to give the final y-value and plot it.
This is the full code right now.
%Constants
L = 79.14.*(10.^-3);
N = 58;
I = 0.7;
mu = 4.*pi*(10.^-7);
% Range and steps of x-values to plot:
x = linspace(0, 90.*(10.^-3), 100);
%Variables;
w = 7.5.*(10.^-3);
u = (8.883).*(10.^-3);
%Calculation;
a = (mu.*I.*N)/(2.*L.*(u-w));
b = sqrt((u.^2)+(x.^2));
c = sqrt((u.^2)+((x-L).^2));
d = sqrt((w.^2)+(x.^2));
e = sqrt((w.^2)+((x-L).^2));
f = (b+u)./(d+w);
g = x.*log(f);
h = (c+u)./(e+w);
i = (x-L).*log(h);
j = (mu.*I.*N)/(2.*L.*(u-w));
y = j.*(g-i);
%Plotting;
figure;
plot(x, y);
xlabel(‘Position Along the Center Axis (m)’);
ylabel(‘Magnetic Field (T)’);
title(‘Magnetic Field Along the Central Axis of the Solenoid’);
What I want to do is repeat the code for 3 sets of values of u and w, add the three results at the end to give the final value of y, and then plot it against x.
I thought about using for loops but I can’t figure out how to sum the values at the end.I have a piece of code I want to execute 3 times, changing a set of two variables to different values each time, and then add the results of the three at the end to give the final y-value and plot it.
This is the full code right now.
%Constants
L = 79.14.*(10.^-3);
N = 58;
I = 0.7;
mu = 4.*pi*(10.^-7);
% Range and steps of x-values to plot:
x = linspace(0, 90.*(10.^-3), 100);
%Variables;
w = 7.5.*(10.^-3);
u = (8.883).*(10.^-3);
%Calculation;
a = (mu.*I.*N)/(2.*L.*(u-w));
b = sqrt((u.^2)+(x.^2));
c = sqrt((u.^2)+((x-L).^2));
d = sqrt((w.^2)+(x.^2));
e = sqrt((w.^2)+((x-L).^2));
f = (b+u)./(d+w);
g = x.*log(f);
h = (c+u)./(e+w);
i = (x-L).*log(h);
j = (mu.*I.*N)/(2.*L.*(u-w));
y = j.*(g-i);
%Plotting;
figure;
plot(x, y);
xlabel(‘Position Along the Center Axis (m)’);
ylabel(‘Magnetic Field (T)’);
title(‘Magnetic Field Along the Central Axis of the Solenoid’);
What I want to do is repeat the code for 3 sets of values of u and w, add the three results at the end to give the final value of y, and then plot it against x.
I thought about using for loops but I can’t figure out how to sum the values at the end. I have a piece of code I want to execute 3 times, changing a set of two variables to different values each time, and then add the results of the three at the end to give the final y-value and plot it.
This is the full code right now.
%Constants
L = 79.14.*(10.^-3);
N = 58;
I = 0.7;
mu = 4.*pi*(10.^-7);
% Range and steps of x-values to plot:
x = linspace(0, 90.*(10.^-3), 100);
%Variables;
w = 7.5.*(10.^-3);
u = (8.883).*(10.^-3);
%Calculation;
a = (mu.*I.*N)/(2.*L.*(u-w));
b = sqrt((u.^2)+(x.^2));
c = sqrt((u.^2)+((x-L).^2));
d = sqrt((w.^2)+(x.^2));
e = sqrt((w.^2)+((x-L).^2));
f = (b+u)./(d+w);
g = x.*log(f);
h = (c+u)./(e+w);
i = (x-L).*log(h);
j = (mu.*I.*N)/(2.*L.*(u-w));
y = j.*(g-i);
%Plotting;
figure;
plot(x, y);
xlabel(‘Position Along the Center Axis (m)’);
ylabel(‘Magnetic Field (T)’);
title(‘Magnetic Field Along the Central Axis of the Solenoid’);
What I want to do is repeat the code for 3 sets of values of u and w, add the three results at the end to give the final value of y, and then plot it against x.
I thought about using for loops but I can’t figure out how to sum the values at the end. solenoid, physics, matlab, for loop, iteration MATLAB Answers — New Questions
fitcdiscr: Why does “ClassNames” have to be provided in alphabetical order or accuracy is terrible?
I couldn’t work out why I was getting terrible results (as if completely random) with fitcdiscr() and I’ve found out that it is because I wasn’t specifying the ClassNames argument in alphabetical order! Is there something I am missing or is this a major pitfall that should be highlighted in the documentation, or a even bug in MATLAB?
Example code that gives terrible accuracy:
load fisheriris
Mdl = fitcdiscr(meas, species, "ClassNames", flip(unique(species, "stable")), "KFold", 10);
validationAccuracy = 1 – kfoldLoss(Mdl)
By simply removing the function flip(), the above code gives the expected accuracy of 0.98, otherwise it gives 0.32.I couldn’t work out why I was getting terrible results (as if completely random) with fitcdiscr() and I’ve found out that it is because I wasn’t specifying the ClassNames argument in alphabetical order! Is there something I am missing or is this a major pitfall that should be highlighted in the documentation, or a even bug in MATLAB?
Example code that gives terrible accuracy:
load fisheriris
Mdl = fitcdiscr(meas, species, "ClassNames", flip(unique(species, "stable")), "KFold", 10);
validationAccuracy = 1 – kfoldLoss(Mdl)
By simply removing the function flip(), the above code gives the expected accuracy of 0.98, otherwise it gives 0.32. I couldn’t work out why I was getting terrible results (as if completely random) with fitcdiscr() and I’ve found out that it is because I wasn’t specifying the ClassNames argument in alphabetical order! Is there something I am missing or is this a major pitfall that should be highlighted in the documentation, or a even bug in MATLAB?
Example code that gives terrible accuracy:
load fisheriris
Mdl = fitcdiscr(meas, species, "ClassNames", flip(unique(species, "stable")), "KFold", 10);
validationAccuracy = 1 – kfoldLoss(Mdl)
By simply removing the function flip(), the above code gives the expected accuracy of 0.98, otherwise it gives 0.32. classnames, fitcdiscr, order, accuracy, unique MATLAB Answers — New Questions
Bayesian Non-linear SSM with time varying coefficients
I have question about using ‘bnlssm’ object when the problem involves time-varying coefs. Per this documentation https://www.mathworks.com/help/econ/bnlssm.html , bnlssm allows the transition matrix to be a nonlinear function of the past values of the state vector, while also being time-varying –, i.e. A_t(x_{t-1}). However the documentation contains an example of the parameter mapping function for a fixed coef. case only:
function [A,B,C,D,Mean0,Cov0,StateType] = paramMap(theta)
A = @(x)blkdiag([theta(1) theta(2); 0 1],[theta(3) theta(4); 0 1])*x;
B = [theta(5) 0; 0 0; 0 theta(6); 0 0];
C = @(x)log(exp(x(1)-theta(2)/(1-theta(1))) + …
exp(x(3)-theta(4)/(1-theta(3))));
D = theta(7);
Mean0 = [theta(2)/(1-theta(1)); 1; theta(4)/(1-theta(3)); 1];
Cov0 = diag([theta(5)^2/(1-theta(1)^2) 0 theta(6)^2/(1-theta(3)^2) 0]);
StateType = [0; 1; 0; 1]; % Stationary state and constant 1 processes
end
Here, A and C are recognised as function handles and show the dimension of 1-by-1 when [A,B,C,D]=paramMap(theta) is run with a given theta. This does not prevent filter/smooth from returning the states estimates, as the function handles become m-by-m matices when the function handles are evaluated.
However, whenever I try replaceing A with a cell array, e.g.:
A=cell(T,1);
for t=1:T
A{t} = @(x)blkdiag([theta(1) theta(2); 0 1],[theta(3) theta(4); 0 1])*x;
end
bnlssm.fixcoeff throws an error that the dimensions of Mean0 are expected to be 1, per validation step in lines 700-707 of ‘bnlssm.m’:
if iscellA; At = A{1};
else; At = A; end
numStates0 = size(At,2);
if
isscalar(mean0)
mean0 = mean0(ones(numStates0,1),:);
elseif
~isempty(mean0)
mean0 = mean0(:);
validateattributes(mean0,{‘numeric’},{‘finite’,’real’,’numel’,numStates0},callerName,’Mean0′);
end
Q: Is it just overlook by MathWorks programmers and creating local verions of bnlssm/filter/smooth with that validation step removed is a viable pathway forward, or is it intentional because ‘bnlssm’ was in fact designed to handle only the time-invariant coef. case correctly?I have question about using ‘bnlssm’ object when the problem involves time-varying coefs. Per this documentation https://www.mathworks.com/help/econ/bnlssm.html , bnlssm allows the transition matrix to be a nonlinear function of the past values of the state vector, while also being time-varying –, i.e. A_t(x_{t-1}). However the documentation contains an example of the parameter mapping function for a fixed coef. case only:
function [A,B,C,D,Mean0,Cov0,StateType] = paramMap(theta)
A = @(x)blkdiag([theta(1) theta(2); 0 1],[theta(3) theta(4); 0 1])*x;
B = [theta(5) 0; 0 0; 0 theta(6); 0 0];
C = @(x)log(exp(x(1)-theta(2)/(1-theta(1))) + …
exp(x(3)-theta(4)/(1-theta(3))));
D = theta(7);
Mean0 = [theta(2)/(1-theta(1)); 1; theta(4)/(1-theta(3)); 1];
Cov0 = diag([theta(5)^2/(1-theta(1)^2) 0 theta(6)^2/(1-theta(3)^2) 0]);
StateType = [0; 1; 0; 1]; % Stationary state and constant 1 processes
end
Here, A and C are recognised as function handles and show the dimension of 1-by-1 when [A,B,C,D]=paramMap(theta) is run with a given theta. This does not prevent filter/smooth from returning the states estimates, as the function handles become m-by-m matices when the function handles are evaluated.
However, whenever I try replaceing A with a cell array, e.g.:
A=cell(T,1);
for t=1:T
A{t} = @(x)blkdiag([theta(1) theta(2); 0 1],[theta(3) theta(4); 0 1])*x;
end
bnlssm.fixcoeff throws an error that the dimensions of Mean0 are expected to be 1, per validation step in lines 700-707 of ‘bnlssm.m’:
if iscellA; At = A{1};
else; At = A; end
numStates0 = size(At,2);
if
isscalar(mean0)
mean0 = mean0(ones(numStates0,1),:);
elseif
~isempty(mean0)
mean0 = mean0(:);
validateattributes(mean0,{‘numeric’},{‘finite’,’real’,’numel’,numStates0},callerName,’Mean0′);
end
Q: Is it just overlook by MathWorks programmers and creating local verions of bnlssm/filter/smooth with that validation step removed is a viable pathway forward, or is it intentional because ‘bnlssm’ was in fact designed to handle only the time-invariant coef. case correctly? I have question about using ‘bnlssm’ object when the problem involves time-varying coefs. Per this documentation https://www.mathworks.com/help/econ/bnlssm.html , bnlssm allows the transition matrix to be a nonlinear function of the past values of the state vector, while also being time-varying –, i.e. A_t(x_{t-1}). However the documentation contains an example of the parameter mapping function for a fixed coef. case only:
function [A,B,C,D,Mean0,Cov0,StateType] = paramMap(theta)
A = @(x)blkdiag([theta(1) theta(2); 0 1],[theta(3) theta(4); 0 1])*x;
B = [theta(5) 0; 0 0; 0 theta(6); 0 0];
C = @(x)log(exp(x(1)-theta(2)/(1-theta(1))) + …
exp(x(3)-theta(4)/(1-theta(3))));
D = theta(7);
Mean0 = [theta(2)/(1-theta(1)); 1; theta(4)/(1-theta(3)); 1];
Cov0 = diag([theta(5)^2/(1-theta(1)^2) 0 theta(6)^2/(1-theta(3)^2) 0]);
StateType = [0; 1; 0; 1]; % Stationary state and constant 1 processes
end
Here, A and C are recognised as function handles and show the dimension of 1-by-1 when [A,B,C,D]=paramMap(theta) is run with a given theta. This does not prevent filter/smooth from returning the states estimates, as the function handles become m-by-m matices when the function handles are evaluated.
However, whenever I try replaceing A with a cell array, e.g.:
A=cell(T,1);
for t=1:T
A{t} = @(x)blkdiag([theta(1) theta(2); 0 1],[theta(3) theta(4); 0 1])*x;
end
bnlssm.fixcoeff throws an error that the dimensions of Mean0 are expected to be 1, per validation step in lines 700-707 of ‘bnlssm.m’:
if iscellA; At = A{1};
else; At = A; end
numStates0 = size(At,2);
if
isscalar(mean0)
mean0 = mean0(ones(numStates0,1),:);
elseif
~isempty(mean0)
mean0 = mean0(:);
validateattributes(mean0,{‘numeric’},{‘finite’,’real’,’numel’,numStates0},callerName,’Mean0′);
end
Q: Is it just overlook by MathWorks programmers and creating local verions of bnlssm/filter/smooth with that validation step removed is a viable pathway forward, or is it intentional because ‘bnlssm’ was in fact designed to handle only the time-invariant coef. case correctly? transferred MATLAB Answers — New Questions
how to fix in imresize error ??
i try to make the function for resizing image, this the code:
>> a=imread(‘1.png’);
>> b=20000;
>> l=sum(sum(a));
>> [m,n]=size(a);
>> m1=fix(m*sqrt(b/l));
>> n1=fix(m*sqrt(b/l));
>> c=imresize(a,[m1 n1]);
but, i get the error :
Error using imresize>scaleOrSize (line 382)
Invalid scale or size input argument.
Error in imresize>parsePreMethodArgs (line 354)
[scale, output_size] = scaleOrSize(next, next_arg);
Error in imresize>parseInputs (line 248)
[params.A, params.map, params.scale, params.output_size] = …
Error in imresize (line 141)
params = parseInputs(varargin{:});
i don’t understand what cause of the error and i try to fix the error but i can’t fix it. Please Help, Thank youi try to make the function for resizing image, this the code:
>> a=imread(‘1.png’);
>> b=20000;
>> l=sum(sum(a));
>> [m,n]=size(a);
>> m1=fix(m*sqrt(b/l));
>> n1=fix(m*sqrt(b/l));
>> c=imresize(a,[m1 n1]);
but, i get the error :
Error using imresize>scaleOrSize (line 382)
Invalid scale or size input argument.
Error in imresize>parsePreMethodArgs (line 354)
[scale, output_size] = scaleOrSize(next, next_arg);
Error in imresize>parseInputs (line 248)
[params.A, params.map, params.scale, params.output_size] = …
Error in imresize (line 141)
params = parseInputs(varargin{:});
i don’t understand what cause of the error and i try to fix the error but i can’t fix it. Please Help, Thank you i try to make the function for resizing image, this the code:
>> a=imread(‘1.png’);
>> b=20000;
>> l=sum(sum(a));
>> [m,n]=size(a);
>> m1=fix(m*sqrt(b/l));
>> n1=fix(m*sqrt(b/l));
>> c=imresize(a,[m1 n1]);
but, i get the error :
Error using imresize>scaleOrSize (line 382)
Invalid scale or size input argument.
Error in imresize>parsePreMethodArgs (line 354)
[scale, output_size] = scaleOrSize(next, next_arg);
Error in imresize>parseInputs (line 248)
[params.A, params.map, params.scale, params.output_size] = …
Error in imresize (line 141)
params = parseInputs(varargin{:});
i don’t understand what cause of the error and i try to fix the error but i can’t fix it. Please Help, Thank you image processing MATLAB Answers — New Questions
Hardware advice on a project
Hello everyone!
I’m developing a project that uses two rotational encoders and one DC motor. An H-bridge drives the DC motor. I need some advice on which acquisition platform to use. I need one compactable with Simulink (something like Arduino support package, easy to develop the control algorithm) that enables small sampling periods in real-time.
Could anyone help me to choose such hardware?
I really appreciate any help you can provide.Hello everyone!
I’m developing a project that uses two rotational encoders and one DC motor. An H-bridge drives the DC motor. I need some advice on which acquisition platform to use. I need one compactable with Simulink (something like Arduino support package, easy to develop the control algorithm) that enables small sampling periods in real-time.
Could anyone help me to choose such hardware?
I really appreciate any help you can provide. Hello everyone!
I’m developing a project that uses two rotational encoders and one DC motor. An H-bridge drives the DC motor. I need some advice on which acquisition platform to use. I need one compactable with Simulink (something like Arduino support package, easy to develop the control algorithm) that enables small sampling periods in real-time.
Could anyone help me to choose such hardware?
I really appreciate any help you can provide. #hardware #simulink #realtime #supportpackage MATLAB Answers — New Questions