{ "metadata": { "name": "", "signature": "sha256:f4049af615a127fe3d3a56cfb6ee4c15193f7fb38b0c9b80194f314ccc094105" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "$$T = \\frac{1}{1+F\\sin^2(\\delta/2)}$$\n", "$$ F = \\frac{4R}{(1-R)^2}$$\n", "$$\\delta = \\frac{2\\pi}{\\lambda}2nl\\cos\\theta$$ is the phase difference between successive transmitted pair" ] }, { "cell_type": "code", "collapsed": false, "input": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "from ipywidgets import StaticInteract, RangeWidget" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "x = np.linspace(-4*np.pi,4*np.pi,501)\n", "def plot(r):\n", " fig, ax = plt.subplots(figsize=(5, 5),\n", " subplot_kw={'axisbg':'#EEEEEE',\n", " 'axisbelow':True})\n", " ax.grid(color='w', linewidth=2, linestyle='solid')\n", " F = 4*r/(1-r)**2\n", " ax.plot(x,1/(1+F*np.sin(x/2)**2), lw=5, alpha=0.4)\n", " #ax.set_xlim(-50, 50)\n", " ax.set_title('plot of transmittance vs phase and reflectance')\n", " ax.set_xlabel('delta')\n", " ax.set_ylabel('transmittance')\n", " ax.set_ylim(0, 1)\n", " return fig" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 10 }, { "cell_type": "code", "collapsed": false, "input": [ "StaticInteract(plot,r=RangeWidget(0.0, 1.0, 0.1))" ], "language": "python", "metadata": {}, "outputs": [ { "html": [ "\n", " \n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", "
\n", " \n", " r: \n", "
\n", " " ], "metadata": {}, "output_type": "pyout", "prompt_number": 11, "text": [ "" ] } ], "prompt_number": 11 }, { "cell_type": "markdown", "metadata": {}, "source": [ "to-do\n", "===\n", "* need to convert xticks to show delta/phase " ] }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }