Þ±cell_dependenciesÞÙ$47690b4e-8df6-4d07-86e2-1399ef75fb3e„´precedence_heuristic §cell_idÙ$47690b4e-8df6-4d07-86e2-1399ef75fb3e´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$3e592d32-9e4a-45e6-b179-b7b72542a60c„´precedence_heuristic §cell_idÙ$3e592d32-9e4a-45e6-b179-b7b72542a60c´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$9457901e-3b05-4c35-ab14-db1399079d79„´precedence_heuristic §cell_idÙ$9457901e-3b05-4c35-ab14-db1399079d79´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$09fe122e-4487-4b1c-9c03-4f160e7e3857„´precedence_heuristic §cell_idÙ$09fe122e-4487-4b1c-9c03-4f160e7e3857´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$d5cb2d7c-a627-4b79-a45e-f11ffe3b15a9„´precedence_heuristic §cell_idÙ$d5cb2d7c-a627-4b79-a45e-f11ffe3b15a9´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$af5521f4-e4a4-4411-972e-c41bc7a85fea„´precedence_heuristic §cell_idÙ$af5521f4-e4a4-4411-972e-c41bc7a85fea´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$901b5296-373c-498b-890a-245b190e68b8„´precedence_heuristic §cell_idÙ$901b5296-373c-498b-890a-245b190e68b8´downstream_cells_map�¬show_outline‘Ù$5259f8a1-2879-45b4-9e74-a2ff0feaff79²upstream_cells_mapФCore�§@md_str�¤Base�·PlutoRunner.create_bond�«PlutoRunner�¨CheckBox�¯Core.applicable�¥@bind�¨Base.get�¨getindex�Ù$a730c823-e0f4-48ed-bf6b-996f8de49460„´precedence_heuristic §cell_idÙ$a730c823-e0f4-48ed-bf6b-996f8de49460´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$be0c5a02-d3f9-447d-83bd-dbc5aba9266d„´precedence_heuristic §cell_idÙ$be0c5a02-d3f9-447d-83bd-dbc5aba9266d´downstream_cells_map€²upstream_cells_map�¯TableOfContents�Ù$f24d7aff-7cc9-4ae4-b757-ed170aaef555„´precedence_heuristic §cell_idÙ$f24d7aff-7cc9-4ae4-b757-ed170aaef555´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$891ad571-27cb-483c-a0a3-4b53ad9d6c17„´precedence_heuristic §cell_idÙ$891ad571-27cb-483c-a0a3-4b53ad9d6c17´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$cf35b757-1783-4316-977c-b182aa10ebb5„´precedence_heuristic §cell_idÙ$cf35b757-1783-4316-977c-b182aa10ebb5´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$9cc8bd7c-7deb-4db9-9175-6d636c6294fe„´precedence_heuristic §cell_idÙ$9cc8bd7c-7deb-4db9-9175-6d636c6294fe´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$976f90d6-eb84-4a45-8d51-8822a283ac35„´precedence_heuristic §cell_idÙ$976f90d6-eb84-4a45-8d51-8822a283ac35´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$877bcf5d-741b-4cfc-b069-ba033eecdcd5„´precedence_heuristic§cell_idÙ$877bcf5d-741b-4cfc-b069-ba033eecdcd5´downstream_cells_mapƒ§PlutoUI�°HypertextLiteral‘Ù$5259f8a1-2879-45b4-9e74-a2ff0feaff79²PlutoTeachingTools�²upstream_cells_map€Ù$008bb52c-c638-11f1-89c7-731f0c90a37d„´precedence_heuristic §cell_idÙ$008bb52c-c638-11f1-89c7-731f0c90a37d´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$450839de-2597-4221-8e68-cd2f6b441e0c„´precedence_heuristic §cell_idÙ$450839de-2597-4221-8e68-cd2f6b441e0c´downstream_cells_map€²upstream_cells_map‚§@md_str�¨getindex�Ù$5259f8a1-2879-45b4-9e74-a2ff0feaff79„´precedence_heuristic §cell_idÙ$5259f8a1-2879-45b4-9e74-a2ff0feaff79´downstream_cells_map�§Sidebar�²upstream_cells_mapŒ¦String�¤@htl�Ù HypertextLiteral.attribute_value�·HypertextLiteral.Result�³RobustLocalResource�·HypertextLiteral.Bypass�¸HypertextLiteral.content�¬show_outline‘Ù$901b5296-373c-498b-890a-245b190e68b8®Markdown.parse�°HypertextLiteral‘Ù$877bcf5d-741b-4cfc-b069-ba033eecdcd5¨Markdown�¤read�´cell_execution_orderÜÙ$877bcf5d-741b-4cfc-b069-ba033eecdcd5Ù$008bb52c-c638-11f1-89c7-731f0c90a37dÙ$f24d7aff-7cc9-4ae4-b757-ed170aaef555Ù$9cc8bd7c-7deb-4db9-9175-6d636c6294feÙ$d5cb2d7c-a627-4b79-a45e-f11ffe3b15a9Ù$891ad571-27cb-483c-a0a3-4b53ad9d6c17Ù$47690b4e-8df6-4d07-86e2-1399ef75fb3eÙ$976f90d6-eb84-4a45-8d51-8822a283ac35Ù$cf35b757-1783-4316-977c-b182aa10ebb5Ù$450839de-2597-4221-8e68-cd2f6b441e0cÙ$3e592d32-9e4a-45e6-b179-b7b72542a60cÙ$af5521f4-e4a4-4411-972e-c41bc7a85feaÙ$9457901e-3b05-4c35-ab14-db1399079d79Ù$a730c823-e0f4-48ed-bf6b-996f8de49460Ù$09fe122e-4487-4b1c-9c03-4f160e7e3857Ù$be0c5a02-d3f9-447d-83bd-dbc5aba9266dÙ$901b5296-373c-498b-890a-245b190e68b8Ù$5259f8a1-2879-45b4-9e74-a2ff0feaff79´last_hot_reload_timeË®process_status¥ready¤pathÙZ/workspace/matmat/teaching_math_aspects_materials/notes/src/90_Revision_matrix_concepts.jl­pluto_version¦v1.1.0ªcell_orderÜÙ$008bb52c-c638-11f1-89c7-731f0c90a37dÙ$877bcf5d-741b-4cfc-b069-ba033eecdcd5Ù$f24d7aff-7cc9-4ae4-b757-ed170aaef555Ù$9cc8bd7c-7deb-4db9-9175-6d636c6294feÙ$d5cb2d7c-a627-4b79-a45e-f11ffe3b15a9Ù$891ad571-27cb-483c-a0a3-4b53ad9d6c17Ù$47690b4e-8df6-4d07-86e2-1399ef75fb3eÙ$976f90d6-eb84-4a45-8d51-8822a283ac35Ù$cf35b757-1783-4316-977c-b182aa10ebb5Ù$450839de-2597-4221-8e68-cd2f6b441e0cÙ$3e592d32-9e4a-45e6-b179-b7b72542a60cÙ$af5521f4-e4a4-4411-972e-c41bc7a85feaÙ$9457901e-3b05-4c35-ab14-db1399079d79Ù$a730c823-e0f4-48ed-bf6b-996f8de49460Ù$09fe122e-4487-4b1c-9c03-4f160e7e3857Ù$be0c5a02-d3f9-447d-83bd-dbc5aba9266dÙ$901b5296-373c-498b-890a-245b190e68b8Ù$5259f8a1-2879-45b4-9e74-a2ff0feaff79­julia_version§v1.12.7±published_objects€¥nbpkgжwaiting_for_permissionÂÙ,waiting_for_permission_but_probably_disabled²installed_versions…§PlutoUI¦0.7.83Ù!__internal_julia_manifest_version¦1.12.7°HypertextLiteral¥1.0.0¸__internal_julia_version¦1.12.7²PlutoTeachingTools¥0.4.7°terminal_outputs„ªnbpkg_syncÚ' Waiting for other notebooks to finish Pkg operations... === Resolving... ===  Project No packages added to or removed from `~/.julia/scratchspaces/c3e4b0f8-55cb-11ea-2926-15256bba5781/pkg_envs/env_kxrynwfjrp/Project.toml`  Manifest No packages added to or removed from `~/.julia/scratchspaces/c3e4b0f8-55cb-11ea-2926-15256bba5781/pkg_envs/env_kxrynwfjrp/Manifest.toml` Instantiating... === Precompiling... === §PlutoUIÚ' Waiting for other notebooks to finish Pkg operations... === Resolving... ===  Project No packages added to or removed from `~/.julia/scratchspaces/c3e4b0f8-55cb-11ea-2926-15256bba5781/pkg_envs/env_kxrynwfjrp/Project.toml`  Manifest No packages added to or removed from `~/.julia/scratchspaces/c3e4b0f8-55cb-11ea-2926-15256bba5781/pkg_envs/env_kxrynwfjrp/Manifest.toml` Instantiating... === Precompiling... === °HypertextLiteralÚ' Waiting for other notebooks to finish Pkg operations... === Resolving... ===  Project No packages added to or removed from `~/.julia/scratchspaces/c3e4b0f8-55cb-11ea-2926-15256bba5781/pkg_envs/env_kxrynwfjrp/Project.toml`  Manifest No packages added to or removed from `~/.julia/scratchspaces/c3e4b0f8-55cb-11ea-2926-15256bba5781/pkg_envs/env_kxrynwfjrp/Manifest.toml` Instantiating... === Precompiling... === ²PlutoTeachingToolsÚ' Waiting for other notebooks to finish Pkg operations... === Resolving... ===  Project No packages added to or removed from `~/.julia/scratchspaces/c3e4b0f8-55cb-11ea-2926-15256bba5781/pkg_envs/env_kxrynwfjrp/Project.toml`  Manifest No packages added to or removed from `~/.julia/scratchspaces/c3e4b0f8-55cb-11ea-2926-15256bba5781/pkg_envs/env_kxrynwfjrp/Manifest.toml` Instantiating... === Precompiling... === §enabledìinstantiated÷restart_recommended_msgÀ´restart_required_msgÀ¯install_time_nsÏ&Võ­busy_packages�«cell_inputsÞÙ$47690b4e-8df6-4d07-86e2-1399ef75fb3e„§cell_idÙ$47690b4e-8df6-4d07-86e2-1399ef75fb3e¤codeÚšmd""" !!! note "Theorem 1" Let $A \in \mathbb C^{n \times n}$. $A$ is injective if and only if $\ker (A) = \{ 0 \}$. > *Proof.* > > $\boxed \Rightarrow$ : Let $v \in \ker (A)$. Then > ```math > \begin{align} > Av = 0 &\Rightarrow A(-v) = - Av = 0 > \\ &\hspace{-0.8em} \stackrel{\text{injective}}{\Rightarrow} v = -v > \\ &\Rightarrow v = 0 > \end{align} > ``` > > $\boxed \Leftarrow$ : Let $A v_1 = A v_2$. Then > ```math > \begin{align} > A(v_1 - v_2) = A v_1 - A v_2 = 0 &\Rightarrow v_1 - v_2 \in \ker (A) > \\&\hspace{-1.4em} \stackrel{\ker (A) = \{ 0\}}{\Rightarrow} v_1 - v_2 = 0 > \\ &\Rightarrow v_1 = v_2 > \\ & && & \square > \end{align} > ``` """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$3e592d32-9e4a-45e6-b179-b7b72542a60c„§cell_idÙ$3e592d32-9e4a-45e6-b179-b7b72542a60c¤codeÚ~md""" !!! tip "Remark" $\| A \|_2$ is also equal to the largest singular value of $A$. In addition, since $A^HA$ (and $A A ^H$) is a Hermitian positive semi-definite matrix, the spectral radius of $A ^H A$ (and $A A ^H$) is simply the largest eigenvalue of $A^HA$, and we can write ```math \| A \|_2 = \sqrt{\lambda_{max} ( A ^H A )} = \sqrt{\lambda_{max} ( A A ^H)} ``` """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$9457901e-3b05-4c35-ab14-db1399079d79„§cell_idÙ$9457901e-3b05-4c35-ab14-db1399079d79¤codeÚ3md""" ## Equivalence of norms !!! note "Definition (Norm equivalence)" The norms $\| \cdot \|_p$ and $\| \cdot \|_q$ over a space $V$ are equivalent if there exists positive constants $\alpha, \beta > 0$ such that for all $A \in V$ ```math \alpha \| A \|_p \leq \| A \|_q \leq \beta \| A \|_p ``` """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$09fe122e-4487-4b1c-9c03-4f160e7e3857„§cell_idÙ$09fe122e-4487-4b1c-9c03-4f160e7e3857¤codeÚºmd""" ## Similarity !!! note "Definition (Similarity and diagonalisability)" $A \in \mathbb C^{n \times n}$ is *similar* to $B \in \mathbb C^{n \times n}$ if there exists an invertible $S \in \mathbb C^{n \times n}$ such that $S^{-1} A S = B$. If $S$ is *unitary*, i.e. $S^{-1} = S^H$, then $A$ and $B$ are *unitarily similar*. $A \in \mathbb C^{n \times n}$ is *diagonalisable* if it is similar to a diagonal matrix. !!! tip "Remark (Eigenvalues of similar matrices)" Similar matrices have equal eigenvalues of equal multiplicities. If $(\lambda, v)$ is an eigenpair of $A$ and $B = S^{-1} A S$, then $(\lambda, S^{-1} v)$ is an eigenpair of $B$. Indeed, since ```math B S^{-1} v = S^{-1} A S S^{-1} v = S^{-1} A v = \lambda S^{-1} v. ``` Thus, the eigenvalues of $A$ and $B$ are equal and their multiplicities the same. !!! tip "Remark" If $A$ and $B$ are unitarily similar and $A$ is Hermitian, so is $B$, since ```math B = U^{-1} A U = U^H A U. ``` !!! tip "Remark" $A \in \mathbb C^{n \times n}$ is diagonalizable if and only is there exists $n$ linearly independent eigenvectors $v_i, i = 1, \dots n$. To see this, we collect all eigenvectors as columns ```math S = \begin{bmatrix} v_1 & v_2 & \dots & v_n \\ \downarrow & \downarrow & & \downarrow \end{bmatrix} ``` implying $AS = S \Lambda$ with $\Lambda = \text{diag} (\lambda_1, \dots, \lambda_n)$. Since $v_i$ are linearly independent, $S$ is invertible, thus $S^{-1} A S = \Lambda$. !!! tip "Remark (Invariance of trace)" The trace is invariant under a similarity transform. Indeed, noticing $\text{tr} (ABC) = \text{tr} (BCA)$, we obtain ```math \text{tr} ( S ^{-1} A S ) = \text{tr} (A S S^{-1}) = \text{tr} (A) ``` """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$d5cb2d7c-a627-4b79-a45e-f11ffe3b15a9„§cell_idÙ$d5cb2d7c-a627-4b79-a45e-f11ffe3b15a9¤codeÚ¼md""" !!! note "Definitions (Trace, kernel, image)" Let $A \in \mathbb C^{n \times n}$. We further define the following quantities ```math \begin{align} &\text{Trace} & \text{tr} (A) &\coloneqq \sum_i A_{ii} \\ &\text{Kernel} & \ker (A) &\coloneqq \{ v \in \mathbb C^n \mid Av = 0\} \\ &\text{Image} & \operatorname{im} (A) &\coloneqq \{ w \in \mathbb C^n \mid \exists v \in \mathbb C^n \text{ s.t. } Av = w \} \end{align} ``` """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$af5521f4-e4a4-4411-972e-c41bc7a85fea„§cell_idÙ$af5521f4-e4a4-4411-972e-c41bc7a85fea¤codeÚmd""" Another important matrix norm is the *Frobenius norm* ```math \| A \|_F = \sqrt{\mathrm{tr} ( A A ^H)} = \sqrt{\sum_{i,j = 1}^n | A _{ij}^2|} ``` Its relevance will become clear once we consider the equivalence of matrix norms, which we treat in the next section. """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$901b5296-373c-498b-890a-245b190e68b8„§cell_idÙ$901b5296-373c-498b-890a-245b190e68b8¤codeÙ=md"Show outline $(@bind show_outline CheckBox(default=true))"¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$a730c823-e0f4-48ed-bf6b-996f8de49460„§cell_idÙ$a730c823-e0f4-48ed-bf6b-996f8de49460¤codeÚömd""" !!! note "Proposition 2 (Equivalence of Frobenius norm and 2-norm)" We consider the vector space of matrices $A \in \mathbb C^{n \times n}$. The Frobenius norm $\| \cdot \|_F$ and the spectral norm $\| \cdot \|_2$ are equivalent. Moreover, ```math \| A \|_2 \stackrel{(a)}{\leq} \| A \|_F \stackrel{(b)}{\leq} \sqrt{n} \| A \|_2 ``` This relationship is relevant as $\| A \|_2$ proves useful in theory but is hard to compute, unlike $\| A \|_F$ for which computation is quite straightforward. > *Proof.* > (a) is left as an exercise. For (b), we have > ```math > \begin{align} > \| A \|_F^2 = \text{tr} ( A ^H A ) = \sum_{i=1}^n \lambda_i( A ^H A) \leq \sum_{i=1}^n \lambda_{max} ( A ^H A ) = n \| A \|_2^2 > && \square > \end{align} > ``` """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$be0c5a02-d3f9-447d-83bd-dbc5aba9266d„§cell_idÙ$be0c5a02-d3f9-447d-83bd-dbc5aba9266d¤code±TableOfContents()¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$f24d7aff-7cc9-4ae4-b757-ed170aaef555„§cell_idÙ$f24d7aff-7cc9-4ae4-b757-ed170aaef555¤codeÙímd""" Some matrix concepts that may be useful to have handy for the course. For more details, see [Youssef Saad. *Numerical Methods for Large Eigenvalue Problems*, SIAM (2011)](https://www-users.cse.umn.edu/~saad/eig_book_2ndEd.pdf). """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$891ad571-27cb-483c-a0a3-4b53ad9d6c17„§cell_idÙ$891ad571-27cb-483c-a0a3-4b53ad9d6c17¤codeÙ´md""" !!! tip "Remark" $\ker(A)$ and $\operatorname{im} (A)$ are subspaces of $\mathbb C^n$, and ```math \dim ( \ker(A)) + \dim (\operatorname{im}(A)) = n. ``` """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$cf35b757-1783-4316-977c-b182aa10ebb5„§cell_idÙ$cf35b757-1783-4316-977c-b182aa10ebb5¤codeÚ€md""" - Recall the conditions a norm should satisfy, for all $A \in \mathbb C^{n\times n}$ : ```math \begin{align} \| A \| &\geq 0 && \| A \| = 0 \iff A = 0 \\ \| \alpha A \| &= | \alpha| \| A \| && \forall \alpha \in \mathbb C \\ \| A + B \| &\leq \| A \| + \| B \| && \forall B \in \mathbb C^{n \times n} \tag{$\triangle$} \end{align} ``` The third condition is known as the triangle inequality and sometimes denoted as $\triangle$. - The $p$-norms of $\mathbb C^n$ ($\| v \|_p = (\sum_{i=1}^n |v_i|^p)^{1/p} \quad \forall v \in \mathbb C^n$) inherit a matrix norm for $p \geq 1$ or $p = \infty$ as ```math \| A \|_p = \max_{0 \neq v \in \mathbb C^n} \frac{\| A v \|_p}{\|v\|_p}, ``` the so-called *induced matrix norm*. - They are always consistent with the associated vector norm, i.e. ```math \| A v \|_p = \| v \|_p \frac{\| A v \|_p }{\|v\|_p} \leq \| v\|_p \max_{0 \neq v \in \mathbb C^n} \frac{\| A v \|_p}{\| v \|_p} = \| A \|_p \| v \|_p. ``` Further, we obtain ```math \| A B \|_p = \max_{x \neq v \in \mathbb C^n} \frac{\| A B v \|_p}{\| v \|_p} \leq \| A \|_p \max_{0 \neq v \in \mathbb C^n } \frac{\| B v \|_p}{\| v\|_p} = \| A \|_p \| B \|_p ``` Thus, ```math \begin{align} \| A ^2 \|_p &\leq \| A \|_p^2. \end{align} ``` - More generally, ```math \begin{align} \| A ^k \|_p &\leq \| A \|_p^k. \end{align} ``` """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$9cc8bd7c-7deb-4db9-9175-6d636c6294fe„§cell_idÙ$9cc8bd7c-7deb-4db9-9175-6d636c6294fe¤code½md""" ## Kernel and Image """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$976f90d6-eb84-4a45-8d51-8822a283ac35„§cell_idÙ$976f90d6-eb84-4a45-8d51-8822a283ac35¤codeÙtmd""" ## Matrix norms Clearly, $A \in \mathbb C^{n \times n}$ is a vector space. What are good norms for it ? """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$877bcf5d-741b-4cfc-b069-ba033eecdcd5„§cell_idÙ$877bcf5d-741b-4cfc-b069-ba033eecdcd5¤codeÙJbegin using PlutoUI using HypertextLiteral using PlutoTeachingTools end¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$008bb52c-c638-11f1-89c7-731f0c90a37d„§cell_idÙ$008bb52c-c638-11f1-89c7-731f0c90a37d¤codeÙ'md""" # Revision of matrix concepts """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$450839de-2597-4221-8e68-cd2f6b441e0c„§cell_idÙ$450839de-2597-4221-8e68-cd2f6b441e0c¤codeÚ�md""" Here are some explicit definitions for induced matrix norms ```math \begin{align} \| A \|_1 &= \max_{j = 1, \dots, n} \sum_{i=1}^n | A_{ij} | \tag{Column-sum norm} \\ \| A \|_\infty &= \max_{i = 1, \dots, n} \sum_{j=1}^n |A_{ij}| \tag{Row-sum norm} \\ \| A\|_2 &= \sqrt{ϱ (A A^H)} = \sqrt{ϱ(A^H A)} \end{align} ``` Here, $A^H$ denotes the Hermitian conjugate, i.e. $A^H = \overline{ A ^T} = (\bar A ) ^T$ and $ϱ$ is the spectral radius of $A$. For those more familiar with linear algebra over the reals, Hermitian matrices ($A^H = A$) generalize many properties of real symmetric matrices, which are themselves all Hermitian. """¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedÃÙ$5259f8a1-2879-45b4-9e74-a2ff0feaff79„§cell_idÙ$5259f8a1-2879-45b4-9e74-a2ff0feaff79¤codeÚHif show_outline RobustLocalResource("https://teaching.matmat.org/math-materials/sidebar.md", "sidebar.md") Sidebar(toc, ypos) = @htl("""""") Sidebar(Markdown.parse(read("sidebar.md", String)), 265) end¨metadataƒ©show_logsèdisabled®skip_as_script«code_foldedënotebook_idÙ$36e2f216-d814-11f1-9cad-d19f17168ef1¥bonds€¬cell_resultsÞÙ$47690b4e-8df6-4d07-86e2-1399ef75fb3eЦqueued¤logs�§running¦output†¤bodyÚÎ

