Política de uso de dados

Ao navegar por este site, você concorda com a política de uso de dados.

Política de uso de dados

Ao navegar por este site, você concorda com a política de uso de dados.

| | | |

Matemática Numérica I

3 Métodos para Sistemas Lineares 3.1 Método da Decomposição LU 3.3 Métodos de Jacobi e de Gauss-Seidel

Ajude a manter o site livre, gratuito e sem propagandas. Colabore!

3.2 Norma e Número de Condicionamento

Nesta seção, fazemos uma rápida discussão sobre normas de vetores e matrizes e sobre o condicionamento de uma matriz.

3.2.1 Norma $L^{2}$

Norma de Vetores

A norma $L^{2}$ de um dado vetor $\boldsymbol{v}=(v_{1},v_{2},\dotsc,v_{n})\in\mathbb{R}^{n}$ é definida por

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|\boldsymbol% {v}\|:=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}.

(3.136)

Lembrando que o produto interno de dois vetores $\boldsymbol{u},\boldsymbol{v}\in\mathbb{R}^{2}$ é definido por

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\boldsymbol{u% }\cdot\boldsymbol{v}:=u_{1}v_{1}+u_{2}v_{2}+\cdots+u_{n}v_{n},

(3.137)

temos que

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|\boldsymbol% {v}\|=\sqrt{\boldsymbol{v}\cdot\boldsymbol{v}}.

(3.138)

Exemplo 3.2.1.

Sejam os vetores

	$\displaystyle\boldsymbol{u}$	$\displaystyle=(1,-2,3,-4),$		(3.139)
	$\displaystyle\boldsymbol{v}$	$\displaystyle=(-1,2,0,1).$		(3.140)

a)

$\displaystyle\boldsymbol{u}\cdot\boldsymbol{v}$ $\displaystyle=u_{1}v_{1}+u_{2}v_{2}+\cdots+u_{n}v_{n}$ (3.141)

$\displaystyle=1\cdot(-1)+(-2)\cdot 2+3\cdot 0+(-4)\cdot 1$ (3.142)

$\displaystyle=-1-4+0-4$ (3.143)

$\displaystyle=-9.$ (3.144)

⬇

1import numpy as np

2u = np.array([1., -2., 3., -4.])

3v = np.array([-1., 2., 0., 1.])

4udv = np.dot(u,v)

5print(f'\nu.v = {udv}')
b)

$\displaystyle\|\boldsymbol{v}\|$ $\displaystyle=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$ (3.145)

$\displaystyle=\sqrt{(-1)^{2}+2^{2}+0^{2}+1^{2}}$ (3.146)

$\displaystyle=\sqrt{1+4+1}$ (3.147)

$\displaystyle=\sqrt{6}=2.4495\mathrm{e}\!+\!0.$ (3.148)

⬇

1import numpy as np

2import numpy.linalg as npla

3v = np.array([-1., 2., 0., 1.])

4norm_v = npla.norm(v)

5print(f'\n||v|| = {norm_v}')

Proposição 3.2.1.

(Propriedades da Norma para Vetores.) Dados os vetores $\boldsymbol{u},\boldsymbol{v}\in\mathbb{R}^{n}$ e um escalar $\lambda\in\mathbb{R}$ , temos:

a)

Positividade

$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \|\boldsymbol{v}\|\geq 0.$ (3.149)

$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \|\boldsymbol{v}\|=0\Leftrightarrow\boldsymbol{v}=\boldsymbol{0}.$ (3.150)
b)

Multiplicação por escalar

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|\lambda% \boldsymbol{v}\|=|\lambda|\cdot\|\boldsymbol{v}\|.$ (3.151)
c)

Desigualdade de Cauchy³²³²endnote: ³²Augustin-Louis Cauchy, 1789-1857, matemático francês. Fonte: Wikipédia: Augustin-Louis Cauchy.-Schwarz³³³³endnote: ³³Karl Hermann Amandus Schwarz, 1843-1921, matemático alemão. Fonte: Wikipédia: Hermann Amandus Schwarz.

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\boldsymbol{u% }\cdot\boldsymbol{v}\leq\|\boldsymbol{u}\|\cdot\|\boldsymbol{v}\|$ (3.152)
d)

Desigualdade triangular

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|\boldsymbol% {u}+\boldsymbol{v}\|\leq\|\boldsymbol{u}\|+\|\boldsymbol{v}\|$ (3.153)

Demonstração.

Sejam dados $\lambda\in\mathbb{R}$ e $\boldsymbol{u},\boldsymbol{v}\in\mathbb{R}^{n}$ .

a)

Positividade.

Observa-se diretamente que $v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}\geq 0$ . Então, como a raiz quadrada é uma função não-negativa, concluímos que

$\displaystyle\|\boldsymbol{v}\|$ $\displaystyle:=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$ (3.154)

$\displaystyle\geq 0$

No caso de $\boldsymbol{v}=\boldsymbol{0}$ , temos $v_{i}=0$ , $i=1,2,\dotsc,n$ , donde

$\displaystyle\|\boldsymbol{v}\|$ $\displaystyle=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$ (3.155)

$\displaystyle\|\boldsymbol{0}\|$ $\displaystyle=\sqrt{0}=0.$

Ou seja, se $\boldsymbol{v}=\boldsymbol{0}$ , então $\|\boldsymbol{v}\|=0$ . Agora, se $\boldsymbol{v}\neq\boldsymbol{0}$ , então tem-se $v_{i}\neq 0$ para algum $i=1,2,\dotsc,n$ . Logo, pela monotonicidade da função raiz quadrada, temos

$\displaystyle v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}>0$ (3.156)

$\displaystyle\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}>\sqrt{0}$ (3.157)

$\displaystyle\|\boldsymbol{v}\|>0.$ (3.158)

Portanto, concluímos que $\|\boldsymbol{v}\|=0$ se, e somente se, $\boldsymbol{v}=\boldsymbol{0}$ .
b)

Multiplicação por escalar.

Observamos que

$\lambda\boldsymbol{v}=(\lambda v_{1},\lambda v_{2},\dotsc,\lambda v_{n}).$ (3.159)

Então, segue por cálculo direto que

$\displaystyle\|\lambda\boldsymbol{v}\|$ $\displaystyle=\sqrt{(\lambda v_{1})^{2}+(\lambda v_{2})^{2}+\cdots+(\lambda v_% {n})^{2}}$ (3.160)

$\displaystyle=\sqrt{\lambda^{2}(v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2})}$ (3.161)

$\displaystyle=\sqrt{\lambda^{2}}\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$ (3.162)

$\displaystyle=|\lambda|\|\boldsymbol{v}\|.$ (3.163)

Desigualdade de Cauchy-Schwarz.

Sem perda de generalidade, se $\boldsymbol{v}=\boldsymbol{0}$ , então a desigualdade é imediatamente satisfeita. Suponhamos, agora, que $\boldsymbol{v}\neq\boldsymbol{0}$ . Para qualquer $\alpha\in\mathbb{R}$ , temos

