Answered

How do I find a D matrix that satisfies [K]*inv(M)*[D]==[D]*inv(M)*[K] ?

Hi Ikram, Using N in place of inv(M) for simplicity, you are looking for a D such that K*N*D = D*N*K Let K*N and N*K have eig...

How do I find a D matrix that satisfies [K]*inv(M)*[D]==[D]*inv(M)*[K] ?

Hi Ikram, Using N in place of inv(M) for simplicity, you are looking for a D such that K*N*D = D*N*K Let K*N and N*K have eig...

2 days ago | 0

| accepted

Answered

Is there anyway to calculate the position in coordinates of points using a n*n distance matrix?

Hi Meme, first of all, if all you have are the distances between points, you can translate the all the points together in the 2...

Is there anyway to calculate the position in coordinates of points using a n*n distance matrix?

Hi Meme, first of all, if all you have are the distances between points, you can translate the all the points together in the 2...

6 days ago | 0

Answered

Summation subroutine giving result as NaN (not a number)!

Hi Louis, the basic problem is that in the product calculation you have for l=0:n if (l~=j) && (l~=i) && (l~=m) prod=p...

Summation subroutine giving result as NaN (not a number)!

Hi Louis, the basic problem is that in the product calculation you have for l=0:n if (l~=j) && (l~=i) && (l~=m) prod=p...

8 days ago | 0

| accepted

Answered

Replace diagonals in a matrix

Hi Hasan, here is one way r = rand(7,7) n = size(r,1); m = (n-1)/2; a = (-m:m)-(-m:m)'; r(abs(a)>1)=0 assuming the main...

Replace diagonals in a matrix

Hi Hasan, here is one way r = rand(7,7) n = size(r,1); m = (n-1)/2; a = (-m:m)-(-m:m)'; r(abs(a)>1)=0 assuming the main...

8 days ago | 0

Answered

unitConvert isn't converting to Hz

