Bonjour à tous,
Je recherche les racines d'un polynômes fractionnels de la forme:
sum(lambda./(lambda-mu)) où mu est l'inconnu et lambda un vecteur de taille 100.
J'ai donc pensé à la fonction solve et je l'ai utilisée de la manière suivante:
Elle fonctionne bien mais prend énormément de temps, or je dois trouver les racines de ce polynôme avec plusieurs expériences aléatoires différentes.
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2 syms mu real; mu=double(solve(sum(lambda./(lambda-mu))/M-1/c==0,mu,'Real', true));
Connaissez une méthode similaire mais plus rapide? J'ai pensé à fzeros() ou d'autres dans le même style, mais je ne peux pas savoir le voisinage des racines.
Je vous remercie d'avance!
Partager