$\displaystyle(\boldsymbol{u}+\alpha\boldsymbol{v})\cdot(\boldsymbol{u}+\alpha% \boldsymbol{v})$	$\displaystyle\geq 0$	(3.164)
$\displaystyle\boldsymbol{u}\cdot\boldsymbol{u}+2\alpha\boldsymbol{u}\cdot% \boldsymbol{v}+\alpha^{2}\boldsymbol{v}\cdot\boldsymbol{v}$	$\displaystyle\geq 0$	(3.165)
$\displaystyle\\|\boldsymbol{u}\\|^{2}+2\alpha\boldsymbol{u}\cdot\boldsymbol{v}+% \alpha^{2}\\|\boldsymbol{v}\\|^{2}$	$\displaystyle\geq 0$	(3.166)

O lado direito desta desigualdade é um polinômio quadrático em $\alpha$ . Para ser não-negativo para todo $\alpha$ , seu discriminante precisa ser não-positivo, i.e.

	$\displaystyle(2\boldsymbol{u}\cdot\boldsymbol{v})^{2}-4\\|\boldsymbol{u}\\|^{2}% \\|\boldsymbol{v}\\|^{2}\leq 0$		(3.167)
	$\displaystyle(\boldsymbol{u}\cdot\boldsymbol{v})^{2}\leq\\|\boldsymbol{u}\\|^{2}% \\|\boldsymbol{v}\\|^{2}$		(3.168)
	$\displaystyle\sqrt{(\boldsymbol{u}\cdot\boldsymbol{v})^{2}}\leq\sqrt{\\|% \boldsymbol{u}\\|^{2}\\|\boldsymbol{v}\\|^{2}}$		(3.169)
	$\displaystyle\|\boldsymbol{u}\cdot\boldsymbol{v}\|\leq\\|\boldsymbol{u}\\|\\|% \boldsymbol{v}\\|$		(3.170)

Donde, concluímos que

\boldsymbol{u}\cdot\boldsymbol{v}\leq\|\boldsymbol{u}\|\|\boldsymbol{v}\|.

(3.171)

d)

Desigualdade triangular

Consulte o Exercício 3.2.6 para a demonstração da desigualdade triangular.

∎

Exemplo 3.2.2.

Vamos verificar a Desigualdade Triangular (3.153) para os vetores

	$\displaystyle\boldsymbol{u}$	$\displaystyle=(1,-2,3,-4),$		(3.172)
	$\displaystyle\boldsymbol{v}$	$\displaystyle=(-1,2,0,1).$		(3.173)

De fato, temos

$\displaystyle\\|\boldsymbol{u}+\boldsymbol{v}\\|$	$\displaystyle=\sqrt{(u_{1}+v_{1})^{2}+(u_{2}+v_{2})^{2}+(u_{3}+v_{3})^{2}+(u_{% n}+v_{n})^{2}}$	(3.174)
	$\displaystyle=\sqrt{0^{2}+0^{2}+3^{2}+3^{2}}$	(3.175)
	$\displaystyle=\sqrt{18}$	(3.176)
	$\displaystyle=4.2426\mathrm{e}\!+\!0$	(3.177)

$\displaystyle\\|\boldsymbol{u}\\|+\\|\boldsymbol{v}\\|$	$\displaystyle=\sqrt{u_{1}^{2}+u_{2}^{2}+u_{3}^{2}+u_{4}^{2}}$
	$\displaystyle+\sqrt{v_{1}^{2}+v_{2}^{2}+v_{3}^{2}+v_{4}^{2}}$	(3.178)
	$\displaystyle=\sqrt{1^{2}+(-2)^{2}+3^{2}+(-4)^{2}}$
	$\displaystyle+\sqrt{(-1)^{2}+2^{2}+0^{2}+1^{2}}$	(3.179)
	$\displaystyle=\sqrt{30}+\sqrt{6}$	(3.180)
	$\displaystyle=7.9267\mathrm{e}\!+\!0$	(3.181)

Ou seja, temos que

\|\boldsymbol{u}+\boldsymbol{v}\|\leq\|\boldsymbol{u}\|+\|\boldsymbol{v}\|

(3.182)

como esperado.

⬇

1import numpy as np

2u = np.array([1., -2., 3., -4.])

3v = np.array([-1., 2., 0., 1.])

4nupv = npla.norm(u+v)

5nupnv = npla.norm(u) + npla.norm(v)

6print('\nu.v <= ||u||.||v||?')

7print(nupv <= nupnv)

Norma de Matrizes

A norma $L^{2}$ induzida de uma dada matriz real $A=[a_{ij}]_{i,j=1}^{n}$ é definida por

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|A\|_{2}:=% \sup_{\boldsymbol{x}\in\mathbb{R}^{n},\|\boldsymbol{x}\|=1}\|A\boldsymbol{x}\|.

(3.183)

Pode-se mostrar que³⁴³⁴endnote: ³⁴Consulte [7, Section 1.3-3] para informações sobre a demonstração.

\|A\|_{2}=\sqrt{\lambda_{max}(A^{T}A)},

(3.184)

onde $\lambda_{max}(A^{T}A):=\max\{|\lambda|;\leavevmode\nobreak\ \lambda\text{ é % autovalor de }A^{T}A\}$ .

Tendo em vista o grande custo computacional em se calcular a norma induzida, vamos trabalhar com a norma de Frobenius (ou norma Euclidiana³⁵³⁵endnote: ³⁵Euclides de Alexandria, 300 a.C., matemático grego. Fonte: Wikipédia: Euclides.)

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|A\|=\left(% \sum_{i=1}^{n}\sum_{j=1}^{n}|a_{i,j}|^{2}\right)^{\frac{1}{2}}

(3.185)

Esta é uma generalização da norma euclidiana para vetores e é equivalente a norma $L^{2}$ induzida, i.e.

\|A\|_{2}\leq\|A\|\leq\sqrt{n}\|A\|_{2}.

(3.186)

Exemplo 3.2.3.

A matriz

A=\begin{bmatrix}1&-1&2\\ -2&\pi&4\\ 7&-5&\sqrt{2}\end{bmatrix}

(3.187)

tem norma

\|A\|=10.5768\mathrm{e}\!+\!0.

(3.188)

⬇

1import numpy as np

2import numpy.linalg as npla

3A = np.array([[1., -1., 2.],

4 [-2., np.pi, 4.],

5 [7., -5., np.sqrt(2)]])

6norm_A = npla.norm(A)

7print(norm_A)

Proposição 3.2.2.

(Propriedades da Norma para Matrizes.) Dadas as matrizes reais $A, B$ $n\times n$ , um vetor $v\in\mathbb{R}^{n}$ e um escalar $\lambda$ , temos

a)

Positividade

$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \|A\|\geq 0$ (3.189)

$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \|A\|=0\Leftrightarrow A=\boldsymbol{0}$ (3.190)
b)

Multiplicação por escalar

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|\lambda A\|% =|\lambda|\cdot\|A\|$ (3.191)
c)

Desigualdade triangular

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|A+B\|\leq\|% A\|+\|B\|$ (3.192)
d)

Sub-multiplicatividade

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|AB\|\leq\|A% \|\|B\|$ (3.193)
e)

Norma da aplicação

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|Av\|\leq\|A% \|\|v\|$ (3.194)

Demonstração.

Consulte o Exercício 3.2.7. ∎

