Answered

Error With Least Squares Approximation

Few other recommendations beside what Torsen and John altrady told you: Lower the polynomial order, do not go above 7, But you ...

Error With Least Squares Approximation

Few other recommendations beside what Torsen and John altrady told you: Lower the polynomial order, do not go above 7, But you ...

20 days ago | 0

Submitted

Free-knot spline approximation

Least squares approximation of 1D data using free-knots spline

20 days ago | 31 downloads |

Question

ANY and ALL with NaN array

According to any documentation page, it " determines if any element is a nonzero number". OK to me NaN is nonzero number so I do...

21 days ago | 1 answer | 0

Answered

Make histogram from multiple cells.

% Generate dummy test data C=arrayfun(@(varargin)randn(5+randi(20),24), 1:30, 'unif', 0); A=cell2mat(C(:)); % m x 24, like c...

Make histogram from multiple cells.

% Generate dummy test data C=arrayfun(@(varargin)randn(5+randi(20),24), 1:30, 'unif', 0); A=cell2mat(C(:)); % m x 24, like c...

22 days ago | 0

| accepted

Answered

My First Derivative is not correctly calculated in matlab

syms yApprox(X1) a0 a1 a2 a3 a4 a5 a6 yApprox = a0 + a1.*X1 + a2.*(X1).^2 + a3.*(X1).^3+ a4.*(X1).^4 + a5.*(X1).^5 + a6.*(X1).^...

My First Derivative is not correctly calculated in matlab

syms yApprox(X1) a0 a1 a2 a3 a4 a5 a6 yApprox = a0 + a1.*X1 + a2.*(X1).^2 + a3.*(X1).^3+ a4.*(X1).^4 + a5.*(X1).^5 + a6.*(X1).^...

23 days ago | 0

| accepted

Answered

for loop vectorization of De Boor algorithm

Old question, but I'll reply it for future readers that wonder how to tackle the same problem. I did a quick and dirdty vectori...

for loop vectorization of De Boor algorithm

Old question, but I'll reply it for future readers that wonder how to tackle the same problem. I did a quick and dirdty vectori...

23 days ago | 0

Answered

Increase plot speed in MATLAB

Try to look at animatedline command Also look at drawnow with or without limitedrate option, it somewhat has certain effect, ev...

Increase plot speed in MATLAB

Try to look at animatedline command Also look at drawnow with or without limitedrate option, it somewhat has certain effect, ev...

26 days ago | 0

Answered

Error calculating inverse for matrix block of LU decomposition of a full row rank matrix