Theorem 1

Let $A \in \mathbb C^{n \times n}$. $A$ is injective if and only if $\ker (A) = \{ 0 \}$.

Proof.

$\boxed \Rightarrow$ : Let $v \in \ker (A)$. Then

$$\begin{align} Av = 0 &\Rightarrow A(-v) = - Av = 0 \\ &\hspace{-0.8em} \stackrel{\text{injective}}{\Rightarrow} v = -v \\ &\Rightarrow v = 0 \end{align}$$

$\boxed \Leftarrow$ : Let $A v_1 = A v_2$. Then

$$\begin{align} A(v_1 - v_2) = A v_1 - A v_2 = 0 &\Rightarrow v_1 - v_2 \in \ker (A) \\&\hspace{-1.4em} \stackrel{\ker (A) = \{ 0\}}{\Rightarrow} v_1 - v_2 = 0 \\ &\Rightarrow v_1 = v_2 \\ & && & \square \end{align}$$

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²))0�·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$47690b4e-8df6-4d07-86e2-1399ef75fb3e¹depends_on_disabled_cells§runtimeΗµpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$3e592d32-9e4a-45e6-b179-b7b72542a60cЦqueued¤logs�§running¦output†¤bodyÚÔ

Remark

$\| A \|_2$ is also equal to the largest singular value of $A$. In addition, since $A^HA$ (and $A A ^H$) is a Hermitian positive semi-definite matrix, the spectral radius of $A ^H A$ (and $A A ^H$) is simply the largest eigenvalue of $A^HA$, and we can write

$$\| A \|_2 = \sqrt{\lambda_{max} ( A ^H A )} = \sqrt{\lambda_{max} ( A A ^H)}$$

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)*ß·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$3e592d32-9e4a-45e6-b179-b7b72542a60c¹depends_on_disabled_cells§runtimeÎȵpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$9457901e-3b05-4c35-ab14-db1399079d79Цqueued¤logs�§running¦output†¤bodyÚ+

Equivalence of norms

Definition (Norm equivalence)

The norms $\| \cdot \|_p$ and $\| \cdot \|_q$ over a space $V$ are equivalent if there exists positive constants $\alpha, \beta > 0$ such that for all $A \in V$

$$ \alpha \| A \|_p \leq \| A \|_q \leq \beta \| A \|_p$$

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)+Æ�·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$9457901e-3b05-4c35-ab14-db1399079d79¹depends_on_disabled_cells§runtimeÎ}[µpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$09fe122e-4487-4b1c-9c03-4f160e7e3857Цqueued¤logs�§running¦output†¤bodyÚ Ï

Similarity

Definition (Similarity and diagonalisability)

$A \in \mathbb C^{n \times n}$ is similar to $B \in \mathbb C^{n \times n}$ if there exists an invertible $S \in \mathbb C^{n \times n}$ such that $S^{-1} A S = B$.

If $S$ is unitary, i.e. $S^{-1} = S^H$, then $A$ and $B$ are unitarily similar.

$A \in \mathbb C^{n \times n}$ is diagonalisable if it is similar to a diagonal matrix.

Remark (Eigenvalues of similar matrices)

Similar matrices have equal eigenvalues of equal multiplicities. If $(\lambda, v)$ is an eigenpair of $A$ and $B = S^{-1} A S$, then $(\lambda, S^{-1} v)$ is an eigenpair of $B$. Indeed, since

$$ B S^{-1} v = S^{-1} A S S^{-1} v = S^{-1} A v = \lambda S^{-1} v.$$

Thus, the eigenvalues of $A$ and $B$ are equal and their multiplicities the same.

Remark

If $A$ and $B$ are unitarily similar and $A$ is Hermitian, so is $B$, since

$$ B = U^{-1} A U = U^H A U.$$

Remark

$A \in \mathbb C^{n \times n}$ is diagonalizable if and only is there exists $n$ linearly independent eigenvectors $v_i, i = 1, \dots n$. To see this, we collect all eigenvectors as columns

$$ S = \begin{bmatrix} v_1 & v_2 & \dots & v_n \\ \downarrow & \downarrow & & \downarrow \end{bmatrix}$$

implying $AS = S \Lambda$ with $\Lambda = \text{diag} (\lambda_1, \dots, \lambda_n)$. Since $v_i$ are linearly independent, $S$ is invertible, thus $S^{-1} A S = \Lambda$.

Remark (Invariance of trace)

The trace is invariant under a similarity transform. Indeed, noticing $\text{tr} (ABC) = \text{tr} (BCA)$, we obtain

$$ \text{tr} ( S ^{-1} A S ) = \text{tr} (A S S^{-1}) = \text{tr} (A)$$

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)-q·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$09fe122e-4487-4b1c-9c03-4f160e7e3857¹depends_on_disabled_cells§runtimeÎÛ¶µpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$d5cb2d7c-a627-4b79-a45e-f11ffe3b15a9Цqueued¤logs�§running¦output†¤bodyÚá