Exemplo 3.2.4.

Vamos ver exemplos das propriedades d) e e) da norma de matriz.

a)

Sub-multiplicatividade

$\displaystyle\|I_{n}^{2}\|=\|I_{n}\|=\sqrt{n}$ (3.195)

$\displaystyle\|I\|\cdot\|I\|=\|I\|^{2}=n$ (3.196)
b)

Norma da aplicação

$\displaystyle\|I_{n}\boldsymbol{v}\|=\|\boldsymbol{v}\|$ (3.197)

$\displaystyle\|I_{n}\|\|v\|=\sqrt{n}\|\boldsymbol{v}\|$ (3.198)

3.2.2 Número de Condicionamento

Em revisão

O número de condicionamento de uma matriz é uma medida referente a propagação de erros que ocorre da sua aplicação. Mais especificamente, assumamos que seja dada uma matriz invertível $A=[a_{ij}]_{i,j=1}^{n,n}$ , um vetor $\boldsymbol{x}\in\mathbb{R}^{n}$ e uma perturbação $\boldsymbol{\delta_{x}}\in\mathbb{R}^{n}$ . Além disso, sejam

	$\displaystyle\boldsymbol{y}$	$\displaystyle=A\boldsymbol{x}$		(3.199)
	$\displaystyle\boldsymbol{y}+\boldsymbol{\delta_{y}}$	$\displaystyle=A(\boldsymbol{x}+\boldsymbol{\delta_{x}}).$		(3.200)

Ou seja, $\boldsymbol{\delta_{y}}$ é a perturbação em $\boldsymbol{y}$ propagada da aplicação de $A$ em $\boldsymbol{x}$ com perturbação $\boldsymbol{\delta_{x}}$ .

Agora, podemos estimar a razão entre os erros relativos

	$\displaystyle e_{rel}(\boldsymbol{y}):=\\|\boldsymbol{\delta_{y}}\\|/\\|% \boldsymbol{y}\\|$		(3.201)
	$\displaystyle e_{rel}(\boldsymbol{x}):=\\|\boldsymbol{\delta_{x}}\\|/\\|% \boldsymbol{x}\\|$		(3.202)

da seguinte forma

$\displaystyle\frac{\\|\boldsymbol{y}\\|/\\|\boldsymbol{\delta_{y}}\\|}{\\|% \boldsymbol{x}\\|/\\|\boldsymbol{\delta_{x}}\\|}$	$\displaystyle=\frac{\\|\boldsymbol{y}\\|}{\\|\boldsymbol{\delta_{y}}\\|}\frac{\\|% \boldsymbol{\delta_{x}}\\|}{\\|\boldsymbol{x}\\|}$	(3.203)
	$\displaystyle=\frac{\\|A\boldsymbol{x}\\|}{\\|\boldsymbol{\delta_{y}}\\|}\frac{\\|A% ^{-1}\boldsymbol{\delta_{y}}\\|}{\\|\boldsymbol{x}\\|}$	(3.204)
	$\displaystyle\leq\frac{\\|A\\|\\|\boldsymbol{x}\\|\\|A^{-1}\\|\\|\boldsymbol{\delta_{% y}}\\|}{\\|\boldsymbol{\delta_{y}}\\|\\|\boldsymbol{x}\\|}$	(3.205)
	$\displaystyle=\\|A\\|\|A^{-1}\\|.$	(3.206)

Logo, temos a seguinte estimativa de propagação de erro

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}e_{rel}(% \boldsymbol{y})\leq\|A\|\|A^{-1}\|e_{rel}(\boldsymbol{x}).

(3.207)

Isto nos motiva a definir o número de condicionamento da matriz $A$ por

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\kappa(A):=\|% A\|\|A^{-1}\|.

(3.208)

A matriz identidade tem o menor número de condicionamento que é

\kappa(I)=1.

(3.209)

temos que $\kappa(A)\geq 1$ e quanto maior, mais mal condicionada é a matriz. Ou seja, quando maior $\kappa(A)$ , maior é a propagação dos erros ao se computar $A\boldsymbol{x}$ .

Exemplo 3.2.5.

Estudamos os seguintes exemplos:

a)

Matriz bem condicionada.

$A=\begin{bmatrix}1&-1&2\\ -2&\pi&4\\ 7&-5&\sqrt{2}\end{bmatrix},$ (3.210)

cujo número de condicionamento é $\kappa(A)=13.997$ .

⬇

1import numpy as np

2import numpy.linalg as npla

3A = np.array([[1., -1., 2.],

4 [-2., np.pi, 4.],

5 [7., -5., np.sqrt(2)]])

6cond_A = npla.cond(A, 'fro')

7print(f'cond(A) = {cond_A}')
b)

Matriz mal condicionada.

$B=\begin{bmatrix}1&0&0\\ 0&0&-2\\ 10^{5}&10^{-4}&10^{5}\end{bmatrix},$ (3.211)

tem número de condicionamento

$\kappa(B)=1.5811\times 10^{14},$ (3.212)

o que indica que $B$ é uma matriz mal condicionada.

⬇

1import numpy as np

2import numpy.linalg as npla

3A = np.array([[1., 0., 0.],

4 [0., 0., -2.],

5 [1e5, 1e-4, 1e5]])

6cond_A = npla.cond(A, 'fro')

7print(f'cond(A) = {cond_A:.4e}')

3.2.3 Exercícios

E. 3.2.1.

Compute a norma de cada um dos seguintes vetores:

a)

$\displaystyle\boldsymbol{u}=(1,0,-1)$
b)

$\displaystyle\boldsymbol{v}=(1,-2,-1,2)$
c)

$\displaystyle\boldsymbol{w}=(1,0,-1,2,-2)$

Resposta 0.

a) $\|u\|=\sqrt{2}$ ; b) $\|v\|=3.1623\mathrm{e}\!+\!0$ ; c) $\|w\|=3.1623\mathrm{e}\!+\!0$ ;

E. 3.2.2.

Compute a norma de cada uma das seguintes matrizes:

a)

$A=\begin{bmatrix}0&1\\ -1&0\end{bmatrix}$ (3.213)
b)

$B=\begin{bmatrix}2&1&-1\\ -1&0&2\\ 3&-2&0\end{bmatrix}$ (3.214)
c)

$C=\begin{bmatrix}2&1&-1&0\\ -1&0&2&0\\ 3&-2&0&-1\\ 0&1&0&2\end{bmatrix}$ (3.215)

Resposta 0.

a) $\|A\|=\sqrt{2}$ ; b) $\|B\|=4.8990\mathrm{e}\!+\!0$ ; c) $\|C\|=5.4772\mathrm{e}\!+\!0$

E. 3.2.3.

Compute o número de condicionamento de cada uma das seguintes matrizes:

a)

$A=\begin{bmatrix}0&1\\ -1&0\end{bmatrix}$ (3.216)
b)

$B=\begin{bmatrix}2&1&-1\\ -1&0&2\\ 3&-2&0\end{bmatrix}$ (3.217)
c)

$C=\begin{bmatrix}2&1&-1&0\\ -1&0&2&0\\ 3&-2&0&-1\\ 0&1&0&2\end{bmatrix}$ (3.218)
d)

