course_planning:183_projects:f21_week_1_problem_voyager_satellite_solution

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Next revision
Previous revision
course_planning:183_projects:f21_week_1_problem_voyager_satellite_solution [2021/09/07 13:34] – created pwirvingcourse_planning:183_projects:f21_week_1_problem_voyager_satellite_solution [2022/09/06 13:10] (current) valen176
Line 1: Line 1:
 ====== Project 1 Solution: Part B: Voyager collision course ====== ====== Project 1 Solution: Part B: Voyager collision course ======
 +
 +Link to starter code: https://www.glowscript.org/#/user/pcubed/folder/incompleteprograms/program/SatteliteCode
  
 <WRAP tip> <WRAP tip>
Line 24: Line 26:
  Asteroid.pos = Asteroid.pos + vAsteroid*dt  Asteroid.pos = Asteroid.pos + vAsteroid*dt
 </code> </code>
 +
  
 <WRAP tip> <WRAP tip>
Line 62: Line 65:
  
 <code python voyager.sol.py> <code python voyager.sol.py>
-from __future__ import division +GlowScript 2.9 VPython 
-from visual import * + 
-from visual.graph import * +get_library('https://cdn.rawgit.com/PERLMSU/physutil/master/js/physutil.js')
-from physutil import *+
  
 #Objects #Objects
 Voyager = sphere(pos=vector(-5e10,-2e10,-3e9), radius=3e9, color=color.green, make_trail=True) Voyager = sphere(pos=vector(-5e10,-2e10,-3e9), radius=3e9, color=color.green, make_trail=True)
-Asteroid = sphere(pos = Voyager.pos + vector(-3e10,8e10,3e9), radius=5e9)+Asteroid = sphere(pos = Voyager.pos + vector(-3e10,8e10,3e9), radius=5e9, color=color.red, make_trail=True) 
 + 
  
 #Parameters and Initial Conditions #Parameters and Initial Conditions
Line 82: Line 85:
 #MotionMap/Graph #MotionMap/Graph
 trackVoyager = MotionMap(Voyager, tf, 5, markerScale=1e5) trackVoyager = MotionMap(Voyager, tf, 5, markerScale=1e5)
 +trackAsteroid = MotionMap(Asteroid, tf, 5, markerSacale=1e5)
 separation = PhysGraph(numPlots=1) separation = PhysGraph(numPlots=1)
 + 
 #Calculation Loop #Calculation Loop
 while t < tf: while t < tf:
- rate(10)+    rate(10) 
 +     
 +    Voyager.pos = Voyager.pos + vVoyager*dt 
 +    Asteroid.pos = Asteroid.pos + vAsteroid*dt 
 +     
 +    trackVoyager.update(t,vVoyager) 
 +    trackAsteroid.update(t,vAsteroid) 
 +    separation.plot(t,mag(Voyager.pos-Asteroid.pos)) 
 +     
 +    t = t + dt 
 +</code>
  
- Voyager.pos = Voyager.pos + vVoyager*dt +Link to solution: https://www.glowscript.org/#/user/pcubed/folder/solutions/program/SatteliteCodeSolution
- Asteroid.pos = Asteroid.pos + vAsteroid*dt +
- +
- trackVoyager.update(t,vVoyager) +
- separation.plot(t,mag(Voyager.pos-Asteroid.pos)) +
- +
- t = t + dt +
-</code>+
  
 <WRAP tip> <WRAP tip>
  • course_planning/183_projects/f21_week_1_problem_voyager_satellite_solution.1631021651.txt.gz
  • Last modified: 2021/09/07 13:34
  • by pwirving