prova simulação de sistemas

Prévia do material em texto

1 – 
- V1 – R1*I1 – R3(I1 – I2) – R2(I1 – I3) = 0
V2 – R6(I2 – I4) – R4(I2 – I3) – R3(I1 – I2) = 0
- R4(I2 – I3) – R5(I3 – I4) – R2(I1 – I3) = 0 
- V3 – R8*14 – R7(I4 – I5) – R5(I4 – I3) – R6(14 – 12) = 0
- V1 – R7(I4 – I5) – R9*I5 = 0 
>> V1 = 50; V2 = 30; V3 = 40;
>> R1 = 20; R2 = 18; R3 = 12; R4 = 9;
>> R5 = 5; R6= 15; R7 = 8; R8 = 13; R9 = 3;
 
>> A = [-(R1+R2+R3) R3 R2 0 0
R3 -(R3+R4+R6) R4 R6 0
R2 R4 -(R2+R4+R5) R5 0
0 R6 R5 -(R8+R7+R5+R6) R7
0 0 0 R7 -(R7+R9)]
A =
 -50 12 18 0 0
 12 -36 9 15 0
 18 9 -32 5 0
 0 15 5 -41 8
 0 0 0 8 -11
>> B = [-V1; V2; 0; -V3; -V1]
B =
 -50
 30
 0
 -40
 -50
>> I= A\B
I =
 2.3017 I1
 1.9021 I2
 2.3477 I3
 3.3152 I4
 6.9565 I5 
Corrente no R1 = 2.3A; Corrente no R2 (I3 – I1) = 0.04A; Corrente no R3 (I1 – I2) = 0.4A; Corrente no R4(I3 – I2) = 0.44A; Corrente no R5(I4 – I3) = 0.97A; Corrente no R6(I4 – I2) = 1.41A; Corrente no R7(I5 – I4) = 3.64A; Corrente no R8(I4) = 3.31A; Corrente no R9(I5) = 6.95ª.
2 – >> V = 250
V =
 250
>> h = [5 6 7 8 9]
h =
 5 6 7 8 9
>> r = sqrt((3*V)/(pi*5))
r =
 6.9099
>> S = pi*6.9099*sqrt(6.9099.^(2)+5.^(2))
S =
 185.1520
>> r = sqrt((3*V)/(pi*6))
r =
 6.3078
>> S = pi*6.3078*sqrt(6.3078.^(2)+6.^(2))
S =
 172.5158
>> r = sqrt((3*V)/(pi*7))
r =
 5.8399
>> S = pi*5.8399*sqrt(5.8399.^(2)+7.^(2))
S =
 167.2505
>> r = sqrt((3*V)/(pi*8))
r =
 5.4627
>> S = pi*5.4627*sqrt(5.4627.^(2)+8.^(2))
S =
 166.2470
>> r = sqrt((3*V)/(pi*9))
r =
 5.1503
>> S = pi*5.1503*sqrt(5.1503.^(2)+9.^(2))
S =
 167.7793
3 – >> g = 9.8
g = 9.8000
>> t = [1; 2; 3; 4; 5; 6; 7; 8; 9; 10]
t =
 1
 2
 3
 4
 5
 6
 7
 8
 9
 10
>> d = (1/2)*g*t.^(2)
d =
 4.9000
 19.6000
 44.1000
 78.4000
 122.5000
 176.4000
 240.1000
 313.6000
 396.9000
 490.0000
4 – >> x = [2 4 6 8 10]
x =
 2 4 6 8 10
>> y = [3 6 9 12 15]
y =
 3 6 9 12 15
>> z = (x(1)*y(1)+y(1)/x(1))/(x(1)+y(1)).^(y(1)-x(1)) +12.^(x/y)
z =
 6.7415
>>z = (x(2)*y(2)+y(2)/x(2))/(x(2)+y(2)).^(y(2)-x(2)) +12.^(x(2)/y(2))
z =
 5.4965
>> z = (x(3)*y(3)+y(3)/x(3))/(x(3)+y(3)).^(y(3)-x(3)) +12.^(x(3)/y(3))
z =
 5.2579
>> z = (x(4)*y(4)+y(4)/x(4))/(x(4)+y(4)).^(y(4)-x(4)) +12.^(x(4)/y(4))
z =
 5.2421
>> z = (x(5)*y(5)+y(5)/x(5))/(x(5)+y(5)).^(y(5)-x(5)) +12.^(x(5)/y(5))
z =
 5.2415
5 – >> fplot('(x.^2-5*x+10)/(x.^2-2*x-3)',[-10,10])
6 –
>> x = [-2*pi:0.1:2*pi];
>> y = (3*x)*sin(x) – 2*x;
>> yd = 3*x*cos(x) + 3*sin(x) – 2;
>> plot(x,y,’-‘,x,yd’—‘)