Computing the Pseudoinverse of Specific Toeplitz Matrices Using Rank-One Updates

<div>ISNR versus <svg height="11.927pt" id="M345" style="vertical-align:-3.291101pt" version="1.1" viewbox="-0.0498162 -8.6359 9.83233 11.927" width="9.83233pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M442 569C442 628 407 660 362 660C328 660 292 648 255 616C200 568 165 508 123 334L106 264C80 242 52 222 23 202L30 176C54 191 76 205 100 222L87 166C73 103 79 43 95 19C108 0 139 -12 169 -12C242 -12 318 38 369 106L354 129C322 93 260 49 215 49C173 49 147 81 170 187L189 274C293 344 442 452 442 569ZM380 565C380 478 267 377 196 322L216 406C254 556 282 623 332 623C364 623 380 596 380 565Z" id="g113-127"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="465" vert-adv-y="465"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,5.46,3.132)"><path d="M387 400C387 425 348 451 303 451C247 451 176 414 132 376C69 322 24 228 24 148C24 43 74 -12 147 -12C211 -12 301 33 363 103L346 128C319 99 249 51 193 51C148 51 112 84 112 165C112 230 130 287 154 330C170 359 199 400 243 400C277 400 304 383 326 354C333 345 343 343 354 348C378 360 387 382 387 400Z" id="g50-100"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="411" vert-adv-y="411"></glyph.data></g></svg> for the removal of Gaussian white noise with zero mean and variance 0.05 and Gaussian blur (<svg height="11.927pt" id="M346" style="vertical-align:-3.291101pt" version="1.1" viewbox="-0.0498162 -8.6359 37.4286 11.927" width="37.4286pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M442 569C442 628 407 660 362 660C328 660 292 648 255 616C200 568 165 508 123 334L106 264C80 242 52 222 23 202L30 176C54 191 76 205 100 222L87 166C73 103 79 43 95 19C108 0 139 -12 169 -12C242 -12 318 38 369 106L354 129C322 93 260 49 215 49C173 49 147 81 170 187L189 274C293 344 442 452 442 569ZM380 565C380 478 267 377 196 322L216 406C254 556 282 623 332 623C364 623 380 596 380 565Z" id="g113-127"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="465" vert-adv-y="465"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,5.46,3.132)"><path d="M397 380C406 395 404 411 396 425S369 451 350 451C302 451 239 372 192 294H189L199 338C214 407 207 451 180 451C152 451 83 405 30 345L48 318C87 354 117 372 125 372S135 362 127 324L55 -5L61 -12C87 -5 117 3 139 5C154 87 168 162 179 207C198 250 240 310 260 332C281 355 297 366 307 366C321 366 333 360 347 346C351 342 360 342 369 348C379 355 389 366 397 380Z" id="g50-115"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="427" vert-adv-y="427"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,13.49,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,24.753,0)"><path d="M285 378C315 398 338 416 353 432C373 451 384 474 384 503C384 579 325 635 236 635H235C182 635 136 610 108 579L65 516L85 496C110 533 150 575 205 575C258 575 300 543 300 481C300 407 232 369 141 339L147 310C163 315 188 321 211 321C268 321 338 284 338 192C338 94 288 40 217 40C160 40 119 68 93 91C85 98 77 97 69 91C60 84 47 71 46 58C44 46 48 35 62 22C75 10 116 -12 162 -12C234 -12 424 62 424 224C424 297 373 359 285 376V378Z" id="g113-52"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,30.993,0)"><path d="M153 550H386L412 615L406 623H120L82 318C104 327 142 338 184 338C294 338 347 275 347 187C347 112 305 39 221 39C160 39 119 71 97 89C88 97 80 96 71 90C59 80 50 67 49 57C48 45 52 36 66 23C80 9 123 -12 169 -12C221 -11 288 15 342 59C403 109 431 165 431 225C431 308 366 395 238 395C212 395 165 379 127 364L153 550Z" id="g113-54"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>, <svg height="12.5794pt" id="M347" style="vertical-align:-3.29107pt" version="1.1" viewbox="-0.0498162 -9.28833 41.4346 12.5794" width="41.4346pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M352 391C352 416 319 448 267 448C236 448 173 423 147 400C107 364 96 332 96 304C96 248 143 210 193 181C241 153 258 124 258 100C258 72 232 38 184 38C151 38 107 66 81 108C77 114 64 116 55 111C34 99 23 84 23 65C23 29 81 -12 134 -12C220 -12 325 61 325 141C325 184 297 215 234 256C194 282 161 309 161 346C161 380 188 401 217 401C255 401 279 380 301 353C308 344 313 341 325 347C341 355 352 371 352 391Z" id="g113-116"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="375" vert-adv-y="375"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,8.507,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,19.77,0)"><path d="M442 569C442 628 407 660 362 660C328 660 292 648 255 616C200 568 165 508 123 334L106 264C80 242 52 222 23 202L30 176C54 191 76 205 100 222L87 166C73 103 79 43 95 19C108 0 139 -12 169 -12C242 -12 318 38 369 106L354 129C322 93 260 49 215 49C173 49 147 81 170 187L189 274C293 344 442 452 442 569ZM380 565C380 478 267 377 196 322L216 406C254 556 282 623 332 623C364 623 380 596 380 565Z" id="g113-127"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="465" vert-adv-y="465"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,25.23,3.132)"><path d="M397 380C406 395 404 411 396 425S369 451 350 451C302 451 239 372 192 294H189L199 338C214 407 207 451 180 451C152 451 83 405 30 345L48 318C87 354 117 372 125 372S135 362 127 324L55 -5L61 -12C87 -5 117 3 139 5C154 87 168 162 179 207C198 250 240 310 260 332C281 355 297 366 307 366C321 366 333 360 347 346C351 342 360 342 369 348C379 355 389 366 397 380Z" id="g50-115"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="427" vert-adv-y="427"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,29.628,0)"><path d="M368 703H309L44 -163H104L368 703Z" id="g113-48"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="412" vert-adv-y="412"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,34.99,0)"><path d="M412 140C382 77 369 73 315 73H129L270 222C362 320 402 379 402 466C402 571 322 635 234 635C177 635 130 609 99 576L42 495L64 475C90 514 133 568 201 568C274 568 318 519 318 435C318 349 255 267 193 193C144 135 87 78 32 23V0H405C417 45 427 89 440 131L412 140Z" id="g113-51"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>, and <svg height="12.5794pt" id="M348" style="vertical-align:-3.29107pt" version="1.1" viewbox="-0.0498162 -9.28833 41.2884 12.5794" width="41.2884pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M352 391C352 416 319 448 267 448C236 448 173 423 147 400C107 364 96 332 96 304C96 248 143 210 193 181C241 153 258 124 258 100C258 72 232 38 184 38C151 38 107 66 81 108C77 114 64 116 55 111C34 99 23 84 23 65C23 29 81 -12 134 -12C220 -12 325 61 325 141C325 184 297 215 234 256C194 282 161 309 161 346C161 380 188 401 217 401C255 401 279 380 301 353C308 344 313 341 325 347C341 355 352 371 352 391Z" id="g113-116"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="375" vert-adv-y="375"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,8.507,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,19.77,0)"><path d="M442 569C442 628 407 660 362 660C328 660 292 648 255 616C200 568 165 508 123 334L106 264C80 242 52 222 23 202L30 176C54 191 76 205 100 222L87 166C73 103 79 43 95 19C108 0 139 -12 169 -12C242 -12 318 38 369 106L354 129C322 93 260 49 215 49C173 49 147 81 170 187L189 274C293 344 442 452 442 569ZM380 565C380 478 267 377 196 322L216 406C254 556 282 623 332 623C364 623 380 596 380 565Z" id="g113-127"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="465" vert-adv-y="465"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,25.23,3.132)"><path d="M387 400C387 425 348 451 303 451C247 451 176 414 132 376C69 322 24 228 24 148C24 43 74 -12 147 -12C211 -12 301 33 363 103L346 128C319 99 249 51 193 51C148 51 112 84 112 165C112 230 130 287 154 330C170 359 199 400 243 400C277 400 304 383 326 354C333 345 343 343 354 348C378 360 387 382 387 400Z" id="g50-100"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="411" vert-adv-y="411"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,29.482,0)"><path d="M368 703H309L44 -163H104L368 703Z" id="g113-48"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="412" vert-adv-y="412"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,34.843,0)"><path d="M412 140C382 77 369 73 315 73H129L270 222C362 320 402 379 402 466C402 571 322 635 234 635C177 635 130 609 99 576L42 495L64 475C90 514 133 568 201 568C274 568 318 519 318 435C318 349 255 267 193 193C144 135 87 78 32 23V0H405C417 45 427 89 440 131L412 140Z" id="g113-51"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>) for Barbara image.</div>

Mathematical Problems in Engineering

fig5

Figure 5

Figure 5: Computing the Pseudoinverse of Specific Toeplitz Matrices Using Rank-One Updates 