Eng. de Computação

[Equações Diferenciais Ordinárias] Lista 3: equações diferenciais separáveis e homogêneas

2012-09-12T20:16:00.002-07:00

1. Resolva as equações diferenciais dadas ou os problemas de valor inicial.
Bem, a maior parte estava integrando tudo do mesmo lado, mas uma coisa que eu gosto particularmente de fazer é isolar x e y um em cada canto ANTES dessa integração, porque aí rola menos trabalho. Costuma mudar apenas na hora do C, mas eu costumo ter um jeito de "igualar" ao do professor, então não se preocupem porque o resultado bate. Veja só:
a) xdx + ydy = 0
$$ydy = -xdx \\
\int y dy = - \int x dx \\
\dfrac{1}{2} y^2 = - \dfrac{1}{2} x^2 + C \\
y^2 = -x^2 + 2C \\
y = \pm \sqrt{2C - x^2} \\
K = 2C \\
y = \pm \sqrt{K - x^2}$$
b) xdx - y³dy = 0
$$y^3 dy = xdx \\
\int y^3 dy = \int xdx \\
\dfrac{1}{4} y^4 = \dfrac{1}{2} x^2 + C \\
y^4 = 2x^2 + 4C \\
y = \pm \sqrt{4C + 2x^2} \\
K = 4C \\
y = \pm \sqrt{K + 2x^2}$$
$$c) dx + \dfrac{1}{y^4} dy = 0 \\
\dfrac{1}{y^4} dy = -dx \\
y^{-4} dy = -dx \\
\int y^{-4} dy = - \int dx \\
- \dfrac{1}{3} y^{-3} = -x + C \\
\dfrac{1}{3} y^{-3} = x - C \\
y^{-3} = 3x - 3C \\
\dfrac{1}{y^3} = 3x - 3C \\
y^3 = \dfrac{1}{3x - 3C} \\
y = \dfrac{\sqrt[3]{1}}{\sqrt[3]{3x - 3C}} = \dfrac{1}{\sqrt[3]{3x - 3C}} \\
K = -3C \\
y = \dfrac{1}{\sqrt[3]{K + 3x}}$$
$$d) (t + 1)dt - \dfrac{1}{y^2} dy = 0 \\
\dfrac{1}{y^2} dy = (t + 1)dt \\
y^{-2} dy = (t + 1)dt \\
\int y^{-2} dy = \int (t + 1)dt \\
-y^{-1} = \dfrac{1}{2} t^2 + t + C \\
- \dfrac{1}{y} = \dfrac{1}{2} t^2 + t + C \\
\dfrac{1}{y} = - \dfrac{1}{2} t^2 - t - C = - (\dfrac{1}{2} t^2 + t + C) \\
y = -(\dfrac{1}{2} t^2 + t + C)^{-1}$$
$$e) \dfrac{4}{t} dt - \dfrac{y - 3}{y} dy = 0 \\
\dfrac{y - 3}{y} dy = 4 (\dfrac{1}{t}) dt \\
1dy - 3 (\dfrac{1}{y}) dy = 4 (\dfrac{1}{t}) \\
\int dy - 3 \int \dfrac{1}{y} dy = 4 \int \dfrac{1}{t} \\
y - 3 \ln{|y|} = 4 \ln{|t|} + C \\
\ln{|y^3|} + \ln{|t^4|} = y - C \\
\ln{|y^3 * t^4|} = y - C \\
y^3 t^4 = e^{y - C} = e^y * e^{-C} \\
K = e^{-C} \\
y^3 t^4 = Ke^y$$
$$f) y' = \dfrac{xe^x}{2y} \\
\dfrac{dy}{dx} = \dfrac{xe^x}{2y} \\
2y dy = xe^x dx \\
\int 2y dy = \int xe^x dx \\
2 \dfrac{1}{2} y^2 = xe^x - \int e^x * 1 dx [1] \\
y^2 = xe^x - (e^x + C) = xe^x - e^x - C \\
y = \pm \sqrt{xe^x - e^x - C}$$
[1] integral por partes: u = x, dv = e elevado a x
$$g) y' = \dfrac{y}{x^2} \\
\dfrac{dy}{dx} = \dfrac{y}{x^2} \\
\dfrac{1}{y} dy = \dfrac{1}{x^2} dx \\
\int \dfrac{1}{y} dy = \int x^{-2} dx \\
\ln{|y|} = - x^{-1} + C = - \dfrac{1}{x} + C \\
y = e^{1/x + C} = e^{1/x} + e^C \\
K = e^C \\
y = Ke^{1/x}$$
h) senx dx + y dy = 0; y(0) = -2
Ok, probleminha de valor inicial. Muda pouca coisa. Primeiro resolvamos:

$$y dy = - \sin{x} dx \\
\int y dy = - \int \sin{x} dx \\
\dfrac{1}{2} y^2 = - (- \cos{x}) + C = \cos{x} + C \\
y^2 = 2 \cos{x} + 2C \\
y = \pm \sqrt{2C + 2 \cos{x}} \\
K = 2C \\
y(0) = -2 \\
y = \pm \sqrt{K + 2 \cos{0}} = -2 \\
y = \pm \sqrt{K + 2} = -2 \\
y^2 = K + 2 = (-2)^2 = 4 \\
K + 2 = 4 \\
K = 4 - 2 = 2 \\
y = - \sqrt{2 + 2 \cos{x}}$$
$$i) xe^{x^2} dx + (y^5 - 1)dy = 0; y(0) = 0 \\
(y^5 - 1)dy = -xe^{x^2} dx \\
\int (y^5 - 1)dy = - \int xe^{x^2} dx \\
\dfrac{1}{6} y^6 - y = - \int xe^u \dfrac{du}{2x} = - \dfrac{1}{2} \int e^u du = - \dfrac{1}{2} e^u = - \dfrac{1}{2} e^{x^2} + C \\
\dfrac{1}{2} e^{x^2} + \dfrac{1}{6} y^6 - y = C \\
y(0) = 0 \\
\dfrac{1}{2} e^0 + \dfrac{1}{6} * 0 - 0 = C \\
C = \dfrac{1}{2} \\
\dfrac{1}{2} e^{x^2} + \dfrac{1}{6} y^6 - y = \dfrac{1}{2}$$
$$j) (1 + x^2)dx + \dfrac{1}{y} dy = 0; y(-1) = 1 \\
\dfrac{1}{y} dy = -(1 + x^2)dx \\
\int \dfrac{1}{y} dy = - \int (1 + x^2)dx \\
\ln{|y|} = - (x + x^3 /3 + C) \\
y = e^{-(x + x^3 /3 + C)} \\
y(-1) = 1 \\
y = e^{-(-1 + (-1)^3 /3 + C)} = 1 \\
e^{-(-1 - 1/3 + C)} = 1 \\
e^{-(- 4/3 + C)} = 1 \\
e^{4/3 - C} = 1 \\
\dfrac{4}{3} - C = \ln{1} = 0 \\
C = \dfrac{4}{3} \\
y = e^{-(x + x^3 /3 + 4/3)} = e^{- (x^3 + 3x + 4)/3}$$

2. Determine se a equação diferencial apresentada é homogênea. Se for, resolva-a.
$$a) y' = \dfrac{y - x}{x}$$

O esquema aqui é o seguinte: vamos lembrar daquele método f(x, y) = f(tx, ty) pra ver se a equação é homogênea. O que acontece é que vamos trocar aonde tem x por tx, e aonde tem y por ty, e terá de ser a mesma coisa. Veja:
$$y' = \dfrac{ty - tx}{tx} = \dfrac{t(y - x)}{tx} = \dfrac{y - x}{x}$$
Os resultados batem, logo ela é homogênea. Sendo homogênea, resolveremos do método que o professor passou na aula, e com a prática entenderemos bem ele: usaremos um V, tal que V = y/x. E por definição, y' = V + xV'. Ok, ótimo, não deu pra entender nada? Vejamos na prática:
$$v = \dfrac{y}{x} \\
y = vx [1] \\
y' = v + xv' \\
y' = \dfrac{y - x}{x} \\
v + xv' = \dfrac{y - x}{x}$$
Lembremos de [1]:
$$v + xv' = \dfrac{vx - x}{x} = v - 1 \\
x \dfrac{dv}{dx} = v - 1 - v = -1 \\
\dfrac{dv}{dx} = - \dfrac{1}{x} \\
dv = - \dfrac{1}{x} dx \\
\int dv = - \int \dfrac{1}{x} dx \\
v = - \ln{|x|} + C$$
Jogando isso na [1]:
$$y = (-\ln{|x|} + C)x = x(C - \ln{|x|}) \\
\ln{|K|} = C \\
y = x(\ln{|K|} - \ln{|x|}) = x (\ln{|K/x|}$$
Honestamente, eu acho essa última parte mais que desnecessária: eu acho que ela atrapalha. Mas só pra chegar na resposta do livro, tá aí.
E sim, equação diferencial homogênea é bem mais chato mesmo. Essa é a mais fácil, inclusive.
$$b) y' = \dfrac{2y + x}{x}$$
Pra descobrir se é homogênea:
$$y' = \dfrac{2(ty) + tx}{tx} = \dfrac{t(2y + x)}{tx} = \dfrac{2y + x}{x}$$
Logo, é homogênea. Logo, bora usar o v = y/x.
$$y = vx \\
y' = v + xv' \\
y' = \dfrac{2y + x}{x} \\
v + xv' = \dfrac{2y + x}{x} = \dfrac{2vx + x}{x} = 2v + 1 \\
x \dfrac{dv}{dx} = 2v + 1 - v = v + 1 \\
\dfrac{dv}{dx} = \dfrac{v + 1}{x} \\
\dfrac{1}{v + 1} dv = \dfrac{1}{x} dx \\
\int \dfrac{1}{v + 1} dv = \int \dfrac{1}{x} dx \\
\ln{|v + 1|} = \ln{|x|} + C \\
C = \ln{|K|} \\
\ln{|v + 1|} = \ln{|x|} + \ln{|K|} = \ln{|Kx|} \\
v + 1 = Kx \\
v = Kx - 1 \\
y = x(Kx - 1)$$
$$c) y' = \dfrac{2y^2 + x^2}{xy}$$
Vendo se é homogênea:
$$y' = \dfrac{2(ty)^2 + (tx)^2}{(tx)(ty)} = \dfrac{t^2 (2y^2 + x^2)}{t^2 (xy)} = \dfrac{2y^2 +

x^2}{xy}$$
Logo, é homogênea. Então vamos resolver:
$$y' = v + xv' = \dfrac{2y^2 + x^2}{xy} \\
y = vx \\
v + xv' = \dfrac{2(vx)^2 + x^2}{x(vx)} = \dfrac{2v^2 x^2 + x^2}{x^2 v} = \dfrac{x^2 (2v^2 +

1)}{x^2 v} = \dfrac{2v^2 + 1}{v} \\
v + xv' = 2v + \dfrac{1}{v} \\
x \dfrac{dv}{dx} = 2v + \dfrac{1}{v} - v = v + \dfrac{1}{v} \\
\dfrac{dv}{dx} = \dfrac{v + 1/v}{x} \\
\dfrac{1}{v + 1/v} dv = \dfrac{1}{x} dx$$
Pra ficar mais fácil, vamos fatorar o primeiro termo: multiplicar em cima e embaixo por v.
$$\dfrac{v}{v^2 + 1} dv = \dfrac{1}{x} dx \\
\int \dfrac{v}{v^2 + 1} dv = \int \dfrac{1}{x} dx$$
Então. A integral com dx é facílima, mas a com dv é um pouquinho mais complicada. Nada realmente difícil: vamos usar regra da substituição, e veremos que vai dar certo. Chamamos 2v² + 1 de u, assim dv = du/(v² + 1)' = du/2v. Substiuímos o u e o du:
$$\int \dfrac{v}{u} \dfrac{du}{2v} = \int \dfrac{1}{x} dx \\
\dfrac{1}{2} \int \dfrac{1}{u} du = \int \dfrac{1}{x} dx \\
\dfrac{1}{2} \ln{|u|} = \ln{|x|} + C \\
\ln{|u^{1/2}|} = \ln{|x|} + C$$
Fazemos o mesmo esquema geral pra unir a constante ao X, mas dessa vez vamos esperar um pouco

pra usar o K. Vocês vão ver porquê. Vamos trocar por A:
$$C = \ln{|A|} \\
\ln{|u^{1/2}|} = \ln{|x|} + \ln{|A|} = \ln{|Ax|} \\
u^{1/2} = Ax \\
u = A^2 x^2$$
Agora sim, vamos chamar A² de K. E lembrar que u = v² + 1:
$$u = Kx^2 \\
v^2 + 1 = Kx^2 \\
v^2 = Kx^2 - 1 \\
v = \sqrt{Kx^2 - 1}$$
Por fim, v = y/x. Logo:
$$\dfrac{y}{x} = \sqrt{Kx^2 - 1} \\
y = \pm x \sqrt{Kx^2 - 1}$$
Que é o exato mesmo resultado da lista.

$$d) y' = \dfrac{2x + y^2}{xy}$$
Verificaremos se é homogênea:
$$y' = \dfrac{2(tx) + (ty)^2}{(tx)(ty)} = \dfrac{2(tx) + t^2 y^2}{t^2 (xy)} \\
\dfrac{t (2x + ty^2)}{t^2 (xy)} = \dfrac{2x + ty^2}{txy} \neq \dfrac{2x + y^2}{xy}$$
Logo, não é homogênea e não nos convém resolvê-la.

$$e) y' = \dfrac{y^2 + x^2}{2xy}$$
Verificaremos se é homogênea:
$$y' = \dfrac{(ty)^2 + (tx)^2}{2(tx)(ty)} = \dfrac{t^2 (y^2 + x^2)}{2t^2 (xy)} = \dfrac{y^2 + x^2}{2xy}$$
Logo, é homogênea. Resolvendo:
$$y' = v + xv' = \dfrac{y^2 + x^2}{2xy} = \dfrac{(vx)^2 + x^2}{2x(vx)} = \dfrac{v^2 x^2 + x^2}{2vx^2} \\
v + xv' = \dfrac{x^2 (v^2 + 1)}{2vx^2} = \dfrac{v^2 + 1}{2v} = \dfrac{v}{2} + \dfrac{1}{2v} \\
x \dfrac{dv}{dx} = \dfrac{v}{2} - v + \dfrac{1}{2v} = \dfrac{1}{2v} - \dfrac{v}{2} \\
x \dfrac{dv}{dx} = 1/2v - 0.5v \\
\dfrac{dv}{dx} = \dfrac{1/2v - 0.5v}{x} \\
\dfrac{1}{1/2v - 0.5v} dv = \dfrac{1}{x} dx$$
Vamos fatorar de maneira igual a c. E fazer a mesma substituição:
$$\dfrac{v}{0.5 - 0.5v^2} dv = \dfrac{1}{x} dx \\
\int \dfrac{v}{0.5 - 0.5v^2} dv = \int \dfrac{1}{x} dx \\
\int \dfrac{v}{u} \dfrac{du}{-1v} = \int \dfrac{1}{x} dx \\
- \int \dfrac{1}{u} du = \int \dfrac{1}{x} dx \\
- \ln{|u|} = \ln{|x|} + C \\
C = \ln{|A|} \\
- \ln{|u|} = \ln{|x|} + \ln{|A|} = \ln{|Ax|}$$
Podemos considerar -ln(u) como ln(1) - ln(u), o que fica:
$$\ln{| \dfrac{1}{u} |} = \ln{|Ax|} \\
\dfrac{1}{u} = Ax \\
\dfrac{1}{0.5 - 0.5v^2} = Ax \\
Ax (0.5 - 0.5v^2) = 1 \\
0.5Ax - 0.5Axv^2 = 1 \\
-0.5Axv^2 = 1 - 0.5Ax \\
Axv^2 = -2 + Ax \\
v^2 = \dfrac{-2}{Ax} + 1 \\
\dfrac{y^2}{x^2} = \dfrac{-2}{Ax} + 1 \\
y^2 = \dfrac{-2x^2}{Ax} + x^2 \\
y^2 = x^2 \\
y = \sqrt{x^2 - \dfrac{2x}{A}} \\
K = \dfrac{2}{A} \\
y = \pm \sqrt{x^2 - Kx}$$
Honestamente, é o que sempre digo: essas manipulações doidas de constante e tudo mais acho extremamente desnecessárias. Foi só pra chegar na resposta final de alguma forma mesmo.

$$f) y' = \dfrac{2xy}{y^2 - x^2}$$
Vendo se é homogênea:
$$y' = \dfrac{2(tx)(ty)}{(ty)^2 - (tx)^2} = \dfrac{2t^2 (xy)}{t^2 (y^2 - x^2)} = \dfrac{2xy}{y^2 - x^2}$$
Logo, é homogênea. Vamos resolver:
$$y' = v + xv' = \dfrac{2xy}{y^2 - x^2} = \dfrac{2x(vx)}{(vx)^2 - x^2} = \dfrac{2x^2 v}{v^2 x^2 - x^2} \\
v + xv' = \dfrac{2x^2 v}{x^2 (v^2 - 1)} = \dfrac{v}{v^2 - 1} \\
x \dfrac{dv}{dx} = \dfrac{v}{v^2 - 1} - v \\
\dfrac{dv}{dx} = ( \dfrac{v}{v^2 - 1} - v)x^{-1} \\
\dfrac{1}{v/(v^2 - 1) - v} dv = \dfrac{1}{x} dx$$
Fatorando a parte debaixo da fração ali (multiplicar tudo por (v² - 1) pra poder cortar aquela fração enxerida):
$$\dfrac{1}{v - v(v^2 - 1)} dv = \dfrac{1}{x} dx \\
\dfrac{1}{v - v^3 + v} dv = \dfrac{1}{x} dx \\
\dfrac{1}{2v - v^3} dv = \dfrac{1}{x} dx \\
\int \dfrac{1}{2v - v^3} dv = \int \dfrac{1}{x} dx$$
Então né, filhão. Pra mim tá bom por hoje, mas já estou postando porque vai ajudar muita gente que vai estudar essa madrugada. Essa não está resolvida e pelo que entendi precisa utilizar de artifícios legais como frações parciais pra resolver, mas quem tiver uma solução melhor (mais fácil, 1000x mais fácil, fração parcial é o inferno) já sou grato desde já.
No mais, boa noite e bons estudos.

[Equações Diferenciais Ordinárias] Lista 2: forma padrão e forma diferencial

2012-09-02T09:57:00.000-07:00

1. Escreva as equações diferenciais na forma padrão. (dy/dx = f(x, y))
a) xy' + y² = 0
$$xy' = -y^2 \\ y' = \dfrac{-y^2}{x} \\ \dfrac{dy}{dx} = \dfrac{-y^2}{x}$$
b) $$e^x y' - x = y'$$
Esse exercício pegou muita gente, mas ele é bem fácil. Veja só, passe y' do lado direito dividindo tudo do lado esquerdo:
$$e^x - \dfrac{x}{y'} = 1$$
Agora tem só um y', só resolver:
$$- \dfrac{x}{y'} = 1 - e^x \\
-x = y'(1 - e^x) \\
y' = \dfrac{-x}{(1 - e^x)} \\
\dfrac{dy}{dx} = \dfrac{x}{(e^x - 1)}$$
c) $$\dfrac{(x+y)}{(x-y)} dx - dy = 0$$
Parece confuso, mas é só juntar dy e dx, de modo que dy fique em cima e dx fique embaixo. Veja só:
$$\dfrac{(x+y)}{(x-y)} dx = dy \\
\dfrac{dy}{dx} = \dfrac{(x+y)}{(x-y)}$$
d) (x-y)dx + y²dy = 0
$$(x-y)dx = -y^2 dy \\
(x-y) = -y^2 \dfrac{dy}{dx} \\
\dfrac{(x-y)}{-y^2} = \dfrac{dy}{dx} \\
\dfrac{dy}{dx} = - \dfrac{(x-y)}{y^2} \\
\dfrac{dy}{dx} = \dfrac{(y-x)}{y^2}$$
e) $$(e^{2x} - y)dx + e^x dy = 0 \\
e^x dy = (y - e^{2x})dx \\
e^x \dfrac{dy}{dx} = (y - e^{2x}) \\
\dfrac{dy}{dx} = \dfrac{(y - e^{2x})}{e^x} \\
\dfrac{dy}{dx} = \dfrac{y}{e^x} - e^x$$

2. Escreva as equações diferenciais na forma diferencial. (f(x, y)dx - f(x, y)dy = 0)
a) y' = xy
$$\dfrac{dy}{dx} = xy \\
dy = (xy)dx \\
(xy)dx - dy = 0$$
b) y' = xy + 1
$$\dfrac{dy}{dx} = xy + 1 \\
dy = (xy + 1)dx \\
(xy+1)dx - dy = 0$$
c) $$y' = \dfrac{x^2}{y^2} \\
\dfrac{dy}{dx} = \dfrac{x^2}{y^2} \\
y^2 \dfrac{dy}{dx} = x^2 \\
y^2 dy = x^2 dx \\
x^2 dx - y^2 dy = 0$$
d) $$y' = \dfrac{-2y}{x} \\
\dfrac{dy}{dx} = \dfrac{-2y}{x} \\
x \dfrac{dy}{dx} = -2y \\
x dy = -2y dx \\
-2y dx - x dy = 0$$
e) $$y' = \dfrac{xy^2}{x^2 y + y^3} \\
\dfrac{dy}{dx} = \dfrac{xy^2}{x^2 y + y^3} \\
(x^2 y + y^3)dy = (xy^2)dx \\
(xy^2)dx - (x^2 y + y^3)dy = 0$$
f) y' = x³y + xy³
$$\dfrac{dy}{dx} = x^3 y + xy^3 \\
dy = (x^3 y + xy^3)dx \\
(x^3 y + xy^3)dx - dy = 0$$

E está resolvida a lista pra quem quiser estudar. Esqueci de postar logo após a entrega, mas de qualquer forma, aí está.

[Cálculo Numérico] Aula 10: sistema hexadecimal

2012-09-01T14:45:00.001-07:00

Sistema Hexadecimal

Composto de 16 valores (0, 1, 2, 3, 4, 5, 6, 7, 8, 9, A, B, C, D, E, F). É um sinal digital com equivalentes em quatro dígitos binários conforme a tabela:

Decimal	Hexadecimal	Binário
0	0	0000
1	1	0001
2	2	0010
3	3	0011
4	4	0100
5	5	0101
6	6	0110
7	7	0111
8	8	1000
9	9	1001
10	A	1010
11	B	1011
12	C	1100
13	D	1101
14	E	1110
15	F	1111

Converter de hexadecimal para decimal

a) FACADA
Só fazer normalmente aquela decomposição de cada algarismo:
Lembrando que F = 15, A = 10, C = 12 e D = 13. Fica:
$$10 * 16^0 + 13 * 16^1 + 10 * 16^2 + 12 * 16^3 + 10 * 16^4 + 15 * 16^5 \\
10 + 208 + 2560 + 49152 + 655360 + 15728640 = 16435930_{10}$$
b) BABACA
Antes: B = 11, A = 10, C = 12.
$$10 * 16^0 + 12 * 16^1 + 10 * 16^2 + 11 * 16^3 + 10 * 16^4 + 11 * 16^5 \\
10 + 192 + 2560 + 45056 + 655360 + 11534336 = 12237514_{10}$$
c) CABECA
Antes: C = 12, A = 10, B = 11, E = 14.
$$10 * 16^0 + 12 * 16^1 + 14 * 16^2 + 11 * 16^3 + 10 * 16^4 + 12 * 16^5 \\
10 + 192 + 3584 + 45056 + 655360 + 12582912 = 13287114_{10}$$
d) 2F4D0
Antes: F = 15, D = 13.
$$0 * 16^0 + 13 * 16^1 + 4 * 16^2 + 15 * 16^3 + 2 * 16^4 \\
0 + 208 + 1024 + 61440 + 131072 = 193744_{10}$$

Converter de hexadecimal para binário

Como ambos são sinais digitais, basta substituir pelos binários tabelados:
2FB10
0010.1111.1011.0001.1100
Logo, 101111101100011100 na base 2.

Converter de decimal para hexadecimal

