Professional Interests: medical image processing, optimization algorithms

PLEASE NOTE: I do not read email sent through my author page. Please post questions about FEX submissions in their respective Comments section.

Answered

exp(-b/x) fit, inf problem when fitting

When x goes to zero, y should go to zero as well since anything power negative infinity is zero. Only if b>=0. g = fittype('ex...

exp(-b/x) fit, inf problem when fitting

When x goes to zero, y should go to zero as well since anything power negative infinity is zero. Only if b>=0. g = fittype('ex...

1 jour ago | 0

Answered

Sum subarrays of a Matrix.

Using sepblockfun https://www.mathworks.com/matlabcentral/fileexchange/48089-separable-block-wise-operations out=sepblockfun(A...

Sum subarrays of a Matrix.

Using sepblockfun https://www.mathworks.com/matlabcentral/fileexchange/48089-separable-block-wise-operations out=sepblockfun(A...

1 jour ago | 0

Question

Recover former state of Matlab Editor

I had quite a few files open in the Matlab Editor, but then inadvertently closed the editor window. Is there any way to re-open ...

3 jours ago | 1 answer | 0

Answered

How to calculate reveral summations

T=triu(ones(600)); I=(1:600).' + (1:600); ur=u(:).'; sigma= (T.*ur(I))*ur(1:600).';

How to calculate reveral summations

T=triu(ones(600)); I=(1:600).' + (1:600); ur=u(:).'; sigma= (T.*ur(I))*ur(1:600).';

3 jours ago | 0

Answered

How do I minimize a particular error function with lsqnonlin?

I don't see much that I think can be improved, except perhaps to formulate it so that your errors are differentiable - the lsqno...

How do I minimize a particular error function with lsqnonlin?

I don't see much that I think can be improved, except perhaps to formulate it so that your errors are differentiable - the lsqno...

3 jours ago | 0

| accepted

Answered

Get "@" back in a table header (R2017b) when saving

For example, T=table(1,'VariableNames',"name_domain") T.Properties.VariableNames = replace(T.Properties.VariableNames,'_','@...

Get "@" back in a table header (R2017b) when saving

For example, T=table(1,'VariableNames',"name_domain") T.Properties.VariableNames = replace(T.Properties.VariableNames,'_','@...

3 jours ago | 0

Answered

How to find splitapply for multiple values

[G, idOcc,idAge] =findgroups(x.occupation,x.Age)

How to find splitapply for multiple values

[G, idOcc,idAge] =findgroups(x.occupation,x.Age)

3 jours ago | 0

| accepted

Answered

Calculation the mean of a portion of a matrix

mean( z(2:end-1,2:end-1) ,'all')

Calculation the mean of a portion of a matrix

mean( z(2:end-1,2:end-1) ,'all')

3 jours ago | 0

| accepted

Answered

Rank one decomposition of a positive semi-definite matrix with inequality trace constraints

Is there any ways I could do the rank one decomposition of matrix X, such that The obvious answer seems to be to test each k to...

Rank one decomposition of a positive semi-definite matrix with inequality trace constraints

Is there any ways I could do the rank one decomposition of matrix X, such that The obvious answer seems to be to test each k to...

4 jours ago | 0

Answered

create matrix of integers with specified sum in column and row

You can try to formulate it as an integer program. M=optimvar('M',[4,4],'type','integer'); Constraints.columns=sum(M,1)<=col...

create matrix of integers with specified sum in column and row

You can try to formulate it as an integer program. M=optimvar('M',[4,4],'type','integer'); Constraints.columns=sum(M,1)<=col...

4 jours ago | 0

Answered

Finding Upper and Lower Indexes

A = [1;2;3;4;5;8]; B = [3.5;6;2]; c=interp1(A,A,B,'previous'); C=[c(:),c(:)+1]

Finding Upper and Lower Indexes

A = [1;2;3;4;5;8]; B = [3.5;6;2]; c=interp1(A,A,B,'previous'); C=[c(:),c(:)+1]

6 jours ago | 1

Answered

Row-normalizing large sparse matrix

The attempt you've posted will only work if vals_W is the same size as sum_vals, which can only occur when there is exactly one ...

Row-normalizing large sparse matrix

The attempt you've posted will only work if vals_W is the same size as sum_vals, which can only occur when there is exactly one ...

7 jours ago | 0

| accepted

Submitted

Object-oriented tools for fitting conics and quadrics

A tool set for fitting various conics and quadric surfaces, e.g., ellipses, cylinders, spheres, planes, and lines.

8 jours ago | 4 downloads |

Answered

Problem with fmincon fval

As Walter suggests, the multiple fvals per iteration are part of finite difference computations needed to evaluate the cost func...

Problem with fmincon fval

As Walter suggests, the multiple fvals per iteration are part of finite difference computations needed to evaluate the cost func...

8 jours ago | 0

Answered

Finding a row with a certain value and the next nth rows after that

Here's an approach that uses logical indexing only: A=[ 0.9058 0.2785 0.9706 0.4218 0.0357 0.1270 0.5469 ...

Finding a row with a certain value and the next nth rows after that

Here's an approach that uses logical indexing only: A=[ 0.9058 0.2785 0.9706 0.4218 0.0357 0.1270 0.5469 ...

8 jours ago | 1

| accepted

Question

Is there a way to interactively move a plot legend in a Live Editor figure?

In a normal figure window, I can reposition a plot legend by clicking and dragging. Is it possible to do something similar in a ...

9 jours ago | 1 answer | 1

Answered

Best way to move multiple objects with a preallocated matrix of coordinates?

I see movement when I do this, xl=xlim; for b=1:1:nStarts*nCarsPerStart for h=nStarts+1:nStarts:nStarts*(t+1) ...

Best way to move multiple objects with a preallocated matrix of coordinates?

I see movement when I do this, xl=xlim; for b=1:1:nStarts*nCarsPerStart for h=nStarts+1:nStarts:nStarts*(t+1) ...

9 jours ago | 0

| accepted

Answered

averaging images and subtracting one image from the other

For example, setA=rand(256,256,10); setB=rand(256,256,10); imageA=mean(setA,3); imageB=mean(setB,3); result = imageB-im...

averaging images and subtracting one image from the other

For example, setA=rand(256,256,10); setB=rand(256,256,10); imageA=mean(setA,3); imageB=mean(setB,3); result = imageB-im...

9 jours ago | 0

| accepted

Answered

Trying to do a curve fitting but having trouble getting it to work.

Do you happen to know why I keep getting the "local minimum possible" Because lsqcurvefit thinks it succeeded in finding a sol...

Trying to do a curve fitting but having trouble getting it to work.

Do you happen to know why I keep getting the "local minimum possible" Because lsqcurvefit thinks it succeeded in finding a sol...

9 jours ago | 0

Answered

Need help creating a variable with different names or a method for saving a variable to different names.

for k=1:K S(k).A(i,j) = readline(s); S(k).B(i,j) = readline(s); S(k).C(i,j) = readline(s); S(k).D(i,j) = readline(s); e...

Need help creating a variable with different names or a method for saving a variable to different names.

for k=1:K S(k).A(i,j) = readline(s); S(k).B(i,j) = readline(s); S(k).C(i,j) = readline(s); S(k).D(i,j) = readline(s); e...

9 jours ago | 0

| accepted

Answered

Index in position 1 exceeds array bounds.

Notice! x=rand(5,500); x=x(1:256,1:256);

Index in position 1 exceeds array bounds.

Notice! x=rand(5,500); x=x(1:256,1:256);

9 jours ago | 0

Answered

I want to combine two plots, one is cone and one is semi-circle

Use hold().

I want to combine two plots, one is cone and one is semi-circle

Use hold().

9 jours ago | 1

Answered

How to find number of elements using a for loop

for i=1:numel(A), A(i)>2, end

How to find number of elements using a for loop

for i=1:numel(A), A(i)>2, end

9 jours ago | 0

Answered

extracting two specific columns from matrix (Indexing)

x=xxx(:,[1,3])

extracting two specific columns from matrix (Indexing)

x=xxx(:,[1,3])

9 jours ago | 0

| accepted

Question

Inheriting help documentation from an Abstract base class method

Suppose I have a base class with many subclasses, each of which will override an abstract method in the base class, classdef Ba...

10 jours ago | 0 answers | 1

Answered

Nonlinear fit instead of setting up a non linear equation?

Your problem appears to be a linear fit. I see no non-linear dependence on x or y. You should therefore use lsqlin, if you also ...

Nonlinear fit instead of setting up a non linear equation?

Your problem appears to be a linear fit. I see no non-linear dependence on x or y. You should therefore use lsqlin, if you also ...

10 jours ago | 1

| accepted

Answered

Would you let me know one of the fastest (or the most efficient) ways to merge data?

final_data = cell2mat(data(:))

Would you let me know one of the fastest (or the most efficient) ways to merge data?

final_data = cell2mat(data(:))

11 jours ago | 0

| accepted

Answered

Inf computed by model function, fitting cannot continue. Try using or tightening upper and lower bounds on coefficients.

the fit created by the function y(x) is very, very poor. Basically, it is straight line The fit looks very good considering you...

Inf computed by model function, fitting cannot continue. Try using or tightening upper and lower bounds on coefficients.

the fit created by the function y(x) is very, very poor. Basically, it is straight line The fit looks very good considering you...

11 jours ago | 1

| accepted