Discussion :

[neon29200]

Bonsoir à tous,

Dans le cadre d'une modélisation d'un problème physique par une méthode de récurrence, je suis amené à trouver les racines une équation du style f(x)=a*tan(x)+b*cos(x)-c avec a,b,c constantes (en réalité elles ne sont pas constantes mais auront été calculées à l'itération précédente).

Par ailleurs, cette même équation peut se mettre sous la forme:
"x=arrcos(e*x)-d". On cherche donc l'abscisse du/des points d'intersection de ces deux courbes.

Nous avons donc tenté:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
clear all
close all
d=1;
e=1.2;
x=0:0.1:pi/2;
% f=a*cos(x)+b*tan(x)-c peut s'écrire 
 
g=acos(e*x)-d;
hold on
plot(x,g,'b-')
hold on
plot(x,x,'g')
grid on
syms x
x0=solve('acos(e*x)-d = x')

mais la réponse de Matlab est:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
Warning: Imaginary parts of complex X and/or Y arguments ignored 
> In Fonction_Test_f at 15
Warning: Explicit solution could not be found. 
> In solve at 81
  In Fonction_Test_f at 20
 
x0 =
 
[ empty sym ]

Voilà si quelqu'un aurait une méthode pour résoudre de telles équations, je suis preneur

Merci

[Invité]

Bonjour,

• Soit tu utilises les variables symboliques sans guillemets
  

  Code :

  Sélectionner tout - Visualiser dans une fenêtre à part

  1
  2
  3
  4
  d=1;
  e=1.2;
  syms x
  x0=solve(acos(e*x)-d - x)

• Soit tu laisses les guillemets mais utilises la fonction subs:
  

  Code :

  Sélectionner tout - Visualiser dans une fenêtre à part

  1
  2
  3
  d=1;
  e=1.2;
  x0 =solve( subs('acos(e*x)-d = x',{'e','d'},{e,d}) )

Mais pas les deux sinon tu ne rentres pas les valeurs de e et d elles restent des inconnues.