Trabalho 2 Matrizes

Prévia do material em texto

Programação para Engenharia 
Trabalho II 
Victor Soares Braz 
 
1) Declare as matrizes A,B e C abaixo: 
A=[ 1 2 3 4 5 6 7 ] B=[ 3 6 9 12 15 18 21 ] 
 C= [0 5 10 15 10 5 0 
 -1 -2 -3 -4 -5 -6 -7] 
Através das matrizes acima, determine a seguir utilizando os comandos já 
mencionados: 
D = [ 4 5 6 7] 
E = [ 7 6 5 4 ] 
F = [ 4 5 6 7 
 7 6 5 4] 
G = [ 4 7 
 5 6 
 6 5 
 7 4] 
H = [12 
 15 
 18 
 21] 
I = [5 
 -6] 
J = [0 5 10 
 -1 -2 -3] 
K = [0 25 100 
 1 4 9] 
L = [0 24 99 
 2 5 10] 
M = [0 -24 -99 
 -2 -5 -10] 
>> A=[1:7] 
B= [3:3:21] 
C= [0 5 10 15 10 5 0;-1 -2 -3 -4 -5 -6 -7] 
A = 1 2 3 4 5 6 7 
B = 3 6 9 12 15 18 21 
C = 0 5 10 15 10 5 0 
 -1 -2 -3 -4 -5 -6 -7 
>> D=A(4:7) 
D = 4 5 6 7 
>> E=sort(D,'descend') 
E =7 6 5 4 
>> F = [D; E] 
F = 4 5 6 7 
 7 6 5 4 
>> G = F' 
G = 4 7 
 5 6 6 5 
 7 4 
>> H = B(4:7) 
H = 12 15 18 21 
>> H = H' 
H =12 
 15 
 18 
 21 
>> C0 = C(1,6) 
C0 = 5 
>> C1 = C(2,6) 
C1 = -6 
>> I = [C0 C1] 
I = 5 -6 
>> I = I’ 
>> I=I' 
I = 5 
 -6 
>> J = C(1:2,1:3) 
J = 0 5 10 
 -1 -2 -3 
>> K = J.*J 
K = 0 25 100 
 1 4 9 
>> Z=[0 -1 -1;1 1 1] 
Z = 0 -1 -1 
 1 1 1 
>> L = K+Z 
L = 0 24 99 
 2 5 10 
>> M = L.*(-1) 
M = 0 -24 -99 
 -2 -5 -10 
2) 
>> A = [2 1.5 1; 1 6 -2; 2 4 0] 
A = 2.0000 1.5000 1.0000 
 1.0000 6.0000 -2.0000 
 2.0000 4.0000 0 
>> B=[13.20; 21.64; 26.62]%vetor de soluções da matriz mat 
B= 13.2000 
 21.6400 
 26.6200 
 
>> X=B\A 
 
X = 
 0.0356 0.1502 0.0565 
3) a) 
>> A=rand(1,50) 
A = Columns 1 through 9 
 0.8147 0.9058 0.1270 0.9134 0.6324 0.0975 0.2785 0.5469 0.9575 
 Columns 10 through 18 
 0.9649 0.1576 0.9706 0.9572 0.4854 0.8003 0.1419 0.4218 0.9157 
 Columns 19 through 27 
 0.7922 0.9595 0.6557 0.0357 0.8491 0.9340 0.6787 0.7577 0.7431 
 Columns 28 through 36 
 0.3922 0.6555 0.1712 0.7060 0.0318 0.2769 0.0462 0.0971 0.8235 
 Columns 37 through 45 
 0.6948 0.3171 0.9502 0.0344 0.4387 0.3816 0.7655 0.7952 0.1869 
 Columns 46 through 50 
 0.4898 0.4456 0.6463 0.7094 0.7547 
b) 
 >> B = A+0.5 
B = Columns 1 through 9 
 1.3147 1.4058 0.6270 1.4134 1.1324 0.5975 0.7785 1.0469 1.4575 
 Columns 10 through 18 
 1.4649 0.6576 1.4706 1.4572 0.9854 1.3003 0.6419 0.9218 1.4157 
 Columns 19 through 27 
 1.2922 1.4595 1.1557 0.5357 1.3491 1.4340 1.1787 1.2577 1.2431 
 Columns 28 through 36 
 0.8922 1.1555 0.6712 1.2060 0.5318 0.7769 0.5462 0.5971 1.3235 
 Columns 37 through 45 
 1.1948 0.8171 1.4502 0.5344 0.9387 0.8816 1.2655 1.2952 0.6869 
 Columns 46 through 50 
 0.9898 0.9456 1.1463 1.2094 1.2547 
