Enzyme Simulation

In our Human Body Systems course, students learn about the digestive system and enzymes that help break down various food macromolecules. As a bonus, most of my HBS students are also taking chemistry, so we have opportunity to share ideas across the disciplines.

In the past, we did a hands-on lab measuring the rate of reaction of catalase to break down hydrogen peroxide. We could vary the temperature, pH, agitation, concentration, or a number of other independent variables to determine their effect on the efficacy of the catalase enzyme. While I love doing hands-on labs with students, this one took several days to collect enough data, and the data we collected was often unreliable–students measured reaction rates using gas pressure sensors. Any leaks in the flasks, stoppers, or hoses would give us bad data, and then using the computer interface to find the slope of the pressure increase (our proxy for reaction rate) was a small hurdle as well.

A couple years ago, my class was running behind and I decided to forego the hands-on portion, but I still needed students to learn the effects of various conditions on enzymes. So, I found a great enzyme simulation at Biology Simulations. Very cool! I love how the collisions between molecules are clearly shown and how we can change so many variables to test their effect on the reaction.

However, for my purposes, this simulation incorporated a bit too much. There are multiple enzymes. There are multiple reactants and products. There are inhibitors! SO many things to play around with. Which is all good, but I needed my students to understand the basics. How do things like pH, temperature, and concentration affect an enzyme…

So, I set out to pare the biology simulations version down to just a few variables and a single enzyme.

Original (Biology Simulations) version

Simplified version

Click the button below to try it out. I also put together a very bare-bones ‘worksheet’ to go along with the simulation. Let me know if this is useful to you.

Basketball Scorebook

For many years, I’ve used the NCAA basketball scorebook for my Jr. High and High School teams, and it’s worked great. However, for the 2023-2024 season, the NFHS changed the rules on tracking team fouls. Teams shoot bonus free throws on the fifth (and subsequent) team foul in each quarter (it used to be after the seventh and tenth foul in each half, but no longer!).

Since the NCAA scorebook only has a place to track team fouls by half, I created my own scorebook that has space to track team fouls by quarter. Check it out, and leave any feedback. I also included an editable Word document below that can be modified to your liking!

Webcam Viewer

I’ve previously posted about my under-$50 “homemade” document camera. I recommended installing a simple program to display the webcam’s output to the screen. However, recent security “improvements” in my school district have prevented me from installing or running such software (not to mention the program I recommended, fscamview, hasn’t been updated in many years and may not be compatible with Windows 10).

As a workaround, I learned that I could create a simple webpage to display the webcam. Webcam images are not sent to my server–your browser simply displays your webcam image so you can project the image to a screen in your classroom or share the browser tab with students on Zoom, Meet, or Teams.

Interactive Measurement Practice

As we begin the school year in the wake of Covid-19, traditional labs and hands-on activities are very limited. Therefore, we need to find ways to do things virtually. I developed this interactive measurement & scale-reading activity using Construct 3. Students can use this to practice reading the scale on a ruler and reading the volume in three different-sized graduated cylinders.

Length Measurement Features

  • Draggable ruler
  • Checks accuracy of measurement and provides feedback.
    • Accepts measurements within ±0.015 cm of the actual length.
    • Provides encouragement to estimate as accurately as possible for measurements within ±0.03 cm of the actual length.
  • Requires input of units (cm) for each measurement.
  • Challenges users to measure five lengths in a row. Any mistakes will cause the count to start over.
  • Button to download a screenshot when finished.

Source Construct 3 file (open with construct.net)

Volume Measurement Features

  • Practice reading scales on 500-ml, 100-ml, and 10-ml graduated cylinders.
  • Checks accuracy of measurement and provides feedback.
    • Accepts a small range of estimates for the final digit in the reading.
    • Requires recording volumes to the level of precision/resolution of the graduated cylinder.
  • Requires input of units (ml) for each measurement.
  • Requires users to correctly measure 6 volumes to complete the activity (2 for each of the cylinder sizes).
  • Button to download a screenshot when finished.

Source Construct 3 file (open with construct.net)

Data Analysis Tool

UPDATE: As of July 2020, I have updated the data analysis tool to use some cool newer JavaScript libraries. My graphing library is Chart.js (it has some fancy animation options and allows copy/pasting of generated graphs) and Handsontable (which uses an excel-like grid for inputting data, which also allows copy/pasting from excel or other spreadsheet programs).

In Modeling Instruction and Advanced Placement science courses, students must be able to analyze data to determine the relationship between two variables. To make this easier, I created a Data Analysis Tool.

I’ve tried a few ways for students to graph data:

