Professional Interests: computational physics, applied mathematics

Answered

How to draw an ellipse

The code discussed at https://www.mathworks.com/matlabcentral/answers/497966-how-to-find-the-vertical-height-of-an-ellipse-y-axi...

How to draw an ellipse

The code discussed at https://www.mathworks.com/matlabcentral/answers/497966-how-to-find-the-vertical-height-of-an-ellipse-y-axi...

8 days ago | 0

Answered

Eigenvalues corresponds to eigenvectors

Absolutely standard: L(k,k) ~ V(:,k). You can check it with the code: for k=1:size(h,1) disp(strcat('k=',num2str(k),'; h*v-la...

Eigenvalues corresponds to eigenvectors

Absolutely standard: L(k,k) ~ V(:,k). You can check it with the code: for k=1:size(h,1) disp(strcat('k=',num2str(k),'; h*v-la...

17 days ago | 0

| accepted

Answered

Subs Command for Multivariable Function

If e^ is the exponent and if you want to estimate the result numerically (from the question it is not quite clear what you want)...

Subs Command for Multivariable Function

If e^ is the exponent and if you want to estimate the result numerically (from the question it is not quite clear what you want)...

18 days ago | 0

Answered

Partial least-squares regression evalutaion

The equation considered in PLS of Matlab is rather Y = X(:,:)*B(2:end,:) + repmat( B(1,:),n,1 ) + E, where: X is an n-by-...

Partial least-squares regression evalutaion

The equation considered in PLS of Matlab is rather Y = X(:,:)*B(2:end,:) + repmat( B(1,:),n,1 ) + E, where: X is an n-by-...

19 days ago | 0

| accepted

Answered

Normalization of complex eigenvector

You probably want V(:,1) = V(:,1)/norm(V(:,1)); Besides, if your theta and phi are supposed to be real, the overall computatio...

Normalization of complex eigenvector

You probably want V(:,1) = V(:,1)/norm(V(:,1)); Besides, if your theta and phi are supposed to be real, the overall computatio...

24 days ago | 0

| accepted

Answered

Eigenvector without calling eigenvalues

This is not an error but a warning that you do not use the eigenvalues, which influences nothing. If you want to avoid it, subst...

Eigenvector without calling eigenvalues

This is not an error but a warning that you do not use the eigenvalues, which influences nothing. If you want to avoid it, subst...

24 days ago | 0

Answered

Can I check if factored and expanded polynomial forms are equal?

simplify((w^2-k^2) - (w+k)*(w-k)) ans = 0 expand((w+k)*(w-k)) ans = w^2 - k^2 factor(w^2-k^2) ans = [ -1, k - w...

Can I check if factored and expanded polynomial forms are equal?

simplify((w^2-k^2) - (w+k)*(w-k)) ans = 0 expand((w+k)*(w-k)) ans = w^2 - k^2 factor(w^2-k^2) ans = [ -1, k - w...

24 days ago | 0

| accepted

Answered

How to call some object handle from opening function of figure?

u0=findobj('Tag','lst1','Style','listbox','Parent',gcf); % search for all the object properties you used to identify it set(u0,...

How to call some object handle from opening function of figure?

u0=findobj('Tag','lst1','Style','listbox','Parent',gcf); % search for all the object properties you used to identify it set(u0,...

25 days ago | 0

| accepted

Answered

Find orthogonal complement for given matrix

B=null(A'); The matrix transposing is important!

Find orthogonal complement for given matrix

B=null(A'); The matrix transposing is important!

26 days ago | 1

| accepted

Answered

How can I write a loop, which can make the parameter "V" change with a cyclic variable "i"?

Use Matlab functions "strcat" and "num2str" within a loop.

How can I write a loop, which can make the parameter "V" change with a cyclic variable "i"?

Use Matlab functions "strcat" and "num2str" within a loop.

27 days ago | 0

| accepted

Answered

How do I calculate acceleration with velocity and the code given?

% sample data t=0:0.1:10; % time x=3+2*t+t.^2; % coordinate [t,ind]=sort(t); % in a case time is not in an ascending order ...

How do I calculate acceleration with velocity and the code given?

% sample data t=0:0.1:10; % time x=3+2*t+t.^2; % coordinate [t,ind]=sort(t); % in a case time is not in an ascending order ...

1 month ago | 0

Answered

How to assign numbers to letters of the alphabet. Ex. A/a=1 B/b=2

The ASCII codes can be found with the function double double('A'); 65 If you want to get its ordinal number in the alphabet, ...

How to assign numbers to letters of the alphabet. Ex. A/a=1 B/b=2

The ASCII codes can be found with the function double double('A'); 65 If you want to get its ordinal number in the alphabet, ...

1 month ago | 0

Answered

How to sum columns of a matrix in MATLAB in specified number of columns?

B=squeeze( sum( reshape(A,size(A,1),10,[]) ,2) );

How to sum columns of a matrix in MATLAB in specified number of columns?

B=squeeze( sum( reshape(A,size(A,1),10,[]) ,2) );

1 month ago | 0

| accepted

Answered

Create a gaussian from x values

L=20; V = linspace(0,L,301); sigma = L/20; T0 = exp(-((V-mean(V))/sigma).^2/2); plot(V,T0);

Create a gaussian from x values

L=20; V = linspace(0,L,301); sigma = L/20; T0 = exp(-((V-mean(V))/sigma).^2/2); plot(V,T0);

1 month ago | 1

| accepted

Answered

extract a double array to multiple arrays

b1=aa(1:4:end); b2=aa(2:4:end); b3=aa(3:4:end); b4=aa(4:4:end);

extract a double array to multiple arrays

b1=aa(1:4:end); b2=aa(2:4:end); b3=aa(3:4:end); b4=aa(4:4:end);

1 month ago | 1

| accepted

Answered

Importing Latex Code into Matlab Scripts

In live script, go to \insert\Equation\Latex equation or just press Cntrl+Shift+L or type your Latex equation embraced by $ symb...

Importing Latex Code into Matlab Scripts

In live script, go to \insert\Equation\Latex equation or just press Cntrl+Shift+L or type your Latex equation embraced by $ symb...

2 months ago | 0

Submitted

DigBit

Image Digitizer: Turn your scanned plots to the numeric form

2 months ago | 0 downloads |

Answered

How can I calculate the integral2 using sum of prod?

symprod does not support the element-wise multiplication, which is needed for integral2. You sholuld probably define your fun v...

How can I calculate the integral2 using sum of prod?

symprod does not support the element-wise multiplication, which is needed for integral2. You sholuld probably define your fun v...

2 months ago | 1

| accepted

Answered

sorting 2 array based on elements of other array

[A,ind]=sort(A); B=B(ind);

sorting 2 array based on elements of other array

[A,ind]=sort(A); B=B(ind);

2 months ago | 0

| accepted

Answered

How to find the maximum distance between two curves?

(1) Specify, what do you mean by the distance between curves. If this is the distance between any points of them, a special (tho...

How to find the maximum distance between two curves?

(1) Specify, what do you mean by the distance between curves. If this is the distance between any points of them, a special (tho...

2 months ago | 2

| accepted

Answered

New to Matlab- How to import a text file for analysis that contains numbers and characters?

If you need it only once (or few times), use the Import utility of Matlab (Home\Import Data - a green downward arrow); chose the...

New to Matlab- How to import a text file for analysis that contains numbers and characters?

If you need it only once (or few times), use the Import utility of Matlab (Home\Import Data - a green downward arrow); chose the...

2 months ago | 0

Answered

help speeding up an incrementing loop

At the first glance, you can replace numerator=0; for k = 1:n_trc numerator = nume...

help speeding up an incrementing loop

At the first glance, you can replace numerator=0; for k = 1:n_trc numerator = nume...

2 months ago | 0

Answered

How to find the vertical height of an ellipse/y axis length of ellipse.

In fact, the height and the width of an ellipse can be easily and more accurately found analytically. Here is Matlab program rea...

How to find the vertical height of an ellipse/y axis length of ellipse.

In fact, the height and the width of an ellipse can be easily and more accurately found analytically. Here is Matlab program rea...

2 months ago | 1

| accepted

Answered

How to expand the Gaussian Kernel into series of eigenfunctions?

Fitst of all, you should specify --nothing can be done with an arbtrary function. Then, you can try to calculate the integrals a...

How to expand the Gaussian Kernel into series of eigenfunctions?

Fitst of all, you should specify --nothing can be done with an arbtrary function. Then, you can try to calculate the integrals a...

2 months ago | 1

| accepted

Answered

Turn function to zero when it reach specific value.

x = linspace (0,180,1800); y = 1. * ((exp (-0.01 .* x)) .* (sin (x)-(x .* cos (x / 2)))); Pu = 30; y(find(y>=Pu,1):end)=0; ...

Turn function to zero when it reach specific value.

x = linspace (0,180,1800); y = 1. * ((exp (-0.01 .* x)) .* (sin (x)-(x .* cos (x / 2)))); Pu = 30; y(find(y>=Pu,1):end)=0; ...

2 months ago | 0

| accepted

Answered

How to show part of a nested structure as a vector?

sn='seine_1'; % sn=fieldnames(all_seines); % sn=sn{1}; nm=fieldnames(all_seines.(sn)); N=numel(nm); A=zeros(N,1); for k=...

How to show part of a nested structure as a vector?

sn='seine_1'; % sn=fieldnames(all_seines); % sn=sn{1}; nm=fieldnames(all_seines.(sn)); N=numel(nm); A=zeros(N,1); for k=...

2 months ago | 0

Answered

Values over 1000 use e

read matlab help on the "format" command. try "format longg".

Values over 1000 use e

read matlab help on the "format" command. try "format longg".

2 months ago | 0

Answered

how to create Chi squared distribution using mean and variance?

see https://www.mathworks.com/help/stats/chi-square-distribution-1.html

how to create Chi squared distribution using mean and variance?

see https://www.mathworks.com/help/stats/chi-square-distribution-1.html

2 months ago | 0

Answered

Not spotting where/how index exceeds number of array elements

Is this what you want? Live script vesion is attached. xmesh = linspace(0.01,0.05,5); solinit = bvpinit(xmesh, @guess); sol...

Not spotting where/how index exceeds number of array elements

Is this what you want? Live script vesion is attached. xmesh = linspace(0.01,0.05,5); solinit = bvpinit(xmesh, @guess); sol...

2 months ago | 0

| accepted