I write documentation for MATLAB mathematical toolboxes, primarily optimization. I have also written documentation for PDE, statistics, symbolic math, and econometrics.

My pre-MathWorks job was with Bell Labs, primarily in mathematical models of data traffic, with a strong interest in parallel computation and in rare events (large deviations).

I do not provide private consulting. If you have a question, please ask on a public forum such as MATLAB Answers.

Answered

someone please help me. I have made a 2 variable code for optimization, I want code for "n" no. of variables.

It is probably a bit easier to write code for NUMERIC minimizaton of an arbitrary-sized expression than a hybrid SYMBOLIC minimi...

someone please help me. I have made a 2 variable code for optimization, I want code for "n" no. of variables.

It is probably a bit easier to write code for NUMERIC minimizaton of an arbitrary-sized expression than a hybrid SYMBOLIC minimi...

1 jour ago | 0

Answered

Writing complicated equations for problem-based optimization

I think that you need to keep clear in your mind (and code) what is an optimization variable and what is an optimization express...

Writing complicated equations for problem-based optimization

I think that you need to keep clear in your mind (and code) what is an optimization variable and what is an optimization express...

2 jours ago | 0

| accepted

Answered

3D Scatterplot and Pareto Front visualization

Did you try the scatteredinterpolant code from the answer? Modify it for your data: F = scatteredInterpolant(numbers(:,1),numbe...

3D Scatterplot and Pareto Front visualization

Did you try the scatteredinterpolant code from the answer? Modify it for your data: F = scatteredInterpolant(numbers(:,1),numbe...

3 jours ago | 0

Answered

Nonlinear Constraints in Multijob optimization

Your constraints can be formulated as nonlinear constraints. Make sure that you read and follow the note on that page. Alan Wei...

Nonlinear Constraints in Multijob optimization

Your constraints can be formulated as nonlinear constraints. Make sure that you read and follow the note on that page. Alan Wei...

7 jours ago | 0

Answered

Error with lsqnonlin : Error in lsqncommon (line 14) if any(~isfinite(initVals.F))

Apparently, some initial values are not finite. You gave a value of x0 so that your objective function evaluated at x0 gives som...

Error with lsqnonlin : Error in lsqncommon (line 14) if any(~isfinite(initVals.F))

Apparently, some initial values are not finite. You gave a value of x0 so that your objective function evaluated at x0 gives som...

9 jours ago | 0

| accepted

Answered

Initial population in genetic algorithm

To get reproducible results, set the random number generator before you call ga. For example, rng default % or rng(1) or whate...

Initial population in genetic algorithm

To get reproducible results, set the random number generator before you call ga. For example, rng default % or rng(1) or whate...

10 jours ago | 0

| accepted

Answered

how to solve this error: objective function is returning undefined values at initial point. fsolve cannot continue.

I suggest that you learn to use the debugger. I put a break point in the solve_nch function, and found that at the first functio...

how to solve this error: objective function is returning undefined values at initial point. fsolve cannot continue.

I suggest that you learn to use the debugger. I put a break point in the solve_nch function, and found that at the first functio...

10 jours ago | 1

Answered

Matrix array in ga

I think that you are making a mistake by choosing the bitstring data type. Doing so loses the possibility of having constraints ...

Matrix array in ga

I think that you are making a mistake by choosing the bitstring data type. Doing so loses the possibility of having constraints ...

15 jours ago | 0

| accepted

Answered

Optimisation Tool GA Custom Plot

Perhaps if you put the names of the functions in curly braces {} you would have success: {@gaplotbestf,@gaplotchange} But come...

Optimisation Tool GA Custom Plot

Perhaps if you put the names of the functions in curly braces {} you would have success: {@gaplotbestf,@gaplotchange} But come...

15 jours ago | 0

| accepted

Answered

How can we insert upper and lower bounds in particle swarm optimization for more than one variable parameters that need to be optimized.

I am not sure that I understand you. You might be asking about bounding some control variables. Generally, all of your control v...

How can we insert upper and lower bounds in particle swarm optimization for more than one variable parameters that need to be optimized.

I am not sure that I understand you. You might be asking about bounding some control variables. Generally, all of your control v...

16 jours ago | 1

Answered

How can we get the surrogate function after the optimisation?

The surrogate model is not available, sorry. Indeed, its quality as a surrogate is unknown, and may be poor. Can you please say...

How can we get the surrogate function after the optimisation?

The surrogate model is not available, sorry. Indeed, its quality as a surrogate is unknown, and may be poor. Can you please say...

21 jours ago | 0

Answered

Nonlinear fitting: how do I split the linear and the nonlinear problems?

You should follow the example more closely. In the example the lambda variables only are declared to be optimization variables; ...

Nonlinear fitting: how do I split the linear and the nonlinear problems?

You should follow the example more closely. In the example the lambda variables only are declared to be optimization variables; ...

21 jours ago | 1

| accepted

Answered

why does this code give an error ( Undefined function or variable 'ObjectiveFunction'. Error in psi_kozhin (line 44) Swarm.Particles(k).O = fobj(currentX); )?

I do not see where you defined ObjectiveFunction. Is it on your MATLAB path? Alan Weiss MATLAB mathematical toolbox documentat...

why does this code give an error ( Undefined function or variable 'ObjectiveFunction'. Error in psi_kozhin (line 44) Swarm.Particles(k).O = fobj(currentX); )?

I do not see where you defined ObjectiveFunction. Is it on your MATLAB path? Alan Weiss MATLAB mathematical toolbox documentat...

28 jours ago | 0

Answered

Find minimum of multi-variable function on fixed interval with additional parameter inputs

You need to pass the parameters correctly. See Passing Extra Parameters. In your case, you need to either nest the function in a...

Find minimum of multi-variable function on fixed interval with additional parameter inputs

You need to pass the parameters correctly. See Passing Extra Parameters. In your case, you need to either nest the function in a...

28 jours ago | 0

| accepted

Answered

Genetic Algorithm: Plot ONLY best fitness value, not mean.

You could try using this slightly modified plot function, gaplotbestf2: (I modified it slightly from gaplotbestf) function stat...

Genetic Algorithm: Plot ONLY best fitness value, not mean.

You could try using this slightly modified plot function, gaplotbestf2: (I modified it slightly from gaplotbestf) function stat...

28 jours ago | 0

Answered

Fitting nonlinear noisy data

Most likely the issue is that there are multiple local minima, as in this example: Nonlinear Data-Fitting Using Several Problem-...

Fitting nonlinear noisy data

Most likely the issue is that there are multiple local minima, as in this example: Nonlinear Data-Fitting Using Several Problem-...

28 jours ago | 1

Answered

Explanation of hyperparameter tuning procedure for regression tree ensembles

You can find all the information later on in that same reference page: Hyperparameter Optimization Options Alan Weiss MATLAB ...

Explanation of hyperparameter tuning procedure for regression tree ensembles

You can find all the information later on in that same reference page: Hyperparameter Optimization Options Alan Weiss MATLAB ...

environ un mois ago | 0

Answered

Fitting procedure using MultiStart - doesn't recognize objective function

I think that you need to have just one input variable, typically called x, and have each of your other named variables be a comp...

Fitting procedure using MultiStart - doesn't recognize objective function

I think that you need to have just one input variable, typically called x, and have each of your other named variables be a comp...

environ un mois ago | 0

Answered

Optimization toolbox with ode45

I guess that I would make a minimization problem for Y(end,2)^2. I don't know what your parameters are, but you can use fminbnd ...

Optimization toolbox with ode45

I guess that I would make a minimization problem for Y(end,2)^2. I don't know what your parameters are, but you can use fminbnd ...

environ un mois ago | 0

Answered

Surrogate Optimization Error: 'Dimensions of arrays being concatenated are not consistent'

You say "Do not mind the AdiabaticReactor function" but I cannot tell what sizes the inputs and outputs are for the objconstrAdi...

Surrogate Optimization Error: 'Dimensions of arrays being concatenated are not consistent'

You say "Do not mind the AdiabaticReactor function" but I cannot tell what sizes the inputs and outputs are for the objconstrAdi...

environ un mois ago | 0

Answered

Error: Complex values are not supported

I don't know where your objective function creates complex values, though I suspect that one or more of your logarithm calls has...

Error: Complex values are not supported

I don't know where your objective function creates complex values, though I suspect that one or more of your logarithm calls has...

environ un mois ago | 0

| accepted

Answered

How to find optimum complex weighting coefficients for antenna beam pattern optimisation? A purely optimisation problem

fmincon requires real values only. Convert your complex-valued problem to twice as many real variables. For an example, see Fit ...

How to find optimum complex weighting coefficients for antenna beam pattern optimisation? A purely optimisation problem

fmincon requires real values only. Convert your complex-valued problem to twice as many real variables. For an example, see Fit ...

environ un mois ago | 0

Answered

system of nonlinear equations

Do you have Optimization Toolbox™ installed? Check by running the ver command. You might need to install or reinstall. Alan Wei...

system of nonlinear equations

Do you have Optimization Toolbox™ installed? Check by running the ver command. You might need to install or reinstall. Alan Wei...

environ un mois ago | 0

| accepted

Answered

Nonlinear fitting depends too heavily on starting points

You might find some relevant information here or here. The point is that fitting problems typically have multiple local optima, ...

Nonlinear fitting depends too heavily on starting points

You might find some relevant information here or here. The point is that fitting problems typically have multiple local optima, ...

environ un mois ago | 0

| accepted

Answered

Genetic Algorithm to Find desired parameters.

You might want to spend a little time reading up on the genetic algorithm: Genetic Algorithm Terminology How the Genetic Algor...

Genetic Algorithm to Find desired parameters.

You might want to spend a little time reading up on the genetic algorithm: Genetic Algorithm Terminology How the Genetic Algor...

environ un mois ago | 0

Answered

I am working on an optimization problem which has 60 optimization variables to get minimum objective function. which solver would be best for this problem? right now I am using multistart with fmincon.

If your problem is smooth and you are looking for a global minimum, then those sound like the right choices, See Table for Choos...

I am working on an optimization problem which has 60 optimization variables to get minimum objective function. which solver would be best for this problem? right now I am using multistart with fmincon.

If your problem is smooth and you are looking for a global minimum, then those sound like the right choices, See Table for Choos...

environ un mois ago | 2

Answered

Find minimum Sum while keeping variable fix

I don't see your equality constraint. If the equality constraint is a nonlinear constraint, you need to pass it in the fmincon n...

Find minimum Sum while keeping variable fix

I don't see your equality constraint. If the equality constraint is a nonlinear constraint, you need to pass it in the fmincon n...

environ un mois ago | 0

Answered

optimization using genetic algorithm using two variable

You might want to try to copy the coding techniques in Optimize an ODE in Parallel. Alan Weiss MATLAB mathematical toolbox doc...

optimization using genetic algorithm using two variable

You might want to try to copy the coding techniques in Optimize an ODE in Parallel. Alan Weiss MATLAB mathematical toolbox doc...

environ un mois ago | 0

Answered

Genetic Algorithm objective and constraint

I suggest that you use the debugger. Put a break point in your nonlinear constraint function (or maybe elsewhere) and find out w...

Genetic Algorithm objective and constraint

I suggest that you use the debugger. Put a break point in your nonlinear constraint function (or maybe elsewhere) and find out w...

environ un mois ago | 0

Answered

lsqnonlin and true parameters value

I can help you with a) and f), but for the rest you will have to look elsewhere. The examples Fit ODE, Problem-Based and Fit an...

lsqnonlin and true parameters value

I can help you with a) and f), but for the rest you will have to look elsewhere. The examples Fit ODE, Problem-Based and Fit an...

environ un mois ago | 0

| accepted