MethodProsCons
Graph Paper* Requires a firm grasp of scales, slope, and y-intercept
* Students gain a better grasp of the meaning of these quantities and where they come from
* Flexible — few limitations
* No technology quirks to learn
* Time-consuming
* Requires multiple iterations to linearize non-linear relationships
* Difficult to test several analysis methods to determine the best fit
* Best-fit lines, slopes, and y-intercepts are not as accurate
TI-84 Calculator* Can be used on the AP Exam & ACT!
* Manual window-setting requires students to consider how they will view their data
* Students can create their own best-fit line by writing a linear equation in Y1
* OR, it has LinReg capabilities
* A bit quicker than graphing by hand
* Data entry, turning stat-plots on, and graphing is not as intuitive as it could/should be
* Regressions are also not very intuitive, and getting the linear fit to display on the graph is another hurdle
* Doesn’t handle units at all — that’s up to the student
* No easy way to get a printout of the graph for inclusion in reports
* Not all students have this calculator (it’s expensive)
Vernier LoggerPro* Quick
* Takes care of units, labels, and titles
* Easily adjust what is graphed on each axis
* Support for force, motion, and other sensors
* Linear fit is available with 1 button click
* Supports graphing multiple data sets, a secondary y-axis scale, and many advanced features
* Linearizing with calculated columns is somewhat cumbersome; requires quite a few steps
* For novices, the large number of features make it difficult to remember the correct steps to get the desired results
* Not likely to use this software in other contexts
* Not widely available (though the site license is generous)
Microsoft Excel or Google Sheets* Powerful and flexible
* Supports graphing multiple data sets, secondary y-axis scale, and many advanced features
* Likely to be used in college, industry, etc.
* Widely available (and for Google Sheets, free)
* Need to set up own data table
* New versions of Excel do not label axes by default
* Google Sheets graphs lack some features and are not the most intuitive (though I haven’t used them in a couple years so perhaps they’ve improved)
* Again, the huge number of features often make it hard to find the desired functions

To start the year, I always have students graph data using graph paper for a couple weeks.  I think students need to be able to do it themselves and understand the basic considerations of choosing scales, deciding what to plot where, and finding slopes and y-intercepts manually before having a device/computer do it for them.

In the past, after students got more comfortable graphing things by hand, I showed them how to use some tools (in the past, this tool has been LoggerPro).  But, after having to guide students through the LoggerPro linearization process time and time again for each lab, I wanted to find a better solution. My current version works pretty well, but I’d welcome your feedback. I also have a version with linearization tools hidden.

Version 1.0

In June of 2016, I wrote a quick online data analysis tool (using the graphing capabilities of CanvasJS).  Check it out here.  It does most of what I need it to, which is to take a set of data, graph it, allow students to linearize it (graph y vs. x2, y vs. 1/x, etc.), and output the best-fit line equation.

I also created a version with the linearization capabilities hidden (to require students to do this manually).

Graphing an inverse squared relationship.
Graphing Y vs. 1/X² gives a linear relationship. The equation of the best-fit line is also available.

Improvements I may try to implement in the future:

  • Data input improvements (arrow keys for navigating) – implemented though still has some quirks
  • Improve linearization interface somehow – still thinking about this
  • Display the squared, cubed, or inverse of a column in the data table (?)
  • Graph multiple data sets simultaneously

Density Lab Simulation

Unit 1 of the Modeling Instruction chemistry curriculum has students develop the ideas of mass and volume and then the relationship between them (density). Consistent with the Modeling Instruction method, students collect data and analyze that data to develop a model to describe a relationship.

This year, my school is beginning virtually. It will be a challenge to transform our chemistry lab activities to the online format until we are able to resume in-person classes (and even then, we will be limited in the types of activities that we can do while “socially distant”). As I thought about facilitating my class online, I discovered Construct.net, which is typically used to produce online games. I found that it can be adapted to simulate labs like this one. I created a couple different versions. Check them out!

Dimensional Measurement Version
Water Displacement version

It simply allows students to measure the mass and volume of several samples of a material (currently, it has steel, aluminum, and wood). The first sample is always the same size for everybody, but subsequent samples are random sizes/masses (so encourage students to include a wide range of sample sizes in their data). Students can then use whatever graphing tools desired to analyze the data. A few years ago, I made a simple data analysis tool that would work well for this.

I like that the water displacement version shows which substances sink and float. Also, you can produce some samples that, because of their size or density (wood) do not completely submerge. This provides a chance to discuss what the water displacement measurement represents.

The water displacement version can also be used to show the relationship between cm3 and ml. Just measure a sample with the ruler and then dunk it in the water.

Leave a comment about any issues you find or requests for features. Here are some things I’m thinking about:

  • More realistic interactions: an on-screen ruler to measure length, width, and height would build measurement skills. But, Contruct.net is really a 2-dimensional tool, so measuring that third dimension would be a challenge.
  • I had a request for cylindrical samples to better match the materials many of us use for the actual lab activity. It’s an idea for the future, especially if this is used to supplement the in-person activity for those who are absent.
  • UPDATE July 27, 2020:
    • Added zinc, copper, and lead to the substances.
    • Changed the button text to say “New Size”
    • Clarified the instructions.

The source files are available for anyone wanting to modify or extend. This was my first time using Construct.net, so keep that in mind before you berate me for my shoddy work 🙂. I do encourage anyone else working on this or similar lab simulations to share what you come up with.

Construct.net Source File – Version 1 (LxWxH Volume measurement)

Construct.net Source File – Version 2 (Water Displacement Volume measurement)

Reaction Time Applet (HBS 2.2.4)

During our PLTW Human Body Systems lesson on the nervous system, we study the speed at which signals travel through neurons, as well as the regions of the brain that control various functions. The curriculum points us to a cool little activity that measures our reaction time when doing tasks of increasing complexity. In Case 4, it becomes quite challenging to have to read a set of instructions and do the opposite as quickly as possible.

The applet used by PLTW is a Shockwave applet by Serendip, but Shockwave hasn’t been well-supported in a long time, and now that my students use Chromebooks, it definitely won’t work. It now appears that the shockwave app has been removed, and the PLTW curriculum now uses a much more basic reaction time app.

So, I set out to recreate this great applet in a more portable and cross-platform language (HTML and JavaScript). The functionality is similar. My attempt is definitely not as refined as the original, but it gets the job done. Also, students are responsible for recording their own times, because my version discards all times when changing between parts of the activity.

Reaction Time Applet

Homemade Wave Driver

In AP Physics 1, students need to learn about the relationship between wave frequency and wavelength (wave speed).  If students are going to investigate this relationship and collect data for analysis themselves, there are a couple options:

  1. Snaky springs – ($13.50/ea from Amazon and comes with slinky and rope, $17/ea from Arbor Sci, $18/ea from PASCO) – Students can set up standing waves in the spring, then measure frequency and wavelength.  It’s very tactile and when it works, it works well — data is reliable.  It is difficult to get a feel for how to oscillate at the correct frequency for the standing waves (and some groups weren’t really successful doing so in the time we had).  Only allows for 4 or 5 harmonics before it becomes very difficult to introduce such high-frequency oscillations.
  2. Computer Simulation – PhET – Allows adjustment of amplitude, frequency, end type, damping, etc. and measurement of wavelength.  A fine option, though not as tactile.  I use it for students to investigate wave reflections, but have not used it for wavelenth & frequency.  I would prefer a real-world, more tactile method if possible — I question whether students always make the connection to real-world situations when their primary experience is a computer simulation.
  3. Wave drivers – (Vernier: $115, but also need a power amplifier, $199; PASCO: $249, but also need a function generator, $259+; CPO $400, there are probably other systems available, but they all seem to be very expensive) – allows students to adjust the frequency of vibrations in a string until a standing wave is set up, then measure wavelength.  Using this method, students are very interested and engaged, and they get to see a real-world physics phenomenon.  The data produced is very reliable.  But this would be a very expensive set-up to provide for 6 lab groups.

I set out to see if I could make my own wave driver for less $$$.

Speaker Version

For my first iteration, I used some Logitech computer speakers ($20 on Amazon).  Essentially, I removed the grille covering the speaker and hot-glued a bent paper clip to the speaker cone, that allows the speaker to drive a string.

I tie the string to a solid object on one end (ring stand) and run it across a table and off the end over a pulley.  Then, hang a weight off the pulley end of the string to introduce a known tension.

Diagram of the speaker version of the apparatus

I connect the speaker system to a laptop, have students go to an online function generator (like this one, though there are many available), and adjust the frequency until a standing wave is set up.  Then they measure the wavelength and record the frequency/wavelength pair.

By graphing wavelength vs. frequency, students can determine that an inverse relationship exists, which they can then linearize (possibly by using my data analysis tool).  The slope then comes out to be a speed, which we determine to be the wave speed in the string.

The speaker works wonderfully.  The wave amplitude isn’t great, but the standing waves are clearly visible.

I also tried this with a cheaper speaker system, but the speaker cone barely moved, so it didn’t produce visible standing waves in the string.  You need a speaker for which the cone moves with a decent amplitude to drive the waves in the string.

Homemade Wave Driver

I wanted to see if I could further cut down on the cost.  I obtained the following supplies:

DescriptionSourceCost/each
TDA7297 Audio Power Amplifier *(There are many similar versions available.  Try to find one with a 5.5 mm OD barrel/coaxial power plug and a 3.5mm aux audio connector.  The one I bought had a 3-pin audio connector, so I had to buy some of those and solder it to a 3.5mm aux cord.  You only need a single channel or just use 1 of a 2-channel amplifier)eBay
(search)
$2.65
(more or less)
DC 12V AC Adapter *(You may be able to find them more cheaply – just need a 5.5 mm OD / 2.1 mm ID plug. I went with a 2-amp output to ensure enough power, but probably 1A is enough.)eBay
or
Amazon
$5.00
or
$4.00
5′ Aux cord (3.5mm headphone plug) *(one cord will make 2 wave drivers – you can use both ends)eBay$0.50
3 pin board connector plug *(This was needed for mine to connect the signal in to the board, but depends on the board you bought)eBay$0.35
2″ or 3″ speaker
A Logitech 3W 4-ohm speaker worked ok.  A 50W car speaker worked ok also, but I’ll probably need to use a heavier string.   Around a 10W speaker seems to be ideal.
eBay
eBay 2
Amazon
Amazon 2
$5.00 – $7.00
Speaker wire – Used scraps I had available
Screws – I have an assortment, but here’s what I used
(3) tiny screws to attach PCB amplifier board to wood board
(2) 2″ drywall screws (or similar) to attach speaker to board
(?) I also screwed down the audio cable with a piece of plastic to keep it from being accidentally pulled off and damaging stuff.
1×4 mounting board (about 6″ long) – From scraps
Total$12.50 – $15.50

