Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
6.3. Opt¸iuni 135<br />
Orice solut¸ie a ecuat¸iei pendulului conservă energia: cantitatea y2(t) 2 −cosy1(t)<br />
este constantă. Vom verifica aceasta prin<br />
>> Ec = 0.5*yc(:,2).ˆ2-cos(yc(:,1));<br />
>> max(abs(Ec(1)-Ec))<br />
ans =<br />
0.0263<br />
y 2 (t)<br />
3<br />
2<br />
1<br />
0<br />
−1<br />
−2<br />
6.3. Opt¸iuni<br />
−3<br />
−5 −4 −3 −2 −1 0<br />
y (t)<br />
1<br />
1 2 3 4 5<br />
Figura 6.2: Grafic de fază pentru ecuat¸ia pendulului<br />
Funct¸ia odeset crează o structură de opt¸iuni care poate fi transmisă unui rezolvitor.<br />
Argumentele lui odeset sunt perechi nume proprietate/valoare proprietate.<br />
Sintaxa este<br />
optiuni=odeset(’nume1’, valoare1, ’nume2’, valoare2, ...)<br />
Aceasta crează o structură de opt¸iuni în care proprietăt¸ile cu numele dat primesc o<br />
valoare specificată. Proprietăt¸ile nespecificate primesc valori implicite. Pentru toate<br />
proprietăt¸ile este suficient să dăm doar caracterele de la început care identifică unic<br />
numele proprietăt¸ii. odeset fără argumente afis¸ează toate numele de proprietăt¸i s¸i<br />
valorile lor posibile; valorile implicite apar între acolade:<br />
>> odeset<br />
AbsTol: [ positive scalar or vector {1e-6} ]<br />
RelTol: [ positive scalar {1e-3} ]<br />
NormControl: [ on | {off} ]<br />
OutputFcn: [ function ]