|
3 | 3 | { |
4 | 4 | "cell_type": "code", |
5 | 5 | "execution_count": null, |
| 6 | + "id": "cdc2b6fb", |
| 7 | + "metadata": {}, |
| 8 | + "outputs": [], |
| 9 | + "source": [ |
| 10 | + "using Pkg\n", |
| 11 | + "Pkg.activate(joinpath(dirname(pwd()), \"conf\", \"01_Ex4b\"))\n", |
| 12 | + "Pkg.instantiate();" |
| 13 | + ] |
| 14 | + }, |
| 15 | + { |
| 16 | + "cell_type": "code", |
| 17 | + "execution_count": 3, |
6 | 18 | "id": "42be4586", |
7 | 19 | "metadata": {}, |
8 | 20 | "outputs": [], |
|
12 | 24 | }, |
13 | 25 | { |
14 | 26 | "cell_type": "code", |
15 | | - "execution_count": null, |
| 27 | + "execution_count": 4, |
16 | 28 | "id": "cc1ff067", |
17 | 29 | "metadata": { |
18 | 30 | "scrolled": false |
19 | 31 | }, |
20 | | - "outputs": [], |
| 32 | + "outputs": [ |
| 33 | + { |
| 34 | + "data": { |
| 35 | + "text/plain": [ |
| 36 | + "3×2 Matrix{Float64}:\n", |
| 37 | + " 1.0 -0.5\n", |
| 38 | + " 0.0 0.0\n", |
| 39 | + " -2.0 1.0" |
| 40 | + ] |
| 41 | + }, |
| 42 | + "metadata": {}, |
| 43 | + "output_type": "display_data" |
| 44 | + } |
| 45 | + ], |
21 | 46 | "source": [ |
22 | 47 | "A = [1 -1/2; 0 0; -2 1]" |
23 | 48 | ] |
|
32 | 57 | }, |
33 | 58 | { |
34 | 59 | "cell_type": "code", |
35 | | - "execution_count": null, |
| 60 | + "execution_count": 5, |
36 | 61 | "id": "938a0cfb", |
37 | 62 | "metadata": {}, |
38 | | - "outputs": [], |
| 63 | + "outputs": [ |
| 64 | + { |
| 65 | + "data": { |
| 66 | + "text/plain": [ |
| 67 | + "SVD{Float64, Float64, Matrix{Float64}, Vector{Float64}}\n", |
| 68 | + "U factor:\n", |
| 69 | + "3×2 Matrix{Float64}:\n", |
| 70 | + " -0.447214 -0.894427\n", |
| 71 | + " 0.0 0.0\n", |
| 72 | + " 0.894427 -0.447214\n", |
| 73 | + "singular values:\n", |
| 74 | + "2-element Vector{Float64}:\n", |
| 75 | + " 2.5000000000000004\n", |
| 76 | + " 9.930136612989092e-17\n", |
| 77 | + "Vt factor:\n", |
| 78 | + "2×2 Matrix{Float64}:\n", |
| 79 | + " -0.894427 0.447214\n", |
| 80 | + " -0.447214 -0.894427" |
| 81 | + ] |
| 82 | + }, |
| 83 | + "metadata": {}, |
| 84 | + "output_type": "display_data" |
| 85 | + } |
| 86 | + ], |
39 | 87 | "source": [ |
40 | 88 | "(U,S,Vt) = svd(A)" |
41 | 89 | ] |
42 | 90 | }, |
43 | 91 | { |
44 | 92 | "cell_type": "code", |
45 | | - "execution_count": null, |
| 93 | + "execution_count": 6, |
46 | 94 | "id": "6ecb1c61", |
47 | 95 | "metadata": {}, |
48 | 96 | "outputs": [], |
|
54 | 102 | }, |
55 | 103 | { |
56 | 104 | "cell_type": "code", |
57 | | - "execution_count": null, |
| 105 | + "execution_count": 7, |
58 | 106 | "id": "da17abd5", |
59 | 107 | "metadata": {}, |
60 | | - "outputs": [], |
| 108 | + "outputs": [ |
| 109 | + { |
| 110 | + "data": { |
| 111 | + "text/plain": [ |
| 112 | + "2×3 Matrix{Float64}:\n", |
| 113 | + " 0.16 0.0 -0.32\n", |
| 114 | + " -0.08 0.0 0.16" |
| 115 | + ] |
| 116 | + }, |
| 117 | + "metadata": {}, |
| 118 | + "output_type": "display_data" |
| 119 | + } |
| 120 | + ], |
61 | 121 | "source": [ |
62 | 122 | "pinvA = pinv(A)" |
63 | 123 | ] |
64 | 124 | }, |
65 | 125 | { |
66 | 126 | "cell_type": "code", |
67 | | - "execution_count": null, |
| 127 | + "execution_count": 8, |
68 | 128 | "id": "f5df81eb", |
69 | 129 | "metadata": {}, |
70 | | - "outputs": [], |
| 130 | + "outputs": [ |
| 131 | + { |
| 132 | + "data": { |
| 133 | + "text/plain": [ |
| 134 | + "2×3 Matrix{Float64}:\n", |
| 135 | + " -4.02814e15 0.0 -2.01407e15\n", |
| 136 | + " 8.05628e15 0.0 4.02814e15" |
| 137 | + ] |
| 138 | + }, |
| 139 | + "metadata": {}, |
| 140 | + "output_type": "display_data" |
| 141 | + } |
| 142 | + ], |
71 | 143 | "source": [ |
72 | 144 | "bad_pinvA = V * Σ^(-1) * Ut" |
73 | 145 | ] |
|
82 | 154 | }, |
83 | 155 | { |
84 | 156 | "cell_type": "code", |
85 | | - "execution_count": null, |
| 157 | + "execution_count": 9, |
86 | 158 | "id": "b6a999b1", |
87 | 159 | "metadata": {}, |
88 | | - "outputs": [], |
| 160 | + "outputs": [ |
| 161 | + { |
| 162 | + "data": { |
| 163 | + "text/plain": [ |
| 164 | + "2×3 Matrix{Float64}:\n", |
| 165 | + " 0.16 0.0 -0.32\n", |
| 166 | + " -0.08 0.0 0.16" |
| 167 | + ] |
| 168 | + }, |
| 169 | + "metadata": {}, |
| 170 | + "output_type": "display_data" |
| 171 | + } |
| 172 | + ], |
89 | 173 | "source": [ |
90 | 174 | "pinv(A, rtol=1e-16)" |
91 | 175 | ] |
92 | 176 | }, |
93 | 177 | { |
94 | 178 | "cell_type": "code", |
95 | | - "execution_count": null, |
| 179 | + "execution_count": 10, |
96 | 180 | "id": "811e8e05", |
97 | 181 | "metadata": {}, |
98 | | - "outputs": [], |
| 182 | + "outputs": [ |
| 183 | + { |
| 184 | + "data": { |
| 185 | + "text/plain": [ |
| 186 | + "2×3 Matrix{Float64}:\n", |
| 187 | + " 4.02814e15 0.0 2.01407e15\n", |
| 188 | + " 8.05628e15 0.0 4.02814e15" |
| 189 | + ] |
| 190 | + }, |
| 191 | + "metadata": {}, |
| 192 | + "output_type": "display_data" |
| 193 | + } |
| 194 | + ], |
99 | 195 | "source": [ |
100 | 196 | "pinv(A, rtol=1e-17)" |
101 | 197 | ] |
102 | 198 | } |
103 | 199 | ], |
104 | 200 | "metadata": { |
105 | 201 | "kernelspec": { |
106 | | - "display_name": "Julia 1.7.2", |
| 202 | + "display_name": "Julia 1.11.4", |
107 | 203 | "language": "julia", |
108 | | - "name": "julia-1.7" |
| 204 | + "name": "julia-1.11" |
109 | 205 | }, |
110 | 206 | "language_info": { |
111 | 207 | "file_extension": ".jl", |
112 | 208 | "mimetype": "application/julia", |
113 | 209 | "name": "julia", |
114 | | - "version": "1.7.2" |
| 210 | + "version": "1.11.4" |
115 | 211 | } |
116 | 212 | }, |
117 | 213 | "nbformat": 4, |
|
0 commit comments