The initial prototype works great.

Prototype. Just needs power to the barrel connector on the PCB and a bent paperclip glued to the speaker cone to interface with the string to produce standing waves.

Assembly:

  1. Cut a piece of 1×4 board in about 6 inch long pieces to serve as a base.
  2. Use 3 tiny screws to secure the amplifier PCB to the 1×4. Arrange it so the power connector hangs off one side of the board.
  3. Cut a 3.5mm aux cord in half.  Strip the insulation from the end you cut.  There are 2 wires inside that carry the left and right channel signal (carefully strip about 1/4″ of insulation from these as well), and bare wires which are neutral.  Of course, not all cables are made the same way, so yours may be different.
  4. Cut a fan extension cable in half so you have the female pin end (shown).  Strip the ends you cut.


  5. Twist together the black fan connector wire with the neutral aux cord wire.  Similarly, connect the red fan connector wire with one channel from the aux cord.  And the yellow fan connector wire with the other channel from the aux cord.  Solder these connections, then cover each one with heat-shrink tubing (or easier: electrical tape).  Sorry, I forgot to take a picture of this…
  6. Plug the female pin connector into the amplifier board.  I also used an electrical cord staple, removed the nails, and screwed it down on top of the aux cord wires to hopefully prevent some jerk from jerking the cord and destroying my wave driver.

  7. Run a wire from the amplifier output (+) to one terminal on the speaker.
  8. Run another wire from the amplifier output (-) to the other terminal of the speaker.  (Doesn’t matter which is which)
  9. Tighten the screw-connectors on the amplifier output terminals.

  10. Solder these other ends on the wires to the speaker.

  11. Use the 2″ screws to attach the speaker to the 1×4.  Depending on the size of your speakers, you may need longer/shorter screws.
  12. Bend a paper clip and hot glue or epoxy it to the speaker cone as shown:

    The cord or string being driven sits in the loop on the paper clip.
  13. Plug it into a computer speaker output, plug in AC adapter, and go to http://www.szynalski.com/tone-generator/ (or any other similar tone/frequency generator).
  14. Run a horizontal string as described in the “Speaker Version” section above.
  15. Adjust the frequency until standing waves are produced….

Any questions?  Comment below.

UPDATE: One of my modeling heroes, Mark Schober, has built something similar (though much nicer looking) – read about it and sign up for the workshop if you’re in the NYC area.

Vernier LoggerPro Function Reference

LoggerPro functions for use in Calculated Columns.  I always Google for this information and can never find it online.  It is hidden in the “Help” for LoggerPro (who would think to look there!?).

Functions

– Trigonometric functions will use degrees or radians as set in the Settings for (file name) in the File menu.

– For more information about Savitsky-Golay methods see Numerical Recipes in C: http://lib-www.lanl.gov/numerical/bookcpdf.html chapter 15

analysis

FunctionDescription
analysis(“X”, startRow, endRow)
Takes all columns named “X” and extracts startRow to endRow for each of those columns, appending the values into a single column.
dataSetsdatasets(“X”)
Appends the dataset names of all datasets that have a column named “X”. Use analysis and datasets together to create a graph (analysis on the vertical axis and datasets on the horizontal). for example, if you had 3 datasets as follows:
DS1  DS2  DS3
X       X      Y
1      11     21
2      12     22and then added analysis(“X”, 1, 1) and datasets(“X”) you would get:datasets    analysis
DS1             1
DS2            11

beats per minuteBeatsPerMinute(“Signal”, “Time”, intervalInSeconds, minPercent, maxPercent, noise)
Returns the number of beats per minute of the values in “Signal” vs. “Time”. This function is similar to the rate function except that the interval given here is always in seconds and the returned value is always in minutes. For example, if “Time” is in seconds then: beatsPerMinute(“Signal”, “Time”, interval, min, max, noise) = 60 * rate(“Signal”, “Time”, interval, min, max, noise)

blood pressure

diastolicThe measured arterial pressure when the heart is at rest. “Pressure” and a “Time” column as inputs and return a single number.
diastolic(“Pressure”, “Time”)
“Pressure”: Pressure values from the BPS
“Time”: Time the pressure values were recorded
Returns the smaller number of blood pressure
meanArterialPressuremeanArterialPressure(“Pressure”, “Time”)“Pressure”: Pressure values from the BPS
“Time”: Time the pressure values were recorded
Returns the pressure value at the max peak used for blood pressure calculations.
OscillationsOscillations(“Pressure”, “Time”)
“Pressure”: Pressure values from the BPS
“Time”: Time the pressure values were recorded
Returns the Oscillations of the peaks and valleys used to calculate systolic and other blood pressure values.
OscillatoryPeaksOscillatoryPeaks(“Pressure”, “Time”)
“Pressure”: Pressure values from the BPS
“Time”: Time the pressure values were recorded
Returns the peaks used to calculate systolic, diastolic, and pulse (the “high” values in “Oscillations”).
pulsepulse(“Pressure”, “Time”)
“Pressure”: Pressure values from the BPS
“Time”: Time the pressure values were recorded
Returns the pulse using the inputs from the Blood Pressure Sensor (similar results, different algorithm as the other beats-per-minute functions)
systolicsystolic(“Pressure”, “Time”)
“Pressure”: Pressure values from the BPS
“Time”: Time the pressure values were recorded
The measured arterial pressure when the heart contracts. Returns the larger number of blood pressure