$D=\begin{bmatrix}1\mathrm{e}\!+\!5&1\mathrm{e}\!-\!1&0&3\mathrm{e}\!-\!3\\ 0&-1\mathrm{e}\!-\!1&0&-1\mathrm{e}\!-\!3\\ 1\mathrm{e}\!+\!5&1\mathrm{e}\!-\!1&1\mathrm{e}\!-\!2&2\mathrm{e}\!-\!3\\ 0&-1\mathrm{e}\!-\!1&0&0\end{bmatrix}$ (3.219)
e)

$E=\begin{bmatrix}1&1\mathrm{e}\!-\!6&0&3\mathrm{e}\!-\!8\\ -1\mathrm{e}\!-\!6&-2\mathrm{e}\!-\!16&-2\mathrm{e}\!-\!1&0\\ 1\mathrm{e}\!+\!7&1\mathrm{e}\!+\!1&1&2\mathrm{e}\!-\!1\\ 0&-1\mathrm{e}\!2\!&0&0\end{bmatrix}$ (3.220)

Qual das matriz acima é a mais mal condicionada? Justifique sua resposta.

Resposta 0.

a) $\kappa(A)=2$ ; b) $\kappa(B)=4.6904\mathrm{e}\!+\!0$ ; c) $\kappa(C)=6.6291\mathrm{e}\!+\!0$ ; d) $\kappa(D)=4.6314\mathrm{e}\!+\!8$ e) $\kappa(E)=1.0000\mathrm{e}\!+\!15$ , a mais mal condicionada.

E. 3.2.4.

Considere o seguinte sistema linear

$\displaystyle 10^{-12}x_{1}+20x_{2}+3x_{3}$	$\displaystyle=-1,$	(3.221)
$\displaystyle 2.001x_{1}+10^{-5}x_{2}+-x_{3}$	$\displaystyle=-2,$	(3.222)
$\displaystyle x_{1}-2x_{2}-0.1x_{3}$	$\displaystyle=0.1.$	(3.223)

a)

Compute a norma do vetor dos termos constantes deste sistema.
b)

Compute a norma matriz dos coeficientes deste sistema.
c)

Compute o número de condicionamento da matriz dos coeficientes deste sistema.

Resposta 0.

a) $2.2383\mathrm{e}\!+\!0$ ; b) $2.0323\mathrm{e}\!+\!1$ ; c) $3.5128\mathrm{e}\!+\!1$

E. 3.2.5.

Considere

A=\begin{bmatrix}1\mathrm{e}\!+\!8&1&2&-1\\ 1\mathrm{e}\!-\!2&1\mathrm{e}\!-\!1&0&1\mathrm{e}\!-\!3\\ 0&1&1\mathrm{e}\!2\!&0\\ 1\mathrm{e}\!-\!5&0&0&1\mathrm{e}\!-\!4\end{bmatrix}

(3.224)

a)

Compute $\kappa(A)$ .
b)

Aloque a matriz diagonal $M$ , cujos elementos da diagonal sejam iguais aos da matriz $A$ .
c)

Verifique que o número de condicionamento de $M^{-1}A$ é melhor que o de $A$ .

Resposta 0.

a) $\kappa(A)=1.0\mathrm{e}\!+\!12$ ; b) M = np.diag(np.diag(A)); c) $\kappa(M^{-1}A)=4.0201\mathrm{e}\!+\!0$

Análise Numérica

E. 3.2.6.

Mostre que a Desigualdade triangular

\|\boldsymbol{u}+\boldsymbol{v}\|\leq\|\boldsymbol{u}\|+\|\boldsymbol{v}\|

(3.225)

vale para quaisquer $\boldsymbol{u},\boldsymbol{v}\in\mathbb{R}^{n}$ .

Resposta 0.

Dica: use a Desigualdade de Cauchy-Schwarz (3.152).

E. 3.2.7.

Mostre a Proposição 3.2.2.

Resposta 0.

Dica: use as propriedades da norma de vetores.

Envie seu comentário

Aproveito para agradecer a todas/os que de forma assídua ou esporádica contribuem enviando correções, sugestões e críticas!

Este texto é disponibilizado nos termos da Licença Creative Commons Atribuição-CompartilhaIgual 4.0 Internacional. Ícones e elementos gráficos podem estar sujeitos a condições adicionais.

Matemática Numérica I

3 Métodos para Sistemas Lineares 3.1 Método da Decomposição LU 3.3 Métodos de Jacobi e de Gauss-Seidel

Ajude a manter o site livre, gratuito e sem propagandas. Colabore!

3.2 Norma e Número de Condicionamento

Nesta seção, fazemos uma rápida discussão sobre normas de vetores e matrizes e sobre o condicionamento de uma matriz.

3.2.1 Norma $L^{2}$

Norma de Vetores

A norma $L^{2}$ de um dado vetor $\boldsymbol{v}=(v_{1},v_{2},\dotsc,v_{n})\in\mathbb{R}^{n}$ é definida por

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|\boldsymbol% {v}\|:=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}.

(3.136)

Lembrando que o produto interno de dois vetores $\boldsymbol{u},\boldsymbol{v}\in\mathbb{R}^{2}$ é definido por

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\boldsymbol{u% }\cdot\boldsymbol{v}:=u_{1}v_{1}+u_{2}v_{2}+\cdots+u_{n}v_{n},

(3.137)

temos que

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|\boldsymbol% {v}\|=\sqrt{\boldsymbol{v}\cdot\boldsymbol{v}}.

(3.138)

Exemplo 3.2.1.

Sejam os vetores

	$\displaystyle\boldsymbol{u}$	$\displaystyle=(1,-2,3,-4),$		(3.139)
	$\displaystyle\boldsymbol{v}$	$\displaystyle=(-1,2,0,1).$		(3.140)

a)

$\displaystyle\boldsymbol{u}\cdot\boldsymbol{v}$ $\displaystyle=u_{1}v_{1}+u_{2}v_{2}+\cdots+u_{n}v_{n}$ (3.141)

$\displaystyle=1\cdot(-1)+(-2)\cdot 2+3\cdot 0+(-4)\cdot 1$ (3.142)

$\displaystyle=-1-4+0-4$ (3.143)

$\displaystyle=-9.$ (3.144)

⬇

1import numpy as np

2u = np.array([1., -2., 3., -4.])

3v = np.array([-1., 2., 0., 1.])

4udv = np.dot(u,v)

5print(f'\nu.v = {udv}')
b)

$\displaystyle\|\boldsymbol{v}\|$ $\displaystyle=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$ (3.145)

$\displaystyle=\sqrt{(-1)^{2}+2^{2}+0^{2}+1^{2}}$ (3.146)

$\displaystyle=\sqrt{1+4+1}$ (3.147)

$\displaystyle=\sqrt{6}=2.4495\mathrm{e}\!+\!0.$ (3.148)

⬇

1import numpy as np

2import numpy.linalg as npla

3v = np.array([-1., 2., 0., 1.])

4norm_v = npla.norm(v)

5print(f'\n||v|| = {norm_v}')

Proposição 3.2.1.