Mesma coisa do resto: divisão por 16 e pegar o resto. Ou fazer aquele esquema de pegar o primeiro algarismo numa multiplicação. Veja bem:
a) 1023
1023/16 = 63, resto 15 (F).
63/16 = 3, resto 15 (F).
3/16 = 0, resto 3.
Logo, 3FF na base 16.

b) 462,24
Parte inteira:
462/16 = 28, resto 14 (E).
28/16 = 1, resto 12 (C).
1/16 = 0, resto 1.
Logo, 1CE na base 16.
Parte fracionária:
0,24 x 16 = 3,84
0,84 x 16 = 13,44 (D)
0,44 x 16 = 7,04
0,04 x 16 = 0,64
E aqui tá bom, né. 3D7.
Assim, podemos considerar 462,24 como 1CE,3D7 na base 16.

c) 16
16/16 = 1, resto 0.
Logo, 10 na base 16.

Exercícios

1. Converta da base 8 para a base 2 ou vice-versa.
a) 101011 na base 8.
NOTE BEM, aqui é a base 8. Não é um binário (não dá pra fazer a base 8 ali bonitinho sem LaTeX e LaTeX que não quebra linha não funciona em todos os navegadores, então só posso pedir para que preste bastante atenção).
Vai ficar assim:
(001)(000)(001)(000)(001)(001)
Ou seja: 1000001000001001 na base 2.

b) 1100 na base 2.
Decompondo: 001.100, que fica 1.4. Logo, 14 na base 8.

c) 11,01 na base 8.
(001)(001),(000)(001)
Logo: 1001,000001.

2. Converta da base 16 para a base 8 ou vice-versa.
Pra isso os srs. terão de converter primeiro pra binário, depois de volta pra base 8 ou pra base 16. Veja só:
a) 10 na base 8.
(001)(000), logo em binário 1000.
1000 na base 16 fica 8.
b) AB na base 16.
(1010)(1011), logo em binário 10101011.
10101011 na base 8 você decompõe da seguinte maneira:
10.101.011, que fica 2.5.3. Que fica 253 na base 8.

3. Converta os números seguintes para as demais bases:
Aqui o esquema é o seguinte: todos os números deverão ser convertidos para todas as bases. Veja só, decimal, binário, octal e hexadecimal: um deles será dado, os outros deverão ser descobertos através dele. Meu conselho é o seguinte:

Descubra sempre o hexadecimal primeiro. O hexadecimal simplifica todos os processos porque ele é o maior, e que vai precisar de menos contas pra resolver. De hexadecimal pra decimal você provavelmente só fará três divisões enquanto octal fará um pouco mais e binário será ridículo, de hexadecimal pra binário é uma tabela maior que octal mas também significa menos decomposição de números, e de hexadecimal pra octal... bem, pra que você vai fazer hexadecimal pra octal primeiro mesmo? "Ah, mas se me derem um octal, eu vou ter que descobrir binário primeiro." então descubra, ué, mas depois faça hexadecimal pra sofrer menos no decimal.
Veja só:

a) 124,5 na base 8.
Pra binário (não dá pra hexadecimal direto):
(001)(010)(100),(101), logo 1010100,101 na base 2.
Pra hexadecimal:

101.0100,1010, não esqueça desse último zero, logo 5.4,A. O que dá 54,A na base 16.

Pra decimal:
$$5 * 16^1 + 4 * 16^0 + \dfrac{10}{16} = 80 + 4 + 0.625 = 84.625_{10}$$

E está respondido o primeiro.

b) 322,21 na base 16.
Pra binário:

(0011)(0010)(0010),(0010)(0001), logo 1100100010,00100001 na base 2.
Pra octal:
1.100.100.010,001.000.010. O que dá 1.4.4.2,1.0.2. Logo, 1442,102 na base 8.
Pra decimal:
$$3 * 16^2 + 2 * 16^1 + 2 * 16^0 + \dfrac{2}{16} + \dfrac{1}{16^2} \\
768 + 32 + 2 + 0.125 + 0.00390625 = 802.12890625_{10}$$
E está respondido o segundo.

c) 100111101 na base 2.
Pra hexadecimal:
1.0011.1101. O que dá 1.3.D. Ou seja, 13D na base 16.
Pra octal:
100.111.101. O que dá 4.7.5. Logo, 475 na base 8.
Pra decimal:
Aqui vamos usar octal porque é mais conveniente que ficar usando letra (tanto octal quanto hexadecimal usam 3 algarismos mesmo). Mas é caso e caso, faça a análise que lhe convier.
5*1 + 7*8 + 4*8² = 5 + 56 + 256 = 317 na base 10.
Resolvido.

d) 7002 na base 10.
Pra hexadecimal:
7002/16 = 437, resto 10 (A).
437/16 = 27, resto 5.
27/16 = 1, resto 11 (B).
1/16 = 0, resto 1.
Ou seja, 1B5A na base 16.
Pra binário:
(0001)(1011)(0101)(1010), ou seja, 1101101011010 na base 2.
Pra octal:
1.101.101.011.010, ou seja, 1.5.5.3.2, logo 15532 na base 8.
E estão resolvidos todos.

Boa noite. \o\

[Cálculo Numérico] Aula 9: sinal digital pra decimal, sistema octal

2012-09-01T13:14:00.002-07:00

1. Converta o sinal digital abaixo (binário) para decimal:
a) 10000001
Se é um sinal digital de 8 bits e está iniciando com 1, temos que é um número negativo. Então vamos convertê-lo pra positivo pra fazer a operação primeiro, fazendo o oposto que fazemos de positivo pra negativo: complemento de um, depois subtrai 1 (soma -1) já que na outra conversão se somava 1.
Complemento de 1
10000001
01111110

Soma (-1)
Primeiro, vale a pena lembrar quem é -1:
1 = 00000001
(Compl. de 1) + 1: (11111110) + 1 = 11111111
Logo, a operação final fica:
01111110 +
11111111
0+1 = 1

1+1 = 10, fica 0 e vai 1
1+1+1 = 11, fica 1 e vai 1
E isso se repete até a última casa, aonde fica 1+1 = 10, fica 0 e vai 1, que some porque sai da última casa do sinal digital.
Sendo assim, temos que esse binário fica: 01111101
Convertendo pra decimal:

$$-n = 1 * 2^0 + 0 * 2^1 + 1 * 2^2 + 1 * 2^3 + 1 * 2^4 + 1 * 2^5 + 1 * 2^6 + 0 * 2^7 \\
1 + 4 + 8 + 16 + 32 + 64 = 125_{10} \\
n = -125_{10}$$
Mas aí uma regra muito doida que eu não entendi se surgiu de algum lugar ou é apenas magia negra diz que você acrescenta -2 ao resultado final, sendo assim, o número é na verdade -127 em decimal. Simplesmente isso.

b) 10010011
Mesma coisa:
Complemento de 1
10010011
01101100

Soma (-1)
01101100 +
11111111
0+1 = 1
0+1 = 1
1+1 = 10, fica 0 e vai 1
1+1+1 = 11, fica 1 e vai 1
1+1 = 10, fica 0 e vai 1
1+1+1 = 11, fica 1 e vai 1
1+1+1 = 11, fica 1 e vai 1
1+1 = 10, fica 0 e vai 1, que é desconsiderado porque explode os 8 bits
Resultado em binário: 01101011
Pra decimal:
$$-n = 1*2^0 + 1*2^1 + 0*2^2 + 1*2^3 + 0*2^4 + 1*2^5 + 1*2^6 + 0*2^7 = \\
1 + 2 + 8 + 32 + 64 = 107 \\
n = -107_{10}$$
Menos 2, -109 decimal.

c) 11111111
Complemento de 1
11111111
00000000

Soma (-1)
00000000 +
11111111
= 11111111
Eu por algum motivo tenho anotado na aula +255 aqui, mas preciso de uma luz nesse. A lógica passada é a seguinte: vai dar -255, menos 2 vai estourar e ir pra +256 e depois +255, mas não lembro se sinal digital vai de -255 a +256 ou de -127 a +128, então para todos os efeitos fiquem com a resposta da aula e perguntemos ao professor depois.

Sistema Octal

Sistema de 8 dígitos: 0, 1, 2, 3, 4, 5, 6, 7, 8.

Converter octal para decimal:

$$123.3_8 = 1 * 8^2 + 2 * 8^1 + 3 * 8^0 + 3 * 8^{-1} \\
123.3_8 = 83.375_{10}$$

Converter decimal em octal:

223
223/8 = 27, resto 7.
27/8 = 3, resto 3.
3/8 = 0, resto 3.
O número então vai dar 0337, ou simplesmente 337 em octal.
495,25
Parte normal
495/8 = 61, resto 7.
61/8 = 7, resto 5.
7/8 = 0, resto 7.
Ou seja, 757.
Parte fracionária: você vai fazer mais ou menos como no binário, multiplicar por 8 o tanto de casas que for preciso e ir pegando o último número. Sobrando 0,00, chega-se à última casa.
0,25 x 8 = 2,00
Pega o 2, deixa o resto, fica 0,00. Então o resultado final é:
757,2

Converter binário em octal

Esse é o mais simples, veja bem:
1011011001
Você vai fragmentar o binário de 3 em 3, contando a partir do primeiro, e vai formar algarismos de 0 a 7 com eles. Simples assim. Veja esse exemplo:
1.011.011.001
1.3.3.1, logo 1331 em octal. (lembrando que o primeiro, isolado, você pode considerar dois zero à esquerda se quiser)

1001110111,100101010
Antes da vírgula, da direita pra esquerda: 1.001.110.111<
1.1.6.7
Depois da vírgula, da esquerda pra direita: 100.101.010
4.5.2
Logo, 1167,452 em octal.

Converter octal em binário

Faça o oposto: pegue cada algarismo do octal e trunque em binários. Veja esse exemplo:
523 em octal
5 = 101
2 = 010
3 = 011
Logo, 101010011. Vale a pena lembrar que a não ser que vá sobrar zero à esquerda ou à direita no algarismo final (ou seja, último algarismo inteiro à esquerda, e fracionário à direita), você obrigatoriamente tem que conservar o 0 do binário.

Exercícios

1. Converta de decimal para octal com 4 dígitos na parte fracionária:
a) 102,4
Parte inteira:
102/8 = 12, sobra 6.
12/8 = 1, sobra 4.
1/8 = 0, sobra 1.
Logo, 146.
Parte fracionária:
0,4 x 8 = 3,2
0,2 x 8 = 1,6
0,6 x 8 = 4,8
0,8 x 8 = 6,4
Logo, 0,3146.
No total, fica: 146,3146 em octal.

b) 1004,5
Parte inteira:
1004/8 = 125, sobra 4.
125/8 = 15, sobra 5.
15/8 = 1, sobra 7.
1/8 = 0, sobra 1.
Logo, 1754.
Parte fracionária:
0,5 x 8 = 4,0
Logo, 0,4.
No total: 1754,4 em octal.

c) 723,45
Parte inteira:

723/8 = 90, sobra 3.
90/8 = 11, sobra 2.
11/8 = 1, sobra 3.
1/8 = 0, sobra 1.
Logo, 1323.
Parte fracionária:
0,45 x 8 = 3,6
0,6 x 8 = 4,8
0,8 x 8 = 6,4
0,4 x 8 = 3,2
Logo, 0,3463.
No total: 1323,3463.

2. Converta de octal para decimal:
a) 732,1

Só fazer a conta tradicional:
$$7 * 8^2 + 3 * 8^1 + 2 * 8^0 + \dfrac{1}{8} = 448 + 24 + 2 + 0.125 = 474.125_{10}$$
b) 423,6
$$4 * 8^2 + 2 * 8^1 + 3 * 8^0 + \dfrac{6}{8} = 256 + 16 + 3 + 0.75 = 275.75_{10}$$

c) 1004,123
$$1 * 8^3 + 0 * 8^2 + 0 * 8^1 + 4 * 8^0 + \dfrac{1}{8} + \dfrac{1}{8^2} + \dfrac{1}{8^3} \\
512 + 4 + 0.125 + 0.03125 + 0.005859375 = 516.162109375_{10}$$

3. Converta de binário para octal:
a) 11100111
011.100.111
3.4.7
Logo, 347.
b) 1100111
001.100.111
1.4.7
Logo, 147.
c) 1111011
001.111.011
1.7.3
Logo, 173.
d) 10,000010001
Parte inteira: 010. Logo, 2.
Parte fracionária:
000.010.001
0.2.1
Logo, 0,021.
Sendo assim, o total é 2,021.

Já vou postando a outra aula e a de equações diferenciais, não confirmo de algoritmo e estrutura de dados ainda porque acho melhor o professor passar o algoritmo dele e tal pra não confundir ninguém. Aliás, estou realmente cogitando postar estrutura de dados porque é uma matéria muito mais complexa que programação comum, mas ainda envolve muito código e acho código extremamente pessoal, daí não sei uma maneira de passar um "código didático" pra vocês que estão interessados.
Discutam aí, se alguém souber essa maneira e quiser passar aqui no blog também, sem problema.

No mais, bom dia a todos.

[Equações Diferenciais Ordinárias] Aula 7: exatas e homogêneas

2012-08-28T19:51:00.000-07:00

Equações Diferenciais Ordinárias de 1ª Ordem

Forma diferencial

$$y' = \dfrac{dy}{dx} = f(x, y) \\
f(x, y)dx - dy = 0 \\
\dfrac{dy}{dx} = \dfrac{M(x, y)}{N(x, y)} \\
M(x, y)dx + N(x, y)dy = 0$$
Exemplo:
$$\dfrac{dy}{dx} = \dfrac{x+y}{y^2} = \dfrac{\dfrac{1}{y^2}}{\dfrac{1}{x+y}}$$
Forma a: (x+y)dx - y²dy = 0
Forma b: $$(\dfrac{x+y}{y^2})dx - 1dy = 0; y \neq 0$$
Forma c: $$1dx - \dfrac{y^2}{x+y} dy = 0; x+y \neq 0$$

Classificação das Equações Diferenciais Ordinárias de 1ª ordem

Equação Diferencial Ordinária Linear
$$y' + a_0 (x)y = f(x) \\ y' + p(x)y = q(x)$$

Equação Diferencial Ordinária Homogênea
y' = f(tx, ty) = f(x, y)

Exemplo: considere: $$f(x, y) = \dfrac{x+y}{y} \\
f(tx, ty) = \dfrac{tx + ty}{tx} = \dfrac{t(x + y)}{tx}$$
Isso nos dá: $$f(tx, ty) = f(x, y)$$
No caso geral: se f(x, y) é homogênea de grau n, então: $$f(tx, ty) = t^n f(x, y)$$
Se n = 0, então f(tx, ty) = f(x, y)

Equação Diferencial Ordinária Separável
M(x)dx + N(y)dy = 0

Equação Diferencial Ordinária Exata
M(x, y)dx + N(x, y)dy = 0

A equação diferencial ordinária acima é exata quando é válida a relação:
$$\dfrac{\partial}{\partial y} M(x, y) = \dfrac{\partial}{\partial x} N(x, y)$$
Exemplo: resolva a equação diferencial ordinária A:
$$A = x^2 dx + xe^y dy = 0$$
a) Verifique se a EDO é separável.
Resolução: dividindo cada termo por x:
$$\dfrac{x^2}{x} dx + \dfrac{x e^y}{x} dy = \dfrac{0}{x} \\
xdx + e^y dy = 0 \\
\int x dx + \int e^y dy = \int 0 \\
\dfrac{x^2}{2} + C_1 + e^y + C_2 = C_3 \\
\dfrac{x^2}{2} + e^y = C_3 - C_1 - C_2$$
Como os três termos são constantes, podemos concatená-los em uma única constante:
$$e^y = C - \dfrac{x^2}{2} \\
y(x) = \ln{|C - \dfrac{x^2}{2}|}$$

Exemplo com equação diferencial ordinária homogênea:
$$y' = f(tx, ty) = f(x, y), n = 0 (f(tx, ty) = t^n f(x, y))$$
Considere: $$z = \dfrac{2y^4 + x^4}{xy^3}$$
Identificando z = f(x, y).
$$z(tx, ty) = \dfrac{2(ty)^4 + (tx)^4}{(tx)(ty)^3} = \dfrac{t^4 2y^4 + t^4 x^4}{tx * t^3 y^3} = z(x, y)$$

Mudança de Variável

Considere y' = f(x, y), onde f(x, y) = f(tx, ty).
Fazendo a seguinte mudança de variável:
$$v = \dfrac{y}{x} \\
y = vx$$
Temos a função:
$$y = vx \\
\dfrac{dy}{dx} = v \dfrac{dx}{dx} + x \dfrac{dv}{dx} = v + x \dfrac{dv}{dx} \\
x \dfrac{dv}{dx} + v = f(v) \\
x \dfrac{dv}{dx} = f(v) - v$$
O que é uma equação separável.
$$\int \dfrac{dv}{f(v) - v} = \int \dfrac{dx}{x} \\
\int \dfrac{dv}{f(v) - v} = \int \dfrac{1}{x} dx + C$$
Observação: a variável v = y/x para o caso em que:
$$y'(x, y) = f(x, x \dfrac{y}{x}) = f(1, v)$$
Já para:
$$y'(x, y) = f(y \dfrac{x}{y}, y) \\
v = \dfrac{x}{y}, f(v, 1)$$

Exemplo: considere a equação diferencial homogênea:
$$y' = \dfrac{x+y}{y}$$
a) Determine a função y(x):
Resolução: definido v = y/x:
$$y = vx \\
y' = v + xv'$$
Substituímos y' pelo valor dado inicialmente ali:
$$\dfrac{x+y}{y} = v + xv'$$
Podemos trocar y por vx, pela igualdade estabelecida logo abaixo da definição de v:
$$\dfrac{x + vx}{x} = v + xv'$$
Colocando x em evidência:
$$\dfrac{x (1+v)}{x} = v + xv' \\
(1+v) = v+xv' \\
1 + v - v = x \dfrac{dv}{dx} \\
1 = x \dfrac{dv}{dx} \\
\dfrac{dx}{x} = dv \\
\int \dfrac{1}{x} dx = dv \\
\ln{|x|} + C = v \\
\dfrac{y}{x} = \ln{|x|} + C \\
y(x) = x(\ln{|x|} + C)$$

Tentei explicar um pouquinho melhor essa última porque muita gente ainda pareceu empacadíssimo nela, mas não sei se fiz um bom trabalho. É tudo muito pré-estabelecido, 1000x mais fácil entender numa resolução de exercício de fato, que provavelmente será a próxima lista.
De qualquer forma, boa noite a todos.

[Cálculo Numérico] Aula 8: subtração de binários

2012-08-25T19:21:00.000-07:00

Subtração de binários

Como os sinais digitais não reconhecem (-) como números negativos o processo de subtração entre A e B deve ser dado por A + (-B) onde -B representa o oposto de B.
Esta determinação de oposto de um binário pode ocorrer de três formas:
1. sinal-magnitude: representado em sinais digitais de 8 bits onde o primeiro bit é o indicativo da magnitude.
Exemplo: (0)0001001 = +9, enquanto (1)0001001 = -9.
Ou seja, todo binário que possui N-1 = 0 é positivo e N-1 = 1 é negativo. Existe assim 2^8 possibilidades de combinação. 2^8 = 256 dígitos, sendo 128 positivos e 128 negativos, mas, observe que 00000000 = +0 e 10000000 = -0 assumem o mesmo valor. Logo:
127 combinações positivas, 127 combinações negativas, e 0. Um total de 255 combinações.

2. complemento de um: troca-se os 0 por 1 e os 1 por zero.
+9 = 00001001
-9 = 11110110
Novamente, 00000000 = +0 e 11111111 = -0. Por isso temos 255 valores de -127 a +127.

3. complemento de dois: realiza-se o complemento de um e soma-se mais um:
+9 = 00001001
-9 = 11110110 (comp. de um de 9) + 1 = 11110111
No caso de 0:
+0 = 00000000
-0 = 11111111 (comp. de um de 0) + 1 = 100000000
Como estamos falando de 8 bits, ignoramos o nono caracter à esquerda. Ou seja: -0 = 00000000. Sendo assim, esse sistema é o mais favorável, pois permite de -128 a +127.

Exercícios:
1. Resolva em binário:
a) (antes dele alterar) 5 - 2

5 + (-2)
5 é muito fácil, né:
4 - 2 - 1
1 - 0 - 1

Logo, 5 = 00000101
Já o -2 é um pouquinho mais complicado, mas nada demais. Lembrando que usaremos complemento de dois:
2 = 00000010
Comp.1 de 2 é 1111101. +1 = 11111110. Logo a operação fica:
00000101 +
11111110

Primeiro algarismo, 1. Segundo algarismo, 1. Terceiro algarismo, 10, fica 0 e vai 1. Aí vai dando 10 (fica 0 e vai 1) até estourar a oitava casa, o que dá: 00000011 (3 no sistema decimal). Mas note que estamos usando um número muito grande, então vamos usar do seguinte aconchego: só vamos até a última casa do número maior, passou disso começaremos a ignorar porque sabemos que vai estourar a oitava casa.
Perceba na b:

b) 3 - 2
3 + (-2)

Não precisamos fazer conta pra saber que 3 = 11.
2 = 10. Com complemento de 1, fica 01. Acrescentando 1, 10.
Logo a operação fica:
11 +
10
1 + 0, fica 1.
1 + 1, fica 0 e vai 1. Mas a última casa é essa, então ignoramos tudo depois dela.

Logo, 3 - 2 dá 01 que, bem, é igual a 1 em decimal.

a) (depois dele alterar) 5 - 1
5 + (-1)

Já sabemos que 5 = 101.
1 é 001. Complemento de 1 fica 110, somando 1 fica 111.
Logo, a conta será:
101 +
111
1+1 dá 10, fica 0 e sobe 1.

0+1+1 dá 10, fica 0 e sobe 1.
1+1+1 dá 11, fica 1 e sobe 1, que ignoramos.
Logo, 5 - 1 fica 100, que é 4.

c) 12 - 6
12 + (-6)

Pra 12:

8 - 4 - 2 - 1
1 - 1 - 0 - 0

Pra 6:
8 - 4 - 2 - 1
0 - 1 - 1 - 0
Complemento de 1: 1001. Soma-se 1: 1010.
Logo a operação fica:

1100 +
1010
0+0 = 0.

0+1 = 1.
1+0 = 1.
1+1 = 10, fica 0 e sobe 1, que ignoramos.
Logo, essa conta dá o resultado 110, que em decimal fica (obviamente) 6.

d) 7 - 5
7 + (-5)

Pra 7:
4 - 2 - 1

1 - 1 - 1

Pra 5:
4 - 2 - 1
1 - 0 - 1
Complemento de 1: 010. Soma-se 1: 011.
Logo a operação fica:

111 +
011
1+1 = 10, fica 0 e sobe 1.

1+1+1 = 11, fica 1 e sobe 1.
1+1 = 10, fica 0 e sobe 1, que ignoramos.
Logo, essa conta dá o resultado de 10, que em decimal fica 2.

E é basicamente isso. Mais uma vez lembro que não sei ao certo qual estrutura de post uso pra essa matéria de binários aqui, porque não sei se está claro dessa forma, ainda mais porque não consigo determinar a base sem usar o LaTeX só que me parece muito inconveniente usar o LaTeX nessa situação.
É vocês que me contatam e decidem. No mais, boa noite.

[Cálculo Numérico] Aula 7: fração decimal pra binária

2012-08-25T18:02:00.002-07:00

Decimal para Binário (Nº fracionário)

Peguemos o número 35,625 na base 10. Vamos convertê-lo para binário.
1º Passo: separamos as partes inteiras das fracionárias:
35 + 0,625
A parte inteira é bem fácil. Como o número é menor que 256, ao invés de fazer divisão, é mais fácil usar aquele método que o professor passou de ordenar decrescentemente as potências de 2:
32 - 16 - 8 - 4 - 2 - 1
É intuitivo, deve ter um método matemático pra descobrir o que usar, mas honestamente... Não é necessário. É só pensar um pouquinho. Queremos chegar ao 35, então 32+2+1 resolve.
Logo:
32 - 16 - 8 - 4 - 2 - 1
1 - 0 - 0 - 0 - 1 - 1