boolean

For the boolean functions a 1 is considered true, 0 false and anything else an invalid input
ANDAND(X, Y) return 1 if and only if X and Y are both 1
NOTNOT(X) return 1 if X is 0; 0 if X is 1
OROR(X, Y) return 1 if X or Y is 1
XORXOR(X, Y) return 1 if X or Y is 1 but not both

calculus

calculus >
derivativederivative(“Y”, “X”)
“Y”: A column of real numbers
“X”: Optional. A column of real numbers
The numerical derivative is the weighted average of the slope of ‘n’ points around each point. You can set ‘n’ in Settings for (Name). If you don’t supply an “X” column, the program will find one.
derivativeSGderivativeSG(“Y”, “X”)
“Y”: A column of real numbers
“X”: Optional. A column of real numbers
Savitsky-Golay derivative. Fits a polynomial to ‘n’ points around each point and computes the derivative of the polynomial at that point. You can set ‘n’ in Settings for (Name). If you don’t supply an “X” column, the program will find one.
derivativeTimeShiftderivativeTimeShift(“Y”, “X”)
Returns the derivative of “Y” with respect to “X”.
This function is specifically designed to be used with photogate and picket fence data. The derivatives returned are adjusted to estimate values at the start of the timing interval, instead of the midpoint. For details see The Physics Teacher, Vol 35, April 1997, p. 220. The article written by William Leonard is entitled “The Dangers of Automated Data Analysis.”
Average velocity during the time interval is equal to the instantaneous velocity at midpoint of the time interval.
Where
integralintegral(“Y”,”X”)
“Y”: A column of real numbers
“X”: Optional. A column of real numbers
The numerical integral is the running sum of the areas of rectangles calculated by the midpoint rule. The i’th rectangle is (Yi – Y(i-1)) / (Xi – X(i-1)). If you don’t supply an “X” column, the program will find one.
secondderivativesecondDerivative(“Y”, “X”)
“Y”: A column of real numbers
“X”: Optional. A column of real numbers
Calculates the numerical second derivative of “Y” with respect to “X”. If you don’t supply an “X” column, the program will find one.
secondderivativeSGsecondDerivativeSG(“Y”, “X”)
“Y”: A column of real numbers
“X”: Optional. A column of real numbers
Savitsky-Golay second derivative. Fits a polynomial to ‘n’ points around each point and computes the second derivative of the polynomial at that point. You can set ‘n’ in Settings for (Name). If you don’t supply an “X” column, the program will find one.
secondderivative
Time Shift
secondDerivativeTimeShift(“Y”, “X”)
“Y”: A column of real numbers
“X”: Optional. A column of real numbers
Numerical time-shifted second derivative. Calculates the second numerical derivative of “Y” with respect to “X”. The values are shifted so that the derivatives are calculated at the midpoints between each two values. If you don’t supply an “X” column, the program will find one.
collapsecollapse(“X”)
Returns a column with all non-numerical cells (blanks and text) removed.
collapseIndirectcollapseIndirect(X, Y)
Returns a column of only the rows in “X” corresponding to rows in “Y” that have valid numerical cells.
constantConstant(x, num)
x: A real number
num: A real number or a column
Generates a constant column filled with the value ‘x’.  The number of values in the returned column is num, or if a column was passed in, the size of the passed-in column.
deltadelta (“X”)
“X”: A column of real numbers
Returns a column of values where the i’th value is the i’th value in “X” minus the (i-1)’th value in “X”.

digital filtering