(Propriedades da Norma para Vetores.) Dados os vetores $\boldsymbol{u},\boldsymbol{v}\in\mathbb{R}^{n}$ e um escalar $\lambda\in\mathbb{R}$ , temos:

a)

Positividade

$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \|\boldsymbol{v}\|\geq 0.$ (3.149)

$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \|\boldsymbol{v}\|=0\Leftrightarrow\boldsymbol{v}=\boldsymbol{0}.$ (3.150)
b)

Multiplicação por escalar

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|\lambda% \boldsymbol{v}\|=|\lambda|\cdot\|\boldsymbol{v}\|.$ (3.151)
c)

Desigualdade de Cauchy³²³²endnote: ³²Augustin-Louis Cauchy, 1789-1857, matemático francês. Fonte: Wikipédia: Augustin-Louis Cauchy.-Schwarz³³³³endnote: ³³Karl Hermann Amandus Schwarz, 1843-1921, matemático alemão. Fonte: Wikipédia: Hermann Amandus Schwarz.

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\boldsymbol{u% }\cdot\boldsymbol{v}\leq\|\boldsymbol{u}\|\cdot\|\boldsymbol{v}\|$ (3.152)
d)

Desigualdade triangular

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|\boldsymbol% {u}+\boldsymbol{v}\|\leq\|\boldsymbol{u}\|+\|\boldsymbol{v}\|$ (3.153)

Demonstração.

Sejam dados $\lambda\in\mathbb{R}$ e $\boldsymbol{u},\boldsymbol{v}\in\mathbb{R}^{n}$ .

a)

Positividade.

Observa-se diretamente que $v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}\geq 0$ . Então, como a raiz quadrada é uma função não-negativa, concluímos que

$\displaystyle\|\boldsymbol{v}\|$ $\displaystyle:=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$ (3.154)

$\displaystyle\geq 0$

No caso de $\boldsymbol{v}=\boldsymbol{0}$ , temos $v_{i}=0$ , $i=1,2,\dotsc,n$ , donde

$\displaystyle\|\boldsymbol{v}\|$ $\displaystyle=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$ (3.155)

$\displaystyle\|\boldsymbol{0}\|$ $\displaystyle=\sqrt{0}=0.$

Ou seja, se $\boldsymbol{v}=\boldsymbol{0}$ , então $\|\boldsymbol{v}\|=0$ . Agora, se $\boldsymbol{v}\neq\boldsymbol{0}$ , então tem-se $v_{i}\neq 0$ para algum $i=1,2,\dotsc,n$ . Logo, pela monotonicidade da função raiz quadrada, temos

$\displaystyle v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}>0$ (3.156)

$\displaystyle\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}>\sqrt{0}$ (3.157)

$\displaystyle\|\boldsymbol{v}\|>0.$ (3.158)

Portanto, concluímos que $\|\boldsymbol{v}\|=0$ se, e somente se, $\boldsymbol{v}=\boldsymbol{0}$ .
b)

Multiplicação por escalar.

Observamos que

$\lambda\boldsymbol{v}=(\lambda v_{1},\lambda v_{2},\dotsc,\lambda v_{n}).$ (3.159)

Então, segue por cálculo direto que

$\displaystyle\|\lambda\boldsymbol{v}\|$ $\displaystyle=\sqrt{(\lambda v_{1})^{2}+(\lambda v_{2})^{2}+\cdots+(\lambda v_% {n})^{2}}$ (3.160)

$\displaystyle=\sqrt{\lambda^{2}(v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2})}$ (3.161)

$\displaystyle=\sqrt{\lambda^{2}}\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$ (3.162)

$\displaystyle=|\lambda|\|\boldsymbol{v}\|.$ (3.163)

Desigualdade de Cauchy-Schwarz.

$\displaystyle(\boldsymbol{u}+\alpha\boldsymbol{v})\cdot(\boldsymbol{u}+\alpha% \boldsymbol{v})$	$\displaystyle\geq 0$	(3.164)
$\displaystyle\boldsymbol{u}\cdot\boldsymbol{u}+2\alpha\boldsymbol{u}\cdot% \boldsymbol{v}+\alpha^{2}\boldsymbol{v}\cdot\boldsymbol{v}$	$\displaystyle\geq 0$	(3.165)
$\displaystyle\\|\boldsymbol{u}\\|^{2}+2\alpha\boldsymbol{u}\cdot\boldsymbol{v}+% \alpha^{2}\\|\boldsymbol{v}\\|^{2}$	$\displaystyle\geq 0$	(3.166)

O lado direito desta desigualdade é um polinômio quadrático em $\alpha$ . Para ser não-negativo para todo $\alpha$ , seu discriminante precisa ser não-positivo, i.e.

	$\displaystyle(2\boldsymbol{u}\cdot\boldsymbol{v})^{2}-4\\|\boldsymbol{u}\\|^{2}% \\|\boldsymbol{v}\\|^{2}\leq 0$		(3.167)
	$\displaystyle(\boldsymbol{u}\cdot\boldsymbol{v})^{2}\leq\\|\boldsymbol{u}\\|^{2}% \\|\boldsymbol{v}\\|^{2}$		(3.168)
	$\displaystyle\sqrt{(\boldsymbol{u}\cdot\boldsymbol{v})^{2}}\leq\sqrt{\\|% \boldsymbol{u}\\|^{2}\\|\boldsymbol{v}\\|^{2}}$		(3.169)
	$\displaystyle\|\boldsymbol{u}\cdot\boldsymbol{v}\|\leq\\|\boldsymbol{u}\\|\\|% \boldsymbol{v}\\|$		(3.170)

Donde, concluímos que

\boldsymbol{u}\cdot\boldsymbol{v}\leq\|\boldsymbol{u}\|\|\boldsymbol{v}\|.

(3.171)

d)

Desigualdade triangular

Consulte o Exercício 3.2.6 para a demonstração da desigualdade triangular.

∎

Exemplo 3.2.2.

Vamos verificar a Desigualdade Triangular (3.153) para os vetores

	$\displaystyle\boldsymbol{u}$	$\displaystyle=(1,-2,3,-4),$		(3.172)
	$\displaystyle\boldsymbol{v}$	$\displaystyle=(-1,2,0,1).$		(3.173)

De fato, temos

$\displaystyle\\|\boldsymbol{u}+\boldsymbol{v}\\|$	$\displaystyle=\sqrt{(u_{1}+v_{1})^{2}+(u_{2}+v_{2})^{2}+(u_{3}+v_{3})^{2}+(u_{% n}+v_{n})^{2}}$	(3.174)
	$\displaystyle=\sqrt{0^{2}+0^{2}+3^{2}+3^{2}}$	(3.175)
	$\displaystyle=\sqrt{18}$	(3.176)
	$\displaystyle=4.2426\mathrm{e}\!+\!0$	(3.177)

$\displaystyle\\|\boldsymbol{u}\\|+\\|\boldsymbol{v}\\|$	$\displaystyle=\sqrt{u_{1}^{2}+u_{2}^{2}+u_{3}^{2}+u_{4}^{2}}$
	$\displaystyle+\sqrt{v_{1}^{2}+v_{2}^{2}+v_{3}^{2}+v_{4}^{2}}$	(3.178)
	$\displaystyle=\sqrt{1^{2}+(-2)^{2}+3^{2}+(-4)^{2}}$
	$\displaystyle+\sqrt{(-1)^{2}+2^{2}+0^{2}+1^{2}}$	(3.179)
	$\displaystyle=\sqrt{30}+\sqrt{6}$	(3.180)
	$\displaystyle=7.9267\mathrm{e}\!+\!0$	(3.181)