Logo:
$$35_{10} = 100011_2$$
Para a fracionada, temos a seguinte metodologia: ao invés de sair dividindo como feito na aula passada, vamos multiplicando por 2, pegando o algarismo depois da vírgula e removendo ele pra colocar no binário. E repetindo até chegar no tanto de casas após a vírgula exigidas. Ao contrário do binário inteiro, o fracionário segue ordem normal.
Explicando fica meio complicado mesmo, então vamos fazer na prática que dá pra entender bem:
$$0.625 \times 2 = 1.250$$
Pegamos esse 1 depois da vírgula, deixamos 0,250 e multiplicamos ele por 2 de novo.
$$0.250 \times 2 = 0.5$$
Pegamos o 0, repetimos o processo:
$$0.5 \times 2 = 1$$
Pegamos o 1 e terminamos o número. Esse caso é bom porque é perfeito: o número tem um fim simples. Muitos deles serão gigantescos e o exercício exigirá um mínimo de casas, então não se preocupe tanto com isso. E o resultado final?
$$0.625_{10} = 0.101_2$$
Logo:
$$35.625_{10} = 100011.101_2$$

Exercícios:
1. Converter para decimal.
a) 11010101011
Binário pra decimal é deveras simples, está na outra aula o método. Só sair somando:
$$1*2^0 + 1*2^1 + 0*2^2 + 1*2^3 + 0*2^4 + 1*2^5 + \\ 0*2^6 + 1*2^7 + 0*2^8 + 1*2^9 + 1*2^{10}$$
$$1 + 2 + 0 + 8 + 0 + 32 + 0 + 128 + 0 + 512 + 1024 \\
1707_{10}$$
b) 1,01010101
Mesma coisa, só que com valores mais feios:
$$(1 * 2^0) + (0 * 2^{-1} + 1 * 2^{-2} + 0 * 2^{-3} + 1 * 2^{-4} + \\ 0 * 2^{-5} + 1 * 2^{-6} + 0 * 2^{-7} + 1 * 2^{-8}$$
$$(1) + (0 + 0.25 + 0 + 0.0625 + 0 + 0.015625 + 0 + 0.00390625) \\
1 + (0.33203125) \\
1.33203125_{10}$$
c) 11011,1101
Agora juntando o útil ao agradável:
$$(1 * 2^4 + 1 * 2^3 + 0 * 2^2 + 1 * 2^1 + 1 * 2^0) + (1 * 2^{-1} + 1 * 2^{-2} + 0 * 2^{-3} + 1 * 2^{-4}) \\
(16 + 8 + 0 + 2 + 1) + (0.5 + 0.25 + 0 + 0.0625) \\
27 + 0.8125 \\
27.8125_{10}$$

2. Converter para binário com 5 dígitos.
a) 234,435
Primeiro, a parte antes da vírgula (234):
128 - 64 - 32 - 16 - 8 - 4 - 2 - 1
128+64 = 192
192+32 = 224
224+16 = 240, então pulamos pro próximo.
224+8 = 232
Sabemos que é somando 2 que temos 234, óbvio, então:
128 - 64 - 32 - 16 - 8 - 4 - 2 - 1
1 - 1 - 1 - 0 - 1 - 0 - 1 - 0

234 = 11101010 em decimal
Pra parte depois da vírgula, 435:
0,435 x 2 = (0),87
0,870 x 2 = (1),74

0,740 x 2 = (1),48
0,480 x 2 = (0),96
0,960 x 2 = (1),92
E tá bom, né, pediu 5 casas decimais depois da vírgula.

Ou seja:
$$0.435 = 0.01101_2 \\
234.435 = 11101010.01101_2$$
b) 12,234
Parte depois da vírgula:
8 - 4 - 2 - 1

1 - 1 - 0 - 0
(acho que é muito "de cabeça" que 12 = 8 + 4)
Pro 0,234:
0,234 x 2 = (0),474
0,474 x 2 = (0),948
0,948 x 2 = (1),896
0,896 x 2 = (1),792
0,792 x 2 = (1),584
Logo:

$$0.234 = 0.00111_2 \\
12.234 = 1100.00111_2$$
c) 43,956
Parte depois da vírgula:

32 - 16 - 8 - 4 - 2 - 1
32 + 16 = 48, passa.

32 + 8 = 40
40 + 4 = 44, passa.
40 + 2 = 42
42 + 1 = 43
Logo:
32 - 16 - 8 - 4 - 2 - 1
1 - 0 - 1 - 0 - 1 - 1
Assim, 43 = 101011.

Agora, pro que vem depois da vírgula:
0,956 x 2 = (1),912
0,912 x 2 = (1),824
0,824 x 2 = (1),648
0,648 x 2 = (1),296
0,296 x 2 = (0),592
E 0,956 fica em 5 casas decimais até 0,11110. E:

$$43.956 = 101011.11110_2$$

Operações com binário

Adição de dois termos: escreva em binário o nº 2.
$$2 + 2 = 4 \\
4_{10} = 10_2 + 10_2 \\
10_2 + 10_2 = 100_2$$
A regra é a seguinte:
0 + 0 = 0
1 + 1 = 10
1 + 0 = 1
0 + 1 = 1

Veja a seguinte conta (com o sinal do outro lado por conveniência):

1101101 +
1001101
Vamos fazer a operação normalmente de trás pra frente, passo a passo, como se fazia no primário. Lembrando que 1+1 = 10. Então fica 0 e sobe 1.
Primeira casa: 1+1 = 10. Fica 0, sobe 1.
Segunda casa: 0+0+1 = 1. Fica 1.
Terceira casa: 1+1 = 10. Fica 0, sobe 1.
Quarta casa: 1+1+1 = 11. Fica 1, sobe 1.
Quinta casa: 0+0+1 = 1. Fica 1.
Sexta casa: 1+0 = 1. Fica 1.
Sétima casa: 1+1 = 10. Como é a última casa, fica 10.
Pega tudo o que está na frente do "fica", obviamente na ordem inversa (da sétima pra primeira):

(10)(1)(1)(1)(0)(1)(0)
Logo, o resultado é: 10111010.

Se tiver qualquer dificuldade nisso, o mais recomendado é pegar algum livro de cálculo numérico ou eletrônica digital e sair queimando todos os exercícios até entender perfeitamente o que está acontecendo. Só dá trabalho, mas é facinho.
Também sei que a estrutura do post ficou meio estranha, não sei ao certo qual seria a melhor pra trabalhar aqui. Se ficou muito pouco claro pra quem tá estudando, sou aberto a sugestões de modelos diferentes.
Gratíssimo, boa noite.

[Física Experimental] Aula 3 e 4: relatório 1 e 2

2012-08-25T17:19:00.000-07:00

Experimento 1: movimento retilíneo uniforme

Cuidados:

não movimentar o cavaleiro com o compressor desligado;
não colocar o compressor em potência máxima;
não alterar a potência do compressor durante o experimento;
não permitir o choque do cavaleiro com o final do trilho;
alinhar os sensores para evitar o choque com o cavaleiro;
deixar os sensores perpendiculares ao trilho;
utilizar para o tempo todos os dígitos do cronômetro.

Equipamentos:

compressor de ar;
trilho de ar;
cavaleiro metálico;
sensor de start;
sensor de stop;
cronômetro digital.

Procedimentos:

sensor start: S0 = 0,2m;
sensor stop: S1 = 0,5m; S2 = 0,6m; S3 = 0,7m; S4 = 0,8m; S5 = 0,9m;
para cada ∆S determinar t1, t2, t3, t4 e t5.

O resto já passou, galera. Não vai adiantar postar metodologia de coisa que já passou e a gente não vai mais usar (no caso, o gráfico). Eu só postei isso aqui, agora, porque não consegui postar antes e porque o experimento 2 que vem logo a seguir vai utilizar:

Experimento 2: movimento retilíneo uniformemente variado (M.U.V.)

Equipamento:

tudo o que foi usado no relatório 1 (só por isso eu postei ele aqui)
barbante
suporte para massa
massa de 50g

Cuidados:

não soltar o cavaleiro com a massa em movimento;
o sensor de start deve ser colocado o mais próximo possível do início do movimento;
não deixar o barbante fora da roldana.

Objetivo geral:

determinar a aceleração do sistema;
determinar a velocidade final.

(obviamente não é assim que tu vai colocar no relatório, né)

Função horária do deslocamento:
Aceleração: $$S = S_0 + V_0 t + \dfrac{at^2}{2}$$
Sendo V0 = 0, temos:
$$S - S_0 = \dfrac{at^2}{2} \\ \Delta S = \dfrac{at^2}{2} \\ 2 \Delta S = at^2 \\ a = \dfrac{2 \Delta S}{t^2}$$
Função horária da velocidade:
$$a = \dfrac{\Delta V}{t} \\ \dfrac{2 \Delta S}{t^2} = \dfrac{\Delta V}{t} \\ \dfrac{2 \Delta S}{t} = \Delta V$$
Lembrando que V0 = 0, e que delta V = V - V0:
$$V = \dfrac{2 \Delta S}{t}$$
Posições que serão/foram usadas no experimento:
S1 = 0,45m; S2 = 0,56m; S3 = 0,67m; S4 = 0,78m; S5 = 0,89m

Utilizar apenas uma massa.

Gráfico: V x t. Usar Microsoft Excel ou Origin 7.0.

[Equações Diferenciais Ordinárias] Aula 6: formas de equações diferenciais

2012-08-25T16:38:00.002-07:00

b) y'' = 2; y(1) = 1; y'(-1) = 1
Primeiramente, sim, estamos utilizando graus.
$$y(1) = 1 \\ c_1 \sin{2} + c_2 \cos{2} = 1$$
$$y'(-1) = 1 \\ 2c_1 \cos{(-2)} - 2c_2 \sin{(-2)} = 1 \\ 2c_1 \cos{(2)} + 2c_2 \sin{(2)}$$
O professor explicou bem a mudança de sinais mas, pra quem não entendeu, cosseno de -90° a 90° é positivo e seno de -180° a 0° é negativo, então você pode igualar cosseno de -2 a cosseno de 2, e inverter o sinal de seno de -2 para obter seno de 2.
Mas bem, agora temos um sistema linear:
$$\begin{array}
c_1 \sin{2} + c_2 \cos{2} = 1 \\
2c_1 \cos{2} + 2c_2 \sin{2} = 1 \end{array}$$
E a solução que o professor nos arranjou foi a seguinte (creio que hajam outras e até mais fáceis, mas por questão de evitar confusão vou usar essa mesmo):
$$-2 \dfrac{\sin{2}}{\cos{2}}$$
Multiplicamos a primeira expressão por esse termo:
$$\begin{array}
c_1 \sin{2} \times (-2) \dfrac{\sin{2}}{\cos{2}} + c_2 \cos{2} \times (-2) \dfrac{\sin{2}}{\cos{2}} = -2 \dfrac{\sin{2}}{\cos{2}} \\
2c_1 \cos{2} + 2c_2 \sin{2} = 1 \end{array}$$
Vamos simplificar o máximo que podemos essa expressão:
$$\begin{array}
-2c_1 \dfrac{\sin^2 {2}}{\cos{2}} - 2c_2 \sin{2} = -2 \dfrac{\sin{2}}{\cos{2}} \\
2c_1 \cos{2} + 2c_2 \sin{2} = 1 \end{array}$$
Somamos o sistema agora:
$$-2c_1 \dfrac{\sin^2 {2}}{\cos{2}} - 2c_2 \sin{2} + 2c_1 \cos{2} + 2c_2 \sin{2} = -2 \dfrac{\sin{2}}{\cos{2}} + 1$$
Vamos deixar c1 de um lado e c2 do outro pra ficar mais fácil de compreender o que vai acontecer. E aproveitar e trocar sen2/cos2 por tangente de 2 do outro lado:
$$(-2c_1 \dfrac{\sin^2 {2}}{\cos{2}} + 2c_1 \cos{2}) + (2c_2 \sin{2} - 2c_2 \sin{2}) = (-2 \tan{2} + 1)$$
Veja bem: podemos cortar os c2, restando só c1 de variável:
$$2c_1 \cos{2} - 2c_1 \dfrac{\sin^2 {2}}{\cos{2}} = -2 \tan{2} + 1$$
Isolando 2c1:
$$2c_1 (\cos{2} - \tan{2} \sin{2}) = 1 - 2 \tan{2} \\ 2c_1 = \dfrac{(1 - 2 \tan{2})}{(\cos{2} - \tan{2} \sin{2})} \\ c_1 = \dfrac{(1 - 2 \tan{2})}{2(\cos{2} - \tan{2} \sin{2})}$$
Jogando isso na nossa maravilhosa calculadora:
$$c_1 \approx \dfrac{0.930}{1.996} = 0.47$$
Substituindo c1 = 0.47 na primeira equação do sistema linear, temos:
$$0.47 \sin{2} + c_2 \cos{2} = 1 \\ c_2 \cos{2} = 1 - 0.47 \sin{2} \\ c_2 = \dfrac{1 - 0.47 \sin{2}}{\cos{2}} = 0.98$$
Portanto para exemplo b) onde a EDO possui as tais condições iniciais, encontramos c1 = 0.47 e c2 = 0.98.

Resolvido isso, tem o slide com a apresentação das formas de equações diferenciais. Aqui estão os exercícios resolvidos:
Passagem para a forma padrão: (y' = f(x, y))
a) $$xy' - y^2 = 0 \\ xy' = y^2 \\ y' = \dfrac{y^2}{x}$$
b) $$e^x y' - e^{2x} y = \sin{x} \\ e^x y' = \sin{x} + e^{2x} y \\ y' = \dfrac{\sin{x} + e^{2x} y}{e^x} = \dfrac{\sin{x}}{e^x} + e^x y$$
c) $$(y - y')^5 = \sin{(\dfrac{y'}{x})}$$
Não creio que haja maneira de colocar isso na forma padrão, inclusive o professor também não. Se alguém conseguir, seria interessantíssimo mostrar. O dilema aqui é o seguinte: dentro do seno tem um y', e dentro do polinômio também. Pra retirar o y' do seno tem que fazer arco seno dos dois lados, daí pra você tirar o y' do polinômio depois você terá de fazer seno dos dois lados, e volta à situação inicial. É paradoxal, não tem como fazer, então deixa ele assim de boa.

Passagem para a forma diferencial: (f(x, y)dx - f(x, y)dy = 0)
a) $$y(yy' - 1) = x \\ y^2 y' - y = x \\ y^2 \dfrac{dy}{dx} = x + y \\ \dfrac{dy}{dx} = \dfrac{x + y}{y^2} \\ y^2 dy = (x+y)dx \\ y^2 dy - (x+y)dx = 0$$

b) $$y' = \dfrac{y}{x} \\ \dfrac{dy}{dx} = \dfrac{y}{x} \\ xdy = ydx \\ xdy - ydx = 0$$

Identificar as equações lineares: (y' + a(x)y = f(x); y'' + a1(x)y' + a2(x)y = f(x); e por aí vai)
a) $$y' = (\sin{x})y + e^x \\ y' - (\sin{x})y = e^x$$
Perceba, y' correto, a(x) = -sen(x) e o y ali, e f(x) = e^x. Perfeito, é linear.
b) $$y' = y^2 + x$$
Pelo y² você já percebe que não é linear.
c) $$y' + xy^5 = 0$$
Mesmo motivo da b, só que com y elevado a quinta.
d) $$xy' + y = \sqrt{y}$$
Perceba que temos f(y) do outro lado, e y é a função incógnita. Pra ser linear, teria de ser f(x), já que x é a variável independente. Não é linear.
e) $$y' + xy = e^x y$$
Do outro lado temos f(x, y), não f(x). Não é linear.
f) $$y' + \dfrac{x}{y} = 0$$
y está elevado a -1, não a 1, então não é linear.
Sendo assim, a letra A é linear.

Essa aula aqui vai ajudar demais vocês na lista 2, já que ela é praticamente só isso. No mais, boa noite a todos.

[Equações Diferenciais Ordinárias] Aula 5: soluções de equações e problemas de valor inicial

2012-08-25T15:29:00.001-07:00

Equações Diferenciais Ordinárias

