Vai all'indice generale

Paragrafo 2.7 risposte ad alcuni degli esercizi proposti

Sottoparagrafo 2.7.1

Risposta all'esercizio Sottoparagrafo 2.6.4. 

sig = @(x) ...
[5*x(2)*x(3) 3*x(2)^2 0;
3*x(2)^2 0 -x(1);
0 -x(1) 0];

n = @(x) ...
[2*x(1)/sqrt((2*x(1))^2+(2*x(2))^2+1);
2*x(2)/sqrt((2*x(1))^2+(2*x(2))^2+1);
1/sqrt((2*x(1))^2+(2*x(2))^2+1)];

t = sig([1/2; sqrt(3)/2; -1])*n([1/2; sqrt(3)/2; -1])
Listato 2.7.1.
Si osserva che, come al solito, le precedenti istruzioni non rappresentano l'unico modo possibile per ottenere il risultato richiesto.

Sottoparagrafo 2.7.2

Risposta all'esercizio Sottoparagrafo 2.6.9. 

x = sym('x', [3 1], 'real');
syms rho g real
sig = [x(2) x(3) 0;
x(3) x(1) 0;
0 0 rho*g*x(3)];
b = [0; 0; -rho*g];

% primo quesito
divergence(sig(1,1:3), x) + b(1)
divergence(sig(2,1:3), x) + b(2)
divergence(sig(3,1:3), x) + b(3)

% secondo quesito
t1 = subs(sig, x(1), 1)*[1; 0; 0]
t2 = subs(sig, x(2), 1)*[0; 1; 0]
t3 = subs(sig, x(3), 1)*[0; 0; 1]
t1_ = subs(sig, x(1), 0)*[-1; 0; 0]
t2_ = subs(sig, x(2), 0)*[0; -1; 0]
t3_ = subs(sig, x(3), 0)*[0; 0; -1]

% terzo quesito
T1 = int(int(t1, x(2), 0, 1), x(3), 0, 1)
T2 = int(int(t2, x(1), 0, 1), x(3), 0, 1)
T3 = int(int(t3, x(1), 0, 1), x(2), 0, 1)
T1_ = int(int(t1_, x(2), 0, 1), x(3), 0, 1)
T2_ = int(int(t2_, x(1), 0, 1), x(3), 0, 1)
T3_ = int(int(t3_, x(1), 0, 1), x(2), 0, 1)
B = int(int(int(b, x(1), 0, 1), x(2), 0, 1), x(3), 0, 1)

R = T1 + T2 + T3 + T1_ + T2_ + T3_ + B
Listato 2.7.2.