Analysis Control Ø Visualization and reconstruction of acoustic field Ø Estimation of source locations and room-acoustic parameters Ø Spatial sound field recording Ø High-fidelity spatial audio reproduction Ø Directivity control and local reproduction Ø Spatial active noise control From theory to application of signal processing and inverse problems for acoustic fields Microphone Loudspeaker
Sound field reconstruction based on harmonic analysis of infinite order [Ueno+ IEEE SPL 2018] • Optimization algorithm for sparse representation of acoustic field [Murata+ IEEE TSP 2018] • Sparse sound field decomposition in reverberant environment [Koyama+ IEEE JSTSP 2019] • Separation of internal and external sound fields [Takida+ EUSIPCO 2018] • Estimation of source parameters based on Reciprocity Gap Functional [Takida+ Elsevier SP 2019] Analysis inside region without sources Analysis inside region including sources
Sound field control based on weighted mode-matching [Ueno+ IEEE/ACM TASLP 2019] • Sound field recording and reproduction in wave-number domain [Koyama+ IEEE(/ACM) TASLP 2013, 2014, JASA 2016] • Super-resolution in recording and reproduction [Koyama+ IEEE JSTSP 2015, JASA 2018] • Optimization of source and sensor placement for sound field control [Koyama+ IEEE/ACM TASLP 2020, IEEE ICASSP 2018] • Spatial active noise control based on kernel interpolation [Ito+ IEEE ICASSP 2019]
Estimation of expansion coeff at arbitrary point March 26, 2020 7 [Ueno+ IEEE SPL 2018] <latexit sha1_base64="TnKjsQx5LCjG+RPMUHDUz7skgmU=">AAAD3HicdZPNjtMwEIDdhp+lwNKFA0hcLCqkVoIqASS4VKqWC5eVFonuLqrbauK6rWnsRLYDLVE4IG6IK+/GG/AYTPoDdNtaSjT6Zr5M4ozDJJLW+f6vUtm7cvXa9YMblZu3bh/eqR7dPbNxarjo8DiKzUUIVkRSi46TLhIXiRGgwkich9PXRf78ozBWxvqdmyeip2Cs5UhycIgG1a9pnYUqM3mDtiizqRpkTKctP+9nTOqRm+drqtLWU0wVCbxTBlEygZU88Bv0w8LM61PTz+p+Ax/4fkn6hZvTOnMT4WCZfMKSiVzVDao1v+kvFt0OglVQI6t1OjgqH7NhzFMltOMRWNsN/MT1MjBO8kjkFZZakQCfwlh0MdSghO1li83K6WMkQzqKDV7a0QX938hAWTtXIVYqcBN7OVfAXblu6kavepnUSeqE5stGozSiLqbFztOhNIK7aI4BcCPxXSmfgAHu8P9sdAkVfoMWn3isFOhhxsCMFcxy3EjsGicZM4oi+1JAFkklnd1hSL3DQLjXsGmybRRwr4Ejsm0UcLeB3RF/FpvOX7pHgtkuaU3/SQqmAvAEONzQCuPtE5idgDNyhsfAtoLn64oYB9Hg1AWXZ2w7OHvWDPxm8PZFrX28mr8D8pA8InUSkJekTd6QU9IhnPwuHZbulx54fe+b9937sSwtl1bOPbKxvJ9/AF68T7Y=</latexit> Spherical wavefunction Expansion coeff at center <latexit sha1_base64="D9HWprNBRjV9vUhnlaU3sm3abQc=">AAADvXicdZLNjtMwEIC9G36W8teFIxeLCqkrRJUAElxWVMsFDisVie5Wqks1cZ3Gqu1EtsO2ROb9eASegivccNLuQretpUjjb+bTJJOJc8GNDcOfe/vBjZu3bh/cady9d//Bw+bhozOTFZqyPs1EpgcxGCa4Yn3LrWCDXDOQsWDn8ex9lT//yrThmfpsFzkbSZgqnnAK1qNxk5IUbEliWRIQeQrO4XZ10+4IH+OaD7i7QkSwxLaXvGe4w8898t0mULOPDhPNp6k9+lK+iFzNjBs3W2EnrA/eDKJV0EKr0xsf7p+QSUYLyZSlAowZRmFuRyVoy6lgrkEKw3KgM5iyoQ8VSGZGZT0Nh595MsFJpv2jLK7p/0YJ0piFjH2lBJua67kKbssNC5u8HZVc5YVlii4bJYXANsPVaPGEa0atWPgAqOb+XTFNQQO1/gesdYml/wbFLmgmJaiJn72eSpi7klRds7wkWmLPvleQCC65NVsMrrYYHu40TJFvGhXcaXCVbBoV3G747h5/Y+vOFd0hwXybdEn/SRJmDPyKWz/QBqHdU5ifgtV87vfcHEevLisymzLtty66vmObwdnLThR2ok+vW92T1f4doCfoKWqjCL1BXfQB9VAfUfQD/UK/0Z/gXcACEahl6f7eynmM1k5w8RcdekYS</latexit> <latexit sha1_base64="8B+ZNBL5UQkvn4J1WUEojaP77iI=">AAADXnicdZJfaxNBEMC3OW1rtPaPL4Ivh0HwKdypoC9CaF98KVQwaSEXytxmLlmzf47duZr0SD+Dr/rNfPOjuJemaprLwMLwm/kxu8ukuRSOoujXViN48HB7Z/dR8/GTvaf7B4dHPWcKy7HLjTT2IgWHUmjskiCJF7lFUKnE83RyUtXPr9A6YfQXmuU4UDDSIhMcyKNekqrSzi8PWlE7WkS4nsTLpMWWcXZ52DhOhoYXCjVxCc714yinQQmWBJc4byaFwxz4BEbY96kGhW5QLq47D195MgwzY/3RFC7o/0YJyrmZSn2nAhq7+7UK1tX6BWUfBqXQeUGo+e2grJAhmbB6ezgUFjnJmU+AW+HvGvIxWODkf2hlSqr8GzR+40Yp0MMyATtSMJ2XSTXV5GViVejZTQUTKZQgV2MIXWN4uNFwRb5uVHCjIXS2blSw3vDTPb7GVecv3SDBtE66o/8kBRMEv4PkP7SZ8M4pTE+BrJj6RXQf47d3HYbGaP3Wxfd3bD3pvWnHUTv+/K7VOV7u3y57wV6y1yxm71mHfWJnrMs4+8q+sx/sZ+N3sB3sBfu3rY2tpfOMrUTw/A+VXSIZ</latexit> Estimates of expansion coeff at <latexit sha1_base64="8B+ZNBL5UQkvn4J1WUEojaP77iI=">AAADXnicdZJfaxNBEMC3OW1rtPaPL4Ivh0HwKdypoC9CaF98KVQwaSEXytxmLlmzf47duZr0SD+Dr/rNfPOjuJemaprLwMLwm/kxu8ukuRSOoujXViN48HB7Z/dR8/GTvaf7B4dHPWcKy7HLjTT2IgWHUmjskiCJF7lFUKnE83RyUtXPr9A6YfQXmuU4UDDSIhMcyKNekqrSzi8PWlE7WkS4nsTLpMWWcXZ52DhOhoYXCjVxCc714yinQQmWBJc4byaFwxz4BEbY96kGhW5QLq47D195MgwzY/3RFC7o/0YJyrmZSn2nAhq7+7UK1tX6BWUfBqXQeUGo+e2grJAhmbB6ezgUFjnJmU+AW+HvGvIxWODkf2hlSqr8GzR+40Yp0MMyATtSMJ2XSTXV5GViVejZTQUTKZQgV2MIXWN4uNFwRb5uVHCjIXS2blSw3vDTPb7GVecv3SDBtE66o/8kBRMEv4PkP7SZ8M4pTE+BrJj6RXQf47d3HYbGaP3Wxfd3bD3pvWnHUTv+/K7VOV7u3y57wV6y1yxm71mHfWJnrMs4+8q+sx/sZ+N3sB3sBfu3rY2tpfOMrUTw/A+VXSIZ</latexit> Matrix product of infinite orders can be analytically derived, and dependency on expansion center is eliminated Observations <latexit sha1_base64="M1+7y5WxOHiRASf8C2J/vkHg7g0=">AAADlXicdZJdixMxFIazHXXX+tXVCy+8CRZhvSkzq+CCFMoqsjeLFexuoVPLmTTThuZjSDLaOoz/xV/jrV76b8xMu2p3poHAy3PykORwooQzY33/917Du3Hz1v7B7eadu/fuP2gdPrwwKtWEDojiSg8jMJQzSQeWWU6HiaYgIk4vo8Wbon75mWrDlPxoVwkdC5hJFjMC1qFJ63UYiSzsG5Z3yzRk+VERdD7xn3/KQgF2buLsLM+r5Umr7Xf8cuFqCDahjTarPzlsnIZTRVJBpSUcjBkFfmLHGWjLCKd5M0wNTYAsYEZHLkoQ1Iyz8pc5fubIFMdKuy0tLun/RgbCmJWI3Mny1ddrBayrjVIbn4wzJpPUUknWF8Upx1bhomV4yjQllq9cAKKZeysmc9BArGvs1i2RcH+Q9AtRQoCcZiHomYBlvu6jSrJQC+zYtwKGnAlmTY3BZI3h4E7DpEnVKOBOg8m4ahSw3nC3O/yVbjt/6Q4JlnXSFf0nCVhQcKNrXUObIemdw/IcrGZLN7+mG7y4OqHsnGo3dcH1GauGi+NO4HeCDy/bvdPN/B2gJ+gpOkIBeoV66Az10QAR9B39QD/RL++x1/Xeeu/WRxt7G+cR2lre+z9lcTeP</latexit> <latexit sha1_base64="PC9pAlu3s2zNT+S7GE8ISzdeo94=">AAADYHicdZJLixNBEIB7Mz5ifOxGb3oZDIKnMOMKehHCevGysILZDWRCqOnUJE36MXTX7CYO8Ud41T/m1V9iTzarZjMpaCi+qo/qbirNpXAURb8OGsGdu/fuNx+0Hj56/OTwqP303JnCcuxzI40dpOBQCo19EiRxkFsElUq8SOcfq/rFJVonjP5CyxxHCqZaZIIDeTRIUlXa1TgaH3WibrSOcDeJN0mHbeJs3G6cJBPDC4WauATnhnGU06gES4JLXLWSwmEOfA5THPpUg0I3KtcXXoWvPJmEmbH+aArX9H+jBOXcUqW+UwHN3O1aBetqw4Ky96NS6Lwg1Px6UFbIkExYvT6cCIuc5NInwK3wdw35DCxw8n+0NSVV/g0ar7hRCvSkTMBOFSxWZVJNNXmZWBV69q2CiRRKkKsxhK4xPNxruCLfNSq41xA62zUqWG/46R5/xW3nL90jwaJOuqH/JAVzBL+F5D+0lfDeKSxOgaxY+FV0H+Ljmw5DM7R+6+LbO7abnL/pxlE3/vy20zvZ7F+TvWAv2WsWs3esxz6xM9ZnnEn2nf1gPxu/g2ZwGLSvWxsHG+cZ24rg+R+OnyK8</latexit> Matrix of translation operators <latexit sha1_base64="PC9pAlu3s2zNT+S7GE8ISzdeo94=">AAADYHicdZJLixNBEIB7Mz5ifOxGb3oZDIKnMOMKehHCevGysILZDWRCqOnUJE36MXTX7CYO8Ud41T/m1V9iTzarZjMpaCi+qo/qbirNpXAURb8OGsGdu/fuNx+0Hj56/OTwqP303JnCcuxzI40dpOBQCo19EiRxkFsElUq8SOcfq/rFJVonjP5CyxxHCqZaZIIDeTRIUlXa1TgaH3WibrSOcDeJN0mHbeJs3G6cJBPDC4WauATnhnGU06gES4JLXLWSwmEOfA5THPpUg0I3KtcXXoWvPJmEmbH+aArX9H+jBOXcUqW+UwHN3O1aBetqw4Ky96NS6Lwg1Px6UFbIkExYvT6cCIuc5NInwK3wdw35DCxw8n+0NSVV/g0ar7hRCvSkTMBOFSxWZVJNNXmZWBV69q2CiRRKkKsxhK4xPNxruCLfNSq41xA62zUqWG/46R5/xW3nL90jwaJOuqH/JAVzBL+F5D+0lfDeKSxOgaxY+FV0H+Ljmw5DM7R+6+LbO7abnL/pxlE3/vy20zvZ7F+TvWAv2WsWs3esxz6xM9ZnnEn2nf1gPxu/g2ZwGLSvWxsHG+cZ24rg+R+OnyK8</latexit>
regression using spherical Bessel function as kernel function March 26, 2020 8 [Ueno+ IEEE SPL 2018, IWAENC 2018] <latexit sha1_base64="Nwb3n7CEyAyxNvvm5JH9L1bryVs=">AAADz3icdZLPb9MwFICdBtjo+NHBkYtFhdQKrUoG0rhUqsYFDpU6iW6T6i56cZ3WNHYi22EtWSau/ImIfwan7YCuraXIz9/z55c4L0xjro3n/XIq7oOHj/b2H1cPnjx99rx2+OJcJ5mirE+TOFGXIWgWc8n6hpuYXaaKgQhjdhFOP5b5i29MaZ7IL2aesqGAseQRp2AsCmoZmYDJs6JBQpGroonbmOhMBLlo+8VVF5OYRaaBF2nS07zAby2z548Al+xz0bzKj/xisdB2Unw8Mc1A4K+BhxtTcrM8+Gg5BYLcBMfNoFb3Wt5i4M3AXwV1tBq94LBySkYJzQSThsag9cD3UjPMQRlOY1ZUSaZZCnQKYzawoQTB9DBf3E+B31gywlGi7CMNXtD/jRyE1nMR2p0CzETfz5VwW26QmejDMOcyzQyTdFkoymJsElxeNh5xxaiJ5zYAqrh9V0wnoIAa+0vWqoTCfoNk1zQRAuQoJ6DGAmZFTsqqSZoTJbBltyUkMRfc6C0Gl1sMC3caOks3jRLuNLiMNo0SbjdsdYu/s3XnL90hwWybdEf/SQKmDGzTG3uhVUI7XZh1wSg+s52v2/67ux2JmTBlu86/32Obwflxy/da/tn7eud01X/76BV6jRrIRyeogz6hHuojin47jlN1Dtwz99q9dX8st1aclfMSrQ335x+c9kl3</latexit> Gram matrix consisting of spherical Bessel functions Observations Position of th mic <latexit sha1_base64="0uStd8PEwCtmctIFyVxU/PU2PFc=">AAADWXicdZLNattAEIDXVtu47l/SHHMRNYWejNQG2kvApJdcAgnEScAyYbQe2Yv3R+yOWrvCfYFc24crfZmuHKetY3lgYfhmPmZ3mTSXwlEU/Wo0g0ePn+y0nrafPX/x8tXu3utLZwrLsc+NNPY6BYdSaOyTIInXuUVQqcSrdPq5ql99QeuE0Rc0z3GoYKxFJjiQR+fqZrcTdaNlhJtJvEo6bBVnN3vN42RkeKFQE5fg3CCOchqWYElwiYt2UjjMgU9hjAOfalDohuXypovwrSejMDPWH03hkv5vlKCcm6vUdyqgiXtYq2BdbVBQ9mlYCp0XhJrfDcoKGZIJq2eHI2GRk5z7BLgV/q4hn4AFTv5z1qakyr9B41dulAI9KhOwYwWzRZlUU01eJlaFnn2vYCKFEuRqDKFrDA+3Gq7IN40KbjWEzjaNCtYbfrrH33Dd+Uu3SDCrk+7pP0nBFMGvH/kPbSe8dwqzUyArZn4H3VH84b7D0ASt37r44Y5tJpfvu3HUjc8PO73j1f612AF7w96xmH1kPXbCzlifcYbslv1gP5u/g0bQCtp3rc3GytlnaxHs/wFu7B+/</latexit> Kernel ridge regression with constraint that interpolation function satisfies Helmholtz eq. <latexit sha1_base64="E0VINkcXCNAwA1Jtzy3Bn53aBRg=">AAAELXicjZNNb9MwGMedhpcR3jo4crGoQONAlQDSxmFSNS5cKhWJbpPqKnIcpzWN7ch2RksWvhIfBnFAQlz5Gjh9GWubSliK9Nf/+f/0OPbjKEuZNr7/02m4N27eur13x7t77/6Dh839R6da5orQPpGpVOcR1jRlgvYNMyk9zxTFPErpWTR5V9XPLqjSTIqPZpbRIccjwRJGsLFW2PyOIl6gnmYlPIYeiuiIiSLi2Cg2Lb1PoX8wgeiyCqkyDF6uBLp8AZ9DRGJptBW7ct0qh5CHLpZBFK/EwrG1Dbb7nz2613t4iIr4atths+W3/fmC2yJYihZYrl643zhBsSQ5p8KQFGs9CPzMDAusDCMpLT2Ua5phMsEjOrBSYE71sJiffgmfWSeGiVT2EwbO3etEgbnWMx7ZpN3gWG/WKrOuNshNcjQsmMhyQwVZNEryFBoJq6uEMVOUmHRmBSaK2b1CMsYKE2MvfK1LxO0/CPqZSM6xPSmE1YjjaVmgqqvMCqQ4tN7XykQp48zoGoKJGsKaOwmdZ9tEZe4kmEi2icqsJ2x3a3+h68yVuwPC0zpo5f6DOJ5QbJ+UsQfqIdLp4ml3PmL2Xenj4PUqIc2YKjt1weaMbYvTV+3Abwcf3rQ6J8v52wNPwFNwAAJwCDrgPeiBPiDOWyd0xg5zv7k/3F/u70W04SyZx2BtuX/+AlX4a+A=</latexit>
Noise cancellation using adaptive filter and loudspeaker Ø Noise control at discrete positions in conventional method March 26, 2020 10 ANC inside continuous region: spatial ANC Loudspeaker Error mic Noise source Control region Reference mic [Ito+ IEEE ICASSP 2019]
regional noise power which is estimated by using kernel interpolation method for acoustic field March 26, 2020 12 <latexit sha1_base64="8sG3MC3lrsBMuNZaKxR7z4Tos08=">AAAEU3icdZPPb9MwFMedtsAoDDo4crGokLpL1WxIcJlUxmUHJopEt0lNV724TmvVdoLtjJYs3LjCX8eBPwVxwUnTbf1lKdLT532/fi/2sx9xpk2r9ccplSv37j/YeVh99Hj3ydPa3rMzHcaK0C4JeagufNCUM0m7hhlOLyJFQficnvuT91n+/IoqzUL52cwi2hcwkixgBIxFg9pfT4AZE+DJhxQfYY9JM0i8j4KOIMXXccPzRaLS/evLA5wrlUiG6RxmchvQ9DLJUzpITtI0Z51VxmlgGiu7Z0JvAlEEmRq0SRfV7qZu2U15XNT3FBuNzX5RsWgGe/RLzK62tvZuIRzU6q1mK194PXCLoI6K1RnslY69YUhiQaUhHLTuua3I9BNQhhFO06oXaxoBmcCI9mwoQVDdT/I7SvErS4Y4CJX9pME5vetIQGg9E75V5h2v5jK4KdeLTfC2nzAZxYZKMi8UxBybEGcXjodMUWL4zAZAFLO9YjIGBcTYsViq4gv7D5J+JaEQIIf2TtRIwDSdn2EYJZ4S2LLvGfQ4E8zoDQ4mNzgs3OrQcbTuyOBWB5PBuiODmx22usXf6LLnhm4xwXSTaUFvTQImFOzDM/ZAqx5pn8L0FIxiU/v69JF7uFCEZkyVnTp3dcbWg7ODpttqup9e19vHxfztoBfoJWogF71BbXSCOqiLiDNwfjg/nV/l3+V/lVKlMpeWnMLzHC2tyu5/lcOAww==</latexit> <latexit sha1_base64="Xjw/jwTwxeAj54v/lgakjUqgUCg=">AAADiXicdZLdbtMwFIC9hp/R8dPBJTcWFdIQoiQwiQkJqdpu4GJSkeg2qS7VieO0Vv2T2c5oicJb8DTcwkPwNjhdB3RNjxTp5PP5cuKjE2eCWxeGv7cawY2bt25v32nu3L13/0Fr9+GJ1bmhrE+10OYsBssEV6zvuBPsLDMMZCzYaTw9qs5PL5ixXKtPbp6xoYSx4imn4DwatV6SWBa9EhN2nvMLvFe9kp7lJX6OifDfSQBX7EP57HPxIipHrXbYCReB15NombTRMnqj3cYhSTTNJVOOCrB2EIWZGxZgHKeClU2SW5YBncKYDXyqQDI7LBY3K/FTTxKcauMf5fCC/m8UIK2dy9hXSnATe/2sgnVng9ylB8OCqyx3TNHLRmkusNO4GhNOuGHUiblPgBru/xXTCRigzg9zpUss/R0U+0K1lKCSgoAZS5iVBam66qwgRmLPvlWQCC65szUGVzWGhxsNm2frRgU3Glyl60YF6w3f3eOvbNX5SzdIMKuTrug/ScKUgV9X5wfaJLR7DLNjcIbP/M7ad9HrqwrtJsz4rYuu79h6cvKqE4Wd6ON+u3u43L9t9Bg9QXsoQm9QF71HPdRHFH1HP9BP9CvYCaLgIHh7WdrYWjqP0EoER38AQoIw7Q==</latexit> Gram matrix for sound field interpolation Kernel function for sound field interpolation Regional noise power Weighted NLMS can be applied Ø Not only on discrete points, but regional noise control. Ø Computational cost for filter update is equivalent to conventional method [Ito+ IEEE ICASSP 2019]
Sound field satisfies inhomogeneous Helmholtz eq. Source Microphone Reflected wave <latexit sha1_base64="VXB56x8cV3M2DCu+b66chBSNCqQ=">AAADi3icdZJfaxNBEMC3Of/UaG2qj74sBiFFDXepoIiFUBF8KbRg2kIuDXObvWS53b1jd04Tj/g1/DS+6mfw27iXpmqay8DC8Jv9MbvDRJkUFn3/91bNu3X7zt3te/X7D3Ye7jb2Hp3ZNDeM91gqU3MRgeVSaN5DgZJfZIaDiiQ/j5L3Zf38MzdWpPoTzjI+UDDWIhYM0KFhI2iFGiIJl53nyWVnn+atUAFOorgw8xfJPj2kL09X0bDR9Nv+Iuh6EiyTJlnGyXCvdhSOUpYrrpFJsLYf+BkOCjAomOTzephbngFLYMz7LtWguB0Ui7/N6TNHRjROjTsa6YL+bxSgrJ2pyN0sX2lv1kpYVevnGL8ZFEJnOXLNrhrFuaSY0nJQdCQMZyhnLgFmhHsrZRMwwNCNc6VLpNwfNP/CUqVAj4oQzFjBdF4s5pZmRWgUdexbCUMplEBbYQhdYTi40bB5tm6UcKMhdLxulLDacN0d/spXnb90gwTTKuma/pMUJBzcwqIbaD1k3WOYHgMaMXVbaw+Dg+sbKU64cVsX3Nyx9eSs0w78dnD6qtk9Wu7fNnlCnpIWCchr0iUfyQnpEUa+kx/kJ/nl7XgH3lvv3dXV2tbSeUxWwvvwB0omMJU=</latexit> Unknown boundary condition outside <latexit sha1_base64="okEXIY/pPSAfenmSSDCRtw4HgOI=">AAADanicdZLdahNREIBPs/7U+NO0XklvFqNQQcKuCnojhHrjTaGCaQvZEGZPZpNDzs9yzmxNXOKTeKvv5Dv4EJ5NUzXNZmBh+OZ8zOwwaS6Foyj6tdMIbt2+c3f3XvP+g4eP9lr7B2fOFJZjjxtp7EUKDqXQ2CNBEi9yi6BSiefp9ENVP79E64TRn2me40DBWItMcCCPhq1WcZQooEmalXbxcvoiHLbaUSdaRriZxKukzVZxOtxvHCcjwwuFmrgE5/pxlNOgBEuCS1w0k8JhDnwKY+z7VINCNyiXoy/C556MwsxY/2kKl/R/owTl3Fyl/mU1prtZq2BdrV9Q9m5QCp0XhJpfNcoKGZIJqz2EI2GRk5z7BLgVftaQT8ACJ7+ttS6p8v+g8Qs3SoEelQnYsYLZolwuzuRlYlXo2bcKJlIoQa7GELrG8HCr4Yp806jgVkPobNOoYL3hu3v8Fdedv3SLBLM66Zr+kxRMEfw9kl9oM+HdE5idAFkx80fp3sevr18YmqD1VxffvLHN5OxVJ4468ac37e7x6v522SF7yo5YzN6yLvvITlmPcXbJvrMf7Gfjd3AQPAkOr542dlbOY7YWwbM/H9wlzg==</latexit> Source distribution <latexit sha1_base64="ff/e1A/38XpRts4D3M/3UmeqwhY=">AAADXnicdZJfaxNBEMC3Oa01WvvHF8GXwyD4FO6soC9CaF/6Uqxg0kIulLnNXLJm/xy7c5p4pJ/BV/1mvvlR3EtTNc1lYGH4zfyY3WXSXApHUfRrqxHcu7/9YOdh89Hj3Sd7+weHPWcKy7HLjTT2MgWHUmjskiCJl7lFUKnEi3RyUtUvvqB1wuhPNMtxoGCkRSY4kEe95IPCEVztt6J2tIhwPYmXSYst4/zqoHGcDA0vFGriEpzrx1FOgxIsCS5x3kwKhznwCYyw71MNCt2gXFx3Hr70ZBhmxvqjKVzQ/40SlHMzlfpOBTR2d2sVrKv1C8reDUqh84JQ85tBWSFDMmH19nAoLHKSM58At8LfNeRjsMDJ/9DKlFT5N2j8yo1SoIdlAnakYDovk2qqycvEqtCz6womUihBrsYQusbwcKPhinzdqOBGQ+hs3ahgveGne/wNV52/dIME0zrplv6TFEwQ/A6S/9BmwjtnMD0DsmLqF9G9j49uOwyN0fqti+/u2HrSe92Oo3b88U2rc7zcvx32nL1gr1jM3rIOO2XnrMs4+8y+sx/sZ+N3sB3sBns3rY2tpfOUrUTw7A+oOyHJ</latexit> This ill-posed inverse problem requires some constraint on source distribution for reconstruction
spatially sparse and reconstruct sound field from observations Ø Group sparsity in time-frequency domain can also be taken into consideration March 26, 2020 16 <latexit sha1_base64="3unqeYcLzapR+kwzDMNoEVC+dbw=">AAADxHicdZJNb9NAEIY3MR8lfKVw5LIiQkJCRHZAopdKUUGIS6UikbZSnEbjzTpZsru2dtfgYNzfyJU/whXG+aBNk4xkafadeTTr2TdKpbDO93/X6t6t23fu7t1r3H/w8NHj5v6TU5tkhvEeS2RiziOwXArNe044yc9Tw0FFkp9F0/dV/ewbN1Yk+oubpXygYKxFLBg4lIbNr6ESWijxgw+LUIGbRHGRlyUNYwOsCMqig/nPVWVW0td0dfhQXgHYM+xcdOgrGkocPoJrUF5iMb1Ih82W3/bnQTeTYJm0yDJOhvv1o3CUsExx7ZgEa/uBn7pBAcYJJnnZCDPLU2BTGPM+phoUt4NivpSSvkBlROPE4KcdnavXiQKUtTMVYWd1UXuzVonbav3MxQeDQug0c1yzxaA4k9QltNowHQnDmZMzTIAZgXelbAK4TofvsDYlUvgPmn9niVKgR0UIZqwgLxcvkaRFaBRF7bISQ4nP5OwWQugtBIo7CZulm0Ql7iSEjjeJStxOrDy1zvxXd0CQb4NW6hWkYMoBne5woY2QdY8hPwZnRI52t4fBm1VH4ibcoOuCmx7bTE477cBvB5/ftrpHS//tkWfkOXlJAvKOdMknckJ6hJFf5A/5WyPeR0961ssWrfXaknlK1sK7/AfrOko1</latexit> <latexit sha1_base64="KdqboG5ErTF8Tsp8xTcQLhLEdoY=">AAADZXicdZLNbtNAEIC3MT8lUNoC4sKBFRFSuUQ2IMEBpKhcuFQqEmkrxVY13oyTVfaP3TUkWOE1uMJj8QS8Bus0BdI4I1kafTOfZj2a3AjufBz/2mpF167fuLl9q337zs7d3b39eydOl5Zhn2mh7VkODgVX2PfcCzwzFkHmAk/zybu6fvoZreNaffQzg5mEkeIFZ+ADyg5i+oYamgqkybPzvU7cjRdB15NkmXTIMo7P91uH6VCzUqLyTIBzgyQ2PqvAes4Ezttp6dAAm8AIByFVINFl1eLVc/o0kCEttA2f8nRB/zcqkM7NZB46Jfixu1qrYVNtUPridVZxZUqPil0MKkpBvab1CuiQW2RezEICzPLwVsrGYIH5sKiVKbkM/6DwC9NSghpWKdiRhOm8Suup2lSplTSwbzVMBZfcuwaDqwYjwI2GK826UcONBlfFulHDZiNMD/grrjp/6QYJpk3SJf0nSZgghFP0YaHtlPWOYHoE3vJpuEf3Nnlx2aH9GG24uuTqja0nJ8+7SdxNPrzs9A6X97dNHpEn5IAk5BXpkffkmPQJI5/Id/KD/Gz9jnaiB9HDi9bW1tK5T1YievwHCnAi1Q==</latexit> Sparsity-inducing penalty term Non-zero grid point <latexit sha1_base64="CMp0bdMsbHJXrQaEDWLvT3oAJJo=">AAADWXicdZLNattAEIDXVtu47l/SHHMRNYWejNQG2kvApJdcAgnEScAyYbQe2Yv3R+yOWrvCfYFc24crfZmuHKetY3lgYfhmPmZ3mTSXwlEU/Wo0g0ePn+y0nrafPX/x8tXu3utLZwrLsc+NNPY6BYdSaOyTIInXuUVQqcSrdPq5ql99QeuE0Rc0z3GoYKxFJjiQR+f6ZrcTdaNlhJtJvEo6bBVnN3vN42RkeKFQE5fg3CCOchqWYElwiYt2UjjMgU9hjAOfalDohuXypovwrSejMDPWH03hkv5vlKCcm6vUdyqgiXtYq2BdbVBQ9mlYCp0XhJrfDcoKGZIJq2eHI2GRk5z7BLgV/q4hn4AFTv5z1qakyr9B41dulAI9KhOwYwWzRZlUU01eJlaFnn2vYCKFEuRqDKFrDA+3Gq7IN40KbjWEzjaNCtYbfrrH33Dd+Uu3SDCrk+7pP0nBFMGvH/kPbSe8dwqzUyArZn4H3VH84b7D0ASt37r44Y5tJpfvu3HUjc8PO73j1f612AF7w96xmH1kPXbCzlifcYbslv1gP5u/g0bQCtp3rc3GytlnaxHs/wFx4R/A</latexit> <latexit sha1_base64="f0AM1dFjbF0IvWOeeDQ4Zux15bE=">AAADWXicdZLNattAEIDXVtu47l/SHHMRNYWejNQG2kvApJdcAgnEScAyYbQe2Yv3R+yOWrvCfYFc24crfZmuHKetY3lgYfhmPmZ3mTSXwlEU/Wo0g0ePn+y0nrafPX/x8tXu3utLZwrLsc+NNPY6BYdSaOyTIInXuUVQqcSrdPq5ql99QeuE0Rc0z3GoYKxFJjiQR+d0s9uJutEyws0kXiUdtoqzm73mcTIyvFCoiUtwbhBHOQ1LsCS4xEU7KRzmwKcwxoFPNSh0w3J500X41pNRmBnrj6ZwSf83SlDOzVXqOxXQxD2sVbCuNigo+zQshc4LQs3vBmWFDMmE1bPDkbDISc59AtwKf9eQT8ACJ/85a1NS5d+g8Ss3SoEelQnYsYLZokyqqSYvE6tCz75XMJFCCXI1htA1hodbDVfkm0YFtxpCZ5tGBesNP93jb7ju/KVbJJjVSff0n6RgiuDXj/yHthPeO4XZKZAVM7+D7ij+cN9haILWb138cMc2k8v33TjqxueHnd7xav9a7IC9Ye9YzD6yHjthZ6zPOEN2y36wn83fQSNoBe271mZj5eyztQj2/wCDnx/G</latexit> <latexit sha1_base64="HW0NkVNdxtG8JSmuoj2JEWTs8bY=">AAADWXicdZLNattAEIDXVtu47l/SHHMRNYWejNQG2kvApJdcAgnEScAyYbQe2Yv3R+yOWrvCfYFc24crfZmuHKetY3lgYfhmPmZ3mTSXwlEU/Wo0g0ePn+y0nrafPX/x8tXu3utLZwrLsc+NNPY6BYdSaOyTIInXuUVQqcSrdPq5ql99QeuE0Rc0z3GoYKxFJjiQR+fZzW4n6kbLCDeTeJV02CrObvaax8nI8EKhJi7BuUEc5TQswZLgEhftpHCYA5/CGAc+1aDQDcvlTRfhW09GYWasP5rCJf3fKEE5N1ep71RAE/ewVsG62qCg7NOwFDovCDW/G5QVMiQTVs8OR8IiJzn3CXAr/F1DPgELnPznrE1JlX+Dxq/cKAV6VCZgxwpmizKpppq8TKwKPftewUQKJcjVGELXGB5uNVyRbxoV3GoInW0aFaw3/HSPv+G685dukWBWJ93Tf5KCKYJfP/If2k547xRmp0BWzPwOuqP4w32HoQlav3Xxwx3bTC7fd+OoG58fdnrHq/1rsQP2hr1jMfvIeuyEnbE+4wzZLfvBfjZ/B42gFbTvWpuNlbPP1iLY/wNaOR+4</latexit> <latexit sha1_base64="0uStd8PEwCtmctIFyVxU/PU2PFc=">AAADWXicdZLNattAEIDXVtu47l/SHHMRNYWejNQG2kvApJdcAgnEScAyYbQe2Yv3R+yOWrvCfYFc24crfZmuHKetY3lgYfhmPmZ3mTSXwlEU/Wo0g0ePn+y0nrafPX/x8tXu3utLZwrLsc+NNPY6BYdSaOyTIInXuUVQqcSrdPq5ql99QeuE0Rc0z3GoYKxFJjiQR+fqZrcTdaNlhJtJvEo6bBVnN3vN42RkeKFQE5fg3CCOchqWYElwiYt2UjjMgU9hjAOfalDohuXypovwrSejMDPWH03hkv5vlKCcm6vUdyqgiXtYq2BdbVBQ9mlYCp0XhJrfDcoKGZIJq2eHI2GRk5z7BLgV/q4hn4AFTv5z1qakyr9B41dulAI9KhOwYwWzRZlUU01eJlaFnn2vYCKFEuRqDKFrDA+3Gq7IN40KbjWEzjaNCtYbfrrH33Dd+Uu3SDCrk+7pP0nBFMGvH/kPbSe8dwqzUyArZn4H3VH84b7D0ASt37r44Y5tJpfvu3HUjc8PO73j1f612AF7w96xmH1kPXbCzlifcYbslv1gP5u/g0bQCtp3rc3GytlnaxHs/wFu7B+/</latexit> [Koyama+ ICASSP 2014] <latexit sha1_base64="TcKdDousfow9G7Auxdzcf+KrT68=">AAADfnicdZLdahNBFICnWX9q/Ev10pvBqFSocdcWLEIh1BtvChVMG8iGcHYymwyZmV1mzurGZX0An8ZbfRPfxtk0UdNsDgwcvjMfZ+ZwolQKi77/e6fh3bh56/bunebde/cfPGztPbqwSWYY77FEJqYfgeVSaN5DgZL3U8NBRZJfRrP3Vf3yMzdWJPoTzlM+VDDRIhYM0KFR61moAKdRXPTL/XcHeBC/pCd0xfJyVDhUjlptv+Mvgm4mwTJpk2Wcj/Yap+E4YZniGpkEaweBn+KwAIOCSV42w8zyFNgMJnzgUg2K22Gx+E5JnzsypnFi3NFIF/R/owBl7VxF7mb1UHu9VsG62iDD+HhYCJ1myDW7ahRnkmJCq9nQsTCcoZy7BJgR7q2UTcEAQzfBtS6Rcn/Q/AtLlAI9LkIwEwV5WSxGl6RFaBR17FsFQymUQFtjCF1jOLjVsFm6aVRwqyF0vGlUsN5w3R3+ytedv3SLBHmdtKL/JAUzDm5H0Q20GbLuGeRngEbkblHtSXC4upHglBu3dcH1HdtMLt50Ar8TfDxqd0+X+7dLnpCnZJ8E5C3pkg/knPQII9/JD/KT/PKI98J75b2+utrYWTqPyVp4x38AGPEs1Q==</latexit>
sparse in time-frequency domain Ø Sparse optimization using mixed-norm penalty March 26, 2020 17 <latexit sha1_base64="jt8UZmLAo4CM+kOHoDXBXXv/2js=">AAAD3HicdZPfihMxFMbTjn/WqmtXLxS8CRZhxW6ZqYLeCGX1QoSFFexuodMdzqSZNnSSmU0y2hrHC/FOvPXdfAMfw0zbce22DQx8+XJ+nDMnJ2EaM6Vd93el6ly5eu36zo3azVu3d+/U9+6eqCSThHZJEieyF4KiMRO0q5mOaS+VFHgY09Nw8ro4P/1IpWKJ+KBnKR1wGAkWMQLaWkH9m8+ZYJx9poHxOehxGJlenmM/kkCMl5u21SrjgdHNyMovZdAsxwe43LzJSzUtYoL2WRs/xX5s6xjCIopAbN7lgUmb583zfP8i15Og3nBb7nzhdeEtRQMt13GwVz30hwnJOBWaxKBU33NTPTAgNSMxzWt+pmgKZAIj2rdSAKdqYObNyvFj6wxxlEj7CY3n7v+EAa7UjIc2sqhRXT4rzE1n/UxHLweGiTTTVJBFoiiLsU5w0Xk8ZJISHc+sACKZrRWTMdg2a3s/K1lCbv9B0E8k4RzE0PggRxym+eKGktT4kmPrfS1MP7bXp9UGgokNhDW3EipL14nC3EowEa0ThbmZKGdtlfnnboFgugkq3QuIw4SCfQHaNrTmk84RTI9ASza1z0C98p6VEYkeU2mnzrs8Y+vipN3y3Jb3/nmjc7icvx30ED1C+8hDL1AHvUXHqIsI+lPZrdyvPHDOnO/OD+fnIrRaWTL30Mpyfv0FTtFRkg==</latexit> <latexit sha1_base64="XMTa05C+9ds+LrBKIVtxIUVcJmU=">AAADbXicdZLNihNBEIB7M/6s8S/r4kmRxiCulzCjgh4UwnrxsrCC2V3IhFDTqUma9J/dPZo4xFfxqo/kU/gK9mSzajaThhmKr+qjuovKjODOx/GvnUZ05eq167s3mjdv3b5zt7V378TpwjLsMS20PcvAoeAKe557gWfGIshM4Gk2fVflTz+jdVyrj35ucCBhrHjOGfiAhq39g5i+oYamAumn5T95Nmy14068PHQzSFZBm6zO8XCvcZiONCskKs8EONdPYuMHJVjPmcBFMy0cGmBTGGM/hAokukG5vP2CPglkRHNtw6c8XdL/jRKkc3OZhUoJfuIu5ypYl+sXPn89KLkyhUfFzhvlhaBe02oUdMQtMi/mIQBmebgrZROwwHwY2FqXTIY3KPzCtJSgRmUKdixhtijTqqs2ZWolDexbBVPBJfeuxuCqxghwq+EKs2lUcKvBVb5pVLDeCN0D/orrzl+6RYJZnXRB/0kSpghhJX0YaDNl3SOYHYG3fBb20r1NXlxUaD9BG7Yuubxjm8HJ804Sd5IPL9vdw9X+7ZIH5DE5IAl5RbrkPTkmPcLInHwnP8jPxu/ofvQwenRe2thZOftk7URP/wA0wyUg</latexit> <latexit sha1_base64="jNNevYceU70kmgc9X62mIHeNJOc=">AAADwHicdZLdbtMwFIC9hp9R/jq45MaiQmqlqkoACW4mVeUGIU0aEt0qNV114jqtVdtJbAfaeeEJeQFeg1u4wGk7RtfWUqST75xPJz45UcqZNr7/86Di3bl77/7hg+rDR4+fPK0dPTvTSa4I7ZGEJ6ofgaacSdozzHDaTxUFEXF6Hs0+lPnzr1RplsgvZpHSoYCJZDEjYBwa1eJQgJkS4PZTMbJpK2tlRWPJotj2iyY+xqHOxUjikNPYNFZv1rTiAl/d1DVky6Hm1UWGQ8UmU9O8sGGsgNi0sFlRjGp1v+0vD94OgnVQR+tzOjqqdMNxQnJBpSEctB4EfmqGFpRhhNOiGuaapkBmMKEDF0oQVA/tciAFfuXIGMeJco80eEn/NywIrRcicpXlFfTtXAl35Qa5id8PLZNpbqgkq0ZxzrFJcDldPGaKEsMXLgCimPtWTKbg5mDcP9joEgl3B0m/kUQIkGMbgpoImBd2OdQktaES2LHvJQw5E8zoHQaTOwwH9xo6T7eNEu41mIy3jRLuNlx3hy/ppvOP7pFgvku6pjeSgBkFt+XGDbQaks4JzE/AKDZ3q66PgzfXFYmZUuW2Lri9Y9vB2et24LeDz2/rne56/w7RC/QSNVCA3qEO+ohOUQ8R9AP9Qr/RH6/rTb3Ey1allYO18xxtHO/yL3fVSK4=</latexit> <latexit sha1_base64="UO/CtIGxHZlCwhDll94dyV4o5hI=">AAADZHicdZJfaxNBEMC3Of/UWLW1+CTIYhB8KOGuCvVFCPXFl0IF0xZyIcxt5pIl+6+7ezXxiB/DV/1afgE/h3tpqqa5DBwMv5kfszdMZgR3Po5/bTWiO3fv3d9+0Hy48+jxk929p2dOF5Zhl2mh7UUGDgVX2PXcC7wwFkFmAs+zyYeqfn6F1nGtPvuZwb6EkeI5Z+ADSlMUYlCag8uDy/lgtxW340XQ9SRZJi2yjNPBXuM4HWpWSFSeCXCul8TG90uwnjOB82ZaODTAJjDCXkgVSHT9cvHoOX0VyJDm2oZPebqg/xslSOdmMgudEvzY3a5VsK7WK3z+rl9yZQqPil0PygtBvabVBuiQW2RezEICzPLwVsrGYIH5sKeVKZkM/6DwC9NSghqWKdiRhOm8TKup2pSplTSwbxVMBZfcuxqDqxojwI2GK8y6UcGNBlf5ulHBeiNMD/grrjp/6QYJpnXSDf0nSZgghEv0YaHNlHVOYHoC3vJpOEf3Pnlz06H9GG24uuT2ja0nZ4ftJG4nn962OsfL+9smz8lL8pok5Ih0yEdySrqEEUO+kx/kZ+N3tBPtR8+uWxtbS2efrET04g+IbSRa</latexit> Optimization algorithm based on Majorization minimization method [Murata+ IEEE TSP 2018] Acoustic signal in time- frequency domain
data – Reproduction of speech signal from loudspeaker March 26, 2020 18 Spatial aliasing artifacts originating from microphone intervals can be alleviated Proposed Method based on plane- wave decomposition [Koyama+ JASA 2018]
field inside using distributed loudspeakers March 26, 2020 20 <latexit sha1_base64="ff/e1A/38XpRts4D3M/3UmeqwhY=">AAADXnicdZJfaxNBEMC3Oa01WvvHF8GXwyD4FO6soC9CaF/6Uqxg0kIulLnNXLJm/xy7c5p4pJ/BV/1mvvlR3EtTNc1lYGH4zfyY3WXSXApHUfRrqxHcu7/9YOdh89Hj3Sd7+weHPWcKy7HLjTT2MgWHUmjskiCJl7lFUKnEi3RyUtUvvqB1wuhPNMtxoGCkRSY4kEe95IPCEVztt6J2tIhwPYmXSYst4/zqoHGcDA0vFGriEpzrx1FOgxIsCS5x3kwKhznwCYyw71MNCt2gXFx3Hr70ZBhmxvqjKVzQ/40SlHMzlfpOBTR2d2sVrKv1C8reDUqh84JQ85tBWSFDMmH19nAoLHKSM58At8LfNeRjsMDJ/9DKlFT5N2j8yo1SoIdlAnakYDovk2qqycvEqtCz6womUihBrsYQusbwcKPhinzdqOBGQ+hs3ahgveGne/wNV52/dIME0zrplv6TFEwQ/A6S/9BmwjtnMD0DsmLqF9G9j49uOwyN0fqti+/u2HrSe92Oo3b88U2rc7zcvx32nL1gr1jM3rIOO2XnrMs4+8y+sx/sZ+N3sB3sBns3rY2tpfOUrUTw7A+oOyHJ</latexit> Desired sound field Synthesized sound field Driving signal of th loudspeaker multiplied by transfer function Difficult to directly optimize because integration w.r.t. is included <latexit sha1_base64="ff/e1A/38XpRts4D3M/3UmeqwhY=">AAADXnicdZJfaxNBEMC3Oa01WvvHF8GXwyD4FO6soC9CaF/6Uqxg0kIulLnNXLJm/xy7c5p4pJ/BV/1mvvlR3EtTNc1lYGH4zfyY3WXSXApHUfRrqxHcu7/9YOdh89Hj3Sd7+weHPWcKy7HLjTT2MgWHUmjskiCJl7lFUKnEi3RyUtUvvqB1wuhPNMtxoGCkRSY4kEe95IPCEVztt6J2tIhwPYmXSYst4/zqoHGcDA0vFGriEpzrx1FOgxIsCS5x3kwKhznwCYyw71MNCt2gXFx3Hr70ZBhmxvqjKVzQ/40SlHMzlfpOBTR2d2sVrKv1C8reDUqh84JQ85tBWSFDMmH19nAoLHKSM58At8LfNeRjsMDJ/9DKlFT5N2j8yo1SoIdlAnakYDovk2qqycvEqtCz6womUihBrsYQusbwcKPhinzdqOBGQ+hs3ahgveGne/wNV52/dIME0zrplv6TFEwQ/A6S/9BmwjtnMD0DsmLqF9G9j49uOwyN0fqti+/u2HrSe92Oo3b88U2rc7zcvx32nL1gr1jM3rIOO2XnrMs4+8y+sx/sZ+N3sB3sBns3rY2tpfOUrUTw7A+oOyHJ</latexit>
loudspeakers March 26, 2020 21 Desired pressure distribution Transfer function matrix Loudspeaker driving signals - Discretize target region and obtain driving signals so that pressure on control points corresponds to desired one <latexit sha1_base64="ff/e1A/38XpRts4D3M/3UmeqwhY=">AAADXnicdZJfaxNBEMC3Oa01WvvHF8GXwyD4FO6soC9CaF/6Uqxg0kIulLnNXLJm/xy7c5p4pJ/BV/1mvvlR3EtTNc1lYGH4zfyY3WXSXApHUfRrqxHcu7/9YOdh89Hj3Sd7+weHPWcKy7HLjTT2MgWHUmjskiCJl7lFUKnEi3RyUtUvvqB1wuhPNMtxoGCkRSY4kEe95IPCEVztt6J2tIhwPYmXSYst4/zqoHGcDA0vFGriEpzrx1FOgxIsCS5x3kwKhznwCYyw71MNCt2gXFx3Hr70ZBhmxvqjKVzQ/40SlHMzlfpOBTR2d2sVrKv1C8reDUqh84JQ85tBWSFDMmH19nAoLHKSM58At8LfNeRjsMDJ/9DKlFT5N2j8yo1SoIdlAnakYDovk2qqycvEqtCz6womUihBrsYQusbwcKPhinzdqOBGQ+hs3ahgveGne/wNV52/dIME0zrplv6TFEwQ/A6S/9BmwjtnMD0DsmLqF9G9j49uOwyN0fqti+/u2HrSe92Oo3b88U2rc7zcvx32nL1gr1jM3rIOO2XnrMs4+8y+sx/sZ+N3sB3sBns3rY2tpfOUrUTw7A+oOyHJ</latexit> Difficult to densely measure transfer functions inside region <latexit sha1_base64="ff/e1A/38XpRts4D3M/3UmeqwhY=">AAADXnicdZJfaxNBEMC3Oa01WvvHF8GXwyD4FO6soC9CaF/6Uqxg0kIulLnNXLJm/xy7c5p4pJ/BV/1mvvlR3EtTNc1lYGH4zfyY3WXSXApHUfRrqxHcu7/9YOdh89Hj3Sd7+weHPWcKy7HLjTT2MgWHUmjskiCJl7lFUKnEi3RyUtUvvqB1wuhPNMtxoGCkRSY4kEe95IPCEVztt6J2tIhwPYmXSYst4/zqoHGcDA0vFGriEpzrx1FOgxIsCS5x3kwKhznwCYyw71MNCt2gXFx3Hr70ZBhmxvqjKVzQ/40SlHMzlfpOBTR2d2sVrKv1C8reDUqh84JQ85tBWSFDMmH19nAoLHKSM58At8LfNeRjsMDJ/9DKlFT5N2j8yo1SoIdlAnakYDovk2qqycvEqtCz6womUihBrsYQusbwcKPhinzdqOBGQ+hs3ahgveGne/wNV52/dIME0zrplv6TFEwQ/A6S/9BmwjtnMD0DsmLqF9G9j49uOwyN0fqti+/u2HrSe92Oo3b88U2rc7zcvx32nL1gr1jM3rIOO2XnrMs4+8y+sx/sZ+N3sB3sBns3rY2tpfOUrUTw7A+oOyHJ</latexit>
loudspeakers March 26, 2020 22 - Control points only on boundary of ? <latexit sha1_base64="ff/e1A/38XpRts4D3M/3UmeqwhY=">AAADXnicdZJfaxNBEMC3Oa01WvvHF8GXwyD4FO6soC9CaF/6Uqxg0kIulLnNXLJm/xy7c5p4pJ/BV/1mvvlR3EtTNc1lYGH4zfyY3WXSXApHUfRrqxHcu7/9YOdh89Hj3Sd7+weHPWcKy7HLjTT2MgWHUmjskiCJl7lFUKnEi3RyUtUvvqB1wuhPNMtxoGCkRSY4kEe95IPCEVztt6J2tIhwPYmXSYst4/zqoHGcDA0vFGriEpzrx1FOgxIsCS5x3kwKhznwCYyw71MNCt2gXFx3Hr70ZBhmxvqjKVzQ/40SlHMzlfpOBTR2d2sVrKv1C8reDUqh84JQ85tBWSFDMmH19nAoLHKSM58At8LfNeRjsMDJ/9DKlFT5N2j8yo1SoIdlAnakYDovk2qqycvEqtCz6womUihBrsYQusbwcKPhinzdqOBGQ+hs3ahgveGne/wNV52/dIME0zrplv6TFEwQ/A6S/9BmwjtnMD0DsmLqF9G9j49uOwyN0fqti+/u2HrSe92Oo3b88U2rc7zcvx32nL1gr1jM3rIOO2XnrMs4+8y+sx/sZ+N3sB3sBns3rY2tpfOUrUTw7A+oOyHJ</latexit> Significant degradation at several frequencies (Forbidden frequency problem) Can be alleviated by placing several control points inside target region, but to determine positions and number of these points is very difficult. <latexit sha1_base64="ff/e1A/38XpRts4D3M/3UmeqwhY=">AAADXnicdZJfaxNBEMC3Oa01WvvHF8GXwyD4FO6soC9CaF/6Uqxg0kIulLnNXLJm/xy7c5p4pJ/BV/1mvvlR3EtTNc1lYGH4zfyY3WXSXApHUfRrqxHcu7/9YOdh89Hj3Sd7+weHPWcKy7HLjTT2MgWHUmjskiCJl7lFUKnEi3RyUtUvvqB1wuhPNMtxoGCkRSY4kEe95IPCEVztt6J2tIhwPYmXSYst4/zqoHGcDA0vFGriEpzrx1FOgxIsCS5x3kwKhznwCYyw71MNCt2gXFx3Hr70ZBhmxvqjKVzQ/40SlHMzlfpOBTR2d2sVrKv1C8reDUqh84JQ85tBWSFDMmH19nAoLHKSM58At8LfNeRjsMDJ/9DKlFT5N2j8yo1SoIdlAnakYDovk2qqycvEqtCz6womUihBrsYQusbwcKPhinzdqOBGQ+hs3ahgveGne/wNV52/dIME0zrplv6TFEwQ/A6S/9BmwjtnMD0DsmLqF9G9j49uOwyN0fqti+/u2HrSe92Oo3b88U2rc7zcvx32nL1gr1jM3rIOO2XnrMs4+8y+sx/sZ+N3sB3sBns3rY2tpfOUrUTw7A+oOyHJ</latexit>
sound field March 26, 2020 23 [Ueno+ IEEE/ACM TASLP 2019] - Optimization problem including regional integral results in minimization problem w.r.t. expansion coefficients Expansion coeff of desired field Weighting matrix determined by regional integral Expansion coeff of transfer func
sound field March 26, 2020 24 [Ueno+ IEEE/ACM TASLP 2019] ü Truncation of expansion order can be avoided by weighting matrix ü Regional priority of control can be taken into consideration - Optimization problem including regional integral results in minimization problem w.r.t. expansion coefficients
loudspeakers and control points in pressure matching March 26, 2020 27 Select optimal placement from candidates w.r.t. control accuracy and filter stability
loudspeakers and control points in pressure matching March 26, 2020 28 - Optimal placement method based on Empirical Interpolation Method [Koyama+ IEEE ICASSP 2018] Sound field control problem is regarded as function interpolation problem: approximate desired sound field by using transfer functions as interpolation functions, control points as sampling points. Applied Empirical Interpolation Method proposed in numerical analysis of PDEs.
analysis of infinite orders and its application to spatial ANC, Reconstruction based on sparse decomposition – Control: Weighted mode-matching, Optimization of source and sensor placement for sound field control Ø Recent topics – Optimization of sensor placement for field estimation, Interpolation of region-to-region acoustic transfer function, Binaural reproduction from measurements of distributed microphones, Spatial ANC Ø Keywords – Kernel method, Gaussian process, Reproducing kernel Hilbert space, Sparse modeling, Adaptive filter, Convex optimization, Physical acoustics, Partial differential equation March 26, 2020 30