Introdução: uma equação diferencial é uma equação que envolve uma função incógnita e suas derivadas.
Ex.: $$y' = 2x + 3; y'' = 2 \\ \dfrac{d^2 y}{dy^2} + (\dfrac{dy}{dx})^2 = 1 \\ y''' + \sin{(x)} y'' + 5xy = 0 \\ y''' - 2y'' + 2y = \sin{(x)}$$

Equações Diferenciais Ordinárias Lineares

São do tipo:
$$y' + a_0 (x) y = f(x) \\ y'' + a_1 (x)y' + a_0 (x)y = f(x)$$
Ex.: $$y'' + \sin{(x)} y' + xy = \cos{(x)}$$
É dada linear, pois y¹.
Caso não-linear:
$$(1 + y^2) \dfrac{d^2 y}{dy^2} + t \dfrac{dy}{dx} + y = e^t$$
É dada não-linear, pois y².

De modo geral, temos a seguinte condição:
$$\sum_{i = 1}^{n} a_1 (x) \dfrac{d^{(i)} y}{dx^{(i)}} = f(x) \\ a_n (x) \neq 0$$
E x pertencendo aos números reais.
Notação:
$$y(x) \\ y' = \dfrac{dy}{dx}; y'' = \dfrac{d^2 y}{dx^2}; y^{[n]} = \dfrac{d^n y}{dx^n}$$
$$x(t) \\ \dot{x} = \dfrac{dx}{dt} \\ \ddot{x} = \dfrac{d^2 x}{dt^2}$$

Exemplo de Solução Verdadeira

y'' + 4y = 0. Verifique se y = sen(2x) e y = cos(2x) são soluções da EDO.
a) Solução para y = sen(2x).
$$y = \sin{(2x)} \\ y' = 2 \cos{(2x)} \\ y'' = -4 \sin{(2x)}$$
Substituindo na EDO, temos:
$$-4 \sin{(2x)} + 4 \sin{(2x)} = 0$$
Logo, y = sen(2x) é solução.
Se procedermos de maneira análoga, verificaremos que y = cos(2x) também é solução da EDO em questão. Portanto a soma das duas soluções também será uma solução. De modo geral, a solução linear das soluções:
$$y_{geral} = y_1 + y_2$$
Onde, nesse caso:
$$y_1 = c_1 \sin{(2x)} \\ y_2 = c_2 \cos{(2x)}$$
Portanto,
$$y_{geral} = c_1 \sin{(2x)} + c_2 \cos{(2x)}$$
E a solução particular se dá por:
$$\begin{array}
y_1 = c_1 \sin{(2x)} \\
y_2 = c_2 \cos{(2x)} \end{array}$$
Onde c1 e c2 pertencem aos números reais.

A determinação dos coeficientes constantes c1 e c2 se dá por meio de:
1. condições iniciais (PVI - problema de valor inicial)
2. condições de contorno (PVC - problema de valor de contorno)

Inserindo uma condição inicial qualquer no y geral que usamos no outro exercício:
y'' = 2; y'(0) = 0; y(0) = 0
Aplicamos as condições e temos:
$$y(x) = c_1 \sin{(2x)} + c_2 \cos{(2x)}$$
$$y(0) = 0 \\ 0 = c_1 \sin{0} + c_2 \cos{0} = 0 + c_2 * 1 \\ c_2 = 0$$
$$y'(x) = 2c_1 \cos{(2x)} - 2c_2 \sin{(2x)} \\ y(0) = 0 \\ 0 = 2c_1 \cos{0} - 2c_2 \sin{0} = 2c_1 * 1 - 0 = 2c_1 \\ c_1 = 0$$
A segunda eu vou deixar todinha pra outra aula, embora tenha se iniciado nessa. É melhor pra todos.

[Equações Diferenciais Ordinárias] Lista 1: classificações, soluções e valores iniciais e de contorno

2012-08-23T17:08:00.000-07:00

1. Determine a ordem, a função incógnita e a variável independente para cada uma das equações diferenciais dadas.

a) (y'') - 3yy' + xy = 0
Ordem: o termo mais derivado é derivado duas vezes, logo segunda ordem.
Função incógnita: a função que está sendo derivada é y, logo, y.
Variável independente: y está sendo derivada com relação a x (não tem outra variável na equação), então, x.

b) t²s'' - ts' = 1 - sen(t)
Ordem: o termo mais derivado é derivado duas vezes, logo segunda ordem.
Função incógnita: a função que está sendo derivada é s.
Variável independente: s está sendo derivada com relação a t.

c) $$\dfrac{d^2 x}{dy^2} = y^2 + 1$$
Ordem: o termo mais derivado é derivado duas vezes, logo segunda ordem.
Função incógnita: a função que está sendo derivada é x.
Variável independente: x está sendo derivada com relação a y.

d) $$(\dfrac{d^2 y}{dx^2})^{3/2} + y = x$$
Ordem: o termo mais derivado é derivado duas vezes, logo segunda ordem.
Função incógnita: a função que está sendo derivada é y.
Variável independente: y está sendo derivada com relação a x.

e) $$(\dfrac{db}{dp})^7 = 3p$$
Ordem: o termo mais derivado é derivado uma vez, logo primeira ordem.
Função incógnita: a função que está sendo derivada é b.
Variável independente: b está sendo derivada com relação a p.

2. Quais das seguintes funções são soluções da equação diferencial y' - 5y = 0?
Só sair testando, o resultado tem que dar 0:

a) y = 5

(5)' - 5(5) = 0 - 25 = -25

Logo, não é solução.
b) y = 5x
(5x)' - 5(5x) = 5 - 25x
Logo, não é solução.
c)
$$y = x^5 \\ (x^5)' - 5(x^5) = 5x^4 - 5x^5 \neq 0$$
Logo, não é solução.
d)
$$y = e^{5x} \\ (e^{5x})' - 5(e^{5x}) = 5e^{5x} - 5e^{5x} = 0$$

Logo, é solução.
e)
$$y = 2e^{5x} \\ (2e^{5x})' - 5(2e^{5x}) = 2(5e^{5x}) - 10e^{5x} = 10e^{5x} - 10e^{5x} = 0$$

Logo, é solução.
f)
$$y = 5e^{5x} \\ (5e^{5x})' - 5(5e^{5x}) = 5(5e^{5x}) - 25e^{5x} = 25e^{5x} - 25e^{5x} = 0$$
Logo, é solução.

3. Quais das seguintes funções são soluções da equação diferencial y' - 3y = 6?
a) y = -2

(-2)' - 3(-2) = 0 + 6 = 6

Logo, é solução.
b) y = 0
0' - 3(0) = 0 - 0 = 0
Logo, não é solução.
c)
$$y = e^{3x} - 2 \\ (e^{3x} - 2)' - 3(e^{3x} - 2) = 3e^{3x} - 3e^{3x} + 6 = 0 + 6 = 6$$

Logo, é solução.
d)
$$y = e^{2x} - 3 \\ (e^{2x} - 3)' - 3(e^{2x} - 3) = 2e^{2x} - 3e^{2x} + 9 = -e^{2x} + 9$$
Logo, não é solução.
e)
$$y = 4e^{3x} - 2 \\ (4e^{3x} - 2)' - 3(4e^{3x} - 2) = 12e^{3x} - 12e^{3x} + 6 = 0 + 6 = 6$$
Logo, é solução.

4. Quais das seguintes funções são soluções da equação diferencial A?
$$A = \dfrac{dy}{dx} = \dfrac{2y^4 + x^4}{xy^3}$$

Essa é complicada. Vamos fazer o seguinte: está permitido passar os valores de um lado pro outro, se os dois lados se igualarem no final está certo. Então:
a) y = x
$$\dfrac{dx}{dx} = \dfrac{2x^4 + x^4}{x(x^3)} \\ 1 = \dfrac{3x^4}{x^4} = 3 \neq 1$$

Logo, não é solução.
b) Não sei se recomendo estudarem esse
$$y = x^8 - x^4 \\ \dfrac{d}{dx} [x^8 - x^4] = \dfrac{2[x^8 - x^4]^4 + x^4}{x[x^8 - x^4]^3} \\ 8x^7 - 4x^3 = \dfrac{2(x^8 - x^4)^4 + x^4}{x(x^8 - x^4)^3} \\ (8x^8 - 4x^4)(x^8 - x^4)^3 = 2(x^8 - x^4)^4 + x^4 \\ 4(2x^8)(x^8 - x^4)^3 - 4(x^4)(x^8 - x^4)^3 = 2(x^8 - x^4)^4 + x^4 \\ 4(2x^8) - 4(x^4) = 2(x^8 - x^4) + \dfrac{x^4}{(x^8 - x^4)^3} \\ 2x^8 - x^4 = \dfrac{(x^8 - x^4)}{2} + \dfrac{x^4}{4(x^8 - x^4)^3} \\ 2 - 1 = \dfrac{1}{2} + \dfrac{x^4}{4(x^8 - x^4)^3} = 1 \\ 2 = 1 + \dfrac{x^4}{4(x^8 - x^4)^3} \\ 1 = \dfrac{x^4}{4(x^8 - x^4)^3} \\ 4 = \dfrac{x^4}{(x^8 - x^4)^3} \\ 4(x^8 - x^4)^3 = x^4 \\ 4(x^8 - x^4)^3 - x^4 = 0$$
E se isso não te provar que é ridículo isso ser solução geral pra equação, sugiro trocar o curso pra bacharel em matemática. De qualquer forma, um breve resumo informal do porque não é solução geral é que x a oitava vai ser sempre positivo então mesmo se por algum motivo mágico os x a quarta se anulassem, sobraria x^24 (x^8 garantidamente vai se multiplicar 3 vezes) positivo que, pra dar 0, só sendo x = 0.
c) Outro horrível
$$y = \sqrt{x^8 - x^4} \\ \dfrac{d}{dx} [x^8 - x^4]^{1/2} = \dfrac{2(x^8 - x^4)^2 + x^4}{x(x^8 - x^4)^{3/2}} \\ \dfrac{1}{2} (x^8 - x^4)^{-1/2} \times (8x^7 - 4x^3) = \dfrac{2(x^8 - x^4)^2 + x^4}{x(x^8 - x^4)^{3/2}} \\ \dfrac{1}{2} (8x^7 - 4x^3) = \dfrac{2(x^8 - x^4)^2 + x^4}{x(x^8 - x^4)^2} \\ (4x^7 - 2x^3) = 2 + \dfrac{x^3}{(x^8 - x^4)^2} \\ \dfrac{(4x^7 - 2x^3)}{(x^8 - x^4)} = \dfrac{2(x^8 - x^4)^2 + x^3}{(x^8 - x^4)^3} \\ \dfrac{4}{x} - \dfrac{2}{x} = \dfrac{2(x^8 - x^4)^2 + x^3}{(x^8 - x^4)^3} \\ 4 - 2 = \dfrac{2x(x^8 - x^4)^2 + x^4}{(x^8 - x^4)^3} = 2 \\ \dfrac{2x}{(x^8 - x^4)} + \dfrac{x^4}{(x^8 - x^4)^3} = 2 \\
2x + x^4/(x^8 - x^4)^2 = 2x^8 - 2x^4 \\
x^4/(x^8 - x^4)^2 = 2x^8 - 2x^4 - 2x \\
x^4 = (2x^8 - 2x^4 - 2x)(x^{16} - 2x^{12} + x^8) \\
x^4 = (2x^{24} - 4x^{20} + 2x^{16} - 2x^{20} + 4x^{16} - 2x^{12} - 2x^{17} + 4x^{13} - 2x^9) \\
x^4 = 2x^{24} - 6x^{20} - 2x^{17} + 6x^{16} + 4x^{13} - 2x^9$$

Desnecessário dizer, a igualdade é falsa, logo não é solução.
d) A esperança
$$y = \sqrt[4]{(x^8 - x^4)} \\ \dfrac{d}{dx} [(x^8 - x^4)^{1/4}] = \dfrac{2(x^8 - x^4) + x^4}{x(x^8 - x^4)^{3/4}} = \dfrac{1}{4} (x^8 - x^4)^{-3/4} (8x^7 - 4x^3) \\ 2x^7 - 4x^3 = \dfrac{2(x^8 - x^4) + x^4}{x} = 2x^7 - 2x^3 + x^3 = 2x^7 - 4x^3$$
E se as igualdades batem, a letra d é solução da equação. Ufa?

5. Quais das seguintes funções são soluções da equação diferencial y'' - xy' + y = 0?
a) y = x²
Felizmente muito mais fácil, né. Vejamos:
(x²)'' - x(x²)' + x² = (2x)' - x(2x) + x² = 2 - 2x² + x² = 2 - x²
O que não configura solução da equação.
b) y = x
(x)'' - x(x)' + x = (1)' - x(1) + x = 0 - x + x = 0
Logo, é solução.
c) y = 1 - x²
(1 - x²)'' - x(1 - x²)' + (1 - x²) = (-2x)' - x(-2x)' + 1 - x² = -2 + 2x² + 1 - x² = -1 + x²

O que não é solução.
d) y = 2x² - 2
(2x² - 2)'' - x(2x² - 2)' + (2x² - 2) = (4x)' - x(4x) + 2x² - 2 = 4 - 4x² + 2x² - 2 = 2 - 2x²

O que não é solução.
e) y = 0
(0)'' - x(0) + 0 = 0 - 0 + 0 = 0

É solução.

6. Solução de A.
$$A = x'' + 4x' + 4x = e^t$$

Essa é chatíssima, não recomendo. Mas é fácil, veja:
$$a) x = e^t \\ (e^t)'' + 4(e^t)' + 4e^t = (e^t)' + 4e^t + 4e^t = e^t + 4e^t + 4e^t = 9e^t \neq e^t$$
Logo, não é solução.
$$b) x = e^{2t} \\ (e^{2t})'' + 4(e^{2t})' + 4e^{2t} = 2(e^{2t})' + 8e^{2t} + 4e^{2t} = 4e^{2t} + 8e^{2t} + 4e^{2t} = 16e^{2t} \neq e^{2t}$$
Logo, não é solução.
$$c) x = e^{2t} + e^t \\ (e^{2t} + e^t)'' + 4(e^{2t} + e^t)' + 4(e^{2t} + e^t) = (2e^{2t} + e^t)' + 4(2e^{2t} + e^t) + 4e^{2t} + 4e^t = \\ (4e^{2t} + e^t) + 8e^{2t} + 4e^t + 4e^{2t} + 4e^t = 12e^{2t} + 9e^t \neq e^t$$
$$d) x = te^{2t} + e^t \\ (te^{2t} + e^t)'' + 4(te^{2t} + e^t)' + 4(te^{2t} + e^t) = \\ (e^{2t} + 2te^{2t} + e^t)' + 4(e^{2t} + 2te^{2t} + e^t) + 4te^{2t} + 4e^t = \\ (2e^{2t} + 2e^{2t} + 4te^{2t} + e^t) + 4e^{2t} + 8te^{2t} + 4e^t + 4te^{2t} + 4e^t = \\ 8e^{2t} + 16te^{2t} + 9e^t \neq e^t$$
Logo, não é solução.
$$e) x = e^{2t} + te^t \\ (e^{2t} + te^t)'' + 4(e^{2t} + te^t)' + 4(e^{2t} + te^t) = \\ (2e^{2t} + e^t + te^t)' + 4(2e^{2t} + e^t + te^t) + 4e^{2t} + 4te^t = \\ (4e^{2t} + e^t + e^t + te^t) + 8e^{2t} + 4e^t + 4te^t + 4e^{2t} + 4te^t = \\ 16e^{2t} + 6e^t + 9te^t \neq e^t$$
Logo, não é solução.

7. Determine c de modo que y(x) = c(1-x²) satisfaça a condição:
a) y(0) = 1

É só substituir y e x. y, no caso, por 1; x, no caso, por 0. Não tem segredo:
1 = c(1-0²) = c, logo c = 1
b) y(1) = 0
$$0 = c(1-1^2) \\ c = \dfrac{0}{0}$$

Logo c é indeterminado.

8. Especifique c1 e c2 de modo que y(x) = c1 sen(x) + c2 cos(x) satisfaça as condições indicadas. Determine se tais condições são condições iniciais ou de contorno.
a) y(0) = 1; y'(0) = 2

É só substituir, também.
$$y(0) = c_1 \sin{0} + c_2 \cos{0} = 1 \\ y(0) = c_2 = 1$$
$$y'(0) = c_1 \cos{0} - c_2 \sin{0} = 2 \\ y'(0) = c_1 = 2$$
Todos os valores de x usados são o mesmo, as condições são iniciais.
$$b) y(\dfrac{\pi}{2}) = 1; y'(\dfrac{\pi}{2}) = 2 \\ y = c_1 \sin{\dfrac{\pi}{2}} + c_2 \cos{\dfrac{\pi}{2}} = 1 \\ y = c_1 = 1 \\ y'(\dfrac{\pi}{2}) = c_1 \cos{\dfrac{\pi}{2}} - c_2 \sin{\dfrac{\pi}{2}} = 2 \\ y'(\dfrac{\pi}{2}) = -c_2 = 2 \\ c_2 = -2$$
Condições iniciais pelo mesmo motivo.
$$c) y(0) = 1, y(\dfrac{\pi}{2}) = 1$$
Pro y(0) é só você repetir o processo que fez na letra A. Pro y(pi/2) só repetir o que fez na B.
Aí você encontra que c2 = 1 e c1 = 1.
Quanto a condição: os valores de x diferem, então é condição de contorno.
$$d) y'(0) = 1; y'(\dfrac{\pi}{2}) = 1 \\ y'(0) = c_1 \cos{0} - c_2 \sin{0} = c_1 = 1 \\ y'(\dfrac{\pi}{2}) = c_1 \cos{\dfrac{\pi}{2}} - c_2 \sin{\dfrac{\pi}{2}} = -c_2 = 1 \\ c_2 = -1$$
Pelo mesmo motivo da anterior ele usa condições de contorno.

[Cálculo Numérico] Aula 6: sistemas de numeração

2012-08-19T12:50:00.000-07:00

Sistemas de Numeração

Surge com a necessidade de contar e tem seus primeiros sistemas relacionados com o homem primitivo.
As primeiras contagens são relacionadas aos dedos da mão (sistema fingers). Com a necessidade de aumento é incorporado o sistema toes (dedos do pé) que leva ao surgimento dos sistemas decimais e vigesimais.
Os sistemas são utilizados para descrever o que chamamos de Sinais (eventos que ocorrem em relação a um determinado tempo).

Sinal Analógico

São sinais contínuos no tempo, neste tipo de sinal ocorre a passagem de um evento de forma suave e podem assumir qualquer valor.

Sinal Digital

São sinais descontínuos com valores discretos (não assumem qualquer valor).

Exemplo: imagine um termômetro digital e um analógico. Enquanto o digital retorna valores discretos, pode-se, com a utilização de uma lupa, devolver valores suaves para um termômetro analógico.
Os fenômenos físicos ocorrem todos com sinal analógico.
Contudo, os sinais digitais apresentam algumas vantagens:

são mais fáceis de projetar;
facilidade em seu armazenamento;
menos suscetível a ruídos.

Pois pequenas variações de amplitude não afetam seu significado. Como na física os sinais são analógicos e na eletrônica o sinal muitas vezes é digital devemos trabalhar com essa transformação.

Descrição de Números

Todo sistema de numeração, seja ele digital ou analógico, segue um padrão de representação padrão. Este é dado por:
$$... + XB^y + ...$$
Onde:
X = dado do sistema de numeração
B = base do sistema de numeração
Y = posição do dígito com relação a vírgula
Obs.: y = 0 para o primeiro dígito à esquerda da vírgula, y = -1 para o primeiro dígito à direita da vírgula.
$$\epsilon = 323841.52_{10} \\ 3*10^5 + 2*10^4 + 3*10^3 + 8*10^2 + 4*10^1 + 1*10^0 + 5*10^{-1} + 2*10^{-2}$$
Quando um número não apresenta indicativo de base a direita, sua base é 10.

Sistema Binário

São representados pela base 2 e apresentam apenas dois dígitos: 1 e 0. O sistema binário possui uma correspondência direta entre os sinais analógicos e digitais, uma vez que sua conversão ocorre de forma direta.
Exemplo:
$$1001101_2 \\ 1*2^6 + 0*2^5 + 0*2^4 + 1*2^3 + 1*2^2 + 0*2^1 + 1*2^0 = \\ 64+0+0+8+4+0+1 = 77_{10}$$
Ocorre assim a conversão de binário em decimal. Como realizar o processo inverso?
Duas maneiras:
1. Vá dividindo por 2 e coletando os restos. Exemplo: para 23 na base 10:
23/2 = 11, resto 1
11/2 = 5, resto 1
5/2 = 2, resto 1
2/2 = 1, resto 0
Agora veja bem: você pega o resultado da primeira divisão, 1, e vai acrescentando ao número os restos em ordem inversa (ou seja, nesse caso, de baixo pra cima). 23 em binário fica 10111.

2. Para até 8 bits, faça a seguinte fila:
1 - 2 - 4 - 8 - 16 - 32 - 64 - 128 - 256
Ou seja, 2 elevado a um número de 0 a 8. Aí é só localizar o número entre esses e analisar qual soma é necessária pra chegar a ele. Peguemos 23 novamente:

ele está antes do 32 e depois do 16, então já convém pegar o 16;
16+8 = 24, passou de 23, não compensa pegar o 8. 16+4 = 20;
20+2 = 22;
22+1, 23.

Agora é só pegar novamente a sequência e inserir 1 abaixo dos números usados, e 0 abaixo dos não-usados:
1 - 2 - 4 - 8 - 16
1 - 1 - 1 - 0 - 1
Invertendo, temos 10111. O mesmo resultado usando a outra maneira. Essa é meio pior de explicar porque depende bastante da prática, mas é só ir pegado.

Números fracionários
$$1101.111_2$$
Agora você vai imaginar o seguinte: 2 elevado a -1 é 0.5, a -2 é 0.25, a -3 0.125, e a cada expoente n+1 você vai dividindo o resultado por 2 (lógico).
Mas a conversão em si de binário pra decimal é bem fácil. É só decompor os números normalmente:
$$1*2^3 + 1*2^2 + 0*2^1 + 1*2^0 + 1*2^{-1} + 1*2^{-2} + 1*2^{-3} \\ 8 + 4 + 0 + 1 + 0.5 + 0.25 + 0.125 = 13 + 0.875 = 13.875$$

E essa foi a aula de sexta. Boa tarde a todos.

[Equações Diferenciais Ordinárias] Aula 4: solução e classificação de equações diferenciais

2012-08-14T12:31:00.000-07:00

Bom, eu sei que tivemos uma infinidade de exercícios passados em uma lista, mas a entrega é segunda então a partir de terça-feira estará quentinho aqui. Mas deixarei em registro os exercícios que fizemos em sala, possivelmente ajudarão bastante.

1. A equação A é solução geral de y'' + 4y = 0?
A: $$y = c_1 \sin{2x} - c_2 \cos{2x}$$
Ok, simples pra caramba. Vamos pegar o y fornecido em A e substituir pelos y da equação passada, e ver se os lados se igualam.
$$(c_1 \sin{2x} - c_2 \cos{2x})'' + 4(c_1 \sin{2x} - c_2 \cos{2x}) = 0 \\ (2c_1 \cos{2x} + 2c_2 \sin{2x})' + 4c_1 \sin{2x} - 4c_2 \cos{2x} = 0 \\ (-4c_1 \sin{2x} + 4c_2 \cos{2x}) + 4c_1 \sin{2x} - 4c_2 \cos{2x} = 0$$
Como vemos, os senos se anulam e os cossenos também. Logo, a expressão é verdadeira e a equação A é solução de y'' + 4y = 0.

2. A equação y = x² - 1 é solução geral de B?
B: $$(y')^4 + y^2 = -1$$

(vale a pena lembrar que estou usando isso de equação A, B, etc. por conveniência de formatação apenas)
Vamos fazer a mesma coisa: jogar x² - 1 em tudo que é y e ver se os resultados batem.
$$((x^2 - 1)')^4 + (x^2 - 1)^2 = -1 \\ (2x)^4 + (x^4 - 2x^2 + 1) = -1 \\ 16x^4 + x^4 - 2x^2 + 1 = -1 \\ 17x^4 - 2x^2 + 2 = -1$$
Só de bater o olho já dá pra perceber que x² - 1 não é solução geral de B. Dá pra passar -1 pro outro lado como +1 e fazer uma resolução de equação do segundo grau, mas será pra uma solução particular com valores particulares, não uma solução geral. Exercício resolvido.

3. Determine a ordem, função incógnita e variável independente das equações diferenciais abaixo: (vou postar uma forma alternativa sempre que achar necessário, pra ficar mais fácil de enxergar)
a) $$y''' - 5xy' = e^x + 1 \\ \dfrac{d^3 y}{dx^3} - 5x \dfrac{dy}{dx} = e^x + 1$$
O termo mais derivado é derivado três vezes, logo terceira ordem. Veja que a função derivada é y com relação a x (dy em cima, dx embaixo), logo a função incógnita é y e a variável independente é x.
b) $$ty'' + t^2 y' - \sin{t} \sqrt{y} = t^2 - t + 1 \\ t \dfrac{d^2 y}{dt^2} + t^2 \dfrac{dy}{dt} - \sin{t} \sqrt{y} = t^2 - t + 1$$
O termo mais derivado é derivado duas vezes, logo segunda ordem. A função derivada é y com relação a t, então função incógnita y e variável independente t.
c) $$s^2 \dfrac{d^2 t}{ds^2} + st \dfrac{dt}{ds} = s$$
(não creio ser necessário a forma alternativa porque acredito ser muito mais visível essa forma que a outra, com a função incógnita em cima e a variável independente embaixo e tal, apesar de nem sempre ser assim)
O termo mais derivado é derivado duas vezes, segunda ordem. A função derivada é t com relação a s, função incógnita t e variável independente s.

E foi só isso mesmo. Já está pronta a lista aqui, e não pretendo demorar a deixá-la pronta no blog também, postarei o mais cedo possível pra ter bastante conteúdo pra estudar pras provas... Já que essa é uma das três? matérias com provas que podem ser bem pesadas.
Vale a pena lembrar também que estamos só vendo como matéria específica algo que vimos generalizado no cálculo III, então pra esse começo o conteúdo desse blog de Cálculo Diferencial e Integral III ajudará muito.

[Equações Diferenciais Ordinárias] Aula 3: revisão de integrais

2012-08-09T20:35:00.000-07:00

Integrais simples
a)
$$\int x^8 dx = (\dfrac{1}{9} x^9) + C = \dfrac{x^9}{9} + C$$
b)
$$\int \dfrac{1}{x^3} dx = \int x^{-3} dx = (\dfrac{1}{-2} x^{-2}) = \dfrac{-1}{2x^2} + C$$
c)
$$\int \sqrt[3]{x^2} dx = \int x^{2/3} dx = (\dfrac{3}{5} x^{5/3}) = \dfrac{3 \sqrt[3]{x^5}}{5} + C$$
d)
$$\int (5x^3 + 2 \cos{x}) dx = 5 \int x^3 dx + 2 \int \cos{x} dx = 5 (\dfrac{1}{4}) x^4 + 2 \sin{x} = \\ \dfrac{5}{4} x^4 + 2 \sin{x} + C$$
e)
$$\int (8t^2 + 6e^t + \dfrac{1}{t})dt = 8 \int t^2 dt + 6 \int e^t dt + \int \dfrac{1}{t} dt = 8 (\dfrac{1}{3} t^3) + 6e^t + \ln{t} = \\ \dfrac{8}{3} t^3 + 6e^t + \ln{t} + C$$
f)
$$\int \dfrac{(x-1)^2}{x^2} dx = \int \dfrac{(x^2 - 2x + 1)}{x^2} dx = \int \dfrac{x^2}{x^2} dx - 2 \int \dfrac{x}{x^2} dx + \int \dfrac{1}{x^2} dx = \\ \int 1 dx - 2 \int \dfrac{1}{x} dx + \int x^{-2} dx = x - 2 \ln{x} + (-x^{-1}) = x - 2 \ln{x} - \dfrac{1}{x} + C$$
g)
$$\int \dfrac{\sin{x}}{\tan{x}} dx = \int \dfrac{\sin{x}}{\dfrac{\sin{x}}{\cos{x}}} dx = \int \dfrac{\sin{x} \cos{x}}{\sin{x}} = \int \cos{x} dx = \sin{x} + C$$

Integrais por substituição
a)
$$\int (x^2 + 3x)^3 (2x + 3)dx \\ u = x^2 + 3x, du = (2x + 3)dx \\ \int u^3 du = \dfrac{1}{4} u^4 = \dfrac{1}{4} (x^2 + 3x)^4 + C$$
b)
$$\int \dfrac{5x}{(x^2 + 1)} dx \\ u = x^2 + 1, du = (2x)dx \\ \int \dfrac{5x}{u} \dfrac{du}{2x} = \dfrac{5}{2} \int \dfrac{1}{u} du = \dfrac{5}{2} \ln{|u|} = \dfrac{5}{2} \ln{|(x^2 + 1)|} + C$$
c)
$$\int \dfrac{x^2}{\sqrt[3]{1 - 2x^3}} dx = \int x^2 ((1 - 2x^3)^{-1/3} dx \\ u = 1 - 2x^3, du = (-6x^2)dx \\ \int x^2 u^{-1/3} \dfrac{du}{-6x^2} = \dfrac{-1}{6} \int u^{-1/3} du = \\ \dfrac{-1}{6} (\dfrac{3}{2} u^{2/3}) = - \dfrac{3}{12} (1 - 2x^3)^{2/3} = - \dfrac{1}{4} \sqrt[3]{(1 - 2x^3)^2} + C$$
d)
$$\int 2e^{-x} dx \\ u = -x, du = (-1)dx \\ 2 \int e^u \dfrac{du}{-1} = -2e^u = -2e^{-x} + C$$
e)
$$\int \sin{3x} dx \\ u = 3x, du = 3dx \\ \sin{u} \dfrac{du}{3} = - \dfrac{1}{3} \cos{3x} + C$$
f)
$$\int 2x \cos{x^2} dx \\ u = x^2, du = (2x)dx \\ \int \cos{u} du = \sin{u} = \sin{x^2} + C$$
g)
$$\int \sin{x} \cos{x} dx \\ u = \sin{x}, du = \cos{x} dx \\ \int u du = \dfrac{1}{2} u^2 = \dfrac{1}{2} \sin^2 {x} + C$$

Integrais por partes
a)
$$\int x \ln{x} dx \\ u = \ln{x}, du = \dfrac{1}{x}, dv = (x)dx, v = \dfrac{x^2}{2} \\ \int x \ln{x} dx = \dfrac{x^2}{2} \ln{x} - \int \dfrac{1}{2} x^2 \dfrac{1}{x} dx = \dfrac{x^2}{2} \ln{x} - \dfrac{1}{2} \int x dx = \dfrac{x^2}{2} \ln{x} - \dfrac{1}{2} (\dfrac{1}{2} x^2) = \\ \dfrac{x^2}{2} \ln{x} - \dfrac{x^2}{4} = \dfrac{x^2}{2} (\ln{x} - \dfrac{1}{2})$$
b)
$$\int x \sin{x} dx \\ u = x, du = 1, dv = \sin{x} dx, v = - \cos{x} \\ \int x \sin{x} dx = x (- \cos{x}) - \int (- \cos{x})*1 dx = \\ -x \cos{x} + \int \cos{x} dx = -x \cos{x} + \sin{x} + C$$
c)
$$\int xe^{2x} dx \\ u = x, du = 1, dv = e^{2x} dx, v = \dfrac{e^{2x}}{2} \\ \int xe^{2x} dx = x (\dfrac{e^{2x}}{2}) - \int \dfrac{e^{2x}}{2} dx = \dfrac{xe^{2x}}{2} - \dfrac{1}{2} \dfrac{e^{2x}}{2} = \dfrac{xe^{2x}}{2} - \dfrac{e^{2x}}{4} = \dfrac{e^{2x}}{2} (x - \dfrac{1}{4}) + C$$

