[MATLAB] itération de model

Fermé
Newenda Messages postés 75 Date d'inscription lundi 1 décembre 2008 Statut Membre Dernière intervention 26 avril 2021 - 13 déc. 2008 à 13:17
Newenda Messages postés 75 Date d'inscription lundi 1 décembre 2008 Statut Membre Dernière intervention 26 avril 2021 - 13 déc. 2008 à 16:04
Bonjour,

Ma question est assez simple

Je veux calculer des modèles par itération : i.e. modele(n+1) = f ( model(n) ) et placé tous ces modèles dans une matrice.


**********************exemple****************

model1 = model0 - 2*model0
model2 = model1 - 2*model1
model3 = model2 - 2*model2
model4 = ....etc...

*********************************************


Merci

1 réponse

Newenda Messages postés 75 Date d'inscription lundi 1 décembre 2008 Statut Membre Dernière intervention 26 avril 2021 23
13 déc. 2008 à 16:04
plus exactement j'aimerais itérer ça

model1= ( mprior' - ( (inv(jacobi(mprior)'*invCd*jacobi(mprior)+invCm))*(jacobi(mprior)'*invCd*(calculA(mprior)'-Aobs)'+invCm*(mprior - mprior)' ) ))'
model2= ( model1' - ( (inv(jacobi(model1)'*invCd*jacobi(model1)+invCm))*(jacobi(model1)'*invCd*(calculA(model1)'-Aobs)'+invCm*(model1 - mprior)' ) ))'
model3= ( model2' - ( (inv(jacobi(model2)'*invCd*jacobi(model2)+invCm))*(jacobi(model2)'*invCd*(calculA(model2)'-Aobs)'+invCm*(model2 - mprior)' ) ))'
model4= ( model3' - ( (inv(jacobi(model3)'*invCd*jacobi(model3)+invCm))*(jacobi(model3)'*invCd*(calculA(model3)'-Aobs)'+invCm*(model3 - mprior)' ) ))'
model5= ( model4' - ( (inv(jacobi(model4)'*invCd*jacobi(model4)+invCm))*(jacobi(model4)'*invCd*(calculA(model4)'-Aobs)'+invCm*(model4 - mprior)' ) ))'
model6= ( model5' - ( (inv(jacobi(model5)'*invCd*jacobi(model5)+invCm))*(jacobi(model5)'*invCd*(calculA(model5)'-Aobs)'+invCm*(model5 - mprior)' ) ))'
model7= ( model6' - ( (inv(jacobi(model6)'*invCd*jacobi(model6)+invCm))*(jacobi(model6)'*invCd*(calculA(model6)'-Aobs)'+invCm*(model6 - mprior)' ) ))'
model8= ( model7' - ( (inv(jacobi(model7)'*invCd*jacobi(model7)+invCm))*(jacobi(model7)'*invCd*(calculA(model7)'-Aobs)'+invCm*(model7 - mprior)' ) ))'
model9= ( model8' - ( (inv(jacobi(model8)'*invCd*jacobi(model8)+invCm))*(jacobi(model8)'*invCd*(calculA(model8)'-Aobs)'+invCm*(model8 - mprior)' ) ))'
0