Iniciar sesiónCrear cuenta

# Matlab - Error calculos cinematicos

Vista:

## Error calculos cinematicos

Publicado por iker (1 intervención) el 25/05/2017 09:48:14
Buenas, el programa mostrado a continuación muestra una serie de errores si alguno es capaz de resolverlas sería de gran utilidad. Gracias de antemano

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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
`clear all`
`close all`
` `
```%aldagaiak
```
```syms x1 x2 x3 x4 x5 x6 %x1=F x2=Na x3=Nb x4=Nc x5=Ne x6=Nd
```
` `
```%parametroak
```
`syms nuk m1 g b h L m2 tita VerlA A w`
`syms sinBeta cosBeta Ro atanB`
`sinBeta = (A*cos(tita)*w)/(sqrt(1+A*A*(cos(tita))*(cos(tita))*w*w));`
`cosBeta = 1/(sqrt(1+A*A*(cos(tita))*(cos(tita))*w*w));`
`Ro = (1+A*A*A*(cos(tita))*(cos(tita))*(cos(tita))*w*w*w)/(A*sin(tita)*w*w);`
`atanB = -((VerlA/cosBeta)*(VerlA/cosBeta)*sinBeta)/(Ro*cosBeta);`
` `
`eqn1 = x1-nuk*x2-nuk*x3-sinBeta*x4-nuk*x4*cosBeta;`
`eqn2 = x3-x2+sinBeta*nuk*x4+cosBeta*x4-m1*g;`
`eqn3 = x3*b-nuk*x3*h-m1*g*(b/2);`
`eqn4 = x5-x6-cosBeta*nuk*x4+sinBeta*x4;`
`eqn5 = nuk*x5+nuk*x6-cosBeta*x4-sinBeta*nuk*x4-m2*g-m2*atanB;`
`eqn6 = -L*sinBeta*x4+L*cosBeta*nuk*x4-x5*(L-(L*(h/2)*cosBeta));`
`[x1, x2, x3, x4, x5, x6] = solve(eqn1,eqn2,eqn3,eqn4,eqn5,eqn6);`
` `
`soluzio=subs([x1 x2 x3 x4 x5 x6],[nuk,m1,g,b,h,L,m2,A,w],[0.57,0.25,9.8,2*pi,0.5,0.2,0.5,20,742.6]);`
` `
`figure`
`ezplot(soluzio(1),[0,2*pi])`
`grid on`
`title('V_b');`
`xlabel('tita (rad)');`
`ylabel('V_b (m/s)');`
` `
`figure`
`ezplot(soluzio(2),[0,2*pi])`
`grid on`
`title('w_2');`
`xlabel('tita (rad)');`
`ylabel('w_2 (rad/s)');`
` `
` `
`figure`
`ezplot(soluzio(3),[0,2*pi])`
`grid on`
`title('a_b');`
`xlabel('tita (rad)');`
`ylabel('a_b (m/s^2)');`
` `
`figure`
`ezplot(soluzio(4),[0,2*pi])`
`grid on`
`title('alfa_2');`
`xlabel('tita (rad)');`
`ylabel('x1 (rad/s^2)');`
` `
`figure`
`ezplot(soluzio(5),[0,2*pi])`
`grid on`
`title('V_b');`
`xlabel('tita (rad)');`
`ylabel('V_b (m/s)');`
` `
`figure`
`ezplot(soluzio(6),[0,2*pi])`
`grid on`
`title('V_b');`
`xlabel('tita (rad)');`
`ylabel('V_b (m/s)');`
Valora esta pregunta
0
Revisar política de publicidad