R语言丨求解矩阵的最小奇异值

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")