\begin{bmatrix} v_1 & \dot v_1 \\ v_2 & \dot v_2 \\ v_3 & \dot v_3 \\ \vdots & \vdots \end{bmatrix} \begin{bmatrix} m \\ b \end{bmatrix} = \begin{bmatrix} f_1\\f_2\\f_3 \\ \vdots \end{bmatrix} こんな感じになります.サンプリング数を$n_s$とすれば,上の式の行列の行数が$n_s$になります.この式を単純に $$Y\Theta=F$$ とまとめましょう.ここで$Y\in\mathbb{R}^{n_s\times 2}$は観測値行列,$\Theta\in\mathbb{R^2}$が求めたいパラメータの行列,$F\in\mathbb{R}^{n_s}$が入力値行列になります.$Y$と$F$は既知ですので,
![パラメータ推定(1)最小二乗法 - Qiita](https://cdn-ak-scissors.b.st-hatena.com/image/square/700b3d3ee739452b3e653a9af58d5d3b6e0c3c3d/height=288;version=1;width=512/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn.qiita.com%252Fassets%252Fpublic%252Farticle-ogp-background-412672c5f0600ab9a64263b751f1bc81.png%3Fixlib%3Drb-4.0.0%26w%3D1200%26mark64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTk3MiZoPTM3OCZ0eHQ9JUUzJTgzJTkxJUUzJTgzJUE5JUUzJTgzJUExJUUzJTgzJUJDJUUzJTgyJUJGJUU2JThFJUE4JUU1JUFFJTlBJUVGJUJDJTg4JUVGJUJDJTkxJUVGJUJDJTg5JUU2JTlDJTgwJUU1JUIwJThGJUU0JUJBJThDJUU0JUI5JTk3JUU2JUIzJTk1JnR4dC1jb2xvcj0lMjMyMTIxMjEmdHh0LWZvbnQ9SGlyYWdpbm8lMjBTYW5zJTIwVzYmdHh0LXNpemU9NTYmdHh0LWFsaWduPWxlZnQlMkN0b3Amcz0zZGE2ZjA4NGJlMTMyZjk4MmI0ZTNkMTA0ZjVmNjZkOA%26mark-x%3D142%26mark-y%3D57%26blend64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZoPTc2Jnc9NzcwJnR4dD0lNDB0YWthX2hvcmliZSZ0eHQtY29sb3I9JTIzMjEyMTIxJnR4dC1mb250PUhpcmFnaW5vJTIwU2FucyUyMFc2JnR4dC1zaXplPTM2JnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnM9OTc2YjNjZTcwYjFkODU2NGU3NWJhOTZhZTEzNzhiZGU%26blend-x%3D142%26blend-y%3D486%26blend-mode%3Dnormal%26s%3Dd824149d6513b34e44e5cf52f18b0473)