d) Integral cíclica (vou sempre pegar os números de Euler como u e as trigonométricas como v)
$$\int e^{x} \cos{x} dx = e^{x} \sin{x} - \int e^{x} \sin{x} dx = e^{x} \sin{x} - [e^{x} (- \cos{x}) - \int (- \cos{x}) e^{x} dx] = \\ e^{x} \sin{x} - [-e^x \cos{x} + \int e^{x} \cos{x} dx] = e^{x} \sin{x} + e^x \cos{x} - \int e^{x} \cos{x} dx \\ 2 \int e^{x} \cos{x} dx = e^{x} \sin{x} + e^x \cos{x} \\ \int e^{x} \cos{x} dx = \dfrac{1}{2} (e^{x} \sin{x} + e^{x} \cos{x}) = \dfrac{e^x}{2} (\sin{x} + \cos{x}) + C$$
e)
$$\int x^2 e^{-x} dx \\ u = x^2, du = 2x, dv = e^{-x} dx, v = -e^{-x} \\ \int x^2 e^{-x} dx = x^2 (-e^{-x}) - \int (-e^{-x}) 2x dx = -x^2 e^{-x} + \int e^{-x} 2x dx \\ u = 2x, du = 2, dv = e^{-x} dx, v = -e^{-x} \\ -x^2 e^{-x} + [2x (-e^{-x}) - \int (-e^{-x}) 2 dx] = -x^2 e^{-x} + [-2xe^{-x} + 2 \int e^{-x} dx] = \\ -x^2 e^{-x} + [-2xe^{-x} - 2e^{-x}] = -x^2 e^{-x} - 2xe^{-x} - 2e^{-x} = \\ -e^{-x} (x^2 + 2x + 2)$$

É isso aí, galera. Acho que agora não tem mais como enrolar e vamos voltar a ver EDOs, mas só vendo pra saber. Boa noite a todos.

[Física Experimental] Aula 2: erros de medição

2012-08-07T20:22:00.002-07:00

Ok, hoje tivemos bastante matéria nova, então trarei aqui uma síntese do que nos foi passado e todos os "exercícios" (de certa forma foram mesmo exercícios) também. Claro, a parte de cálculos, o resto o professor disponibiliza tudo no slide bonitinho que é muito melhor que eu explicando do meu jeito enrolado.

As equações usadas:
1. Valor médio.
$$V_{m} = \dfrac{\sum_{i=1}^{n} x}{n}$$
2. Desvio absoluto.
$$\delta_a = |V_m - V_{ex}|$$

3. Valor medido.
$$V_{mdo} = |V_m \pm \delta_m|$$

4. Desvio médio.
$$\delta_{m} = \dfrac{\sum_{i=1}^{n} \delta_{ai}}{n}$$
5. Erro relativo.
$$\delta_R = \dfrac{\delta_a}{V_m} \times 100\%$$

Não chegamos a usar operações básicas com erros então deixo as fórmulas pra um momento mais apropriado. Todas essas estão nos slides, estou deixando essas porque é uma justificativa para o que usaremos a seguir.
Tivemos que montar uma tabela com os seguintes valores experimentais:

5.84
5.90
5.80
5.82
5.81
5.84
5.85
5.83
5.84
5.83

Ok? Todos os valores experimentais estão aqui, listados. Precisamos agora terminar de preencher a seguinte tabela:

Medida	$$V_{exp} (cm)$$	$$\delta_a (cm)$$	$$\delta_R (\%)$$
1	5.84
2	5.90
3	5.80
4	5.82
5	5.81
6	5.84
7	5.85
8	5.83
9	5.84
10	5.83
Média

Ok? Primeiro, calculemos o valor médio. Feito isso, teremos os valores para jogar na equação dos desvios absolutos. E tendo o desvio absoluto, podemos calcular o desvio médio, que é o que usaremos junto com o valor médio para descobrir aquele desvio relativo no final da tabela. Um leva ao outro, não tem erro - piadas infames de lado.
À média:
$$V_m = \dfrac{5.84 + 5.90 + 5.80 + 5.82 + 5.81 + 5.84 + 5.85 + 5.83 + 5.84 + 5.83}{10} = \dfrac{58.36}{10} \\ V_m = 5.836 \approx 5.84$$
Lembrando que aqui vamos sempre arredondar para ter duas casas decimais e seguir o padrão proposto até agora. Atualizando a tabela:

Medida	$$V_{exp} (cm)$$	$$\delta_a (cm)$$	$$\delta_R$$
1	5.84
2	5.90
3	5.80
4	5.82
5	5.81
6	5.84
7	5.85
8	5.83
9	5.84
10	5.83
Média	5.84

O próximo passo é descobrir os desvios absolutos. São dez desvios, mais o médio, que faremos com facilidade efetuando a soma de todos os desvios e dividindo por 10.
$$\delta_{a1} = |5.84 - 5.84| = 0.00 \\ \delta_{a2} = |5.84 - 5.90| = 0.06 \\ \delta_{a3} = |5.84 - 5.80| = 0.04 \\ \delta_{a4} = |5.84 - 5.82| = 0.02 \\ \delta_{a5} = |5.84 - 5.81| = 0.03 \\ \delta_{a6} = |5.84 - 5.84| = 0.00 \\ \delta_{a7} = |5.84 - 5.85| = 0.01 \\ \delta_{a8} = |5.84 - 5.83| = 0.01 \\ \delta_{a9} = |5.84 - 5.84| = 0.00 \\ \delta_{a10} = |5.84 - 5.83| = 0.01$$
Agora para o desvio médio:
$$\delta_m = \dfrac{0.00 + 0.06 + 0.04 + 0.02 + 0.03 + 0.00 + 0.01 + 0.01 + 0.00 + 0.01}{10} = \dfrac{0.18}{10} = \\ \delta_m = 0.018 \approx 0.02$$

Medida	$$V_{exp} (cm)$$	$$\delta_a (cm)$$
1	5.84	0.00
2	5.90	0.06
3	5.80	0.04
4	5.82	0.02
5	5.81	0.03
6	5.84	0.00
7	5.85	0.01
8	5.83	0.01
9	5.84	0.00
10	5.83	0.01
Média	5.84	0.02

E agora por fim os erros relativos:
$$\delta_{R1} = \delta_{R6} = \delta_{R9} = \dfrac{0}{5.84} \times 100\% = 0 \times 100\% = 0.00\% \\ \delta_{R2} = \dfrac{0.06}{5.84} \times 100\% \approx 0.0103 \times 100\% = 1.03\% \\ \delta_{R3} = \dfrac{0.04}{5.84} \approx 0.0068 \times 100\% \approx 0.68\% \\ \delta_{R4} = \delta_{Rm} = \dfrac{0.02}{5.84} \times 100\% \approx 0.0034 \times 100\% = 0.34\% \\ \delta_{R5} = \dfrac{0.03}{5.84} \times 100\% \approx 0.0051 \times 100\% = 0.51\% \\ \delta_{R7} = \delta_{R8} = \delta_{R10} = \dfrac{0.01}{5.84} \times 100\% \approx 0.0017 \times 100\% = 0.17\%$$
Colocando tudo na tabela:

Medida	$$V_{exp} (cm)$$	$$\delta_a (cm)$$	$$\delta_R$$
1	5.84	0.00	0.00%
2	5.90	0.06	1.03%
3	5.80	0.04	0.68%
4	5.82	0.02	0.34%
5	5.81	0.03	0.51%
6	5.84	0.00	0.00%
7	5.85	0.01	0.17%
8	5.83	0.01	0.17%
9	5.84	0.00	0.00%
10	5.83	0.01	0.17%
Média	5.84	0.02	0.34%

O resto era questão de arredondamento de algarismos significativos. Mais uma vez, coisa simples. Como não vou usar LaTeX e sim texto normal usarei o símbolo de igual mesmo, mas fica claro que estamos aproximando valores. É que será um recurso útil pra deixar em negrito os mais importantes e deixar catalogada aquela regrinha de arredondamento passada:
a) 2.75 = 2.8
b) 2.49 = 2.5
c) 3.95 = 4.0
d) 4.0501 = 4.1
e) 7.95002 = 8.0
f) 6.95 = 7.0
g) 7.849 = 7.8
h) 3.45 = 3.4
Os dois negritados nos levam àquela regra que o professor passou no quadro:
"Se o dígito a ser eliminado for 5, temos os seguintes casos:
- se o antecessor for par, continua par;
- se o antecessor for ímpar, vira par (ou seja, logicamente, é arredondado pra cima);"
E é isso aí. Não sei dizer se vou colocar muita coisa disso de física experimental porque será mais relatório, e relatório é pessoal, mas o que tiver de matéria mais suave assim é bom postar pra ter referência.

Boa noite a todos.

[Equações Diferenciais Ordinárias] Aula 2: revisão de derivadas

2012-08-06T20:03:00.000-07:00

Boa noite, galera. Como tivemos que resolver nossos primeiros exercícios, por mais fáceis que eles sejam, aqui apareço para deixar suas resoluções matemáticas (não creio que haja como ficar explicando passo-a-passo como se deriva, nem mesmo sei como) pra todo mundo.

Parte 1: derivadas simples
$$y = t^5 - t^4 + 6 \\ y' = (t^5)' - (t^4)' + (6)' = 5t^4 - 4t^3 + 0 = 5t^4 - 4t^3$$
$$y = 2t - \dfrac{1}{t^2} + \sin{t} \\ y' = (2t)' - (t^{-2})' + (\sin{t})' = 2 - (-2t^{-3}) + \cos{t} = 2 + \dfrac{2}{t^3} + \cos{t}$$
$$y = \dfrac{\sqrt{t}}{2} + \ln{t} \\ y' = (\dfrac{1}{2} t^{1/2})' + (\ln{t})' = \dfrac{1}{2} \dfrac{1}{2} t^{-1/2} + \dfrac{1}{|t|} = \dfrac{1}{4t^{1/2}} + \dfrac{1}{|t|} = \dfrac{1}{4 \sqrt{t}} + \dfrac{1}{|t|}$$
$$f(t) = e^t \sin{t} \\ f'(t) = (e^t)' \sin{t} + e^t (\sin{t})' = e^t \sin{t} + e^t \cos{t}$$
$$y = \dfrac{x^2 - 1}{x + 3} \\ y' = \dfrac{(x^2 - 1)' (x + 3) - (x^2 - 1) (x+3)'}{(x+3)^2} = \dfrac{(2x)(x + 3) - (x^2 - 1)(1)}{(x+3)^2} \\ = \dfrac{2x^2 + 6x - x^2 + 1}{(x+3)^2} = \dfrac{x^2 + 6x + 1}{(x+3)^2}$$
$$y = \tan{z} \\ y' = (\dfrac{\sin{z}}{\cos{z}})' = (\sin{z} \cos^{-1}{z})' = (\sin{z})' \cos^{-1}{z} + \sin{z} (\cos^{-1}{z})' = \\ \dfrac{\cos{z}}{\cos{z}} + \sin{z} [(-1)(\cos^{-2}{z}) \times \sin{z}] = \\ 1 - (\sin{z}) (\cos^{-2}{z}) (\sin{z}) = 1 - \dfrac{\sin^2 {z}}{\cos^2 {z}} = 1 - \tan^2 {z}$$

Parte 2: regra da cadeia (destaquei entre colchetes o que é a derivada "de fora", e o que é a derivada "de dentro")
$$f(x) = (x^2 - 5x + 1)^4 \\ f'(x) = [2x - 5][4(x^2 -5x + 1)^3] = (8x - 20)(x^2 - 5x + 1)^3$$
$$f(x) = \sqrt{x^2 - 4} \\ f'(x) = ((x^2 - 4)^{1/2})' = [2x][\dfrac{1}{2} (x^2 - 4)^{-1/2}] = x(x^2 - 4)^{-1/2} = \dfrac{x}{\sqrt{x^2 - 4}}$$
$$y = \dfrac{2}{(3x+7)^5} \\ y' = (2(3x+7)^{-5})' = 2[3][(-5)(3x+7)^{-6}] = \dfrac{-30}{(3x+7)^6}$$
$$y = \sin{3x} \\ y' = [3][\cos{3x}] = 3 \cos{3x}$$
$$y = e^{2x} \\ y' = [2][e^{2x}] = 2e^{2x}$$
$$y = \ln{2x} \\ y' = [2][\dfrac{1}{|x|}] = \dfrac{2}{|x|}$$
$$y = \ln{\dfrac{x}{x+1}} \\ y' = (\ln{x} - \ln{x+1})' = (\ln{x})' - (\ln{x+1})' = \dfrac{1}{x} - \dfrac{1}{x+1}$$

E esses foram todos os exercícios passados. Pode ajudar alguns a revisarem porque o professor não chegou a resolver todos os de regra da cadeia, e sei que lááááááá do primeiro semestre ainda tem gente que tem dúvida como funciona.
É uma resolução bem boba e puramente matemática, eu sei, pode não estar muito detalhada de fato. Sintam-se livres pra perguntar qualquer coisa e bem-vindos de volta.

Portfolio I Capítulo 22 - Física Geral II

2012-07-10T17:45:00.001-07:00

Capítulo 22 - Campo elétrico

1. Na figura, as linhas do campo elétrico do lado esquerdo têm uma separação duas vezes maior que as linhas do lado direito.
a) Se o módulo do campo elétrico no ponto A é 40N/C, qual é o módulo da força a que é submetido um próton no ponto A?
b) Qual é o módulo do campo elétrico no ponto B?

Esse exercício deve ser o único realmente fácil do capítulo 22, junto com o 5. Inclusive, esse é realmente fácil. Só precisamos saber de duas informações para resolver a A:
$$F = Eq \\
q_{p} = 1.6*10^{-19}C$$
No caso, que a força é igual ao campo elétrico vezes a carga inserida e que a carga de um próton é de 1.6 vezes 10 a -19 Coulombs positivo. Se temos o módulo do campo elétrico e temos a carga, resolver a força é bem fácil:
$$F = (40N/C) * 1.6*10^{-19}C = 64*10^{-19}N$$
Para a B, é só saber que, no campo elétrico, quanto maior a distância entre as linhas de campo, menor é seu módulo. Como a distância entre as linhas de campo no ponto B é o dobro da distância de no ponto A:
$$E_b = \dfrac{E_a}{2} = \dfrac{40N/C}{2} = 20N/C$$
Resolvido.

5. O núcleo de um átomo de plutônio-239 contém 94 prótons. Suponha que o núcleo é uma esfera com 6.64fm de raio e que a carga dos prótons está distribuída uniformemente nessa esfera. Determine:
a) o módulo;
b) o sentido (para dentro ou para fora) do campo elétrico produzido pelos prótons na superfície do núcleo.
Primeiro, vamos analisar as informações que temos:
$$q = 94*1.6*10^{-19}C = 150.4*10^{-19}C \\
r = 6.64*10^{-15}m$$
No caso, a carga do plutônio é a carga contida em 94 prótons, e o raio da esfera é igual a 6.64 vezes 10 a -15 metros.

Para determinar o módulo do campo elétrico, temos a seguinte equação:
$$E = \dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q}{d^2} = 8.99*10^9 Nm^2/C^2 * \dfrac{150.4*10^{-19}C}{6.64*10^{-15}m}^2 = \dfrac{1352.096*10^{-10} Nm^2/C}{44.0896*10^{-30}m^2} \\
\approx 30.667*10^{20} N/C$$
Bem-respondida está a letra A.

Quanto a letra B, não tem o que equacionar. A carga do próton é positiva, então o campo elétrico produzido pelos prótons é divergente, ou seja, para fora.

6. Duas partículas são mantidas fixas sobre o eixo x: a partícula 1, de carga q1 no ponto x = 6.00cm, e a partícula q2 no ponto x = 21.0cm. Qual é o campo elétrico total a meio caminho entre as partículas, em termos dos vetores unitários?
$$q_1 = -2.00*10^{-7}C \\
q_2 = +2.00*10^{-7}C \\
x_1 = 6.00cm = 0.06m \\
x_2 = 21.0cm = 0.21m$$
Outro exercício bem fácil, mas pega na questão da distância, de maneira um pouco diferente daqueles exercícios do capítulo 21. A distância é a seguinte: você vai fixar o campo elétrico num ponto entre a primeira partícula e a segunda, e a distância será entre esse ponto e o x da partícula.
Primeiro, o ponto entre a primeira e a segunda:
$$x_E = \dfrac{0.21m + 0.06m}{2} = \dfrac{0.27m}{2} = 0.135m$$
Depois, as distâncias (que logicamente tem de ter módulo igual, mas sinais opostos):
$$d_1 = 0.135m - 0.06m = 0.075m\\
d_2 = 0.135m - 0.21m = -0.075m$$
Agora, se o campo elétrico total é a soma do campo com relação à primeira partícula e com relação à segunda:
$$E_T = E_1 + E_2 = K \dfrac{|q_1|}{d_1 ^2} + K \dfrac{|q_2|}{d_2 ^2} = K (\dfrac{|2*10^{-7}C|}{(0.075m)^2} + {|2*10^{-7}C|}{(-0.075m)^2}\\
= K (\dfrac{2*10^{-7}C}{0.005625m^2} + \dfrac{2*10^{-7}C}{0.005625m^2}) = K \dfrac{4*10^{-7}C}{0.005625m^2} \\
\approx 8.99*10^9 Nm^2/C^2 * 711.11*10^{-7}C/m^2 = 6392.8789*10^2 N/C$$

Se for pra expressar em termos dos vetores unitários, bem, o exercício especifica que as partículas são mantidas fixas no eixo x, aonde y = 0. Se o campo trabalha completamente no eixo x, o ângulo trabalhado é 0 - a equação no eixo x será multiplicada por cosseno de 0, e no eixo y pelo seno de 0.
Logo, ele é equivalente ao módulo em x e, por consequência, igual a 0 em y.

10. Na figura, as quatro partículas são mantidas fixas e têm cargas q1 = q2 = +5e, q3 = +3e e q4 = -12e. A distância d = 5.0 micrômetros. Qual é o módulo do campo elétrico no ponto P?

Esse exercício é facílimo de interpretar e dá pra responder sem fazer nenhuma conta, apenas usando deduções, mas se quiser usar as equações fica um pouquinho mais complicado. Não muito, de qualquer forma. Queremos saber o campo resultante dos 4 campos com relação ao ponto P, então o ponto de partida é esse:
$$E = E_1 + E_2 + E_3 + E_4$$
Uma informação que podemos ver pelo gráfico é que, com relação ao ponto P, a primeira carga está oposta à segunda, então vamos estabelecer a seguinte relação:
$$E_1 = -E_2$$
Também é bom lembrar que a distância entre P e q4 não é d, e sim 2d, como podemos notar pelo próprio gráfico. E que, com relação a P, um dos campos elétricos têm de ter carga negativa, já que a carga 3 é positiva e a carga 4 é negativa. Agora é só terminar o trabalho com a equação:
$$E = E_1 - E_1 + E_3 + E_4 = E_3 + E_4 = -K \dfrac{q_3}{d_3 ^2} + K \dfrac{q_4}{d_4 ^2}] = K (- \dfrac{3e}{d^2} + \dfrac{12e}{(2d)^2}) \\
K (-\dfrac{3e}{d^2} + \dfrac{12e}{4d^2}) = K * 0 = 0$$
Logo, o módulo do campo elétrico no ponto P é 0. Isso é facilmente explicado: todas as linhas de campo são na direção de 4, a única carga oposta, e como nenhuma delas é forte o suficiente para tocar o campo P "por acaso", o campo em P é nulo.

13. Na figura, as três partículas são mantidas fixas no lugar e têm cargas q1 = q2 = +e e q3 = +2e. A distância a = 6.00 micrômetros. Determine:
a) o módulo e b) a direção do campo elétrico no ponto P.

Outro exercício semelhante a um que maior parte patinou na prova, e de fato era chatinho, em especial pela questão da distância, mas nem tão problemático assim. Temos um ponto de partida extremamente semelhante ao exercício passado porque, se você perceber bem, com relação ao ponto P a primeira e a segunda carga se anulam. Exatamente porque as duas têm cargas iguais e o ponto P está entre elas.
Visto isso, notamos que o campo elétrico em P está integralmente relacionado ao campo elétrico da terceira carga. Logo:
$$E_P = E_3$$
Temos que a carga de 3 corresponde a dois elétrons, que resulta em:
$$q_3 = 2*1.6*10^{-19}C = 3.2*10^{-19}C$$
E o problema maior é a distância, mas podemos facilmente fazer essa relação notando que P está entre a primeira e a segunda carga, que estão em eixos opostos no plano cartesiano. O P está no meio da hipotenusa de um triângulo fechado pelas cargas q1 e q2. Se cada cateto é a distância a (6 micrômetros), temos que a hipotenusa é:
$$h = \sqrt{6^2 + 6^2} = \sqrt{36 + 36} = \sqrt{72} \mu m$$
Ótimo, temos a hipotenusa... Mas algo que pegou muita gente é que a distância entre q3 e o ponto P não é a hipotenusa. O ponto P está no centro dessa hipotenusa, ou seja:
$$d = \dfrac{\sqrt{72}}{2} \mu m = \dfrac{\sqrt{72}}{2}*10^{-6}m$$
(lembrando que não aproximei o resultado porque elevaremos ele ao quadrado e podemos descobrir algo mais exato)
Com essa distância, temos tudo o que precisamos para descobrir o módulo do campo elétrico:
$$E = K \dfrac{q}{d^2} = 8.99*10^9 Nm^2/C^2 * \dfrac{3.2*10^{-19}C}{(\dfrac{\sqrt{72}}{2} * 10^{-6}m)^2} = \dfrac{28.768*10^{-10} Nm^2/C}{\dfrac{72}{4} * 10^{-12}m^2} \\
\dfrac{28.768*4*10^2 N/C}{72} \approx 1.598*10^2 N/C \approx 160N/C$$
E essa é a resposta. Parece complicadinho, mas depois de deixar de estriar com a distância o exercício fica praticamente pronto.

Quanto ao ângulo (direção), é fácil de deduzir. Ainda no triângulo q1-q2-q3, temos um triângulo cujos dois catetos são iguais, logo têm dois ângulos de 45 graus e um de 90 - esse, no caso, é q3. Mas não estamos tratando desse triângulo, estamos tratando de um triângulo q1-q3-P ou q2-q3-P.
Não seria difícil fazer essa relação geometricamente falando, mas o exercício nos dá uma dica valiosíssima que não nos deixa perder tempo: o ângulo no ponto P do triângulo q2-q3-P é 90°. Se q2 continua tendo 45°, q3 nessa ocasião só pode ter 45° também.
Logo, a direção desse campo elétrico é de 45°.

16. A figura mostra um anel de plástico de raio R = 50.0cm. Duas pequenas contas coloridas estão sobre um anel: a conta 1, de carga +2.00 microC, que é mantida fixa na extremidade esquerda, e a conta 2, de carga +6.00 microC, que pode ser deslocada ao longo do anel. As duas contas, juntas, produzem um campo elétrico de módulo E no centro do anel. Determine:
a) um valor positivo do ângulo para o dado valor de E e b) um valor negativo.
$$E = 2.00*10^5 N/C \\
q_1 = +2.00 \mu C = +2.00*10^{-6} C \\
q_2 = +6.00 \mu C = +6.00*10^{-6} C \\
R = 50cm = 0.5m$$

Esse exercício é outro que precisa de umas sacadinhas de trigonometria pra resolver, mas não é nada realmente preocupante. Veja bem, temos um campo elétrico total que nos foi dado, mas sabemos que esse campo total é o que fica entre a primeira e a segunda conta. Ou seja, ele leva em conta tanto a primeira conta quanto a segunda.
$$E = E_1 + E_2$$
Mas lembremos que não podemos simplesmente somar os campos elétricos da primeira e da segunda conta para descobrir o que está no centro. Perceba que eles são grandezas vetoriais e, como tais, têm seu próprio sentido - o campo elétrico da primeira "aponta" para frente com relação ao centro, enquanto o segundo aponta para trás e para baixo numa direção desconhecida (que, inclusive, é o que queremos saber); sendo assim, seu sentido é negativo.
Vamos montar as equações de cada campo em cada direção pra exemplificar melhor:
$$E_1 = (\dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q_1}{R^2} \cos{0})\hat{i} + (\dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q_1}{R^2} \sin{0})\hat{j} = (\dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q_1}{R^2})\hat{i} \\
E_2 = -(\dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q_2}{R^2} \cos{\theta})\hat{i} - (\dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q_2}{R^2} \sin{\theta})\hat{j}$$
Perceba que já defini o ângulo do campo da primeira conta como 0, afinal, ela está no eixo x.

