Skip to content

Instantly share code, notes, and snippets.

@jzuhone
Created August 7, 2014 17:06
Show Gist options
  • Save jzuhone/a80dff15b1e9686ab838 to your computer and use it in GitHub Desktop.
Save jzuhone/a80dff15b1e9686ab838 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"metadata": {
"name": "",
"signature": "sha256:641ca8140cb7761860f1ea28b7ea7445acddd67decee00c0750a8d398053312f"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"For individual point differences, we have"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"$s = (v_i - v_j)^2$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"so "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"$\\sigma_s^2 = \\left(\\frac{\\partial{s}}{\\partial{v_i}}\\right)^2\\sigma_{v_i}^2 + \\left(\\frac{\\partial{s}}{\\partial{v_j}}\\right)^2\\sigma_{v_j}^2$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"and "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"$\\left(\\frac{\\partial{s}}{\\partial{v_i}}\\right) = 2(v_i - v_j)$\n",
"\n",
"$\\left(\\frac{\\partial{s}}{\\partial{v_j}}\\right) = -2(v_i - v_j)$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"therefore "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"$\\sigma_s^2 = 4(v_i-v_j)^2\\sigma_{v_i}^2 + 4(v_i-v_j)^2\\sigma_{v_j}^2$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"but"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"$\\sigma_{v_i} = \\sigma_{v_j} = \\sigma_v$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"so "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"$\\sigma_s^2 = 8(v_i-v_j)^2\\sigma_v^2 = 8s\\sigma_v^2$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"but the structure function is not $s$, but the *average* of $s$:"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"$SF = \\langle(v_i-v_j)^2\\rangle$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"So I'm not sure how to take it past this point, e.g. how to compute the error on the bin from $\\sigma_s$."
]
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment