Ayuda calculo cinemático
Publicado por unai (1 intervención) el 24/05/2017 21:02:55
Hola buenas, tengo un problema con el programa planteado para realizar un calculo cinemático. Si a alguien se le ocurre como solucionarlo seria de gran ayuda. Gracias!
Programa:
Programa:
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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
%Aldagaiak
syms x1 x2 x3 x4 % x1=VerlA2 x2=VerlB2 x3=AerlA2 x4=AerlB2
%Parametroak
syms r W Ro sinTita tanBeta cosBeta sinBeta cosTita
Beta = cos(W*1.5*cosTita);
Ro =(W^2*1.5*sinTita)/(1+(W*1.5*cosTita)^2)^(3/2);
sinTita=(1-(cosTita)^2)^(1/2);
eqn1 = (W*r*tanBeta)-x1;
eqn2 = ((W*r)/cosBeta)-x2;
eqn3 = x3+((x2^2)/Ro*tanBeta);
eqn4 = x4+((x2^2)/Ro*((cosBeta^2-sinBeta^2)/cosBeta));
[x1, x2, x3, x4] = solve(eqn1, eqn2, eqn3, eqn4);
solution=subs([x1 x2 x3 x4], [W,r,], [7.5,14000,]);
figure
ezplot(solution (1), [0,2*pi]);
grid on
title('VerlA2');
xlabel('Beta(rad)');
ylabel('VerlA2(m/s)');
figure
ezplot(solution(2), [0,2*pi]);
grid on
title('VerlB2');
xlabel('Beta (rad)');
ylabel('VerlB2(m/s)');
figure
ezplot(solution(3), [0,2*pi]);
grid on
title('AerlA2');
xlabel('Beta (rad)');
ylabel('AerlA2 (m/s^2)');
figure
ezplot(solution(4), [0,2*pi]);
grid on
title('AerlB2');
xlabel('Beta (rad)');
ylabel('AerlB2 (rad/s^2)');
Valora esta pregunta


0