Difference between revisions of "EGR 103/Fall 2019/Minilab 3"

From PrattWiki
Jump to navigation Jump to search
Line 12: Line 12:
 
:* [[Python:Interpolation]]
 
:* [[Python:Interpolation]]
 
:* [[Python:Extrema]]
 
:* [[Python:Extrema]]
:* [https://docs.scipy.org/doc/scipy-0.18.1/reference/generated/scipy.interpolate.CubicSpline.html scipy.interpolate.CubicSpline¶] at [http://docs.scipy.org docsscipy.org]
+
:* [https://docs.scipy.org/doc/scipy-0.18.1/reference/generated/scipy.interpolate.CubicSpline.html scipy.interpolate.CubicSpline¶] at [http://docs.scipy.org docs.scipy.org]
  
 
* Problem 2.5.2
 
* Problem 2.5.2
** [https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.simps.html scipy.integrate.simps] at [http://docs.scipy.org docsscipy.org]
+
** [https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.simps.html scipy.integrate.simps] at [http://docs.scipy.org docs.scipy.org]
 
** Note that the denominator in the line of action calculation is $$f_t$$, which you already calculated.
 
** Note that the denominator in the line of action calculation is $$f_t$$, which you already calculated.
 
** If you want to see a graph of the cross section (i.e. Figure 19.9(b)), assuming you call the height above the bottom $$z$$ and the width at that height $$wz$$, you can add the following code:
 
** If you want to see a graph of the cross section (i.e. Figure 19.9(b)), assuming you call the height above the bottom $$z$$ and the width at that height $$wz$$, you can add the following code:
Line 28: Line 28:
  
 
* Problem 2.5.3
 
* Problem 2.5.3
** [https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.trapz.html scipy.integrate.trapz] at [http://docs.scipy.org docsscipy.org]
+
** [https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.trapz.html scipy.integrate.trapz] at [http://docs.scipy.org docs.scipy.org]
** [https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.simps.html scipy.integrate.simps] at [http://docs.scipy.org docsscipy.org]
+
** [https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.simps.html scipy.integrate.simps] at [http://docs.scipy.org docs.scipy.org]
 
** Since $$n$$ goes from 2 to 30, it is not '''exactly''' what you need for an index to an array...think carefully about how to store things.
 
** Since $$n$$ goes from 2 to 30, it is not '''exactly''' what you need for an index to an array...think carefully about how to store things.
  

Revision as of 02:29, 3 December 2019

  • Problem 2.5.1
    • To load data froman Excel file with headers:
import pandas
 
# %% Load data
edata = pandas.read_excel("file.xlsx")
col_1_stuff = edata.values[:, 0].copy()
col_2_stuff = edata.values[:, 1].copy()
  • Problem 2.5.2
    • scipy.integrate.simps at docs.scipy.org
    • Note that the denominator in the line of action calculation is $$f_t$$, which you already calculated.
    • If you want to see a graph of the cross section (i.e. Figure 19.9(b)), assuming you call the height above the bottom $$z$$ and the width at that height $$wz$$, you can add the following code:
fig = plt.figure(num=1, clear=True)
ax = fig.add_subplot(1, 1, 1)
zval = np.block([z[::-1], z[:]])
wval = np.block([-wz[::-1] / 2, wz[:] / 2])
ax.plot(wval, zval, "k-")
ax.axis("equal")