hw_ch2_ch2pr1.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "source": [
    "### Inverse Problems and Parameter Estimation, GEOS 627/427, University of Alaska Fairbanks\n",
    "\n",
    "- script hw_ch2_ch2pr1.ipynb"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "%run lib_header.py"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy.matlib\n",
    "from matplotlib import path\n",
    "from numpy.linalg import inv\n",
    "#import scipy.special\n",
    "from scipy.stats import chi2 as chi_2   # chi2pdf\n",
    "\n",
    "from lib_geos import * # plot_histo, corrcov\n",
    "from lib_inverse import plot_ellipse\n",
    "from lib_peip import chi2cdf\n",
    "from lib_peip import chi2inv"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "ndata = 6     # number of data points\n",
    "nparm = 2     # number of model parameters\n",
    "\n",
    "# vector of distances\n",
    "# start of refraction\n",
    "xstart = 5 + 1/3\n",
    "# end of refraction\n",
    "xmax = 5*xstart\n",
    "x = np.linspace(np.ceil(xstart),xmax,ndata)\n",
    "x = x.reshape((ndata,1))\n",
    "\n",
    "# vector of first arrival times\n",
    "t = np.array([3.4935,4.2853,5.1374,5.8181,6.8632,8.1841])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "\n",
    "\n"
   ]
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "7f341b321178e62fbbfacfb91167ed733f9907f6f0692efdceb8394e8303ba73"
  },
  "kernelspec": {
   "display_name": "inverse [conda env:inverse]",
   "language": "python",
   "name": "conda-env-inverse-py"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}