Definitions (Trace, kernel, image)

Let $A \in \mathbb C^{n \times n}$. We further define the following quantities

$$\begin{align} &\text{Trace} & \text{tr} (A) &\coloneqq \sum_i A_{ii} \\ &\text{Kernel} & \ker (A) &\coloneqq \{ v \in \mathbb C^n \mid Av = 0\} \\ &\text{Image} & \operatorname{im} (A) &\coloneqq \{ w \in \mathbb C^n \mid \exists v \in \mathbb C^n \text{ s.t. } Av = w \} \end{align}$$

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)'ö ·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$d5cb2d7c-a627-4b79-a45e-f11ffe3b15a9¹depends_on_disabled_cells§runtimeÎåµpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$af5521f4-e4a4-4411-972e-c41bc7a85feaЦqueued¤logs�§running¦output†¤bodyÚ‹

Another important matrix norm is the Frobenius norm

$$\| A \|_F = \sqrt{\mathrm{tr} ( A A ^H)} = \sqrt{\sum_{i,j = 1}^n | A _{ij}^2|}$$

Its relevance will become clear once we consider the equivalence of matrix norms, which we treat in the next section.

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)+r÷·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$af5521f4-e4a4-4411-972e-c41bc7a85fea¹depends_on_disabled_cells§runtimeÎôjµpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$901b5296-373c-498b-890a-245b190e68b8Цqueued¤logs�§running¦output†¤bodyÙ‰