Agora o que faremos é uma soma disso:
$$E = E_1 + E_2 = ((\dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q_1}{R^2})\hat{i}) + (-(\dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q_2}{R^2} \cos{\theta})\hat{i} - (\dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q_2}{R^2} \sin{\theta})\hat{j})$$
Vamos separar os i dos j:
$$E = (\dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q_1}{R^2} - \dfrac{q_2}{R^2} \cos{\theta}) - (\dfrac{1}{4 \pi \varepsilon_0} * \dfrac{q_2}{R^2} \sin{\theta})\hat{j})$$
Agora temos que fazer uma análise trigonométrica do caso: queremos o ângulo theta, mas estamos com equações vetoriais e precisamos convertê-las para que fique escalar. Normalmente usaríamos o teorema de Pitágoras, mas o próprio ângulo é uma incógnita e... Inclusive a que queremos saber. Convém-nos então utilizar a lei dos cossenos (ou teorema de Pitágoras level 2 pros que gostam de associar):
$$a^2 = b^2 + c^2 - 2ab \cos{\theta}$$
Substituindo os valores por uma hipotenusa E, uma direção Ei e outra direção Ej:
$$E^2 = (K \dfrac{q_1}{R^2})^2 + (K \dfrac{q_2}{R^2})^2 - 2(K \dfrac{q_1}{R^2})(K \dfrac{q_2}{R^2}) \cos{\theta} \\
E^2 = K^2 \dfrac{q_1 ^2}{(R^2)^2} + K^2 \dfrac{q_2 ^2}{(R^2)^2} - 2 K^2 \dfrac{q_1 * q_2}{(R^2)^2} \cos{\theta} \\
E^2 = K^2 \dfrac{q_1 ^2}{R^4} + K^2 \dfrac{q_2 ^2}{R^4} - 2 K^2 \dfrac{q_1 * q_2}{R^4} \cos{\theta} \\
E^2 = \dfrac{K^2}{R^4} (q_1 ^2 + q_2 ^2 - 2 q_1 q_2 \cos{\theta})$$
Agora vamos isolar quem realmente interessa: quem está com o theta.
$$E + \dfrac{K^2}{R^4}*(2 q_1 q_2 \cos{\theta}) = \dfrac{K^2}{R^4} (q_1 ^2 + q_2 ^2) \\
\dfrac{K^2}{R^4} * (2 q_1 q_2 \cos{\theta}) = \dfrac{K^2}{R^4} (q_1 ^2 + q_2 ^2) - E^2 \\
2 q_1 q_2 \cos{\theta} = \dfrac{K^2 * R^4}{R^4 * K^2} (q_1 ^2 + q_2 ^2) - \dfrac{R^4}{K^2} E^2 \\
2 q_1 q_2 \cos{\theta} = q_1 ^2 + q_2 ^2 - \dfrac{R^4}{K^2} E^2 \\
\cos{\theta} = \dfrac{q_1 ^2 + q_2 ^2 - \dfrac{R^4}{K^2} E^2}{2 q_1 q_2} \\
\theta = \cos^{-1}{\dfrac{q_1 ^2 + q_2 ^2 - \dfrac{R^4}{K^2} E^2}{2 q_1 q_2}}$$
Agora substitua as variáveis dessa equação pelos valores dados pelo exercício na sua companheira calculadora científica. Ou no seu Matlab. A conta é extremamente extensa e cansativa, mas é algo que a calculadora pode fazer sem nenhum trabalho.
A resposta tem que dar aproximadamente 67.8°.

Quanto ao ângulo negativo... Ele será, bem, -67.8° mesmo.

19. A figura mostra um dipolo elétrico. Determine:
a) o módulo;
b) a orientação (em relação ao semi-eixo x positivo) do campo elétrico produzido pelo dipolo em um ponto P situado a uma distância r >> d.

É um exercício de raciocínio simples. O campo elétrico está no centro e, como tal, tem influência do campo gerado pela carga positiva e pela carga negativa. E, claro, como uma carga está divergindo e outra está convergindo, elas não se anulam e sim se acumulam num ponto só. Essa é a sacada do dipolo, que aliás é simplíssimo - o bom é que não é dado valor nenhum, então não há resultado constante a se considerar, é pra chegar numa equação de dipolo genérica.
Mas veja bem, como eu disse: as cargas não se acumulam, elas se juntam. Só com isso já dá pra imaginar o resultado do campo elétrico... Mas façamos uma análise mais minuciosa, pra não pecar em não entender o exercício como um todo.
Note que as duas cargas estão no mesmo eixo y, de ângulo desconhecido, mas possivelmente ignorando o eixo x. Como dito, elas se "complementam", então estão no mesmo sentido. Com isso concluímos que:
$$E_+ = E_{+y}\hat{j} = E \sin{\theta} \hat{j} \\ E_- = E_{-y}\hat{j} = E \sin{\theta} \hat{j} = E \hat{j} \\ E_{+x} = E_{-x} \\ |\vec{E}| = E_+ + E_- = E\hat{j} + E\hat{j} = 2E \sin{\theta}$$
Até aqui ok, vem agora outra parte: definir a equação completa pro campo. Talvez a pior.
Lembremos que a equação básica do campo elétrico é:
$$E = K \dfrac{q}{d^2}$$
Não temos o valor de q então ele ficará no resultado final, enquanto a distância d pode ser descoberta através da seguinte relação:

Ainda não defini a resultante porque se não dá aquele problema do cara querer chegar na resposta final direto e whatever. O fato é: veja que formamos triângulos retângulos entre o eixo x e y e essas linhas, logo essas linhas são hipotenusas. Podemos considerar o eixo x desse campo a variável r, enquanto o eixo y é representado por d/2. Pelo teorema de Pitágoras:
$$d_{E} = \sqrt{(\dfrac{d}{2})^2 + r^2}$$
Jogando na equação já final:
$$|\vec{E}| = 2(K \dfrac{q}{(\sqrt{(\dfrac{d}{2})^2 + r^2})^2}) \sin{\theta} = 2(K \dfrac{q}{(\dfrac{d}{2})^2 + r^2}) \sin{\theta}$$
Quanto ao seno de theta, vamos fazer relações trigonométricas básicas do ensino médio. O básico, seno é igual ao cateto oposto sobre a hipotenusa. De uma série de maneiras, uma delas sendo o fato de saber que d/2 está relacionado com o seno (quanto mais perto de 90° o ângulo, mais perto o vetor chega de d/2 absoluto), podemos definir o cateto oposto como d/2. A hipotenusa acabamos de expressar. Logo, no lugar de seno colocamos:
$$\sin{\theta} = \dfrac{\dfrac{d}{2}}{\sqrt{(\dfrac{d}{2})^2 + r^2}}$$
E jogamos na equação final:
$$|\vec{E}| = 2 \times (K \dfrac{q}{(\dfrac{d}{2})^2 + r^2}) \times \dfrac{\dfrac{d}{2}}{\sqrt{(\dfrac{d}{2})^2 + r^2}} = \dfrac{2Kqd}{2 \sqrt{((\dfrac{d}{2})^2 + r^2)^3}} = \dfrac{Kqd}{ \sqrt{((\dfrac{d}{2})^2 + r^2)^3}}$$
Sabemos que r >> d, logo, acrescentar qualquer valor de d para r se torna razoavelmente insignificante, então podemos simplesmente descartá-lo e deixar a equação dessa maneira:
$$|\vec{E}| = \dfrac{Kqd}{(\sqrt{r^2})^3} = \dfrac{Kqd}{r^3}$$
E está respondido o primeiro.

O segundo acabamos considerando a resposta no meio do raciocínio do primeiro. Em algum momento, dissemos da relação de d/2 com o seno e concluímos que ele tem 90°. Lembrando que as linhas vão em sentido à carga negativa, então o ângulo correto é -90°.

25. Na figura, duas barras curvas de plástico, uma de carga +q e outra de carga -q, formam uma circunferência de raio R = 8.50cm no plano xy. O eixo x passa pelos dois pontos de ligação entre os arcos, e a carga está distribuída uniformemente nos dois arcos. Se q = 15.0 picoC, determine:
a) o módulo;
b) a orientação (em relação ao semi-eixo x positivo) do campo elétrico E no ponto P, situado no centro da circunferência.

Galera, eu tinha uma resposta quase certa mas fracassei num detalhe que compromete toda a conta... No caso, um pi³ no lugar de um pi² no divisor. De qualquer forma, a expressão usada por uma resolução automática é decompor a carga dessa maneira:

$$q = \lambda \pi R$$

E jogar na equação do campo elétrico:

$$|\vec{E}| = E_+ + E_- = 2E \sin{\theta}$$

(até aí, igualzinho ao exercício anterior)

$$|\vec{E}| = 2 \dfrac{1}{4 \pi \varepsilon_0} \times \dfrac{\lambda \pi R}{\pi R^2} \sin{\theta} = \dfrac{\lambda \pi R}{2R^2 \pi^2 \varepsilon_0} (\sin{-90} - \sin{90}) = \\ \dfrac{\lambda}{2R \pi \varepsilon_0} \times (-2) = -\dfrac{q}{\pi^2 \varepsilon_0 R^2}$$

Daí jogando os valores:

$$|\vec{E}| = -\dfrac{15 \times 10^{-12}}{\pi^2 \varepsilon_0 (0.085)^2} = -\dfrac{15 \times 10^{-12}C}{8.85 \times 10^{-12} \times \pi^2 \times 0.007225} \\ = -\dfrac{15}{0.06394125 \pi^2} \approx \dfrac{234.59}{\pi^2} \approx 23.77N/C$$

E é isso aí. Deve ter reconhecido um problema no piR², ou está certo e eu estou confuso na matemática básica, sei lá. Mas o resultado é esse.

Quanto ao ângulo, perceba que a carga negativa está na parte inferior do círculo. É a mesma associação do exercício passado, -90°. Ao menos isso eu entendi bem, rs. Perdão pelo exercício, galera. Próximo.

30. Uma barra fina não-condutora, com uma distribuição uniforme de carga positiva Q, tem a forma de um círculo de raio R. O eixo central do anel é o eixo z, com a origem no centro do anel. Determine o módulo do campo elétrico:

a) no ponto z = 0;

b) no ponto z tendendo a infinito;

c) em termos de R, para que valor positivo de z o módulo do campo é máximo?

d) se R = 2.00cm e Q = 4.00 microC, qual é o valor máximo do campo?

Ok, essa questão é chata... Acho que caiu na primeira prova, até. Aqui o esquema é o seguinte: podemos fazer uma análise por equações prontas que o livro do Halliday nos ensina. Podemos fazer aquela análise por integral também. Vou fazer o completo (aproveitar que o livro tem tudo rs) mas a parte importante mesmo é só a equação final.

Veja bem:

$$\cos{\theta} = \dfrac{z}{r} = \dfrac{z}{(z^2 + R^2)^{1/2}}$$

(o r menor virar aquilo é devido ao fato do r menor ser a hipotenusa de um triângulo retângulo formado pelos catetos z e R)

Sendo:

$$dE = K \dfrac{\lambda ds}{(z^2 + R^2)}$$

Juntamos dE com o cosseno:

$$dE \cos{\theta} = K \lambda \dfrac{z}{(z^2 + R^2)^{3/2}} ds$$

Integramos os dois lados (s sendo o arco do anel, indo de 0 a 2piR):

$$E = \int dE \cos{\theta} = \int K \lambda \dfrac{z}{(z^2 + R^2)^{3/2}} ds = K \lambda \dfrac{z}{(z^2 + R^2)^{3/2}} 2 \pi R$$

O lambda com o 2piR forma uma carga q, o resto permanece igual:

$$E = \dfrac{Kqz}{(z^2 + R^2)^{3/2}}$$

Essa é a equação fixa do anel. Para o exercício A, precisamos substituir z por 0.

$$E = \dfrac{Kq}{(0 + R^2)^{3/2}} \times 0 = 0$$

Para o exercício B, jogamos um limite de z tendendo a infinito após trabalhar a equação:

$$E = \lim_{z\to\infty} Kq \dfrac{z}{(z^2 + R^2)^{3/2}}$$

Dá pra escapar daqui por um monte de método, mas no blog é bom economizar espaço. Vamos então aplicar o teorema de L'Hopital (aquele do cálculo I, lembram?), derivar em cima e embaixo com relação a z pra facilitar as coisas:

$$E = \lim_{z\to\infty} Kq \dfrac{1}{2z \times \dfrac{3}{2} (z^2 + R^2)^{1/2}} = \lim_{z\to\infty} Kq \dfrac{1}{3z \times (z^2 + R^2)^{1/2}}$$

Jogue os limites, o resultado tende a 0. Grato ao arroba @armiquilino (Alexandre Miquilino) por me lembrar do L'Hopital. (começando as participações especiais no blog!)

Para o exercício C, queremos o valor de z para o campo máximo. A conclusão a se chegar é fácil, o único chato é fazer mesmo. Quando queremos o valor máximo de alguma coisa usando alguma variável, já nos vem de cara setar a derivada da função incógnita com relação à variável a 0. Assim, em prática:

$$\dfrac{df(x)}{dx} = 0$$

No caso:

$$\dfrac{df(z)}{dz} = \dfrac{d}{dz} [Kq \dfrac{z}{(z^2 + R^2)^{3/2}}] = Kq \dfrac{d}{dz} [z \times (z^2 + R^2)^{-3/2}] = \\ Kq [(1 \times (z^2 + R^2)^{-3/2}) + (z \times (2z \times \dfrac{-3}{2} (z^2 + R^2)^{-5/2}))] = \\ Kq [\dfrac{1}{(z^2 + R^2)^{3/2}} - \dfrac{3z^2}{(z^2 + R^2)^{5/2}}] = \\ Kq [\dfrac{(z^2 + R^2) - 3z^2}{(z^2 + R^2)^{5/2}}] = Kq \dfrac{R^2 - 2z^2}{(z^2 + R^2)^{5/2}} = 0$$

Agora queremos isolar o z, não é? Vamos fazer o seguinte então. Note que a expressão acima pode ser vista também como:

$$Kq \dfrac{R^2}{(z^2 + R^2)^{5/2}} - Kq \dfrac{2z^2}{(z^2 + R^2)^{5/2}} = 0$$

Passando a parte do -2z pro outro lado como positivo:

$$Kq \dfrac{2z^2}{(z^2 + R^2)^{5/2}} = Kq \dfrac{R^2}{(z^2 + R^2)^{5/2}}$$

Cortando as partes iguais dos dois lados:

$$2z^2 = R^2 \\ z^2 = \dfrac{R^2}{2} \\ z = \sqrt{\dfrac{R^2}{2}} = \dfrac{R}{\sqrt{2}} \approx 0.707R$$

Isso, claro, porque queremos o valor positivo de z. Poderia ser -0.707R também por causa da raiz.

Bem trabalhoso, não é?

Para o exercício D, temos dois valores dados, a carga C (4 microC) e o raio R (2cm). O valor de z é o máximo, logo vamos substituir z pelo valor descoberto no exercício passado. Só jogar tudo na equação do anel:

$$E = Kq\dfrac{z}{(z^2 + R^2)^{3/2}} = (8.99 \times 10^9 Nm^2/C^2) \times (4 \times 10^{-6} C) \times \dfrac{\dfrac{R}{\sqrt{2}}}{((\dfrac{R}{\sqrt{2}})^2 + (0.02m)^2)^{3/2}} = \\ 35.96 \times 10^3 Nm^2/C \times \dfrac{\dfrac{0.02m}{\sqrt{2}}}{((\dfrac{0.02m}{\sqrt{2}})^2 + 0.0004m^2)^{3/2}} = \\ 35.96 \times 10^3 Nm^2/C \times \dfrac{0.02m}{\sqrt{2} \times ((\dfrac{0.0004m^2}{2}) + 0.0004m^2)^{3/2}} = \\ 35.96 \times 10^3 Nm^2/C \times \dfrac{0.02m}{\sqrt{2} \times (0.0002m^2 + 0.0004m^2)^{3/2}} = \\ 35.96 \times 10^3 Nm^2/C \times \dfrac{0.02m}{\sqrt{2} \times (0.0006m^2)^{3/2}} \approx \\ 35.96 \times 10^3 Nm^2/C \times \dfrac{0.02m}{\sqrt{2} \times (1.47 \times 10^{-5})} = \\ 35.96 \times 10^3 Nm^2/C \dfrac{1360.83}{\sqrt{2} m^2} = \\ \dfrac{48935 \times 10^3 Nm^2/C}{\sqrt{2} m^2} \approx 34602.59 \times 10^3 N/C \approx 3.5 \times 10^7 N/C$$

Enfim está encerrado. Acho que esse 30 é o mais escroto do capítulo 22. Depois disso é tudo mais pacífico.

O resto eu não vou fazer por motivos de: já passou o semestre e eu nem tenho mais as coisas aqui, e o exame já foi e etc. Então fiquem com esses exercícios mesmo, talvez sejam úteis. :Ç

Queda Livre para Exame - Cálculo Diferencial e Integral III

2012-07-10T13:51:00.000-07:00

Ok, me ofereceram uma bolada de 150 dólares pra esquecer um pouco as férias e curtir cálculo III com a galera que ficou de exame, então aqui estou eu. O pedido foi queda livre, então peguei os... três exercícios do livro sobre queda livre pra resolver aqui detalhadamente pra vocês.
Como de costume, ressaltarei no final. Mas estudem bastante o primeiro, o segundo é deveras desnecessário do ponto de vista de uma prova de cálculo, e o terceiro será apenas no caso de ele querer usar queda livre como exam killer.

1. Deixa-se cair um corpo de 10 "slugs" de massa de uma altura de 1000 pés sem velocidade inicial. A resistência do ar é proporcional à velocidade do corpo. Se a velocidade limite é 320 pés/s, determine:
a) uma expressão para a velocidade no instante t;
b) o tempo necessário para o corpo atingir a velocidade de 160 pés/s.
Ok, galera. Queda livre é um negócio menos complicado do que parece. Vamos começar lá do princípio, puxar a equação pelo esqueleto dela, e aí começamos a trabalhar porque aí dá pra entender bem. A equação usada aqui é:

$$F_g = ma + Kv$$
Sendo Fg a força da gravidade, m a massa do objeto, a a aceleração, K uma pequena proporcionalidade da resistência do ar e v a velocidade. Podemos igualar F a mg, ficando assim:
$$mg = ma + Kv$$
Passamos m dividindo tudo:
$$g = a + K \dfrac{v}{m}$$
Ok, mais uma consideração: quando a resistência do ar é proporcional à velocidade do corpo, temos que aquela constante K é definida pela razão entre o peso (massa vezes aceleração da gravidade) e a velocidade inicial. Difícil de entender? Seguinte:
$$K = \dfrac{mg}{v_0}$$
Sempre quando o exercício disser "A resistência do ar é proporcional à velocidade do corpo".
Substituindo isso na equação temos:
$$g = a + \dfrac{mg}{v_0} \times \dfrac{v}{m} = a + \dfrac{gv}{v_0}$$
Sendo assim, a equação está quase pronta pra gente trabalhar nela. Falta uma única consideração antes de jogar os valores constantes: note que a aceleração a é a derivada da velocidade com relação ao tempo, logo:
$$\dfrac{dv}{dt} + \dfrac{gv}{v_0} = g$$
Agora sim temos uma equação diferencial bonita e gostosa de se trabalhar. Ao jogar os valores, só precisamos fazer uma consideração: no sistema de medidas utilizado no exercício (slugs, pés, etc.) a gravidade se dá por 32.2 pés/s², que o exercício nos dá a colher de chá pra arrendondar pra 32 pés/s². Logo, jogamos tudo o que o exercício nos dá:
$$\dfrac{dv}{dt} + \dfrac{32}{320} v = 32 \\ \dfrac{dv}{dt} + 0.1v = 32$$
Resolvendo a equação diferencial usando o método do fator integrante:
$$I(t, v) = e^{\int 0.1 dt} = e^{0.1t}$$
$$e^{0.1t} [\dfrac{dv}{dt} + 0.1v] = 32e^{0.1t} \\ \dfrac{d}{dt} [ve^{0.1t}] = 32e^{0.1t}$$
Integrando dos dois lados:
$$ve^{0.1t} = 32 \dfrac{e^{0.1t}}{0.1} + C = 320e^{0.1t} + C$$
Isolando v:
$$v = 320 + Ce^{-0.1t}$$

Agora vamos analisar as informações que o exercício nos dá pra determinar a constante C: ele nos diz que não há velocidade inicial, o que indica que a velocidade inicial é 0. Logo, quando t = 0, v = 0. Podemos usar isso na equação:
$$v(0) = 320 + Ce^0 = 0 \\ 320 + C = 0 \\ C = -320$$
Sendo assim:
$$v = 320 - 320e^{-0.1t}$$

Essa é a resposta da letra a, que pede uma expressão para a velocidade no instante t. Como vê, a única variável da equação é t, o que indica que chegamos no melhor estado possível pra equação.
Tendo essa equação, no entanto, fica fácil resolver a letra b, que pede o tempo necessário para o corpo atingir 160 pés/s. Só substituir v por 160 e isolar t. Veja:
$$V(t) = 320 - 320^{-0.1t} = 160 \\ 320e^{-0.1t} = 320 - 160 = 160 \\ e^{-0.1t} = 0.5 \\ -0.1t = \ln{0.5} \approx -0.693 \\ t = \dfrac{-0.693}{-0.1} = 6.93$$
Ou seja, o tempo necessário é de 6.93 segundos. Resolvido o primeiro exercício.

2. Lança-se um corpo de massa m verticalmente para cima com velocidade inicial v0. Supondo nula resistência do ar, determine:
a) a equação do movimento no sistema de coordenadas da figura;
b) uma expressão para a velocidade do corpo no instante t;
c) o instante tm em que o corpo atinge altura máxima;
d) uma expressão para a posição do corpo no instante t;
e) a altura máxima atingida pelo corpo.
Vamos reconsiderar a equação que elaboramos anteriormente:

$$\dfrac{dv}{dt} + \dfrac{Kv}{m} = g$$
Veja bem, o exercício nos dá uma ajudinha aqui: "supondo nula resistência do ar". Resistência do ar nula indica que a constante K é 0. E outra consideração relevante é o fato do corpo estar sendo lançado verticalmente para cima, sendo assim a gravidade se opõe a ele e numa análise vetorial mais complexa chegaríamos à inversão de sinais de g.
Mas olha, vamos apenas usar a lógica básica: no outro exercício o objeto tava descendo então g era positivo, nesse ele está subindo então g é negativo. Simples assim.
Logo:
$$\dfrac{dv}{dt} = -g$$

Essa é a resposta pra letra a. Para as outras vamos chegar através de integrais a todas as equações que nos jogaram no ensino médio. Veja bem, a letra b pede a velocidade do corpo no instante t. Se temos dv/dt, só integrar:
$$v = \int \dfrac{dv}{dt} dt = \int -g \\ v = -gt + C$$
Temos que a velocidade inicial (t = 0) é v0, então façamos essa substituição:
$$v(0) = C = v_0$$
Sendo assim:
$$v = v_0 - gt$$

Para determinar o instante tm aonde a altura é máxima, temos que usar aquele princípio básico de equação de segundo grau. Lembra? Máximos e mínimos. O valor máximo de uma parábola com concavidade pra baixo é descoberto ao igualar a derivada da função a 0. A função, no caso, é a posição. Ou seja:
$$\dfrac{ds}{dt} = 0$$
Sabemos pelas leis fundamentais da mecânica que a velocidade é a derivada da posição, logo usamos a equação que descobrimos anteriormente:
$$\dfrac{ds}{dt} = v_0 - gt \\ v_0 - gt = 0$$
Agora vamos fazer o seguinte: já que não temos valores constantes pra colocar aí, vamos isolar t e criar uma equação pra isso. É suficiente.
$$gt = v_0 \\ t = \dfrac{v_0}{g}$$