Ou seja, temos que

\|\boldsymbol{u}+\boldsymbol{v}\|\leq\|\boldsymbol{u}\|+\|\boldsymbol{v}\|

(3.182)

como esperado.

⬇

1import numpy as np

2u = np.array([1., -2., 3., -4.])

3v = np.array([-1., 2., 0., 1.])

4nupv = npla.norm(u+v)

5nupnv = npla.norm(u) + npla.norm(v)

6print('\nu.v <= ||u||.||v||?')

7print(nupv <= nupnv)

Norma de Matrizes

A norma $L^{2}$ induzida de uma dada matriz real $A=[a_{ij}]_{i,j=1}^{n}$ é definida por

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|A\|_{2}:=% \sup_{\boldsymbol{x}\in\mathbb{R}^{n},\|\boldsymbol{x}\|=1}\|A\boldsymbol{x}\|.

(3.183)

Pode-se mostrar que³⁴³⁴endnote: ³⁴Consulte [7, Section 1.3-3] para informações sobre a demonstração.

\|A\|_{2}=\sqrt{\lambda_{max}(A^{T}A)},

(3.184)

onde $\lambda_{max}(A^{T}A):=\max\{|\lambda|;\leavevmode\nobreak\ \lambda\text{ é % autovalor de }A^{T}A\}$ .

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|A\|=\left(% \sum_{i=1}^{n}\sum_{j=1}^{n}|a_{i,j}|^{2}\right)^{\frac{1}{2}}

(3.185)

Esta é uma generalização da norma euclidiana para vetores e é equivalente a norma $L^{2}$ induzida, i.e.

\|A\|_{2}\leq\|A\|\leq\sqrt{n}\|A\|_{2}.

(3.186)

Exemplo 3.2.3.

A matriz

A=\begin{bmatrix}1&-1&2\\ -2&\pi&4\\ 7&-5&\sqrt{2}\end{bmatrix}

(3.187)

tem norma

\|A\|=10.5768\mathrm{e}\!+\!0.

(3.188)

⬇

1import numpy as np

2import numpy.linalg as npla

3A = np.array([[1., -1., 2.],

4 [-2., np.pi, 4.],

5 [7., -5., np.sqrt(2)]])

6norm_A = npla.norm(A)

7print(norm_A)

Proposição 3.2.2.

(Propriedades da Norma para Matrizes.) Dadas as matrizes reais $A, B$ $n\times n$ , um vetor $v\in\mathbb{R}^{n}$ e um escalar $\lambda$ , temos

a)

Positividade

$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \|A\|\geq 0$ (3.189)

$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \|A\|=0\Leftrightarrow A=\boldsymbol{0}$ (3.190)
b)

Multiplicação por escalar

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|\lambda A\|% =|\lambda|\cdot\|A\|$ (3.191)
c)

Desigualdade triangular

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|A+B\|\leq\|% A\|+\|B\|$ (3.192)
d)

Sub-multiplicatividade

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|AB\|\leq\|A% \|\|B\|$ (3.193)
e)

Norma da aplicação

$\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\|Av\|\leq\|A% \|\|v\|$ (3.194)

Demonstração.

Consulte o Exercício 3.2.7. ∎

Exemplo 3.2.4.

Vamos ver exemplos das propriedades d) e e) da norma de matriz.

a)

Sub-multiplicatividade

$\displaystyle\|I_{n}^{2}\|=\|I_{n}\|=\sqrt{n}$ (3.195)

$\displaystyle\|I\|\cdot\|I\|=\|I\|^{2}=n$ (3.196)
b)

Norma da aplicação

$\displaystyle\|I_{n}\boldsymbol{v}\|=\|\boldsymbol{v}\|$ (3.197)

$\displaystyle\|I_{n}\|\|v\|=\sqrt{n}\|\boldsymbol{v}\|$ (3.198)

3.2.2 Número de Condicionamento

Em revisão

	$\displaystyle\boldsymbol{y}$	$\displaystyle=A\boldsymbol{x}$		(3.199)
	$\displaystyle\boldsymbol{y}+\boldsymbol{\delta_{y}}$	$\displaystyle=A(\boldsymbol{x}+\boldsymbol{\delta_{x}}).$		(3.200)

Ou seja, $\boldsymbol{\delta_{y}}$ é a perturbação em $\boldsymbol{y}$ propagada da aplicação de $A$ em $\boldsymbol{x}$ com perturbação $\boldsymbol{\delta_{x}}$ .

Agora, podemos estimar a razão entre os erros relativos

	$\displaystyle e_{rel}(\boldsymbol{y}):=\\|\boldsymbol{\delta_{y}}\\|/\\|% \boldsymbol{y}\\|$		(3.201)
	$\displaystyle e_{rel}(\boldsymbol{x}):=\\|\boldsymbol{\delta_{x}}\\|/\\|% \boldsymbol{x}\\|$		(3.202)

da seguinte forma

$\displaystyle\frac{\\|\boldsymbol{y}\\|/\\|\boldsymbol{\delta_{y}}\\|}{\\|% \boldsymbol{x}\\|/\\|\boldsymbol{\delta_{x}}\\|}$	$\displaystyle=\frac{\\|\boldsymbol{y}\\|}{\\|\boldsymbol{\delta_{y}}\\|}\frac{\\|% \boldsymbol{\delta_{x}}\\|}{\\|\boldsymbol{x}\\|}$	(3.203)
	$\displaystyle=\frac{\\|A\boldsymbol{x}\\|}{\\|\boldsymbol{\delta_{y}}\\|}\frac{\\|A% ^{-1}\boldsymbol{\delta_{y}}\\|}{\\|\boldsymbol{x}\\|}$	(3.204)
	$\displaystyle\leq\frac{\\|A\\|\\|\boldsymbol{x}\\|\\|A^{-1}\\|\\|\boldsymbol{\delta_{% y}}\\|}{\\|\boldsymbol{\delta_{y}}\\|\\|\boldsymbol{x}\\|}$	(3.205)
	$\displaystyle=\\|A\\|\|A^{-1}\\|.$	(3.206)

Logo, temos a seguinte estimativa de propagação de erro

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}e_{rel}(% \boldsymbol{y})\leq\|A\|\|A^{-1}\|e_{rel}(\boldsymbol{x}).

(3.207)

Isto nos motiva a definir o número de condicionamento da matriz $A$ por

\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}\kappa(A):=\|% A\|\|A^{-1}\|.

(3.208)

A matriz identidade tem o menor número de condicionamento que é

\kappa(I)=1.

(3.209)

temos que $\kappa(A)\geq 1$ e quanto maior, mais mal condicionada é a matriz. Ou seja, quando maior $\kappa(A)$ , maior é a propagação dos erros ao se computar $A\boldsymbol{x}$ .

Exemplo 3.2.5.

