https://mahurin.us/robtasm/enContents © 2026 <a href="mailto:robert.s.mahurin@gmail.com">rob mahurin</a> Sat, 31 Jan 2026 03:07:36 GMTNikola (getnikola.com)http://blogs.law.harvard.edu/tech/rsshttps://mahurin.us/robtasm/posts/2025/07/31/its-only-fitting/rob mahurin<figure><img src="https://mahurin.us/robtasm/posts/2025/07/31/its-only-fitting/assets/signal_fitted.svg"></figure> <div><p>Here's a quick-n-dirty Gaussian curve fitter, which I seem to reinvent about twice a year. There are a number of canned solutions that I can never remember how to use, but I also frequently find myself wanting to fit weird functions.</p> <p>Below in a minimum working example. But the short-short version is</p> <div class="code"><pre class="code literal-block"><span class="k">def</span><span class="w"> </span><span class="nf">func</span><span class="p">(</span><span class="n">x</span><span class="p">,</span> <span class="o">*</span><span class="n">params</span><span class="p">):</span> <span class="o">...</span> <span class="n">x</span><span class="p">,</span><span class="n">y</span> <span class="o">=</span> <span class="n">read_some_data</span><span class="p">()</span> <span class="n">guess_params</span> <span class="o">=</span> <span class="p">[</span><span class="o">...</span><span class="p">]</span> <span class="kn">from</span><span class="w"> </span><span class="nn">scipy</span><span class="w"> </span><span class="kn">import</span> <span class="n">optimize</span> <span class="n">better_params</span><span class="p">,</span> <span class="n">covariance</span> <span class="o">=</span> <span class="n">optimize</span><span class="o">.</span><span class="n">curve_fit</span><span class="p">(</span> <span class="n">func</span><span class="p">,</span> <span class="n">x</span><span class="p">,</span> <span class="n">y</span><span class="p">,</span> <span class="n">p0</span><span class="o">=</span><span class="n">guess_params</span><span class="p">)</span> </pre></div> <p><a href="https://mahurin.us/robtasm/posts/2025/07/31/its-only-fitting/">Read more…</a> (2 min remaining to read)</p></div>analysiscurvespythonscipyhttps://mahurin.us/robtasm/posts/2025/07/31/its-only-fitting/Thu, 31 Jul 2025 19:19:00 GMT