O outro exercício nos pede uma expressão para posição. Raciocínio lógico puro. Se a velocidade é a derivada da posição, a posição é a integral da velocidade, sendo assim vamos integrar aquela equação que achamos para a letra b.
$$S = \int (v_0 - gt) dt = v_0 t - g(\dfrac{1}{2} t^2) + C = v_0 t - \dfrac{gt^2}{2} + C$$
Temos que a posição inicial é 0. Logo, (t = 0, v = 0). Jogando os valores:
$$S(0) = 0v_0 - \dfrac{0g}{2} + C = 0 \\ C = 0$$
E a equação final é:
$$S = V_0 t - \dfrac{gt^2}{2}$$

E o último pede a altura máxima, não o tempo necessário pra ela. Mas bem, se temos o tempo necessário pra ela, é só substituir t pelo tm que descobrimos na letra c.
$$S_{max} = V_0 \dfrac{V_0}{g} - \dfrac{g}{2} (\dfrac{V_0}{g})^2 = \dfrac{V_0^2}{g} - \dfrac{gV_0^2}{2g^2} = \dfrac{V_0^2}{g} - \dfrac{V_0^2}{2g} = \dfrac{V_0^2}{2g}$$

E enfim está tudo resolvido. Parece meio confuso a parte de máximo e tal, mas duvido bastante que isso vá cair na prova, é só pra ter queda livre dissecada passo a passo.

3. Um corpo de 1 "slug" de massa é solto no espaço com velocidade inicial de 1 pé/s, e encontra uma resistência do ar dada exatamente por -8v². Determine a velocidade no instante t.
Parece o exercício mais fácil, mas é muito mais difícil que o primeiro também. Devido ao fato do corpo ser lançado no espaço, temos que estabelecer uma força com relação ao corpo, e não à gravidade. Assim:
$$F_c = F_g - 8v^2 \\ ma = mg - 8v^2 \\ m \dfrac{dv}{dt} = mg - 8v^2$$
Substituindo os valores que temos:
$$1 \dfrac{dv}{dt} = 1 \times 32 - 8v^2 = 32 - 8v^2$$
Ok, e agora? Fator integrante? Não. Não temos a equação no formato de fator integrante, o que vamos fazer é "brincar" com ela um pouquinho com aquele método distributivo que aprendemos lá no começo.
$$dv = (32 - 8v^2)dt \\ \dfrac{dv}{32 - 8v^2} = dt \\ \dfrac{dv}{32 - 8v^2} - dt = 0$$
Integrando toda a equação temos:
$$\int \dfrac{1}{32 - 8v^2} dv - \int dt = C$$
A segunda é fácil, a primeira vamos usar a ajuda das nossas boas calculadoras porque usa um método horrível. É muito duvidoso que caia algo do tipo na prova, então é só pra passar isso logo e ir pra parte que interessa mais: o resultado.
$$(\dfrac{1}{32} [\ln{|v + 2|} - \ln{|2 - v|}]) - t = C$$
Usando mais propriedades (lnA - lnB = lnA/lnB) pra simplificar o que queremos:
$$\dfrac{1}{32} [\ln{|\dfrac{2 + v}{2 - v}|}] = C + t \\ \ln{|\dfrac{2 + v}{2 - v}|} = 32C + 32t = K + 32t \\ \dfrac{2 + v}{2 - v} = e^{K + 32t}$$
Uma das respostas do livro chega nesse formato. A outra é v isolado bonitinho. Quem quiser isolar completamente v, tem essa saída brincando com álgebra até cansar (vou chamar convenientemente K + 32t de x).
A princípio, jogamos o 2 pro outro lado:
$$\dfrac{v}{2 - v} = e^x - \dfrac{2}{2 - v}$$
Depois multiplicamos os dois lados por 2-v pra deixar v sozinho:
$$v = e^x (2 - v) - \dfrac{2(2 - v)}{2 - v} = (2e^x) - ve^x - 2$$
Então passamos ve^x pro outro lado de novo, porque vamos fazer um jogo de fatoração bem interessante:
$$v + ve^x = 2e^x - 2$$
Fatorando:
$$v(1 + e^x) = 2e^x - 2$$
Isolando v:
$$v = \dfrac{2e^x - 2}{e^x + 1} = 2 \times \dfrac{e^{K + 32t} - 1}{e^{K + 32t} + 1}$$
Agora, claro, vamos descobrir o valor de K. Que aliás usaremos a seguinte propriedade:
$$e^{K + 32t} = e^K \times e^{32t} = xe^{32t}$$
Ou seja: descobriremos x. Muito mais fácil.
Sabemos que a velocidade inicial é 1, logo, quanto t = 0, v = 1. Então:
$$v(0) = \dfrac{2x e^0 - 2}{x e^0 + 1} = 1 \\ \dfrac{2x - 2}{x + 1} = 1 \\ 2x - 2 = x + 1 \\ 2x - x = 1 + 2 \\ x = 3$$
Sendo assim, as equações finais ficam, a seu critério, ou:
$$\dfrac{2 + v}{2 - v} = 3e^{32t}$$
Ou:
$$v = \dfrac{6e^{32t} - 2}{3e^{32t} + 1} = 2 \dfrac{(3e^{32t} - 1)}{(3e^{32t} + 1)}$$
E isso é tudo que o Bronson mais antigo tem sobre queda livre. Recomendo mesmo estudarem o primeiro se for cair novamente, os outros dois são muito mais complexos e duvido bastante que ele vá soltar na prova... Mas deixei aqui porque, se vocês entenderam os três, vocês dominaram a matéria e não passarão necessidade na prova de forma alguma.

Bons estudos, galera. Qualquer coisa é só me dar um toque.

Aula X - Circuitos Elétricos

2012-06-25T12:25:00.001-07:00

Números complexos (imaginários)
$$j = \sqrt{-1}$$

Forma retangular:
$$a + (b)j$$
Forma polar:
$$C \angle{\theta°}$$
Representação gráfica:

Transformação de retangular para polar

$$c^2 = a^2 + b^2 \\ \theta = \tan^{-1}{\dfrac{b}{a}}$$

Exemplo:
$$3 + 4j \\ C = \sqrt{8^2 + 4^2} = \sqrt{9 + 16} = \sqrt{25} = 5 \\ \theta = \tan^{-1}{\dfrac{4}{3}} \approx 53.1°$$
Na calculadora: aperte Pol, insira o valor real, vírgula, depois o valor imaginário, feche parênteses. Dando ok aparecerá o módulo C. Então apertando RCL e tangente, aparecerá o ângulo theta.

Transformação de polar para retangular
$$\sin{\theta} = \dfrac{b}{C} \\ b = C \sin{\theta}$$

$$\cos{\theta} = \dfrac{a}{C} \\ a = C \cos{\theta}$$
Na calculadora: aperte shift, Pol, insira o valor do módulo C, vírgula, ângulo theta, feche parênteses e dê ok. Aparecerá o valor de A. Apertando RCL e tangente, aparecerá o valor de B.

Operações matemáticas
$$A + a_1 + b_1 j \\ B = a_2 + b_2 j$$

Adição:
$$A+B = (a_1 + b_1 j) + (a_2 + b_2 j) = (a_1 + a_2) + (b_1 + b_2)j$$

Subtração:
$$A - B = (a_1 + b_1 j) - (a_2 + b_2 j) = (a_1 - a_2) + (b_1 - b_2)j$$

Multiplicação:
$$D = C_1 \angle{\theta_1} \\ E = C_2 \angle{\theta_2} \\ D \times E = C_1 \angle{\theta_1} \times C_2 \angle{\theta_2} = (C_1 \times C_2) \angle{\theta_1 + \theta_2}$$
Divisão:
$$\dfrac{D}{E} = \dfrac{C_1 \angle{\theta_1}}{C_2 \angle{\theta_2}} = \dfrac{C_1}{C_2} \angle{\theta_1 - \theta_2}$$

Exercícios:
Calcule:
a) $$\dfrac{(2+3j) \times 5\angle{30°}}{2 - 2j + 3 + 5j}$$

Vamos fazer da seguinte forma: é mais fácil multiplicar/dividir usando forma polar e mais fácil somar/subtrair usando forma retangular, então sempre converteremos antes de fazer operações por aí. Lá em cima tem uma multiplicação, então converteremos 2+3j para polar:
$$C_1 = \sqrt{2^2 + 3^2} = \sqrt{4 + 9} = \sqrt{13} \approx 3.605 \\ \theta_1 = \tan^{-1}{\dfrac{3}{2}} \approx 56.31°$$
Lá embaixo está tudo convertido, então vamos continuar:
$$\dfrac{(3.605 \angle{56.31°}) \times (5 \angle{30°})}{2 - 2j + 3 + 5j} = \dfrac{18.025 \angle{86.31°}}{5 + 3j}$$
Convertamos o abaixo novamente:
$$C_2 = \sqrt{5^2 + 3^2} = \sqrt{25 + 9} = \sqrt{34} \approx 5.831 \\ \theta_2 = \tan^{-1}{\dfrac{3}{5}} \approx 30.964°$$
Voltando à conta:
$$\dfrac{18.025 \angle{86.31°}}{5.831 \angle{30.961°}} = 3.091 \angle{55.349°}$$

b) $$\dfrac{(2 + 5j)(3 - 2j)}{(1 + 2j)^2}$$
Ok, às conversões:
$$C_1 = \sqrt{2^2 + 5^2} = \sqrt{4 + 25} = \sqrt{29} \approx 5.385 \\ \theta_1 = \tan^{-1}{\dfrac{5}{2}} \approx 68.199°$$
$$C_2 = \sqrt{3^2 + (-2)^2} = \sqrt{9 + 4} = \sqrt{13} \approx 3.605 \\ \theta_2 = \tan^{-1}{\dfrac{-2}{3}} \approx -33.69007°$$
$$C_3 = \sqrt{1^2 + 2^2} = \sqrt{1 + 4} = \sqrt{5} \approx 2.236 \\ \theta_3 = \tan^{-1}{\dfrac{2}{1}} \approx 63.435°$$
E sim, fizemos C3 separado sem elevar ao quadrado porque fica melhor de enxergar assim. Veja adiante:
$$\dfrac{5.385 \angle{68.199°} \times 3.605 \angle{-33.6907°}}{2.236 \angle{63.435°} \times 2.236 \angle{63.435°}} = \dfrac{19.412925 \angle{34.5083°}}{4.999696 \angle{126.87°}} \approx 3.883 \angle{-92.3617}$$

Calcule:
$$A = 3 \angle{90°} \\ B = 2 \angle{0°} \\ C = 5 \angle{-90°}$$

a) $$\dfrac{A \times B}{A + B} = \dfrac{3 \angle{90°} \times 2 \angle{0°}}{3 \angle{90°} + 2 \angle{0°}} = \dfrac{6 \angle{90°}}{3j + 2}$$
Agora já vamos usar a calculadora pra converter de retangular pra polar, né. Já deu pra entender como funciona.
$$\dfrac{6 \angle{90°}}{3.605 \angle{56.31°}} = 1.664 \angle{33.69°}$$
b) $$\dfrac{C \times B}{C + B} + A \\ \dfrac{5 \angle{-90} \times 2 \angle{0°}}{-5j + 2} + 3j \approx \dfrac{10 \angle{-90}}{5.385 \angle{-68.199°}} + 3j \approx \\ 1.857 \angle{-21.801°} + 3j \approx 1.724 - 0.6897j + 3j = 1.72 + 2.3103j \approx 2.8803 \angle{53.333°}$$

Senóides e Fasores

$$V = A_m \times \sin{\omega t}$$

Fasores
$$V = A_m \times \sin{(\omega t \pm \theta)}$$
Domínio do tempo.
$$v = A \angle{\pm \theta°}$$
Domínio do fasor.

Não muito bem feito porque, ao contrário do previsto, eu não me lembrei de baixar software de plotar gráfico

$$v_1 = 2 \sin{(\omega t + \theta_1)} \\ v_2 = 1.5 \sin{(\omega t - \theta_2)}$$

Converta os seguintes fasores para o domínio do tempo, caso a frequência seja 60Hz.
a) $$v_1 = 10 \angle{20°}$$
Temos que ômega se torna 2pi vezes a frequência 60Hz. O resto é só usar o que está no exercício:
$$v_1 = 10 \sin{(\omega t + \theta)} = 10 \sin{((2\pi \times 60)t + 20°)} = 10 \sin{(120\pi t + 20°)}$$
b) $$v_2 = 12 \angle{-60°} \\ v_2 = 12 \sin{(120\pi t - 60°)}$$

Determine a defasagem entre as ondas e se i está adiantada ou atrasada em relação a v.
a) i = 2sen(377t - 30°), v = 5sen(377t - 40°)
A corrente i está adiantada com relação a v (-30 > -40), a defasagem (diferença de ângulo entre um e outro) é de 10°.

b) $$i = 3 \sin{(\omega t + 45°)} \\ v = 7 \sin{(\omega t - 35°)}$$
Mais uma vez, i está adiantado com relação a v (45 > -35). A defasagem, no caso, é de 80°.

Aula IX - Circuitos Elétricos

2012-06-25T10:35:00.000-07:00

Indutores

Obviamente com muito, mas muito mais espiras

Armazenam energia.

$$V(t) = L \dfrac{dI}{dt} \\ I(t) = \dfrac{1}{L} \int V(t) dt$$

L: indutância. Unidade: henry (H).

Formas de onda

Indutores em série:

$$L_{eq} = L_1 + L_2 + ... + L_n$$

Indutores em paralelo:

$$\dfrac{1}{L_{eq}} = \dfrac{1}{L_1} + \dfrac{1}{L_2} + ... + \dfrac{1}{L_n}$$

Exercícios:
A corrente em um indutor de 50H é mostrada no gráfico abaixo. Determine a tensão sobre o indutor:

Ok, esse exercício é extremamente semelhante ao do capacitor feito anteriormente. Vamos traçar uma série de equações de reta pra cada alteração que houver na corrente. A primeira equação é sempre padrão: não temos b porque a equação da reta começa em 0, e analisando o gráfico vemos que em 2 segundos a corrente atinge -10mA, sendo assim em 1 segundo ela atinge -5mA e podemos traçar a primeira equação como i = -5t. Sendo assim:
$$V_1 = L \dfrac{dI}{dt} = 50 \times \dfrac{d}{dt}[-5t] = 50 \times -5 = -250mV$$
Para a segunda, temos uma variação que já não começa em 0, logo b tem algum valor relevante e temos duas incógnitas na função. Sendo assim, façamos aquela análise de sistemas. Primeiro, pegaremos os valores das extremidades: (t, I): (2, -10), (4, 5). Agora, façamos um sistema linear substituindo I e t na equação I = at + b:
$$\begin{cases} -10 = 2a + b \\ 5 = 4a + b \end{cases}$$
Como de praxe, multipliquemos uma das equações por (-1) pra cortar os b. A segunda. Na soma, resta então -10 = 2a e -5 = -4a, que dá -2a = -15, assim a = 7.5. Substituímos a em uma das equações para descobrir b:
$$2 \times 7.5 + b = -10 \\ b = -10 - 15 = -25$$
E temos a equação final:
$$I = 7.5t - 25$$
Que colocaremos na equação da tensão:
$$V_2 = L \dfrac{dI}{dt} = 50 \times 7.5 = 375mV$$
Para a parte entre 4s e 6s, por favor, nem façamos análise. É uma reta constante, então sua derivada é 0, sendo assim a voltagem nesse ponto é 0. Resta então o ponto entre 6 e 8 segundos, que terá análise semelhante ao ponto entre 2 e 4.
Vejamos, os dois valores que podemos usar são: (t, I): (6, 5), (8,0). Sendo assim, o sistema fica:
$$\begin{cases} 5 = 6a + b \\ 0 = 8a + b \end{cases}$$
Multiplicando a equação de baixo por -1, eliminamos b, e sobra pra somar 5 = 6a e 0 = -8a, cujo resultado real é -2a = 5, e a = -2.5. Jogando a numa das equações para descobrir b:
$$8 \times -2.5 + b = 0 \\ b = 20 \\ I = -2.5t + 20$$
Assim, na equação da voltagem temos:
$$V_4 = L \dfrac{dI}{dt} = 50 \times (-25) = -125mV$$

E temos tudo aí. Recapitulando para plotar o gráfico:
Entre 0s e 2s, a tensão é de -250mV.
Entre 2s e 4s, a tensão é de 375mV.
Entre 4s e 6s, a tensão é 0.
Entre 6s e 8s, a tensão é -125mV.
E o gráfico fica:

Encontre a indutância equivalente.

Exercício extremamente simples, mas desenhei os nós pra não enganar ninguém:

Quando os circuitos "se cruzam" mas tem uma dessas curvas no cruzamento, é porque um fio não encosta no outro, assim aquilo no meio não é um nó; entende mais fácil quem já jogou Zelda ou alguns RPGs 2D. Bem, de qualquer forma, bem simples.
Analisando, temos AB em paralelo, BC em paralelo, e AB em série com BC. Jogando tudo nas equações:
$$L_{AB} = \dfrac{6 \times 12}{6+12} = \dfrac{72}{18} = 4H \\ L_{BC} = 2 || 2 = 1H \\ L_{eq} = L_{AB} + L_{BC} = 4 + 1 = 5H$$

Outro bem simples. Com os nós, vemos que os indutores de 2mH e 1mH que estão em série estão em curto. Sendo assim, sobra uma indutância equivalente entre os três de 6mH. E como sabemos, paralelo de valores iguais, você conta quantos dispositivos tem e divide o valor original por esse número. No caso, L = 6/3 = 2mH.

Viu só? Facílimo.

Aula VIII - Circuitos Elétricos

2012-06-25T08:15:00.000-07:00

Capacitores
- Função: armazenar carga

$$C = \dfrac{Q}{V} \\ Q = C \times V$$
Onde:
C = capacitância em F (Faraday)
Q = carga em C (Coulombs)
V = tensão em V (Volts)

Temos que:
$$i = \dfrac{\Delta Q}{\Delta t} \\ Q = C \times V \\ i = C \dfrac{\Delta V}{\Delta t} \\ i_C = C \dfrac{dV}{dt}$$
Integrando a equação da corrente, temos:
$$\int i_C dt = \int C \dfrac{dV}{dt} \\ \int i dt = C \times V \\ V_C = \dfrac{1}{C} \int i_C dt$$
Tendo:

Capacitância equivalente
Série:

$$V_{eq} = V_1 + V_2 + ... + V_n \\ \dfrac{1}{C_{eq}} \int i dt = \dfrac{1}{C_1} \int i dt + \dfrac{1}{C_2} \int i dt + ... + \dfrac{1}{C_n} \int i dt \\ \dfrac{1}{C_{eq}} \int i dt = (\dfrac{1}{C_1} + \dfrac{1}{C_2} + ... + \dfrac{1}{C_n}) \int i dt \\ \dfrac{1}{C_{eq}} = \dfrac{1}{C_1} + \dfrac{1}{C_2} + ... + \dfrac{1}{C_n}$$
Paralelo:

$$i_T = i_1 + i_2 + ... + i_n \\ C_{eq} \dfrac{dV}{dt} = C_1 \dfrac{dV}{dt} + C_2 \dfrac{dV}{dt} + ... + C_n \dfrac{dV}{dt} \\ C_{eq} \dfrac{dV}{dt} = (C_1 + C_2 + ... + C_n) \dfrac{dV}{dt} \\ C_{eq} = C_1 + C_2 + ... + C_n$$

Exercícios:
1. Encontre a capacitância equivalente do circuito.

Bem simples. Vamos resolver o que está em paralelo primeiro, depois estabeleceremos uma equação geral pra resistência em série.
$$C_1 = 1+1 = 2 \mu F \\ C_2 = 3+3 = 6 \mu F$$
E agora...
$$\dfrac{1}{C_{eq}} = \dfrac{1}{3} + \dfrac{1}{6} + \dfrac{1}{2} + \dfrac{1}{3} + \dfrac{1}{6} + \dfrac{1}{12} + \dfrac{1}{4} \\ \dfrac{1}{C_{eq}} = \dfrac{4+2+6+4+2+1+3}{12} = \dfrac{22}{12} = \dfrac{11}{6} \\ C_{eq} = \dfrac{1}{6} \approx 0.54 \mu F$$

2. Encontre a forma de onda da corrente em um capacitor de 1.4 microF quando temos a forma de onda da tensão abaixo:

Ok, seguinte situação: precisamos determinar a corrente em cada ponto a partir da equação do capacitor inteira. Sendo a equação a seguinte:
$$i = C \dfrac{dV}{dt}$$
Sim, a derivada. É que faremos tal procedimento: traçaremos a equação da reta do início ao fim de todos os segmentos da curva, e ela será o que colocaremos no lugar de V (veja bem, não dV/dt, apenas V). Siga o primeiro exemplo, a primeira reta:
(t, V): (0, 0), (1, 2)
Sendo uma reta, já fica bem claro que a equação da reta é 2t. Sendo assim, o gráfico da corrente entre t = 0 e t = 1 é:
$$i_{C1} = C \dfrac{dV}{dt} = C \dfrac{d}{dt} [2t] = C \times 2 = 1.4 \mu \times 2 = 2.8 \mu A$$
Fácil, não é? Se tem algum problema, ele é apenas traçar a equação da reta. E não é problema, de qualquer forma. Vamos pra segunda corrente, entre o tempo 1ms e o tempo 3ms - veja bem, a voltagem permanece 2 estático, então V = 2. Sendo assim:
$$i_{C2} = C \dfrac{d}{dt} [2] = C \times 0 = 0$$
Ótimo. Para a terceira voltagem, temos realmente que fazer os esquemas básicos pra equação da reta, mas bem simples. A equação da reta é a seguinte:
y = ax+b (y função, x variável)
Nesse caso, temos que:
V = at+b (V função, t variável)
Sendo assim, podemos jogar nos valores de V e de t o que temos no gráfico, que são os pontos. Os dois melhores exemplos sempre são os extremos, nesse caso (t, V): (3, 2), (5, -2). Formamos aí um sistema linear de equações:
$$\begin{cases} 2 = 3a + b \\ -2 = 5a + b \end{cases}$$
Multiplicamos a parte debaixo por (-1) pra poder cortar os b, e temos 2 = 3a e 2 = -5a pra somar. Assim, -2a = 4. Isolando, temos que a = -2. Substituímos a numa equação qualquer para descobrir b:
$$3 \times (-2) + b = 2 \\ b = 2 + 6 = 8$$
E temos que a equação final é V = -2t + 8. Colocamos isso na equação da corrente:
$$i_{C3} = C \dfrac{d}{dt} [-2t + 8] = 1.4 \mu \times (-2) = -2.8A$$
Ótimo, agora resta apenas uma equação da reta pra traçar. Ela é entre t = 5ms e t = 6ms, aonde há um aumento na voltagem. Assim como a última equação, vamos fazer uma análise dos pontos. (t, V): (5, -2), (6, 0).
E agora montar o sistema de equações:
$$\begin{cases} -2 = 5a + b \\ 0 = 6a + b \end{cases}$$
Aonde mais uma vez multiplicamos a segunda equação por (-1) pra anular os b, sobrando 5a = -2 e -6a = 0 na soma, e ficando -a = -2. Só inverter os sinais, temos a = 2. Então só colocar em alguma das duas equações pra descobrir b:
$$5 \times 2 + b = -2 \\ b = -2 - 10 = -12 \\ V = 2t - 12$$
Agora substituindo na equação da corrente:
$$i_{C4} = C \dfrac{d}{dt} [2t - 12] = 1.4\mu \times 2 = 2.8\mu A$$

Com isso dá pra plotar o gráfico da corrente com relação ao tempo. Mas pra ficar bem fácil e entendível, vamos fazer um sumário de tudo o que coletamos aqui.
Entre t = 0ms e t = 1ms, a corrente é 2.8 microA.
Entre t = 1ms e t = 3ms, a corrente é 0.
Entre t = 3ms e t = 5ms, a corrente é -2.8 microA.
Entre t = 5ms e t = 6ms, a corrente é 2.8 microA.

Assim, plotaremos o gráfico:

Viu? Até fácil, só trabalhosinho.

Sinto estar postando algo que talvez nem caia na prova, galera, mas tava vendo aqui e realmente não tem como não explicar capacitor e indutor antes de entrar na corrente alternada. Tem gente que tá entendendo errado as coisas até agora. Vou acelerar o que puder, ver se posto indutor ainda hoje e arriscar números complexos também.
No mais, bom dia a todos, e bons estudos.