>> B=B' 
B = 1.3147 
 1.4058 
 0.6270 
 1.4134 1.1324 
 0.5975 
 0.7785 
 1.0469 
 1.4575 
 1.4649 
 0.6576 
 1.4706 
 1.4572 
 0.9854 
 1.3003 
 0.6419 
 0.9218 
 1.4157 
 1.2922 
 1.4595 
 1.1557 
 0.5357 
 1.3491 
 1.4340 
 1.1787 
 1.2577 
 1.2431 
 0.8922 
 1.1555 
 0.6712 
 1.2060 
 0.5318 
 0.7769 
 0.5462 
 0.5971 
 1.3235 
 1.1948 
 0.8171 
 1.4502 
 0.5344 
 0.9387 
 0.8816 
 1.2655 
 1.2952 
 0.6869 
 0.9898 
 0.9456 
 1.1463 
 1.2094 
 1.2547 
>> B=sort(B,'descend') 
B = 1.4706 
 1.4649 
 1.4595 
 1.4575 
 1.4572 
 1.4502 
 1.4340 
 1.4157 
 1.4134 
 1.4058 
 1.3491 
 1.3235 
 1.3147 
 1.3003 
 1.2952 
 1.2922 
 1.2655 
 1.2577 
 1.2547 
 1.2431 
 1.2094 
 1.2060 
 1.1948 
 1.1787 
 1.1557 
 1.1555 
 1.1463 
 1.1324 
 1.0469 
 0.9898 
 0.9854 
 0.9456 
 0.9387 
 0.9218 
 0.8922 
 0.8816 
 0.8171 
 0.7785 
 0.7769 
 0.6869 
 0.6712 
 0.6576 
 0.6419 
 0.6270 
 0.5975 
 0.5971 
 0.5462 
 0.5357 
 0.5344 
 0.5318 
4) 
>> A=[6 -1 2] 
A = 6 -1 2 
>> B=[-2 3 -4] 
B = -2 3 -4 
>> C=[-3 1 5] 
C = -3 1 5 
>> AB=B-A 
AB = -8 4 -6 
>> AC=C-A 
AC = -9 2 3 
>> prodoutodvetorial=cross(AB, AC) 
prodoutodvetorial = 24 78 20 
>> norm(AB) 
ans = 10.7703 
>> norm(ac) 
ans = 9.2736 
>> produtoescalar = dot(x,y) 
produtoescalar = 104.4222 
>> angulobcanoverticea=acos(produtoescalar/x*y) 
angulobcanoverticea = 0 + 5.2364i 
5) a) 
>> Y=min(X) 
Y = 3 2 1 
b) 
 >> y=sum(Y) 
y = 6 
c) 
 >> Z=max(X) 
Z = 65 32 45 
d) 
 >> z=prod(Z) 
z = 93600 
e) 
 >> x1=[6 2 45] 
X1=norm(x1) 
x1 = 6 2 45 
X1 = 45.4423 
>> a=norm(X) 
a = 74.7098 
>> x2=[65 32 9] 
X2=norm(x2) 
x2 = 65 32 9 
X2 = 73.0068 
>> x3=[3 8 1] 
X3=norm(x3) 
x3 = 3 8 1 
X3 = 8.6023>> t=norm(X) 
t = 74.7098 
6) 
>> delta = [17/6 -1/3; 1/6 -1.2/2] 
delta = 2.8333 -0.3333 
 0.1667 -0.6000 
>> delta1 = [16 -1/3; 1/10 -1.2/2] 
delta1 = 16.0000 -0.3333 
 0.1000 -0.6000 
>> delta2 = [17/6 16; 1/6 1/10] 
delta2 = 2.8333 16.0000 
 0.1667 0.1000 
>> delta = det (delta) 
delta = -1.6444 
>> delta1 = det(delta1) 
delta1 = -9.5667 
>> delta2 = det(delta2) 
delta2 = -2.3833 
>> v1=delta1/delta 
v1 = 5.8176 
>> v2=delta2/delta 
v2 = 1.4493