Estudamos os seguintes exemplos:

a)

Matriz bem condicionada.

$A=\begin{bmatrix}1&-1&2\\ -2&\pi&4\\ 7&-5&\sqrt{2}\end{bmatrix},$ (3.210)

cujo número de condicionamento é $\kappa(A)=13.997$ .

⬇

1import numpy as np

2import numpy.linalg as npla

3A = np.array([[1., -1., 2.],

4 [-2., np.pi, 4.],

5 [7., -5., np.sqrt(2)]])

6cond_A = npla.cond(A, 'fro')

7print(f'cond(A) = {cond_A}')
b)

Matriz mal condicionada.

$B=\begin{bmatrix}1&0&0\\ 0&0&-2\\ 10^{5}&10^{-4}&10^{5}\end{bmatrix},$ (3.211)

tem número de condicionamento

$\kappa(B)=1.5811\times 10^{14},$ (3.212)

o que indica que $B$ é uma matriz mal condicionada.

⬇

1import numpy as np

2import numpy.linalg as npla

3A = np.array([[1., 0., 0.],

4 [0., 0., -2.],

5 [1e5, 1e-4, 1e5]])

6cond_A = npla.cond(A, 'fro')

7print(f'cond(A) = {cond_A:.4e}')

3.2.3 Exercícios

E. 3.2.1.

Compute a norma de cada um dos seguintes vetores:

a)

$\displaystyle\boldsymbol{u}=(1,0,-1)$
b)

$\displaystyle\boldsymbol{v}=(1,-2,-1,2)$
c)

$\displaystyle\boldsymbol{w}=(1,0,-1,2,-2)$

Resposta 0.

a) $\|u\|=\sqrt{2}$ ; b) $\|v\|=3.1623\mathrm{e}\!+\!0$ ; c) $\|w\|=3.1623\mathrm{e}\!+\!0$ ;

E. 3.2.2.

Compute a norma de cada uma das seguintes matrizes:

a)

$A=\begin{bmatrix}0&1\\ -1&0\end{bmatrix}$ (3.213)
b)

$B=\begin{bmatrix}2&1&-1\\ -1&0&2\\ 3&-2&0\end{bmatrix}$ (3.214)
c)

$C=\begin{bmatrix}2&1&-1&0\\ -1&0&2&0\\ 3&-2&0&-1\\ 0&1&0&2\end{bmatrix}$ (3.215)

Resposta 0.

a) $\|A\|=\sqrt{2}$ ; b) $\|B\|=4.8990\mathrm{e}\!+\!0$ ; c) $\|C\|=5.4772\mathrm{e}\!+\!0$

E. 3.2.3.

Compute o número de condicionamento de cada uma das seguintes matrizes:

a)

$A=\begin{bmatrix}0&1\\ -1&0\end{bmatrix}$ (3.216)
b)

$B=\begin{bmatrix}2&1&-1\\ -1&0&2\\ 3&-2&0\end{bmatrix}$ (3.217)
c)

$C=\begin{bmatrix}2&1&-1&0\\ -1&0&2&0\\ 3&-2&0&-1\\ 0&1&0&2\end{bmatrix}$ (3.218)
d)

$D=\begin{bmatrix}1\mathrm{e}\!+\!5&1\mathrm{e}\!-\!1&0&3\mathrm{e}\!-\!3\\ 0&-1\mathrm{e}\!-\!1&0&-1\mathrm{e}\!-\!3\\ 1\mathrm{e}\!+\!5&1\mathrm{e}\!-\!1&1\mathrm{e}\!-\!2&2\mathrm{e}\!-\!3\\ 0&-1\mathrm{e}\!-\!1&0&0\end{bmatrix}$ (3.219)
e)

$E=\begin{bmatrix}1&1\mathrm{e}\!-\!6&0&3\mathrm{e}\!-\!8\\ -1\mathrm{e}\!-\!6&-2\mathrm{e}\!-\!16&-2\mathrm{e}\!-\!1&0\\ 1\mathrm{e}\!+\!7&1\mathrm{e}\!+\!1&1&2\mathrm{e}\!-\!1\\ 0&-1\mathrm{e}\!2\!&0&0\end{bmatrix}$ (3.220)

Qual das matriz acima é a mais mal condicionada? Justifique sua resposta.

Resposta 0.

E. 3.2.4.

Considere o seguinte sistema linear

$\displaystyle 10^{-12}x_{1}+20x_{2}+3x_{3}$	$\displaystyle=-1,$	(3.221)
$\displaystyle 2.001x_{1}+10^{-5}x_{2}+-x_{3}$	$\displaystyle=-2,$	(3.222)
$\displaystyle x_{1}-2x_{2}-0.1x_{3}$	$\displaystyle=0.1.$	(3.223)

a)

Compute a norma do vetor dos termos constantes deste sistema.
b)

Compute a norma matriz dos coeficientes deste sistema.
c)

Compute o número de condicionamento da matriz dos coeficientes deste sistema.

Resposta 0.

a) $2.2383\mathrm{e}\!+\!0$ ; b) $2.0323\mathrm{e}\!+\!1$ ; c) $3.5128\mathrm{e}\!+\!1$

E. 3.2.5.

Considere

A=\begin{bmatrix}1\mathrm{e}\!+\!8&1&2&-1\\ 1\mathrm{e}\!-\!2&1\mathrm{e}\!-\!1&0&1\mathrm{e}\!-\!3\\ 0&1&1\mathrm{e}\!2\!&0\\ 1\mathrm{e}\!-\!5&0&0&1\mathrm{e}\!-\!4\end{bmatrix}

(3.224)

a)

Compute $\kappa(A)$ .
b)

Aloque a matriz diagonal $M$ , cujos elementos da diagonal sejam iguais aos da matriz $A$ .
c)

Verifique que o número de condicionamento de $M^{-1}A$ é melhor que o de $A$ .

Resposta 0.

a) $\kappa(A)=1.0\mathrm{e}\!+\!12$ ; b) M = np.diag(np.diag(A)); c) $\kappa(M^{-1}A)=4.0201\mathrm{e}\!+\!0$

Análise Numérica

E. 3.2.6.

Mostre que a Desigualdade triangular

\|\boldsymbol{u}+\boldsymbol{v}\|\leq\|\boldsymbol{u}\|+\|\boldsymbol{v}\|

(3.225)

vale para quaisquer $\boldsymbol{u},\boldsymbol{v}\in\mathbb{R}^{n}$ .

Resposta 0.

Dica: use a Desigualdade de Cauchy-Schwarz (3.152).

E. 3.2.7.

Mostre a Proposição 3.2.2.

Resposta 0.

Dica: use as propriedades da norma de vetores.

Envie seu comentário

Aproveito para agradecer a todas/os que de forma assídua ou esporádica contribuem enviando correções, sugestões e críticas!

Este texto é disponibilizado nos termos da Licença Creative Commons Atribuição-CompartilhaIgual 4.0 Internacional. Ícones e elementos gráficos podem estar sujeitos a condições adicionais.

Pedro H A Konzen

| | | |