Hi Henrique. you can get there with u = symunit L = 10*u.uH; C = 10*u.fF; unitConvert(simplify(unitConvert(1/sqrt((L*C)), u...

unitConvert isn't converting to Hz

Hi Henrique. you can get there with u = symunit L = 10*u.uH; C = 10*u.fF; unitConvert(simplify(unitConvert(1/sqrt((L*C)), u...

11 days ago | 0

| accepted

Answered

Optimizing calculation of eigenvectors and eigenvalues

Hi Guiseppi, there is nothing guaranteed about the order of the eigenvalues that are produced by eig. And if you do the first ...

Optimizing calculation of eigenvectors and eigenvalues

Hi Guiseppi, there is nothing guaranteed about the order of the eigenvalues that are produced by eig. And if you do the first ...

11 days ago | 2

| accepted

Answered

Creating Polar mesh in matlab

Hi Reza, % make a rectangular grid of r and theta, % then define x and y in the usual way rr = 0:1:20; thth = (0:.05:1)*2*p...

Creating Polar mesh in matlab

Hi Reza, % make a rectangular grid of r and theta, % then define x and y in the usual way rr = 0:1:20; thth = (0:.05:1)*2*p...

13 days ago | 0

| accepted

Answered

Why when calculating impulse response of a filter (system) I get a “tail” at the end

Hi Yurii, It looks like the 'lowpass' function is doing something unexpected, which causes this. In order to get the transfer ...

Why when calculating impulse response of a filter (system) I get a “tail” at the end

Hi Yurii, It looks like the 'lowpass' function is doing something unexpected, which causes this. In order to get the transfer ...

14 days ago | 0

Answered

Peculiar Result with Solve (Symbolic Math Toolbox) with Three Equations and Three Unknowns. Can Anyone Explain?

Hi Paul, certainly in the first case, [eps/5, 4eps/5,eps] is a solution for any nonzero eps, and if you let eps--> 0 and invo...

Peculiar Result with Solve (Symbolic Math Toolbox) with Three Equations and Three Unknowns. Can Anyone Explain?

Hi Paul, certainly in the first case, [eps/5, 4eps/5,eps] is a solution for any nonzero eps, and if you let eps--> 0 and invo...

18 days ago | 0

Answered

Find the common eigenvectors and eigenvalues between 2 matrices

Hi petit, Eigenvectors calculated by Matlab are normalized, but neither (a) the the overall phase of each one or (b) the order ...

Find the common eigenvectors and eigenvalues between 2 matrices

Hi petit, Eigenvectors calculated by Matlab are normalized, but neither (a) the the overall phase of each one or (b) the order ...

22 days ago | 0

| accepted

Answered

means squared error and dB

Hi MatG, As with a lot of things, the answer depends. In your example, which is the first case below, the variation in power i...

means squared error and dB

Hi MatG, As with a lot of things, the answer depends. In your example, which is the first case below, the variation in power i...

26 days ago | 0

| accepted

Answered

Given One Partition of a Matrix, What is the Best Way to Find a Second Partition that Ensures the Matrix is Nonsingular?

Hi Paul, First of all, if C is not close to full rank, there can be numerical problems with most any calculation involving C. ...

Given One Partition of a Matrix, What is the Best Way to Find a Second Partition that Ensures the Matrix is Nonsingular?

Hi Paul, First of all, if C is not close to full rank, there can be numerical problems with most any calculation involving C. ...

26 days ago | 0

| accepted

Answered

convolution of exponential with unit step ... using conv command

Hi Muhammed, since your reply is a copy of the the orignal code, nothing changes. Anyway, the reason for the effect is that wh...

convolution of exponential with unit step ... using conv command

Hi Muhammed, since your reply is a copy of the the orignal code, nothing changes. Anyway, the reason for the effect is that wh...

26 days ago | 0

Answered

Please help me convert spectral integral equation to matlab code.

Hi Emma, Seeing the actual problem helped. It sounds like this is supposed to be a region-1-only problem with region 2 out of ...

Please help me convert spectral integral equation to matlab code.

Hi Emma, Seeing the actual problem helped. It sounds like this is supposed to be a region-1-only problem with region 2 out of ...

29 days ago | 0

Answered

Logical Indexing via multiplication

Hi Inna, not done by multiplication, but: ind = find(A==1) rows_you_want = B(ind,:) the colon means to take every column in w...

Logical Indexing via multiplication

Hi Inna, not done by multiplication, but: ind = find(A==1) rows_you_want = B(ind,:) the colon means to take every column in w...

1 month ago | 0

| accepted

Answered

how to check if a point is inside an rotated ellipsoid?

Hi Amin, Since it appears that you know the rotation matrix, you can rotate the point p in question in the opposite direction t...

how to check if a point is inside an rotated ellipsoid?

Hi Amin, Since it appears that you know the rotation matrix, you can rotate the point p in question in the opposite direction t...

1 month ago | 1

| accepted

Answered

Since the matrix is near to singular I have tried pinv to solve but getting too deviated answers. Can you please let me know if there is any other way to solve the equations

Hi VS, Usually in these cases, finding the determinant is not as useful as finding the condition number. Unless you can show t...

Since the matrix is near to singular I have tried pinv to solve but getting too deviated answers. Can you please let me know if there is any other way to solve the equations

Hi VS, Usually in these cases, finding the determinant is not as useful as finding the condition number. Unless you can show t...

1 month ago | 0

| accepted

Answered

How to find the delta cycles (change in the number of cycles) in two sine signals with nearly identical frequencies?

Hi Jay, Using the hilbert transform on the signal gives the so-called analytic signal. The transform creates an imaginary part...

How to find the delta cycles (change in the number of cycles) in two sine signals with nearly identical frequencies?

Hi Jay, Using the hilbert transform on the signal gives the so-called analytic signal. The transform creates an imaginary part...

1 month ago | 0

Answered

Guess update for Bessel function

Hi JS, For j = 1, figure 1 shows the general behavior of the function. Since it's a plot of log(abs(fun)) there can be no zero...

Guess update for Bessel function

Hi JS, For j = 1, figure 1 shows the general behavior of the function. Since it's a plot of log(abs(fun)) there can be no zero...

2 months ago | 0

| accepted

Answered

What kind of sorting algorithm is this?

Hello Karen, I guess it could depend on what is meant by 'corresponds to', since this sort is O(n^2) as is the bubble sort. So...

What kind of sorting algorithm is this?

Hello Karen, I guess it could depend on what is meant by 'corresponds to', since this sort is O(n^2) as is the bubble sort. So...

2 months ago | 1

| accepted

Answered

sum two variables ignoring NaN

Hi alpedhuez, sum(Var1,Var2,2,'omitnan') Here the '2' indicates that you are summing rows rather than columns (the default is ...

sum two variables ignoring NaN

Hi alpedhuez, sum(Var1,Var2,2,'omitnan') Here the '2' indicates that you are summing rows rather than columns (the default is ...

2 months ago | 0

Answered

How to make a color scale like this?

Hi pf, I copied in the image in your posting (a jpeg) and called it image1. That results in a 3d array that is 180x882x3 uint8...

How to make a color scale like this?

Hi pf, I copied in the image in your posting (a jpeg) and called it image1. That results in a 3d array that is 180x882x3 uint8...

2 months ago | 0

Answered

FFT of coupled oscillator from time domain (picosecond) to frequency domain (wavenumber, cm-1)

HI Kunjanni, this occurs because the fft contains both positive and negative frequency (wavenumber) components. Zero frequency...

FFT of coupled oscillator from time domain (picosecond) to frequency domain (wavenumber, cm-1)

HI Kunjanni, this occurs because the fft contains both positive and negative frequency (wavenumber) components. Zero frequency...

2 months ago | 0

| accepted

Answered

Exponential function produces error with integrate command

Hi Ben for large x, exp(x) overflows to infinity before the exp(-(x.^2)/2) factor has a chance to bring in down. Try completin...

Exponential function produces error with integrate command

Hi Ben for large x, exp(x) overflows to infinity before the exp(-(x.^2)/2) factor has a chance to bring in down. Try completin...

2 months ago | 0

| accepted

Answered

Even distribution of coordinates subset

Hi Giovanni, here is one method, which assumes that B is supposed to span basically the entire domain of A. Suppose the set of...

Even distribution of coordinates subset

Hi Giovanni, here is one method, which assumes that B is supposed to span basically the entire domain of A. Suppose the set of...

2 months ago | 0

| accepted

Answered

Why does pcolor require you to transpose the matrix when you give it x and y vectors?

Hi Darcy. pcolor works in the same fashion as meshgrid, and creates a grid with nx = length(x) columns and ny = length(y) rows....

Why does pcolor require you to transpose the matrix when you give it x and y vectors?

Hi Darcy. pcolor works in the same fashion as meshgrid, and creates a grid with nx = length(x) columns and ny = length(y) rows....

2 months ago | 0

Answered

Undefined function 'sinc' for input arguments of type 'double'.

Hi boutaina, Since it's easy in Matlab to make your own functions, you don't need toolboxes that are (probably not intentionall...

Undefined function 'sinc' for input arguments of type 'double'.

Hi boutaina, Since it's easy in Matlab to make your own functions, you don't need toolboxes that are (probably not intentionall...

2 months ago | 0

Answered

Choose Lyapunov function of a linear system

Hi Cristina, The first part of this is probably telling you what you already know, but denoting C = c*[1 -1 0; 0 1 -1; -1 0 1]...

Choose Lyapunov function of a linear system

Hi Cristina, The first part of this is probably telling you what you already know, but denoting C = c*[1 -1 0; 0 1 -1; -1 0 1]...

2 months ago | 0

Answered

orthagonal planes and normal vectors

Hi M, A plane satisfies the dot product equation equation r dot v = const where r is the position vector [x y z] and v = [vx ...

orthagonal planes and normal vectors

Hi M, A plane satisfies the dot product equation equation r dot v = const where r is the position vector [x y z] and v = [vx ...

2 months ago | 2

| accepted

Answered

How to generate a random complex unitary matrix whose columns each sum up to 1

MODIFIED to replace previous random function Hi Michael, here is one way. It's based on the idea that if the unitary matrix U...

How to generate a random complex unitary matrix whose columns each sum up to 1

MODIFIED to replace previous random function Hi Michael, here is one way. It's based on the idea that if the unitary matrix U...

2 months ago | 1

| accepted