lowPassFilter(“Y”, “X”, “ripple”, “freqCutoff”)
“Y”: The data column to be filtered
“X”: The associated time column for “Y”
“ripple”: The ripple allowed in the pass-band
“freqCutoff”: Cut-off frequency (-3dB), in hertz
Applies a Chebyshev low-pass filter. For “ripple”, enter a value that is a percent of the pass-band. To apply a Butterworth low-pass filter, set “ripple” to 0.
highPassFilter(“Y”, “X”, “ripple”, “freqCutoff”)
“Y”: The data column to be filtered
“X”: The associated time column for “Y”
“ripple”: The ripple allowed in the pass-band
“freqCutoff”: Cut-off frequency (-3dB), in hertz
Applies a Chebyshev high-pass filter. For “ripple”, enter a value that is a percent of the pass-band. To apply a Butterworth high-pass filter, set ripple to 0.
bandPassFliter(“Y”, “X”, “lowFreq”, “highFreq”)
“Y”: The data column to be filtered
“X”: The associated time column for “Y”
“lowFreq”: Low frequency cut-off (-3dB), in hertz
“highFreq”: High frequency cut-off (-3dB), in hertz
Ripple is automatically set to zero and is not adjustable. The function returns the signal with the frequencies outside the designated frequency range removed.
bandStopFliter(“Y”, “X”, “lowFreq”, “highFreq”)
“Y”: The data column to be filtered
“X”: The associated time column for “Y”
“lowFreq”: Low frequency cut-off (-3dB), in hertz
“highFreq”: High frequency cut-off (-3dB), in hertz
Ripple is automatically set to zero and is not adjustable. The function returns the signal with the frequencies inside the designated frequency range removed.
timeDecayFilter(“Y”, “X”, “decayConstant”)
“Y”: The data column to be filtered
“X”: The associated time column for “Y”
“decayConstant”: A value in seconds that determines the decay of “Y”
Applies an exponential time decay to the signal.
ElectrophoresisInterpolateElectrophoresisInterpolate(“Std. Dist”, “Std. BP”, “Dist”)
“Std. Dist”: Distances from the standard
“Std. BP”: Base Pair Counts from the standard
“Dist”: Distances to interpolate
Returns a column of base pair counts based on the Electrophoresis curve fit for “Std. Dist” vs. “Std. BP” given “Dist”. Will NOT work if curve fit has been deleted. This function is used automatically when doing a Gel Analysis (Electrophoresis).
expexp(“X”)
“X”: A column of real numbers
Returns the exponent, exp(x) = e^x, where e is the natural log base (2.17…).
integerinteger(“X”)
Extracts the integral part of values in “X”.
interpolateinterpolate(“X”)
Fills in missing values using linear interpolation.
lnln(“X”)
“X”: A column of real numbers larger than 0
Returns the natural logarithm. If b = ln(a) then e^b = a  (Where e is the constant 2.17…).
loglog(“X”)
“X”: A column of real numbers larger than 0
(Log base 10) If b = log(a) then 10^b = a.
modulomodulo(“X”, n)
“X”: A column of integers
n: An integer larger than 0
Returns the remainder of each of the numbers in “X” when divided by n.

photogate

Blocked MidTimesBlockedMidTimes(“Time”, “Gate1”, “Gate2”)
“Time”: Optional. A column of real numbers (the times of events)
“Gate1”: A column of photogate states (1’s and 0’s)
“Gate2”: Optional. A column of photogate states (1’s and 0’s)
Calculate the average times between blocked events from Gate 1 to Gate 2. If you don’t enter a “Time” column, the program will find one. If you don’t enter “Gate2”, “Gate1” will be used.
Blocked to BlockedBlockedToBlocked(“Time”, “Gate1”, “Gate2”)
“Time”: Optional. A column of real numbers (the times of events)
“Gate1”: A column of photogate states (1’s and 0’s)
“Gate2”: Optional. A column of photogate states (1’s and 0’s)
Returns a column of the times between successive blocked events in gate 1 and blocked events in gate 2. If you don’t enter a “Time” column, the program will find one. If you don’t enter “Gate2”, “Gate1” will be used.
Blocked to UnblockedBlockedToUnblocked(“Time”, “Gate1”, “Gate2”)
“Time”: Optional. A column of real numbers (the times of events)
“Gate1”: A column of photogate states (1’s and 0’s)
“Gate2”: Optional. A column of photogate states (1’s and 0’s)
Returns a column of the times between successive blocked events in gate 1 and unblocked events in gate 2. If you don’t enter a “Time” column, the program will find one. If you don’t enter “Gate2”, “Gate1” will be used.
Blocked to Unblocked MidTimesBlocked to Unblocked MidTimes
“Time”: Optional. A column of real numbers (the times of events)
“Gate1”: A column of photogate states (1’s and 0’s)
“Gate2”: Optional. A column of photogate states (1’s and 0’s)
Calculate the average time between the blocked events in gate 1 and unblocked events in gate 2. If you don’t enter a “Time” column, the program will find one. If you don’t enter “Gate2”, “Gate1” will be used.
derivativeTimeShiftDerivativeTimeShift (“Y”, “X”)
Returns the derivative of “Y” with respect to “X”.
This function is specifically designed to be used with photogate and picket fence data. The derivatives returned are adjusted to estimate values at the start of the timing interval, instead of the midpoint. For details see The Physics Teacher, Vol 35, April 1997, p. 220.
Average velocity during the time interval is equal to the instantaneous velocity at midpoint of the time interval.
Where
Pendulum PeriodPendulumPeriod(“Time”, “Gate1”)
“Time”: Optional. A column of real numbers (the times of events)
“Gate1”: A column of photogate states (1’s and 0’s)
Calculate the time between every other blocked event on Gate 1. If you don’t enter a “Time” column, the program will find one.
secondDerivativeTimeShiftsecondDerivativeTimeShift(“Y”, “X”)
“Y”: A column of real numbers
“X”: Optional. A column of real numbers
Numerical time-shifted second derivative. Calculates the second numerical derivative of “Y” with respect to “X”. The values are shifted so that the derivatives are calculated at the midpoints between each two values. If you don’t supply an “X” column, the program will find one.
Unblocked to BlockedUnblockedToBlocked(“Time”, “Gate1”, “Gate2”)
“Time”: Optional. A column of real numbers (the times of events)
“Gate1”: A column of photogate states (1’s and 0’s)
“Gate2”: Optional. A column of photogate states (1’s and 0’s)
Returns a column of the times between successive unblocked events in gate 1 and blocked events in gate 2. If you don’t enter a “Time” column, the program will find one. If you don’t enter “Gate2”, “Gate1” will be used.
Unblocked to UnblockedUnblockedToUnblocked(“Time”, “Gate1”, “Gate2”)
“Time”: Optional. A column of real numbers (the times of events)
“Gate1”: A column of photogate states (1’s and 0’s)
“Gate2”: Optional. A column of photogate states (1’s and 0’s)
Returns a column of the times between successive unblocked events in gate 1 and unblocked events in gate 2. If you don’t enter a “Time” column, the program will find one. If you don’t enter “Gate2”, “Gate1” will be used.
Unblocked to Blocked MidTimesUnblocked to Blocked MidTimes
“Time”: Optional. A column of real numbers (the times of events)
“Gate1”: A column of photogate states (1’s and 0’s)
“Gate2”: Optional. A column of photogate states (1’s and 0’s)
Calculate the average time between unblocked events in gate 1 and blocked events in gate 2. If you don’t enter a “Time” column, the program will find one. If you don’t enter “Gate2”, “Gate1” will be used.
Unblocked MidTimesUnblockedMidTimes(“Time”, “Gate1”, “Gate2”)
“Time”: Optional. A column of real numbers (the times of events)
“Gate1”: A column of photogate states (1’s and 0’s)
“Gate2”: Optional. A column of photogate states (1’s and 0’s)
Calculate the average times between unblocked events from Gate 1 to Gate 2. If you don’t enter a “Time” column, the program will find one. If you don’t enter “Gate2”, “Gate1” will be used.