A respeito do resto do terceiro semestre

2012-06-24T03:56:00.000-07:00

Bom dia,
Gostaria primeiramente de agradecer às pessoas que veem o blog, sejam eles da universidade ou externos. Embora pouco dinheiro (que não era o intuito, mesmo, os ads são só pra falar que tem), o número de views foi muito maior do que eu esperava desde a criação daqui. Pra falar a verdade, eu esperava só uns 5 ou 6 acessos desesperados na semana de provas. Tive quase 2000 views em 3 ou 4 meses de blog... Pra blogueiro profissional e blablabla não é muito, mas esse é um blog exclusivamente pela ciência, não esperava ninguém vendo.

Mas enfim, vamos tratar de coisa séria. Temos uma última semana de provas por aí, a respeito de circuitos e materiais elétricos, além dos trabalhos... E, muito importante pra muita gente, temos a semana de exames. Eu gostaria que me respondessem via Facebook ou comentários daqui - já que não me encontrarei na sociedade até quinta - se desejam que eu continue postando o conteúdo todo para exame. Vou trabalhar intensamente circuitos (talvez pulando algumas aulas desnecessárias, estamos em capacitores de corrente contínua ainda e na prova provavelmente só cai corrente alternada) essa semana pra ajudar quem está pendurado na prova, e gostaria de saber porque, embora quase perfeito com cálculo III, o blog está devendo listas detalhadas ou aulas bem descritivas de todas as outras matérias.
Espero que ainda assim tenha ajudado um pouco em todas quem ficou checando. E foi realmente satisfatório ver que estava ajudando gente com as equações diferenciais. Agora que eu aprendi a lidar bem com o blog, quem sabe próximo semestre consigo estar a par de todas as matérias igualmente.
Se tiver alguém muito afim de ajudar com alguma coisa - próximo semestre tem eletrônica digital, tem gente na sala que trabalha com isso -, ou até mesmo tomar uma matéria pra si, pode conversar comigo e me deixar um recado na internet porque se não eu esqueço rs. Ninguém precisa se sentir compelido a ajudar, é só se for de interesse pessoal mesmo.

Sério mesmo, foi um prazer ajudar e estamos aí pro que precisar. Basta pedir. Integrai-vos.

Aula X - Mecânica dos Sólidos

2012-06-13T06:53:00.000-07:00

Substitua o binário e a força mostrados na figura por uma força única equivalente aplicada à alavanca. Determine a distância do eixo ao ponto de aplicação dessa força equivalente.

I. Cálculo do momento produzido pela força F = 400N

Ok, o cálculo de momento é bem claro. Sendo r vetorial F a fórmula, e tendo a força facilmente decomposta em -400N na direção j (veja o ângulo de 90° pra baixo), só precisamos descobrir o r que nomearemos por questão de praticidade de OB.
Veja bem, temos uma hipotenusa OB de 0.3m e um cateto adjacente (OBx) de 0.15m. Para OBy, então, vamos fazer a relação de cateto oposto:
$$\sin{\theta} = \dfrac{\vec{OB}_y}{|\vec{OB}|} \\ \vec{OB}_y = |\vec{OB}| \sin{\theta} \\ \vec{OB}_y = 0.3m \sin{60°} \approx 0.26m$$
Logo:
$$\vec{OB} = (0.15m)\hat{i} + (0.26m)\hat{j}$$
Sendo o cálculo do momento o vetorial da distância r (no caso OB) com a força F, fazemos:
$$\vec{M} = \vec{OB} \times \vec{F} \\ = [(0.15m)\hat{i} + (0.26m)\hat{j}] \times [(-400N)\hat{j}] \\ = [(0.15m)(-400N)]\hat{i} \times \hat{j} + [(0.26m)(-400N)]\hat{j} \times \hat{j} = \\ (-60Nm)\hat{k}$$

II. Cálculo do binário
Agora, ao cálculo do binário. Perceba o esquema da figura abaixo:

Vamos pegar com relação ao centro. O disco tem 60mm e, bem, temos que multiplicá-lo com as forças, que são iguais, e somar os resultados. O que significa que podemos resumir a expressão em:
$$\vec{M} = \vec{r} \times \vec{F} = \\ 2[(0.06m)\hat{j} \times (200N)\hat{i}] = \\ 2 * (-12Nm)\hat{k} = (-24Nm)\hat{k}$$

III. Soma dos momentos
Bem, autoexplicativo:
$$\vec{M}_R = (-60Nm)\hat{k} + (-24Nm)\hat{k} = (-84Nm)\hat{k}$$

IV. Cálculo do novo vetor posição de F = 400N
Agora é a hora de recalcular o vetor posição com essa consideração da soma dos momentos. Veja:
$$\vec{M}_A = \vec{OC} \times \vec{F}$$
Sendo OC a nova posição. Vamos substituir os valores:
$$(-84Nm)\hat{k} = [(OC \cos{60})\hat{i} + (OC \sin{60})\hat{j}] \times [(-400N)\hat{j}] \\ = [(OC \cos{60})(-400N)\hat{i} \times \hat{j}] + [(OC \sin{60})(-400N)\hat{j} \times \hat{j}] \\ = (-400N)(OC \cos{60})\hat{k}$$
Isolando OC:
$$OC = \dfrac{(-84Nm)}{(400N)\cos{60}} \approx 0.48m$$
E está resolvido o exercício.

Equilíbrio de corpos rígidos

Um guindaste fixo tem uma massa de 1000Kg e é usado para suspender um caixote de 2400Kg. O guindaste é mantido na posição indicada na figura por um pino A e um basculante B. O centro de gravidade do guindaste está localizado em G. Determine as componentes das reações em A e B.

(grato ao professor Renato por já disponibilizar uma imagem do exercício com o diagrama, facilita bastante... E acho que pro aprendizado é melhor assim. Ainda mais pra mostrar a importância do centro de gravidade na estrutura)

I. Cálculo de B

A condição de momento pro equilíbrio do corpo é:

$$\Sigma \vec{M}_i = 0$$

Sendo:

$$M = Fd$$

Somamos então todos os momentos do exercício e igualamos a zero:

$$B(1.5) - (9.81 \times 10^3)(2) - (23.5 \times 10^3)(6) = 0 \\ 1.5B = 19.62 \times 10^3 + 141 \times 10^3 \\ B = \dfrac{160.62}{1.5} \times 10^3 = 107.1kN$$

II. Cálculo de Ax

Outra das condições para que haja equilíbrio é:

$$\Sigma \vec{F} = 0$$

Sendo assim:

$$\Sigma \vec{F}_x = 0 \\ \Sigma \vec{F}_y = 0$$

Podemos usar a equação de x que pode ser resolvida agora com o valor de B:

$$A_x + B = 0 \\ A_x + 107.1kN = 0 \\ A_x = -107.1kN$$

III. Cálculo de Ay

Fazemos o mesmo para Ay agora:

$$A_y - 9.81 \times 10^3 - 23.5 \times 10^3 = 0 \\ A_y = 33.3kN$$

IV. Módulo de A

Teorema de Pitágoras basicão:

$$A^2 = A_x ^2 + A_y ^2 \\ A = \sqrt{(-107.1 \times 10^3)^2 + (33.3 \times 10^3)^2} = 112.2kN$$

V. Direção de A

$$\tan{\theta} = \dfrac{A_y}{A_x} = -0.31 \\ \theta = \tan^{-1}{-0.31} = -17.3°$$

Lista VII - Cálculo Diferencial e Integral III

2012-06-12T05:33:00.000-07:00

21. Coloca-se um corpo à temperatura de 0°F em um quarto mantido à temperatura constante de 100°F. Se após 10 minutos a temperatura do corpo é 25°, determine:
a) o tempo necessário para a temperatura do corpo atingir 50° F;
b) a temperatura do corpo após 20 minutos.
Exercício bem simples, bem começo de lista mesmo. O esquema é lembrar já de cara da equação de Newton... Aquela mesmo:
$$T' + KT = KT_m$$
Vamos precisar trabalhar essa equação a partir do fator integrante pra descobrir uma expressão válida pra T. Veja bem. T' é dT/dt, K é uma constante que não temos, T é a temperatura e Tm é a temperatura do ambiente. Disso tudo, fomos informados que o quarto tem temperatura constante de 100°F, e nada mais; então é a única coisa que podemos substituir.
$$T' + KT = 100K$$
Agora vamos resolver a equação diferencial:
$$I(t, T) = e^{\int Kdt} = e^{Kt} \\ e^{Kt} [T' + KT] = e^{Kt} 100K \\ \dfrac{d}{dt} [e^{Kt}T] = e^{Kt} 100K \\ e^{Kt}T = \int e^{Kt} 100K dt = 100e^{Kt} + C \\ T = 100 + \dfrac{C}{e^{Kt}} = 100 + Ce^{-Kt}$$
(se não entendeu a resolução, é recomendável que revise a lista anterior de cálculo)
Ok, temos essa pré-equação aí, mas temos duas constantes não-resolvidas aí. C e K. Temos também duas informações prévias que podemos usar agora, mas não fazia sentido anteriormente. Primeiro, a temperatura do corpo ao ser colocada no ambiente é 0°F, o que indica que sua temperatura inicial T(0) é 0°F. Outra informação é que, em T(10) é 25°F. A informação mais conveniente no momento é T(0) para determinar C.
Veja:
$$T(0) = 100 + Ce^{-0K} = 0 \\ 100 + Ce^0 = 0 \\ 100 + C = 0 \\ C = -100$$
Substituindo na equação:
$$T = 100 - 100e^{-Kt}$$
E agora convém usar a segunda informação a nosso favor.
$$T(10) = 100 - 100e^{-10K} = 25 \\ -100e^{-10K} = 25-100 = -75 \\ e^{-10K} = \dfrac{-75}{-100} = 0.75 \\ -10K = \ln{0.75} \approx -0.288 \\ K = \dfrac{-0.288}{-10} \approx 0.029$$
Substituindo na equação:
$$T = 100 - 100e^{-0.029t}$$
Agora dá pra resolver o exercício.

O tempo necessário para a temperatura do corpo atingir 50°F. Temos que T(t) = 50°F. Queremos descobrir o valor de t e, colocando o resultado, vemos que t é a única variável que sobra, então podemos resolver.
$$T(t) = 100 - 100e^{-0.029t} = 50 \\ -100e^{-0.029t} = 50 - 100 = -50 \\ e^{-0.029t} = \dfrac{-50}{-100} = 0.5 \\ -0.029t = \ln{0.5} \approx -0.693 \\ t = \dfrac{-0.693}{-0.029} \approx 23.9$$
Ou seja, em aproximadamente 23.9 minutos.

A temperatura do corpo após 20 minutos. Mais fácil que o primeiro, bem mais fácil. Temos o valor de t, não temos a resposta, então é só jogar no lugar de t na equação e resolver.
$$T = 100 - 100e^{-0.029 \times 20} = 100 - \dfrac{100}{e^{0.58}} = 100 - 55.99 \approx 44.01°F$$
Resolvido. Próximo.

22. Um corpo com temperatura desconhecida é colocado em um refrigerador mantido à temperatura constante de 0°F. Se após 20 minutos a temperatura do corpo é 40°F e após 40 minutos é 20, determine a temperatura inicial do corpo.
Ok, temos a seguinte informação inicial pra equação: a temperatura do ambiente é 0°F. Nada mais, só isso. E é suficiente pra resolver bem a equação, quer ver?
$$T' + KT = 0K = 0 \\ I(t, T) = e^{Kt} \\ e^{Kt} [T' + KT] = 0 \\ \dfrac{d}{dt} [e^{Kt}T] = 0 \\ e^{Kt}T = C \\ T = Ce^{-Kt}$$
Vejamos agora. Temos duas informações, a primeira é: T(20) = 40, T(40) = 20. Se você analisar, nenhuma delas vai anular alguma variável pra gente trabalhar com a outra, então vamos ver:
$$Ce^{-20K} = 40 \\ \ln{Ce^{-20K}} = \ln{40} \\ \ln{C} + \ln{e^{-20K}} \approx 3.69 \\ \ln{C} - 20K = 3.69$$
$$Ce^{-40K} = 20 \\ ln{Ce^{-40K}} = \ln{20} \\ ln{C} + \ln{e^{-40K}} \approx 2.99 \\ \ln{C} - 40K = 2.99$$
Temos aqui um sistema linear com duas variáveis e duas expressões. Bem simples, podemos resolvê-lo simplesmente invertendo o sinal de uma das expressões pra cortar o lnC, assim:

$$\begin{cases} \ln{C} - 20K = 3.69 \\
-\ln{C} + 40K = -2.99 \end{cases} \\ 20K = 0.7 \\ K = \dfrac{0.7}{20} = 0.035$$
Ok, agora para C:
$$\ln{C} - 20 \times 0.035 = 3.69 \\ \ln{C} = 3.69 + 0.7 = 4.39 \\ C = e^{4.39} \approx 80.64$$
Logo, a equação final é:
$$T = 80.64e^{-0.035t}$$
Sendo assim, o T inicial T(0) é:
$$T(0) = 80.64e^{-0.035 \times 0} = 80.64e^0 = 80.64$$
E essa é a resposta do exercício.

23. Um corpo à temperatura de 50°F é colocado em um forno cuja temperatura é mantida a 150°F. Se após 10 minutos a temperatura do corpo é 75°F, determine o tempo necessário para o corpo atingir a temperatura de 100°F.
Temperatura ambiente Tm é 150°F. O resto não tem nada. Já basta pra iniciar a equação diferencial:
$$T' + KT = 150K \\ I(t, T) = e^{Kt} \\ e^{Kt} [T' + KT] = e^{Kt} 150K \\ \dfrac{d}{dt} [e^{Kt}T] = e^{Kt} 150K \\ e^{Kt}T = \int e^{Kt} 150K dt = 150e^{Kt} + C \\ T = 150 + Ce^{-Kt}$$
Agora outras informações: T(0) = 50, T(10) = 75. T(0) sempre é conveniente pra descobrir a constante C, então usaremos ele:
$$T(0) = 150 + Ce^0 = 50 \\ C = 50 - 150 = -100 \\ T = 150 - 100e^{-Kt}$$
$$T(10) = 150 - 100e^{-10K} = 75 \\ -100e^{-10K} = 75 - 150 = -75 \\ e^{-10K} = \dfrac{-75}{-100} = 0.75 \\ -10K = \ln{0.75} \approx -0.288 \\ K = \dfrac{-0.288}{-10} = 0.029 \\ T = 150 - 100e^{-0.029t}$$
Equação igualzinha a de outro exercício lá. Agora vamos resolver o exercício T(t) = 100. Veja bem:
$$T = 150 - 100e^{-0.029t} = 100 \\ -100e^{-0.029t} = 100 - 150 = -50 \\ e^{-0.029t} = \dfrac{-50}{-100} = 0.5 \\ -0.029t = \ln{0.5} \approx -0.693 \\ t = \dfrac{-0.693}{-0.029} \approx 23.9$$
Sim. É o mesmo resultado da 21. Próxima.

24. Certa substância radioativa decresce a uma taxa proporcional à quantidade de substância presente. Se, para uma quantidade inicial de uma substância de 100 miligramas, se observa um decréscimo de 5% após dois anos, determine:
a) uma expressão para a quantidade restante no tempo t;
b) o tempo necessário para uma redução de 10% da quantidade inicial.
Crescimento e decrescimento é sempre o mais fácil, porque sempre depende de pouca coisa. Veja, a equação diferencial pra decrescimento é: N' - KN = 0. Sendo a massa da substância N e K uma constante da equação. Vamos resolvê-la:
$$I(t, N) = e^{\int -Kdt} = e^{-Kt} \\ e^{-Kt} [N' - KN] = 0 \\ \dfrac{d}{dt} [e^{-Kt} N] = 0 \\ e^{-Kt} N = C \\ N = Ce^{Kt}$$
Ok, ok. Sabemos que N é a massa, e que a massa inicial N(0) é igual a 100mg. Substituímos os valores na equação:
$$100 = Ce^0 = C \\ N = 100e^{Kt}$$
E temos que essa massa diminui 5% após dois anos. Se temos 100mg e diminuímos 5%, temos 95mg (hurr). O tempo é 2 anos. Substituímos na equação:
$$95 = 100e^{2K} \\ e^{2K} = \dfrac{95}{100} = 0.95 \\ 2K = \ln{0.95} \approx -0.051 \\ K = \dfrac{-0.051}{2} = -0.026 \\ N = 100e^{-0.026t}$$
Precisamos agora determinar o tempo necessário pra uma redução de 10% da quantidade inicial. Retirando 10% de 100mg temos 90mg, então só substituir N:
$$90 = 100e^{-0.026t} \\ e^{-0.026t} = 0.9 \\ -0.026t = \ln{0.9} = -0.105 \\ t = \dfrac{-0.105}{-0.026} \approx 4.05$$
Ou seja, 4.05 anos. Essa é a resposta.

25. Certa substância radioativa decresce a uma taxa proporcional à quantidade presente. Se se observa que, após uma hora, houve uma redução de 10% da quantidade inicial da substância, determine a "meia-vida" (half life) da substância. (Sugestão: designe por N0 a quantidade inicial da substância. Não é preciso conhecer N0 explicitamente)
Nem precisamos trabalhar a equação pra saber que será a mesma coisa da anterior. Não tem valor do outro lado, não tem nem como mudar.
Então, vamos seguir a dica do exercício e determinar um N0 para o momento inicial N(0):
$$N_0 = Ce^{0} = C \\ N = N_0 e^{Kt}$$
Agora temos a informação de que há redução de 10% da quantidade inicial quando se passa uma hora. A massa restante, obviamente, é 90%. Representamos 90% da massa inicial por 0.9N0, assim:
$$0.9N_0 = N_0 e^{1K} \\ e^K = 0.9 \\ K = \ln{0.9} \approx -0.105 \\ N = N_0 e^{-0.105t}$$
A meia-vida da substância é quanto ela tem metade, 50%, da massa inicial. Sendo assim, 0.5N0:
$$0.5N_0 = N_0 e^{-0.105t} \\ e^{-0.105t} = 0.5 \\ -0.105t = \ln{0.5} = -0.693 \\ t = \dfrac{-0.693}{-0.105} \approx 6.6$$
Exato, 6.6 horas.

26. Sabe-se que a população de determinada cidade cresce a uma taxa proporcional ao número de habitantes existente. Se, após 10 anos, a população triplica, e após 20 anos é de 150.000 habitantes, determine a população inicial.
A equação é a mesma: N' - KN = 0. Como tal, vamos usar a solução inicial das outras questões de crescimento/decrescimento. Também podemos nos estender ainda mais e usar aquela solução com N0 do exercício passado. No caso:
$$N = N_0 e^{Kt}$$
Temos duas informações: após 10 anos, a população triplica, ou seja, é de 3N0. Após 20 anos, é de 150k habitantes. Perceba que se usarmos N0 do outro lado, podemos cortar os dois N0 e deixar apenas K de variável, assim podendo determiná-la... Então, bem, compensa mais usar a informação que temos com 10 anos.
$$3N_0 = N_0 e^{10K} \\ e^{10K} = 3 \\ 10K = \ln{3} \approx 1.099 \\ K = \dfrac{1.099}{10} \approx 0.109 \\ N = N_0 e^{0.109t}$$
Agora podemos usar a segunda informação para isolar e descobrir N0, veja:
$$150000 = N_0 e^{0.109 \times 20} = N_0 e^{2.18} \\ N_0 = \dfrac{150000}{e^{2.18}} \approx 16956$$
Ou seja, a população inicial era de aproximadamente 16956 habitantes.

Sendo assim, estão resolvidos todos os exercícios do livro recomendado a respeito de variação de temperatura e crescimento/decrescimento de dados diversos. Os próximos são os de circuitos, e os próximos de queda livre. Posso chegar a resolver os de escoamento e de trajetória ortogonal também, mas como não foi conteúdo em sala, não creio que vá cair na prova e não darei prioridade nenhuma a eles.
No mais, é isso. Bom dia a todos!

Aula IX - Mecânica dos Sólidos

2012-06-06T11:47:00.004-07:00

Exercício 2

O que o exercício nos pede é o momento com relação ao ponto B. Já está tudo expresso no gráfico dessa vez, mas é bom treinar para resolver a prova e saber resolver qualquer caso - a lista provavelmente ajudará bastante, mas isso é matéria de fácil aplicação em casos inúmeros da vida profissional e até mesmo com alguns probleminhas em casa.
De qualquer forma, como eu disse, já está tudo expresso no gráfico, é só fazer uma análise minuciosa de cada detalhe de lá - depois disso, vemos o que usaremos mesmo.

Primeiro, temos um vetor posição que mostra aonde está, numa parede, a força aplicada e aponta positivamente no eixo y e negativamente no eixo x. Com relação a B, a parede tem uma distância de 200mm (0.2m, considerando o sentido, -0.2m). A altura da força com relação a B é de 160mm, 0.16m positivo.
Temos um vetor posição representado assim:
$$\vec{r}_{AB} = -(0.2m)\hat{i} + (0.16m)\hat{j}$$
Quanto à força, é nos dado um módulo de 800N e um ângulo de 60° que se localiza no eixo x. Sendo assim, x está relacionado com o cosseno e y com o seno:
$$F_x = F \cos{60} = 800 \times \dfrac{1}{2} = 400N \\ F_y = F \sin{60} = 800 \times \dfrac{\sqrt{3}}{2} \approx 693N \\ \vec{F} = (400N)\hat{i} + (693N)\hat{j}$$

Lembrando que o exercício pede o momento com relação ao ponto B. A fórmula para o momento de um ponto é o produto vetorial entre a posição e a força, então temos que a resolução se dá simplesmente por:
$$\vec{M}_B = \vec{r}_{AB} \times \vec{F} = [(-0.2m)\hat{i} + (0.16m)\hat{j}] \times [(400N)\hat{i} + (693N)\hat{j}] \\ = [-0.2m \times 400N]\hat{i} \times \hat{i} + [-0.2m \times 693N]\hat{i} \times \hat{j} + [0.16m \times 400N]\hat{j} \times \hat{i} + [0.16m \times 693N]\hat{j} \times \hat{j} \\ 0 + (-138.6Nm)\hat{k} + (64Nm)-\hat{k} + 0 = (-138.6Nm - 64Nm)\hat{k} = -202.6Nm$$

E está resolvido.

Exercício 3

Mais uma vez, o gráfico já foi expresso previamente nesse. Temos um módulo de força de 135N sobre um ângulo alfa desconhecido.
Mas veja bem, para esse ângulo alfa, podemos usar uma simples relação do teorema de Tales. A figura te dá um triângulo supostamente análogo ALÉM da força e ainda, de quebra, dá um ângulo de 20°. Veja bem, por associação, podemos ver que o ângulo fecha com o eixo imaginário x 50°... Sabemos que a parte superior tem 20°, então alfa é o que resta, 30°.
Visto isso, prossigamos.

Ok, então decompomos os dois. A força e o vetor posição. O vetor posição podemos referenciar ao ângulo de 50, assim trabalhando x no cosseno e y no seno:
$$\vec{r} = [0.9m \cos{50}]\hat{i} + [0.9m \sin{50}]\hat{j} = (0.58m)\hat{i} + (0.69m)\hat{j}$$
E a força, com o ângulo que descobrimos:
$$\vec{F} = [135N \cos{30}]\hat{i} + [135N \sin{30}]\hat{j} = (116.91N)\hat{i} + (67.5N)\hat{j}$$

Agora lembremos da expressão pro momento, que é o produto vetorial entre a distância e a força. Assim como o professor, resolveremos essa em forma de matriz:
$$\vec{M} = \vec{r} \times \vec{F} = \begin{pmatrix}
\hat{i} & \hat{j} & \hat{k} \\
0.58m & 0.69m & 0 \\
116.91N & 67.5N & 0
\end{pmatrix} = ([0.69m \times 0] - [0 \times 67.5N])\hat{i} \\ - ([0.58m \times 0] - [0 \times 116.91N])\hat{j} + ([0.58m \times 67.5N] - [0.69m \times 116.91])\hat{k} \\ = [39.15Nm - 80.67Nm]\hat{k} = -41.52Nm$$