Trouble in having access to the cvx variable inside the cvx loop.
5 次查看(过去 30 天)
显示 更早的评论
I write this code:
lambda_opt = randn(num_combinations, 1);
cvx_quiet(true);
cvx_begin variable lambda_opt(num_combinations)
minimize(objective_function(lambda_opt, G, rho, optimal_design))
subject to sum(lambda_opt) == 1
lambda_opt >= 0
cvx_end
function obj_param = objective_function(lambda, G, rho, optimal_design)
result_matrix = zeros(size(G,1));
for ind = 1 : length(lambda)
result_matrix = result_matrix + lambda(ind)*G(:,:,ind);
end
% Convert the CVX matrix to a regular MATLAB matrix.
result_matrix = double(result_matrix);
[~,eig_val] = eig(result_matrix);
obj_param = trace_inv(sqrtm(result_matrix)) * max(eig_val);
I did not have access to the matrix "result_matrix" to calculate its eigenvalues. So, I searched about this and find 3 different answers for tackling this and changing this matrix from cvx object to the regular matrix:
- result_matrix = value(result_matrix);
- result_matrix = cvx_optval.result_matrix;
- result_matrix = double(result_matrix);
However, all of these did not work.
I appreciate it if you have any suggestions for this.
2 个评论
Harsha Vardhan
2023-10-3
编辑:Harsha Vardhan
2023-10-3
Hi, please explain what do you mean by 'I did not have access to the matrix "result_matrix" to calculate its eigenvalues.' What is the error you observed if you use 'results_matrix' directly to calculate the eigen values as follows, after the for loop?
[~,eig_val] = eig(result_matrix);
回答(1 个)
arushi
2023-10-17
Hi Alireza,
I understand that you are not able to access the “result_matrix” for calculating the eigenvalues. I see that you are trying to convert the cvx object to MATLAB matrix and then trying to calculate the eigen values.
You can try calculating the eigen values directly from the cvx object as all the three functions value ,cvx_optval.result_matrix and double does not provide direct access to 'result_matrix'.
Code:
% Remove the line: result_matrix = double(result_matrix);
% Calculate the eigenvalues of result_matrix
[~, eig_val] = eig(result_matrix);
eigenvalues = diag(eig_val);
Please let me know if this works or the error that you are facing in directly using the cvx object to calculate the eigen values.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Matrix Computations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!