raterate(“Y”, “X”, t, m1, m2, n)
“Y”: A column of real numbers
“X”: Optional. A column of real numbers
t: Optional. Time interval
m1: Optional. Minimum threshold
m2: Optional. Maximum threshold
n: Optional. Noise threshold
Returns the rate of “Y” with respect to “X”, where t is the time interval measured, m1 is min percentage threshold, m2 is max percentage threshold, and n is noise threshold. “X”, t, m1, m2, and noise are all optional with default values “X” is time column, t = 1/10 the range, m1 = 40%, m2 = 60%, and noise = 0. Details

rotary motion


amplitude(“Data Column”, “Time Column”, “Min Percent”, “Max Percent”, “Time Interval”)
“Data Column”: Data for which you want to calculate amplitude
“Time Column”: Associated time column for “Data Column”
“Min Percent”: Threshold used to detect valleys
“Max Percent”: Threshold used to detect peaks
“Time Interval”: Period of time over which amplitude is calculated (in the time units of the experiment)
Calculates peak to peak amplitude. For Min and Max Percent, enter values between 0 and 100. Smaller values are more sensitive to noise and thus more sensitive to real cycles. Larger values are less sensitive to noise; too large of a value may filter out real cycles. “Time Interval” ends at the row at which the value is calculated (the current time).
period(“Data Column”, “Time Column”, “Min Percent”, “Max Percent”, “Time Interval”)
“Data Column”: Data for which you want to calculate period
“Time Column”: Associated time column for “Data Column”
“Min Percent”: Threshold used to detect valleys
“Max Percent”: Threshold used to detect peaks
“Time Interval”: Period of time over which period is calculated (in the time units of the experiment)
Calculates the period of an oscillating function. For Min and Max Percent, enter values between 0 and 100. Smaller values are more sensitive to noise and thus more sensitive to real cycles. Larger values are less sensitive to noise; too large of a value may filter out real cycles. “Time Interval” ends at the row at which the value is calculated (the current time).
roundround(“X”)
“X”: A column of real numbers
Round. Returns the closest integer to x. If x is equidistant to two integers, round(x) gives the largest of the two (e.g., round(0.5) = 1).
smoothAvesmoothAve(“X”)
“X”: A column of real numbers
Returns a column of moving averages of the values in “X”. The width of the “window” to use when averaging points can be set in Settings for (Name)…

statistics