Show outline

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²W._·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$901b5296-373c-498b-890a-245b190e68b8¹depends_on_disabled_cells§runtimeÎÀPµpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$a730c823-e0f4-48ed-bf6b-996f8de49460Цqueued¤logs�§running¦output†¤bodyÚ¿

Proposition 2 (Equivalence of Frobenius norm and 2-norm)

We consider the vector space of matrices $A \in \mathbb C^{n \times n}$. The Frobenius norm $\| \cdot \|_F$ and the spectral norm $\| \cdot \|_2$ are equivalent. Moreover,

$$ \| A \|_2 \stackrel{(a)}{\leq} \| A \|_F \stackrel{(b)}{\leq} \sqrt{n} \| A \|_2$$

This relationship is relevant as $\| A \|_2$ proves useful in theory but is hard to compute, unlike $\| A \|_F$ for which computation is quite straightforward.

Proof. (a) is left as an exercise. For (b), we have

$$\begin{align} \| A \|_F^2 = \text{tr} ( A ^H A ) = \sum_{i=1}^n \lambda_i( A ^H A) \leq \sum_{i=1}^n \lambda_{max} ( A ^H A ) = n \| A \|_2^2 && \square \end{align}$$

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²),·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$a730c823-e0f4-48ed-bf6b-996f8de49460¹depends_on_disabled_cells§runtimeÎvѵpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$be0c5a02-d3f9-447d-83bd-dbc5aba9266dЦqueued¤logs�§running¦output†¤bodyÚP¾ °persist_js_state¤mime©text/html²last_run_timestampËAÚº²Vê�·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$be0c5a02-d3f9-447d-83bd-dbc5aba9266d¹depends_on_disabled_cells§runtimeÍt„µpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$f24d7aff-7cc9-4ae4-b757-ed170aaef555Цqueued¤logs�§running¦output†¤bodyÚ!

