mathWithPython_harmonograph.ipynb

{"cells":[{"cell_type":"markdown","metadata":{"id":"Brtl6ZgJqgu2"},"source":["# COURSE: Master math by coding in Python\n","# SECTION: Introduction\n","# VIDEO: Create a beautiful harmonograph!\n","\n","\n","### https://www.udemy.com/course/math-with-python/?couponCode=202312\n","#### INSTRUCTOR: Mike X Cohen (http://sincxpress.com)\n","\n","This code roughly matches the code shown in the live recording: variable names, order of lines, and parameter settings may be slightly different."]},{"cell_type":"markdown","source":["<a target=\"_blank\" href=\"https://colab.research.google.com/github/mikexcohen/MathWithPython/blob/main/mathWithPython_harmonograph.ipynb\">\n","  <img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/>\n","</a>"],"metadata":{"id":"JZhUk_36xwBy"}},{"cell_type":"code","execution_count":null,"metadata":{"id":"-ybsgr5FbiZA"},"outputs":[],"source":["import numpy as np\n","import matplotlib.pyplot as plt"]},{"cell_type":"code","execution_count":null,"metadata":{"id":"9voiYnfRbiZC"},"outputs":[],"source":["# leave these parameters as-is\n","n = 1000000\n","t = np.logspace(np.log10(10),np.log10(500),n)\n","\n","# try chaning these parameters\n","A = [  1, 1, 1.5, 1.5 ]\n","d = [ .004, .001, .002, .0015 ]\n","f = [   3, 1, 2, 2.5 ]\n","\n","# generate XY value pairs\n","x = A[0]*np.sin(t*f[0])*np.exp(-d[0]*t) + A[1]*np.sin(t*f[1])*np.exp(-d[1]*t)\n","y = A[2]*np.sin(t*f[2])*np.exp(-d[2]*t) + A[3]*np.sin(t*f[3])*np.exp(-d[3]*t)\n","\n","# plot them!\n","plt.plot(x,y,'k',linewidth=.1)\n","plt.axis('off')\n","plt.show()"]},{"cell_type":"code","execution_count":null,"metadata":{"id":"R2v9uaWebiZC"},"outputs":[],"source":["# what are the elementary components of the plots above? Let's have a look!\n","plt.plot(np.sin(t*f[0])*np.exp(-d[0]*5*t),linewidth=.5)\n","plt.show()"]}],"metadata":{"kernelspec":{"display_name":"Python 3","language":"python","name":"python3"},"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.7.3"},"colab":{"provenance":[]}},"nbformat":4,"nbformat_minor":0}