absabs(“X”)
“X”: A column of real numbers
Absolute value. If x less than 0, then abs(x) = -x.  Otherwise, abs(x) = x.
ceilingceiling(“X”)
“X”: A column of real numbers
Returns the smallest integer larger than or equal to x.
floorfloor(“X”)
“X”: A column of real numbers
Returns the largest integer smaller than or equal to x.
maxmax(“X”)
“X”: A column of real numbers
Compares all the values in a single column and returns a single number-the largest number in the column.
max2max2(“X”, “Y”)
“X”: column of real numbers
“Y”: A column of real numbers or a single number.
Compares all the values in a column against a real number (e.g max2(“X”, 5.1))
meanmean(“X”)
“X”: A column of real numbers
Arithmetic mean. Returns the sum of all the values in “X” divided by the number of values.
medianMedian(“X”).
“X”: A column of real numbers
If m = median(“X”), then half the numbers in “X” are greater than (or equal) to m, and half are less than or equal.
minmin(“X”)
“X”: A column of real numbers
Compares all the values in a single column ,and returns a single number-the smallest number in the column.
min2min(“X”, “Y”)
“X”: A column of real numbers
“Y”: A column of real numbers or a single number
Compares all the values in a column against a real number (e.g min2(“X”, 5.1))
numRowsNumRows(“X”)
“X”: A column of real numbers
Returns a single value-the number of rows in the column.
randIntrandInt(min, max, num):
min: A real number
max: A real number
num: A real number or a column
Random Integer. Returns a column of random integers between min and max (inclusive). The size of the returned column is num. If num is a column, then the size will be the number of rows in that column.
randRealrandReal(min, max, num)
min: A real number
max: A real number
num: A real number or a column
Random Real. Returns a column of random real numbers between min and max (inclusive). The size of the returned column is num. If num is a column, then the size will be the number of rows in that column.
stddevstddev(“X”)
“X”: A column of real numbers
Standard Deviation. Returns a column representing the standard deviations of each of the numbers in a column.
stepstep(start, increment, num, first, skip)
start: Start value
increment: Increment value
num: Number of values to generate
first: Optional. First non-empty row
skip: Optional. Rows to skip between each value
Generates a column “num” rows long starting with “start” and incrementing by “increment”. “num” can be a positive integer or a column name. Optional parameters: “first” is the first non-empty row and “skip” is the number of rows to skip between each value.
StepColumnBasestepColumnBased(“X”, start, increment, first, skip)
start: Start value
increment: increment value
first: Optional. First non-empty row
skip: Optional. Row to skip between each value
Generates a column based on non-empty values in column “X” starting with “start” and incrementing by “increment.” “First” is the first non-empty row and “skip” is the number of rows to skip between each value.
subsetsubset(“X”, startRow, step)
“X”: A column of real numbers
startRow: An integer larger than 0
step: An integer larger than 0
Extract a subset. Returns a column extracted from “X” starting with ‘startRow’ by ‘step’. For example, subset(“X”, 1, 2) will get every second row of “X” starting with row 1.
sumSum(“X”)
“X”: A column of real numbers
Returns a column whose n’th value is the sum of the values in “X” from row 1 to n.
sqrtSquare root. “X”: A column of non-negative real numbers.
If x is the square root of y, then x*x = y.

trigonometric


sinsin(“X”)
“X”: A column of real numbers
In a right triangle with angle between two sides ‘x’, sin(x) is the length of the opposite side divided by the hypotenuse.
coscos(“X”)
“X”: A column of real numbers
In a right triangle with angle between two sides ‘x’, cos(x) is the length of the adjacent side divided by the hypotenuse.
tantan(“X”)
“X”: A column of real numbers
In a right triangle with angle between two sides ‘x’, tan(x) is the length of the opposite side divided by the adjacent side.
asinasin(“X”)
“X”: A column of real numbers between -1 and 1
Arcsine function. asin(x) = the angle whose sine is x.
acosacos(“X”)
“X”: A column of real numbers between -1 and 1
Arccosine function. acos(x) = the angle whose cosine is x.
atanatan(“X”)
“X”: A column of real numbers
Arctangent function. atan(x) = the angle whose tangent is x. The result will be between -pi/2 and pi/2.
sinhsinh(“X”)
“X”: A column of real numbers
Hyperbolic sine.
coshcosh(“X”)
“X”: A column of real numbers
Hyperbolic cosine.
tanhtanh(“X”)
“X”: A column of real numbers
Hyperbolic tangent.

ValueValue(n, “X”)
n: Number of rows backwards (when n < 0) or forwards (n >0) in column “X” to extract a value from.
“X”: Column from which to extract values .
Create a new column from another column by extracting offset values.

If data are imported from an experiment file, you may want to specify the independent column. For example, if the imported data included “time” in the first column but you wanted to calculate the derivative of pH with respect to volume, you have to define the derivative as derivative(“pH”,”Volume”).

Coulomb’s Law Simulation

In the Modeling Instruction materials for Physics: E&M, there are videos provided that can be analyzed to develop quantitative relationships between electrostatic force, separation distance, and charge quantity.  However, this process is rather complex and requires using distance as a proxy for force.  While the inverse-square relationship comes out (though even that can be difficult to determine if students aren’t very careful about the video analysis), it does not allow for a determination of the Coulomb constant, since the charge of the objects in the video is not known.

Using actual charged pith balls and measuring repulsion.

To try to remove some of the difficulties while still using the data analysis to develop model for electrostatic force interactions, I developed a simulation that allows a more direct “measurement” of charge, distance, and force.

Simulation 1: Using a Spring as a Measure of Force

Simulation 2: Direct Measurement of Force

Screen shot of Simulation 2 – direct measurement of force

This simulation is great to use in combination with my data analysis tool to analyze the relationship between distance, charge, and forces.

Close