{
  "nbformat": 4,
  "nbformat_minor": 0,
  "metadata": {
    "colab": {
      "provenance": []
    },
    "kernelspec": {
      "name": "python3",
      "display_name": "Python 3"
    },
    "language_info": {
      "name": "python"
    }
  },
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "7au-B_pAh4pL"
      },
      "source": [
        "**Documente o seu programa**\n",
        "\n",
        "\n",
        "\n"
      ]
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "1I1YN3Uv23Fi"
      },
      "source": [
        "import matplotlib.pyplot as plt\n",
        "import numpy as np\n",
        "\n",
        "from google.colab import files"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "IxVX_Xqf3C1X"
      },
      "source": [
        "#faça o upload dos arquivos com os valores de tempo-distância das ondas direta e refratadas 1 e 2:\n",
        "# files.upload() "
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "VqJtfuXd3SYi"
      },
      "source": [
        "# Complete os comandos a seguir:\n",
        " = np.loadtxt(\" \", unpack=True)\n",
        " = np.loadtxt(\" \", unpack=True)\n",
        " = np.loadtxt(\" \", unpack=True)\n"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "Vf21JxBkf1fO"
      },
      "source": [
        "# Apresente graficamente os dados dos arquivos importados acima\n",
        "\n",
        "\n"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "code",
      "metadata": {
        "id": "-z6sv3B95Se9"
      },
      "source": [
        "### Implemente a \"Inversão\" dos tempos de percurso:\n",
        "#1. Utilize a função \"polyfit\" para ajustar uma reta aos tempos dos eventos sísmicos  \n",
        "#2. A partir dos coeficientes (linear e angular) das retas, implemente as fórmulas para obtenção das velocidades e espessuras das camadas. \n",
        "#3. Faça um print dos resultados, com a descrição do significado dos valores calculados"
      ],
      "execution_count": null,
      "outputs": []
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "wJsod13y6lNy"
      },
      "source": [
        "**Para ajustar uma reta (a + bx) pelo método dos \"mínimos quadrados\"**\n",
        "\n",
        "A função **polyfit** do **numpy** efetua o ajuste de polinômios por mínimos quadrados  (c0 + c1x + c2x^2 + c3x^3 + ... +cnx^n): uma reta é um polinômio de grau 1 \n",
        "\n",
        "A função **polyfit** retorna os coeficientes do polinômio na ordem do maior (grau do polinômio) para o menor (a constante, ou o coeficiente de índice zero) índice.\n",
        "\n",
        "`b,a = np.polyfit(x, y, 1)`"
      ]
    }
  ]
}