回答済み
Sorting eigenvectors using symbolic toolbox for PageRank algorithm
There's an unknown variable in the value you pass to sort, so this won't be sorted by magnitude as the magnitude isn't known. H...

4年弱前 | 0

| 採用済み

回答済み
eigenvalue and eigenvector of free vibration while stiffness matrix is not invertable
If eig has been working well for the size of your problem you could consider using [U, D] = eig(K, M); %this solves K*U = M*U*D...

4年弱前 | 1

| 採用済み

回答済み
Finding a node in graph with most mutually adjacent nodes
Thanks for adding the tag, Steve! Here's another idea for how to do this. So you want to find all pairs of edges that go a->b a...

4年弱前 | 1

| 採用済み

回答済み
eigen value what i am getting from matlab software that is not exact when i am compairing with hand calculated value and the mode shape i am getting is reverse
The results look correct to me, can you say what you are expecting instead? Two points that are maybe relevant: MATLAB does num...

4年弱前 | 0

回答済み
The node number in the digraph does not match the node number used in shortestpath().
Try putting the node numbers in quotes: [path,d, edgepath] = shortestpath(G,"228","231"); You have given each node a name, and...

約4年前 | 1

| 採用済み

回答済み
Unitary matrix with non-orthogonal eigenvectors?
The short answer is that it's possible to compute an orthonormal basis of eigenvectors for an orthogonal matrix, but that MATLAB...

約4年前 | 2

| 採用済み

回答済み
Eigenvalue problem on non-Hermitian matrix
You can think of it like this: If there is an invertible V such that inv(V)*H*V=D, then you can expect that the matrix myW = inv...

約4年前 | 0

| 採用済み

回答済み
Seasonal Breakdown from a Time Series data
Perhaps the trenddecomp function would be useful? It has an optional input period where a vector of several periods can also be ...

約4年前 | 1

回答済み
How to find all the downstream nodes from a node in a graph?
Looking at this again due to the recent comment added, it might be simpler to use a directed graph instead of an undirected one:...

約4年前 | 1

回答済み
Which solver should I use to solve a square and symmetric, linear system of equations?
The first thing to check is the size and density of your matrix. The SYMMLQ function you mentioned is one of a group of methods ...

約4年前 | 1

| 採用済み

回答済み
Can we compute the graph Laplacian matrix for a directed graph?
It depends how you want to define it, there is no one consistent definition of what the graph laplacian of a directed graph is. ...

約4年前 | 1

| 採用済み

回答済み
eigs using 'smallestabs' vs scalar
Hi Jack, I had initially misunderstood that you were getting different results when passing in 'smallestabs' vs. passing in the...

約4年前 | 0

| 採用済み

回答済み
missing eigenvector when using eig on defective matrix
The short answer is that this is a case where a small change in the input matrix has a large effect on the eigenvectors being co...

約4年前 | 1

| 採用済み

回答済み
THE GRAPH FUNCTION DOESN'T CONSIDER WEIGHTS
Hard to say what's going wrong without seeing your code. Here's an example where you can see the weights being used to decide on...

約4年前 | 1

| 採用済み

回答済み
mldivide (backslash): advantage by using sparse matrix instead of a full matrix with respect to the quality of the results?
The condition estimate used in mldivide to decide if a warning should be given is much less accurate for the sparse case. This i...

約4年前 | 2

| 採用済み

回答済み
LDL does not support complex symmetric matrices
Largely the reason LDL doesn't support this is that it's less commonly requested for complex symmetric matrices than for complex...

4年以上前 | 0

| 採用済み

回答済み
Inconsistency in behavior of SVD
Firstly, I agree completely with John's great explanations above: Any singular value below eps*first singular value should be tr...

4年以上前 | 3

| 採用済み

回答済み
Accuracy Problem when solving linear equation system using: lu(S) vs. decomposition(S,'lu')
The version in decomposition does some optional steps of iterative refinement: It uses the same solution you have above based on...

4年以上前 | 0

| 採用済み

回答済み
Need help fixing errors on code written in old MATLAB syntax
I'd also recommend replacing the calls svd(A) and svd(A, 0) with calling svd(A, 'econ'), since the code here is only using the f...

4年以上前 | 0

回答済み
Different results on different computers, Matlab 2021b - 64 bit, windows 64 bit, both Intel chips
First, about "which value to trust?" Both values are equally trustworthy, the differences in results come down to applying the ...

4年以上前 | 1

| 採用済み

回答済み
eigs with Extended Capabilities
The first of these would only be useful if you need to apply EIGS to many problems in parallel, in which case each of these coul...

4年以上前 | 1

回答済み
How to parallelize many small SVD's on a GPU
As of R2021b, there is a function pagesvd which is supported on both CPU and distributed arrays. It's not supported for gpuArray...

4年以上前 | 2

| 採用済み

回答済み
How does shortestpath function work?
Every edge has a number, which is the order in which they appear in the Edges table (try displaying g.Edges). The edgepath conta...

4年以上前 | 1

| 採用済み

回答済み
Eigs passes the wrong StartVector
The standard method used in EIGS is only efficient for matrices that are quite large. When eigs detects that the input matrix is...

4年以上前 | 0

| 採用済み

回答済み
How to enhance the performance of Matalb's function sqrtm?
Not in general. You might be able to find a cheaper but less accurate way to compute this, but if that's acceptable would depend...

4年以上前 | 0

回答済み
fastest path between more than two nodes
To compute simply the shortest-path distance, you can use the distances function and pass in a graph object you've constructed u...

4年以上前 | 1

回答済み
Follow a complex eigenvalue in a spectrum at different times
You could try using the matchpairs function: rng default; d = randn(20, 1) + 1i*randn(20, 1); d2 = d + 0.1*(randn(20, 1) + 1i...

4年以上前 | 1

| 採用済み

回答済み
Error using eig in Matlab 2021b by using Parallels 17 of Macbook Pro (M1 )
Since MATLAB is trying to load a .dll file, I'm assuming you're running a Windows VM on this Macbook Pro. MATLAB doesn't support...

4年以上前 | 0

| 採用済み

回答済み
use of svds with function handles
The svds function needs to be able to apply both A*x and A'*x, so your function handle should accept a second input which is eit...

4年以上前 | 2

回答済み
How to compute cholesky to all slice of a tensor?
There isn't another way to do this right now. We have functions that do other linear algebra operations to each page of an ND-ar...

4年以上前 | 1