It's a little bit perturbant, but reading the lu doc I don't see why U1 must have the same rank than U (or A, meaning 100 on you...

Error calculating inverse for matrix block of LU decomposition of a full row rank matrix

It's a little bit perturbant, but reading the lu doc I don't see why U1 must have the same rank than U (or A, meaning 100 on you...

26 days ago | 0

| accepted

Answered

de boor bspline algorithm

A very late answer, but as I deal myself with de Boor algorithm I look what exists out there. If the question concerns the 3D a...

de boor bspline algorithm

A very late answer, but as I deal myself with de Boor algorithm I look what exists out there. If the question concerns the 3D a...

27 days ago | 0

Answered

For a repeated eigenvalue only one eigenvctor is being returned

Indeed in case of eigenvalue with multiplicity > 1; the problemie is numerical challenging and MATLAB might fail to find the cor...

For a repeated eigenvalue only one eigenvctor is being returned

Indeed in case of eigenvalue with multiplicity > 1; the problemie is numerical challenging and MATLAB might fail to find the cor...

27 days ago | 0

Answered

How "chol" and "qz" MATLAB algorithms are utilised in the "eig" MATLAB function?

See stage 3 described in https://netlib.org/lapack/lug/node56.html or somewhere in this file

How "chol" and "qz" MATLAB algorithms are utilised in the "eig" MATLAB function?

See stage 3 described in https://netlib.org/lapack/lug/node56.html or somewhere in this file

1 month ago | 0

| accepted

Answered

steps to convert spline from B-form to pp-form (fn2fm)

The pp-form stores the polynomial of each each subiterval, the variable is x := (t-ti) where ti is the left knot of interval #i....

steps to convert spline from B-form to pp-form (fn2fm)

The pp-form stores the polynomial of each each subiterval, the variable is x := (t-ti) where ti is the left knot of interval #i....

1 month ago | 0

| accepted

Answered

How to extract a matrix of values from cell array of cell arrays of structs

Generate data (Thanks Voss) n1 = 3; n2 = 4; mse = cell(1,n1); for ii = 1:n1 mse{ii} = cell(1,n2); % row vector for...

How to extract a matrix of values from cell array of cell arrays of structs

Generate data (Thanks Voss) n1 = 3; n2 = 4; mse = cell(1,n1); for ii = 1:n1 mse{ii} = cell(1,n2); % row vector for...

1 month ago | 0

| accepted

Answered

Where is dlfeval tracing information stored?

I don't kwow the dlarray class but IMO it looks like the creating function handle f is stored in some sort of private property i...

Where is dlfeval tracing information stored?

I don't kwow the dlarray class but IMO it looks like the creating function handle f is stored in some sort of private property i...

2 months ago | 0

Question

Default rank revealing tolerance

In few MATLAB functions such as RANK, PINV, ORTH there is a parameter TOL. MATLAB uses the default value as % tol = max(size(A...

2 months ago | 0 answers | 1

Answered

plot() and fplot() give different results for the same function

Your coding is wrong: Replace "/" with "./"

plot() and fplot() give different results for the same function

Your coding is wrong: Replace "/" with "./"

2 months ago | 0

| accepted

Answered

getting warning while using eigs on the matrices obtained from freefem++

Try to follow the suggestion you get " Consider specifying a perturbed numeric sigma value to improve the condition of the matri...

getting warning while using eigs on the matrices obtained from freefem++

Try to follow the suggestion you get " Consider specifying a perturbed numeric sigma value to improve the condition of the matri...

2 months ago | 0

Answered

cubic spline interpolation - mixed boundary conditions possible?

In principle yes. The spline on the second segment (b,c) are piecewise cubic poynomials; starting from b sequentially to each in...

cubic spline interpolation - mixed boundary conditions possible?

In principle yes. The spline on the second segment (b,c) are piecewise cubic poynomials; starting from b sequentially to each in...

2 months ago | 0

| accepted

Answered

How to convert /map uint16 to int16 maintaining the dynamic range so that [0] and [2^16-1] in uint16 beomes [-2^15] and [2^15-1] without going through double

If you have C compiler this simple mex file will do the job: /*****************************************************************...

How to convert /map uint16 to int16 maintaining the dynamic range so that [0] and [2^16-1] in uint16 beomes [-2^15] and [2^15-1] without going through double

If you have C compiler this simple mex file will do the job: /*****************************************************************...

2 months ago | 1

Answered

Faster/better conversion than using swapbytes typecast on udp byte datastream?

Try this: %% Exampel data data =[0 5 187 8 1 136 61 48 0 0 0 0 1 67 58 246 255 127 176 76 10 231 55 137 255 240 253 91 255 22...

Faster/better conversion than using swapbytes typecast on udp byte datastream?

Try this: %% Exampel data data =[0 5 187 8 1 136 61 48 0 0 0 0 1 67 58 246 255 127 176 76 10 231 55 137 255 240 253 91 255 22...

2 months ago | 0

| accepted

Answered

How to add extra columns at specified locations

After computing dataset2 dataset = [1, 2, 3, nan, 4, 5, nan, 6, 7]; % ... dataset2 = [2, 5, 7, 8, 4, 8, 10]; do this: da...

How to add extra columns at specified locations

After computing dataset2 dataset = [1, 2, 3, nan, 4, 5, nan, 6, 7]; % ... dataset2 = [2, 5, 7, 8, 4, 8, 10]; do this: da...

2 months ago | 0

| accepted

Answered

scrit can't error related to Sortrows and matrix

Replace initialization statements sol=[]; % Start with an empty matrix of feasible solutions data = []; with sol=zeros(...

scrit can't error related to Sortrows and matrix

Replace initialization statements sol=[]; % Start with an empty matrix of feasible solutions data = []; with sol=zeros(...

2 months ago | 0

Answered

Max Recursion even when using clear

Don't name your script function clear.m. In command windows type which clear then rename or remove that file to restore the bu...

Max Recursion even when using clear

Don't name your script function clear.m. In command windows type which clear then rename or remove that file to restore the bu...

2 months ago | 0

Answered

Upsampling a matrix with zero elements

x = [1 2 3 4; 5 6 7 8] kron(x, [1 0; 0 0])

Upsampling a matrix with zero elements

x = [1 2 3 4; 5 6 7 8] kron(x, [1 0; 0 0])

2 months ago | 1

Answered

warning in Polyval - Warning: Matrix is singular, close to singular or badly scaled. Results may be inaccurate. RCOND = NaN.

Check your data with this command all(isfinite(rawData(:))) if it returns 0/false then your data is not valid

warning in Polyval - Warning: Matrix is singular, close to singular or badly scaled. Results may be inaccurate. RCOND = NaN.

Check your data with this command all(isfinite(rawData(:))) if it returns 0/false then your data is not valid

2 months ago | 0

Answered

How to detect whether a figure is created by uifigure()

This command sreturns true for uifigure handle fig matlab.ui.internal.isUIFigure(fig)

How to detect whether a figure is created by uifigure()

This command sreturns true for uifigure handle fig matlab.ui.internal.isUIFigure(fig)

3 months ago | 0

| accepted

Answered

Trying to optimise this loop

function WriteOutValuesToAFile(filename,M) fileID = fopen(filename,'wt'); for k=1:size(M,3); s = sprintf('%f,', M(:,:,k))...

Trying to optimise this loop

function WriteOutValuesToAFile(filename,M) fileID = fopen(filename,'wt'); for k=1:size(M,3); s = sprintf('%f,', M(:,:,k))...

3 months ago | 0

Answered

How to do a layered matrix--vector multiplication without a for loop?

NNN = 10; RR = rand(NNN, 3); TT = rand(3,3,NNN); for iii = 1:1:NNN TRR(iii,:) = TT(:,:,iii) * RR(iii, :).'; % Bug fix ...

How to do a layered matrix--vector multiplication without a for loop?

NNN = 10; RR = rand(NNN, 3); TT = rand(3,3,NNN); for iii = 1:1:NNN TRR(iii,:) = TT(:,:,iii) * RR(iii, :).'; % Bug fix ...

3 months ago | 2

| accepted

Answered

Convert Quaternion to Euler angle extrinsically

Here is the answer by MATLAB code: extrinsic angles is flip intrinsic angles % Generate random unit quaternion q = quaternion(...

Convert Quaternion to Euler angle extrinsically

Here is the answer by MATLAB code: extrinsic angles is flip intrinsic angles % Generate random unit quaternion q = quaternion(...

3 months ago | 1