Some matrix concepts that may be useful to have handy for the course. For more details, see Youssef Saad. Numerical Methods for Large Eigenvalue Problems, SIAM (2011).

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)' ·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$f24d7aff-7cc9-4ae4-b757-ed170aaef555¹depends_on_disabled_cells§runtimeÎާµpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$891ad571-27cb-483c-a0a3-4b53ad9d6c17Цqueued¤logs�§running¦output†¤bodyÚ—

Remark

$\ker(A)$ and $\operatorname{im} (A)$ are subspaces of $\mathbb C^n$, and

$$\dim ( \ker(A)) + \dim (\operatorname{im}(A)) = n.$$

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)(ÿ°·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$891ad571-27cb-483c-a0a3-4b53ad9d6c17¹depends_on_disabled_cells§runtimeθµpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$cf35b757-1783-4316-977c-b182aa10ebb5Цqueued¤logs�§running¦output†¤bodyÚ'
°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)*b··has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$cf35b757-1783-4316-977c-b182aa10ebb5¹depends_on_disabled_cells§runtimeÎ ±µpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$9cc8bd7c-7deb-4db9-9175-6d636c6294feЦqueued¤logs�§running¦output†¤bodyÙL

Kernel and Image

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)'Ú­·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$9cc8bd7c-7deb-4db9-9175-6d636c6294fe¹depends_on_disabled_cells§runtimeÎ@µpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$976f90d6-eb84-4a45-8d51-8822a283ac35Цqueued¤logs�§running¦output†¤bodyÙÆ