$\displaystyle\frac{\\|\boldsymbol{y}\\|/\\|\boldsymbol{\delta_{y}}\\|}{\\|% \boldsymbol{x}\\|/\\|\boldsymbol{\delta_{x}}\\|}$	$\displaystyle=\frac{\\|\boldsymbol{y}\\|}{\\|\boldsymbol{\delta_{y}}\\|}\frac{\\|% \boldsymbol{\delta_{x}}\\|}{\\|\boldsymbol{x}\\|}$	(3.203)
	$\displaystyle=\frac{\\|A\boldsymbol{x}\\|}{\\|\boldsymbol{\delta_{y}}\\|}\frac{\\|A% ^{-1}\boldsymbol{\delta_{y}}\\|}{\\|\boldsymbol{x}\\|}$	(3.204)
	$\displaystyle\leq\frac{\\|A\\|\\|\boldsymbol{x}\\|\\|A^{-1}\\|\\|\boldsymbol{\delta_{% y}}\\|}{\\|\boldsymbol{\delta_{y}}\\|\\|\boldsymbol{x}\\|}$	(3.205)
	$\displaystyle=\\|A\\|\|A^{-1}\\|.$	(3.206)

$\displaystyle\frac{\\|\boldsymbol{y}\\|/\\|\boldsymbol{\delta_{y}}\\|}{\\|% \boldsymbol{x}\\|/\\|\boldsymbol{\delta_{x}}\\|}$	$\displaystyle=\frac{\\|\boldsymbol{y}\\|}{\\|\boldsymbol{\delta_{y}}\\|}\frac{\\|% \boldsymbol{\delta_{x}}\\|}{\\|\boldsymbol{x}\\|}$	(3.203)
	$\displaystyle=\frac{\\|A\boldsymbol{x}\\|}{\\|\boldsymbol{\delta_{y}}\\|}\frac{\\|A% ^{-1}\boldsymbol{\delta_{y}}\\|}{\\|\boldsymbol{x}\\|}$	(3.204)
	$\displaystyle\leq\frac{\\|A\\|\\|\boldsymbol{x}\\|\\|A^{-1}\\|\\|\boldsymbol{\delta_{% y}}\\|}{\\|\boldsymbol{\delta_{y}}\\|\\|\boldsymbol{x}\\|}$	(3.205)
	$\displaystyle=\\|A\\|\|A^{-1}\\|.$	(3.206)

$\displaystyle\boldsymbol{u}\cdot\boldsymbol{v}$	$\displaystyle=u_{1}v_{1}+u_{2}v_{2}+\cdots+u_{n}v_{n}$	(3.141)
	$\displaystyle=1\cdot(-1)+(-2)\cdot 2+3\cdot 0+(-4)\cdot 1$	(3.142)
	$\displaystyle=-1-4+0-4$	(3.143)
	$\displaystyle=-9.$	(3.144)

$\displaystyle\\|\boldsymbol{v}\\|$	$\displaystyle=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$	(3.145)
	$\displaystyle=\sqrt{(-1)^{2}+2^{2}+0^{2}+1^{2}}$	(3.146)
	$\displaystyle=\sqrt{1+4+1}$	(3.147)
	$\displaystyle=\sqrt{6}=2.4495\mathrm{e}\!+\!0.$	(3.148)

	$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \\|\boldsymbol{v}\\|\geq 0.$		(3.149)
	$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \\|\boldsymbol{v}\\|=0\Leftrightarrow\boldsymbol{v}=\boldsymbol{0}.$		(3.150)

	$\displaystyle\\|\boldsymbol{v}\\|$	$\displaystyle:=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$		(3.154)
		$\displaystyle\geq 0$		(3.154)

	$\displaystyle\\|\boldsymbol{v}\\|$	$\displaystyle=\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$		(3.155)
	$\displaystyle\\|\boldsymbol{0}\\|$	$\displaystyle=\sqrt{0}=0.$		(3.155)

	$\displaystyle v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}>0$		(3.156)
	$\displaystyle\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}>\sqrt{0}$		(3.157)
	$\displaystyle\\|\boldsymbol{v}\\|>0.$		(3.158)

$\displaystyle\\|\lambda\boldsymbol{v}\\|$	$\displaystyle=\sqrt{(\lambda v_{1})^{2}+(\lambda v_{2})^{2}+\cdots+(\lambda v_% {n})^{2}}$	(3.160)
	$\displaystyle=\sqrt{\lambda^{2}(v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2})}$	(3.161)
	$\displaystyle=\sqrt{\lambda^{2}}\sqrt{v_{1}^{2}+v_{2}^{2}+\cdots+v_{n}^{2}}$	(3.162)
	$\displaystyle=\|\lambda\|\\|\boldsymbol{v}\\|.$	(3.163)

	$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \\|A\\|\geq 0$		(3.189)
	$\displaystyle\color[rgb]{0,0,1}\definecolor[named]{pgfstrokecolor}{rgb}{0,0,1}% \\|A\\|=0\Leftrightarrow A=\boldsymbol{0}$		(3.190)

	$\displaystyle\\|I_{n}^{2}\\|=\\|I_{n}\\|=\sqrt{n}$		(3.195)
	$\displaystyle\\|I\\|\cdot\\|I\\|=\\|I\\|^{2}=n$		(3.196)

	$\displaystyle\\|I_{n}\boldsymbol{v}\\|=\\|\boldsymbol{v}\\|$		(3.197)
	$\displaystyle\\|I_{n}\\|\\|v\\|=\sqrt{n}\\|\boldsymbol{v}\\|$		(3.198)

Política de uso de dados

Política de uso de dados

Matemática Numérica I

3.2 Norma e Número de Condicionamento

3.2.1 Norma L2

Norma de Vetores

Exemplo 3.2.1.

Proposição 3.2.1.

Demonstração.

Exemplo 3.2.2.

Norma de Matrizes

Exemplo 3.2.3.

Proposição 3.2.2.

Demonstração.

Exemplo 3.2.4.

3.2.2 Número de Condicionamento

Exemplo 3.2.5.

3.2.3 Exercícios

E. 3.2.1.

Resposta 0.

E. 3.2.2.

Resposta 0.

E. 3.2.3.

Resposta 0.

E. 3.2.4.

Resposta 0.

E. 3.2.5.

Resposta 0.

Análise Numérica

E. 3.2.6.

Resposta 0.

E. 3.2.7.

Resposta 0.

Envie seu comentário

Matemática Numérica I

3.2 Norma e Número de Condicionamento

3.2.1 Norma L2

Norma de Vetores

Exemplo 3.2.1.

Proposição 3.2.1.

Demonstração.

Exemplo 3.2.2.

Norma de Matrizes

Exemplo 3.2.3.

Proposição 3.2.2.

Demonstração.

Exemplo 3.2.4.

3.2.2 Número de Condicionamento

Exemplo 3.2.5.

3.2.3 Exercícios

E. 3.2.1.

Resposta 0.

E. 3.2.2.

Resposta 0.

E. 3.2.3.

Resposta 0.

E. 3.2.4.

Resposta 0.

E. 3.2.5.

Resposta 0.

Análise Numérica

E. 3.2.6.

Resposta 0.

E. 3.2.7.

Resposta 0.

Envie seu comentário

3.2.1 Norma $L^{2}$

3.2.1 Norma $L^{2}$