Effacer les filtres
Effacer les filtres

How to write a Matlab code for a sequence of functionals and their minimums

1 vue (au cours des 30 derniers jours)
Let A(α) be a matrix (N ,N) depending on a real parameter α ( α ∈ I_α ) and (f_k )_(k≥1) a sequence of vectors in R^N.
We note by u a vector in R^N. We define the following functional sequence
J_k (α)=‖A(α) f_k - u‖^2 (norme square)
and the sequence of minimums
α_(k+1)=min (J_k) for α ∈ I_α
The vector f_1 and the number α_1 are given. We suppose that k varies from 1 to L with L≥1
Question∶How to write the Matlab code corresponding to the J_k and α_(k+1) sequences ?
  6 commentaires
Torsten le 3 Avr 2024
Modifié(e) : Torsten le 3 Avr 2024
And on what does this recurrence relation depend ? On f_(k-1),...,f_(k-m) ? On alpha_k ? On the optimal alpha-values of the preceding problems solved alpha_opt_(k-1),...,alpha_opt_(k-m) ?
Djamel le 3 Avr 2024
In this relation, f_k+1 depends on f_k and the matrix A(α_k). By giving f_1 and α_1, the vectors f_k are all well defined.

Connectez-vous pour commenter.

Réponse acceptée

Torsten le 3 Avr 2024
Modifié(e) : Torsten le 3 Avr 2024
Your code:
number_of_iterations = 3;
A = @(alpha) ...; % Your NxN matrix A depending on alpha
u = ...; % Your Nx1 vector u
f(:,1) = ...; % Your initial Nx1 vector f
alpha0 = ...; % Your initial guess for alpha_opt(1)
for i = 1:number_of_iterations
fun = @(alpha) (A(alpha)*f(:,i)-u).'*(A(alpha)*f(:,i)-u);
[alpha_opt(i),~,exitflag,output] = fminunc(fun,alpha0);
f(:,i+1) = ...; function of A(alpha_opt(i)) and f(:,i) - your recurrence relation
alpha0 = alpha_opt(i);
  5 commentaires
Torsten le 4 Avr 2024
Seems fminbnd does not need an initial guess. You only need to specify start and end point of I (assuming I is finite).
Djamel le 6 Avr 2024
I used fminbnd and it worked. Thank you so much.

Connectez-vous pour commenter.

Plus de réponses (0)


En savoir plus sur Particle & Nuclear Physics dans Help Center et File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by