Bonjour je ne sais pas si je trouverai quelqu'un sur ce site qui puisse m'aider car je souhaite coder sous MATLAB...
Le titre est plutôt clair je souhaite résoudre un systéme de 4 équations non linéaires. J'ai un bon début de code mais j'ai un message d'erreur que je ne comprends pas du tout...
Voici le code :
Et voici le message d'erreur:
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
16
17
18
19
20
21
22
23
24
25 function F=myfun(x) V1=400000 %V2=x(2) %P1=x(1) P2=100000000 %Q1=x(4) Q2=100000000 R=1.575 X=22.275 Z=sqrt(R*R+X*X) teta1=0 %teta2=x(3) d=atan(R/X) dprime=((pi/2)-d) Y=0.00014361 F=[(x(1))-(V1*x(2)/Z)*(sin(teta1-(x(3))-dprime))+(V1*V1/Z)*(sin(dprime)); P2-(V1*x(2)/Z)*(sin(x(3)-teta1-dprime))+(x(2)*x(2)/Z)*(sin(dprime)); x(4)-(V1*x(2)/Z)*(cos(teta1-x(3)-dprime))+(V1*V1)*((cos(dprime)/Z)-Y); Q2-(V1*x(2)/Z)*(cos(x(3)-teta1-dprime))+(x(2)*x(2))*((cos(dprime)/Z)-Y)] x0=[-5;-5;-5;-5] options=optimset('Display','iter'); [x,fval]=fsolve(@myfun,x0,options,V1,P2,Q2,R,X,Z,teta1,d,dprime,Y)
HELP !!??? Input argument "x" is undefined. Error in ==> myfun at 23 F=[(x(1))-(V1*x(2)/Z)*(sin(teta1-(x(3))-dprime))+(V1*V1/Z)*(sin(dprime));
Merci d'avance.
Partager