library(ggplot2)
compute_smallest_singular_SVD <- function(n)
{
# Create the upper-triangular matrix
A <- diag(0.1, n, n)
A[upper.tri(A)] <- 1
# Compute the smallest singular value using svd
svd_values <- svd(A)$d #只取奇异值,忽略左奇异和右奇异向量!
smallest_singular_value <- min(svd_values)
return(smallest_singular_value)
}
compute_min_singular_AtA<-function(n){
# Create the upper-triangular matrix
A <- diag(0.1, n, n)
A[upper.tri(A)] <- 1
# Compute the smallest eigenvalue of A^2
eigenvalues_A <- eigen(t(A) %*% A)$values
smallest_eigenvalue_A <- sqrt(min(eigenvalues_A))
return(smallest_eigenvalue_A)
}
n_values <- 1:30
results <- sapply(n_values,compute_smallest_singular_SVD)
results
outcomes<-sapply(n_values,compute_min_singular_AtA)
results_df <- data.frame(Size = n_values, results)
outcomes_df<-data.frame(size=n_values,outcomes)
outcomes_df
results_df[,1]
qplot(x=results_df[,1], y=log10(results), data=results_df, geom="auto")
qplot(x=outcomes_df[,1],y=log10(outcomes),data = outcomes_df,geom="point")
R语言丨求解矩阵的最小奇异值
最新推荐文章于 2024-08-17 08:45:02 发布