Matrix norms

Clearly, $A \in \mathbb C^{n \times n}$ is a vector space. What are good norms for it ?

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²))Ñþ·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$976f90d6-eb84-4a45-8d51-8822a283ac35¹depends_on_disabled_cells§runtimeÎö,µpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$877bcf5d-741b-4cfc-b069-ba033eecdcd5Цqueued¤logs�§running¦output†¤body °persist_js_state¤mimeªtext/plain²last_run_timestampËAÚº²Veš¦·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$877bcf5d-741b-4cfc-b069-ba033eecdcd5¹depends_on_disabled_cells§runtimeΞ±µpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$008bb52c-c638-11f1-89c7-731f0c90a37dЦqueued¤logs�§running¦output†¤bodyÙb

Revision of matrix concepts

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)'¥Î·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$008bb52c-c638-11f1-89c7-731f0c90a37d¹depends_on_disabled_cells§runtimeÎ,Ÿµpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$450839de-2597-4221-8e68-cd2f6b441e0cЦqueued¤logs�§running¦output†¤bodyÚ

Here are some explicit definitions for induced matrix norms

$$\begin{align} \| A \|_1 &= \max_{j = 1, \dots, n} \sum_{i=1}^n | A_{ij} | \tag{Column-sum norm} \\ \| A \|_\infty &= \max_{i = 1, \dots, n} \sum_{j=1}^n |A_{ij}| \tag{Row-sum norm} \\ \| A\|_2 &= \sqrt{ϱ (A A^H)} = \sqrt{ϱ(A^H A)} \end{align}$$

Here, $A^H$ denotes the Hermitian conjugate, i.e. $A^H = \overline{ A ^T} = (\bar A ) ^T$ and $ϱ$ is the spectral radius of $A$. For those more familiar with linear algebra over the reals, Hermitian matrices ($A^H = A$) generalize many properties of real symmetric matrices, which are themselves all Hermitian.

°persist_js_state¤mime©text/html²last_run_timestampËAÚº²)*�4·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$450839de-2597-4221-8e68-cd2f6b441e0c¹depends_on_disabled_cells§runtimeξIµpublished_object_keys�¸depends_on_skipped_cells§erroredÂÙ$5259f8a1-2879-45b4-9e74-a2ff0feaff79Цqueued¤logs�§running¦output†¤bodyÚO°persist_js_state¤mime©text/html²last_run_timestampËAÚº²WV¢¶·has_pluto_hook_features¬rootassigneeÀ§cell_idÙ$5259f8a1-2879-45b4-9e74-a2ff0feaff79¹depends_on_disabled_cells§runtimeÎ9/’µpublished_object_keys�¸depends_on_skipped_cells§errored©shortpath¾90_Revision_matrix_concepts.jl®last_save_timeËAÚº²)&7‚«in_temp_dir¨metadata€