course_planning:183_projects:s23_week_3_geostationary_orbit

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revisionBoth sides next revision
course_planning:183_projects:s23_week_3_geostationary_orbit [2023/01/25 19:34] – hallsteincourse_planning:183_projects:s23_week_3_geostationary_orbit [2023/01/25 19:51] – hallstein
Line 330: Line 330:
 SatelliteMotionMap = MotionMap(Satellite, tf, 20, markerScale=2000, labelMarkerOrder=False) SatelliteMotionMap = MotionMap(Satellite, tf, 20, markerScale=2000, labelMarkerOrder=False)
 FnetMotionMap = MotionMap(Satellite, tf, 20, markerScale=2000, labelMarkerOrder=False)   # add Fnet force arrows - also need to add in loop) FnetMotionMap = MotionMap(Satellite, tf, 20, markerScale=2000, labelMarkerOrder=False)   # add Fnet force arrows - also need to add in loop)
-#separationGraph = PhysGraph(numPlots=1) 
-graphSatellite = PhysGraph(numPlot=1) 
-graphFnet = PhysGraph(numPlots=1) 
  
-graphpF = PhysGraph(numPlots=2)+graphSeparation = PhysGraph(numPlot=1)   # initialize the magnitude of the separation vector vs time plot 
 +graphFnet = PhysGraph(numPlots=1)   # initialize the magnitude of the net force vs time plot 
 + 
 +graphpF = PhysGraph(numPlots=2)   # initialize the x componets of p and Fnet vs time plot
  
    
Line 352: Line 352:
     FnetMotionMap.update(t, Fnet)  # Add Fnet vector arrows     FnetMotionMap.update(t, Fnet)  # Add Fnet vector arrows
     #separationGraph.plot(t,mag(mSatellite))     #separationGraph.plot(t,mag(mSatellite))
-    graphSatellite.plot(t,mag(Satellite.pos)) +    graphSeparation.plot(t,mag(Satellite.pos))  # update magnitude of sep vector vs time plot 
-    graphFnet.plot(t,mag(Fnet))+    graphFnet.plot(t,mag(Fnet))    # update magnitude of Fnet vs time plot
     scale=1.0E4   # add scale to make fluctuations in Fx visible (units of 10^-4 N)     scale=1.0E4   # add scale to make fluctuations in Fx visible (units of 10^-4 N)
-    graphpF.plot(t,pSatellite.x,scale*Fnet.x)+    graphpF.plot(t,pSatellite.x,scale*Fnet.x)   # update x-components of Fnet and p vs time
    
     t = t + dt     t = t + dt
Line 362: Line 362:
     theta = 7.29e-5*dt     theta = 7.29e-5*dt
     Earth.rotate(angle=theta, axis=vector(0,0,1), origin=Earth.pos)     Earth.rotate(angle=theta, axis=vector(0,0,1), origin=Earth.pos)
-</code>+</code>    
  
 Modification of solution to part B to get |r| vs t, |Fnet| vs t and Fnet,x and P,x vs t: Modification of solution to part B to get |r| vs t, |Fnet| vs t and Fnet,x and P,x vs t:
  
  
  • course_planning/183_projects/s23_week_3_geostationary_orbit.txt
  • Last modified: 2023/10/18 01:20
  • by hallstein