5,585 total contributions since 2009

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

eigs bug for 0 as lowest eigenvalue in parts of the matlab versions

Increasing the SubspaceDimension also works, so it's hard to be certain if this is actually a bug. sp_matrix = sparse( diag(0:2...

eigs bug for 0 as lowest eigenvalue in parts of the matlab versions

Increasing the SubspaceDimension also works, so it's hard to be certain if this is actually a bug. sp_matrix = sparse( diag(0:2...

1 day ago | 1

| accepted

Answered

eigs bug for 0 as lowest eigenvalue in parts of the matlab versions

Here is another possible workaround based on svds [V,d] = svds(sp_matrix.'*sp_matrix,1,'smallest'); Eg_sp= V'*sp_matrix*V;

eigs bug for 0 as lowest eigenvalue in parts of the matlab versions

Here is another possible workaround based on svds [V,d] = svds(sp_matrix.'*sp_matrix,1,'smallest'); Eg_sp= V'*sp_matrix*V;

1 day ago | 0

Answered

eigs bug for 0 as lowest eigenvalue in parts of the matlab versions

If the matrix is known a priori to be diagonal, it makes more sense to just do Eg_sp=min(abs(diag(sp_matrix)))

eigs bug for 0 as lowest eigenvalue in parts of the matlab versions

If the matrix is known a priori to be diagonal, it makes more sense to just do Eg_sp=min(abs(diag(sp_matrix)))

2 days ago | 0

Answered

Subscripted assignment dimension mismatch.

sxint(3,:) is not a scalar, but the left hand side S(i,1) is a scalar location.

Subscripted assignment dimension mismatch.

sxint(3,:) is not a scalar, but the left hand side S(i,1) is a scalar location.

2 days ago | 0

Answered

eigs bug for 0 as lowest eigenvalue in parts of the matlab versions

That does seem like an interesting bug. This seems to work though Eg_sp = eigs(sp_matrix,1,1e-12);

eigs bug for 0 as lowest eigenvalue in parts of the matlab versions

That does seem like an interesting bug. This seems to work though Eg_sp = eigs(sp_matrix,1,1e-12);

2 days ago | 0

Answered

Create a vector of vector exponents

No, a for-loop is fastest, but you want to implement it the right way, with a recursive update, m=repmat(v*A,n,1); %pre-allocat...

Create a vector of vector exponents

No, a for-loop is fastest, but you want to implement it the right way, with a recursive update, m=repmat(v*A,n,1); %pre-allocat...

2 days ago | 0

| accepted

Answered

subtract each pixel from average pixel value

It seems deceptively simple... Image-meanValue

subtract each pixel from average pixel value

It seems deceptively simple... Image-meanValue

2 days ago | 0

| accepted

Answered

fmincon should return a scalar value error

AC Power reetuns 8760*1 value If you know that your objective function code is returning a 8760x1 value, it should be pretty cl...

fmincon should return a scalar value error

AC Power reetuns 8760*1 value If you know that your objective function code is returning a 8760x1 value, it should be pretty cl...

3 days ago | 0

| accepted

Answered

Help in understanding optimization problems and solving them in Matlab (choise of appropriate solver)

You have to reformulate the problem with a convex objective to be certain that the Lagrangian minimization will behave as you're...

Help in understanding optimization problems and solving them in Matlab (choise of appropriate solver)

You have to reformulate the problem with a convex objective to be certain that the Lagrangian minimization will behave as you're...

3 days ago | 0

Answered

How to delete repeating column or duplicate column in MATLAB? of type Double

output_array=unique(input_array.','rows').'

How to delete repeating column or duplicate column in MATLAB? of type Double

output_array=unique(input_array.','rows').'

3 days ago | 2

| accepted

Answered

imregister translation mode, suppress registration in second direction

One hack that I can think of (and I emphasize that it is only a hack) is to deliberately corrupt the images with bright horizont...

imregister translation mode, suppress registration in second direction

One hack that I can think of (and I emphasize that it is only a hack) is to deliberately corrupt the images with bright horizont...

10 days ago | 0

| accepted

Answered

fsolve error using trustnleqn (line 28)

The initial point is [0;0;0;0;0;0;0;0;0]. Calling your function at this point gives something with NaNs >> equilibrium([0;0;0;0...

fsolve error using trustnleqn (line 28)

The initial point is [0;0;0;0;0;0;0;0;0]. Calling your function at this point gives something with NaNs >> equilibrium([0;0;0;0...

10 days ago | 0

| accepted

Answered

how to full fill the diagonal of a matrix by a vector

You can use interpMatrix from the File Exchange >> full(interpMatrix([1,2,3],1,5,1,'circ')) ans = 1 0 0 ...

how to full fill the diagonal of a matrix by a vector

You can use interpMatrix from the File Exchange >> full(interpMatrix([1,2,3],1,5,1,'circ')) ans = 1 0 0 ...

10 days ago | 1

Answered

For loop not running properly

17.99 does not have an exact binary floating point representation >> format long >> 17.99 ans = 17.989999999999998 So...

For loop not running properly

17.99 does not have an exact binary floating point representation >> format long >> 17.99 ans = 17.989999999999998 So...

10 days ago | 0

| accepted

Answered

Efficiently calculating the trace of a matrix product

Bt=B.'; traceProduct = A(:).'*Bt(:);

Efficiently calculating the trace of a matrix product

Bt=B.'; traceProduct = A(:).'*Bt(:);

11 days ago | 2

| accepted

Answered

Solving a linear equation using least-squares (Calibration Matrix)

Are these equations for projective transformations? If so, they are not really linear equations. They are accurate only up to so...

Solving a linear equation using least-squares (Calibration Matrix)

Are these equations for projective transformations? If so, they are not really linear equations. They are accurate only up to so...

11 days ago | 0

Answered

Using the global operator

No, you should not use global. Your code should look like this: function [ll, mm, kk]= ssample(e2,t,kk) ... end

Using the global operator

No, you should not use global. Your code should look like this: function [ll, mm, kk]= ssample(e2,t,kk) ... end

11 days ago | 0

| accepted

Answered

Best Way to Set Private Properties in Class Constructor

There's nothing wrong with the first version of the code you posted (the one with the nargin tests). Also, there is no reason yo...

Best Way to Set Private Properties in Class Constructor

There's nothing wrong with the first version of the code you posted (the one with the nargin tests). Also, there is no reason yo...

11 days ago | 0

| accepted

Answered

unbounded problem in linprog but not in fmincon

If the problem is as below, then I obtain the same solution essentially from both linprog and fmincon FUN= @(x) 150*x(1)+230*x(...

unbounded problem in linprog but not in fmincon

If the problem is as below, then I obtain the same solution essentially from both linprog and fmincon FUN= @(x) 150*x(1)+230*x(...

12 days ago | 0

Answered

Solving a linear equation using least-squares (Calibration Matrix)

W=[ones(1,n);X]; Z=F/W; A=Z(:,1); C=Z(:,2:end);

Solving a linear equation using least-squares (Calibration Matrix)

W=[ones(1,n);X]; Z=F/W; A=Z(:,1); C=Z(:,2:end);

12 days ago | 0

| accepted

Answered

fmincon unidentified variable in objective function

This should fix it, FUN= @(x)150*x(1)+230*x(2)+260*x(3)+238*x(5)-170*x(4)+210*x(7)-150*x(6)-36*x(8)-10*x(9); however, you shou...

fmincon unidentified variable in objective function

This should fix it, FUN= @(x)150*x(1)+230*x(2)+260*x(3)+238*x(5)-170*x(4)+210*x(7)-150*x(6)-36*x(8)-10*x(9); however, you shou...

12 days ago | 1

| accepted

Answered

Dimensions of arrays being concatenated are not consistent.

No need to loop, E=[0,40,65,55,inf]; [binct,~,G]=histcounts(p,E); ypartmean=splitapply(@mean,d,G);

Dimensions of arrays being concatenated are not consistent.

No need to loop, E=[0,40,65,55,inf]; [binct,~,G]=histcounts(p,E); ypartmean=splitapply(@mean,d,G);

12 days ago | 0

| accepted

Answered

Choosing a GPU workstation

Yes, the GPU uses its own local RAM to do computations (or at least that's when it is fast), so you should get a graphics card w...

Choosing a GPU workstation

Yes, the GPU uses its own local RAM to do computations (or at least that's when it is fast), so you should get a graphics card w...

12 days ago | 1

| accepted

Answered

How to give index to each replicated number?

Yes correct! In that case, S=sparse( 1:numel(a), findgroups(a), 1); S(S>0)=nonzeros( cumsum(sort(S)) ); b=full(s...

How to give index to each replicated number?

Yes correct! In that case, S=sparse( 1:numel(a), findgroups(a), 1); S(S>0)=nonzeros( cumsum(sort(S)) ); b=full(s...

12 days ago | 0

Answered

How to give index to each replicated number?

G=cumsum(diff([0,a])~=0); b=cell2mat(splitapply(@(g){1:numel(g)}, G,G));

How to give index to each replicated number?

G=cumsum(diff([0,a])~=0); b=cell2mat(splitapply(@(g){1:numel(g)}, G,G));

12 days ago | 0

Answered

Can I use the Jacobian provided by 'lsqnonlin' to compute the confidence intervals using 'nlparci'?

nlparci probably wouldn't be applicable if you used lsqnonlin with bound constraints, unless the bounds were inactive at the sol...

Can I use the Jacobian provided by 'lsqnonlin' to compute the confidence intervals using 'nlparci'?

nlparci probably wouldn't be applicable if you used lsqnonlin with bound constraints, unless the bounds were inactive at the sol...

15 days ago | 3

Answered

Insert variable name into another variable

You would re-write all of your code to look like this instead Faces=StarFace{Number} or this Faces=StarFace(Number) See also...

Insert variable name into another variable

You would re-write all of your code to look like this instead Faces=StarFace{Number} or this Faces=StarFace(Number) See also...

17 days ago | 1

Answered

Can you constrain the boundary conditions of a smoothing spline fit?

It's possible with this, https://www.mathworks.com/matlabcentral/fileexchange/24443-slm-shape-language-modeling

Can you constrain the boundary conditions of a smoothing spline fit?

It's possible with this, https://www.mathworks.com/matlabcentral/fileexchange/24443-slm-shape-language-modeling

17 days ago | 2

| accepted

Answered

How do you edit the fitoptions of a cfit object generated by the fit function?

No, you have to re-run the fit.

How do you edit the fitoptions of a cfit object generated by the fit function?

No, you have to re-run the fit.

17 days ago | 0

| accepted

Answered

Unable to generate same plot in MATLAB as generated in Curve Fitting tool app

Use the File>Print to Figure menu option.

Unable to generate same plot in MATLAB as generated in Curve Fitting tool app

Use the File>Print to Figure menu option.

17 days ago | 0

| accepted