Difference between revisions of "EGR 103/Fall 2019/Minilab 3"
Jump to navigation
Jump to search
(Created page with "* Problem 2.5.1 ** To load data froman Excel file with headers: ::<syntaxhighlight lang=python> import pandas # %% Load data edata = pandas.read_excel('file.xlsx') col_1_stu...") |
|||
| Line 13: | Line 13: | ||
:* [[Python:Extrema]] | :* [[Python:Extrema]] | ||
:* [https://docs.scipy.org/doc/scipy-0.18.1/reference/generated/scipy.interpolate.CubicSpline.html scipy.interpolate.CubicSpline¶] at [docs.scipy.org docsscipy.org] | :* [https://docs.scipy.org/doc/scipy-0.18.1/reference/generated/scipy.interpolate.CubicSpline.html scipy.interpolate.CubicSpline¶] at [docs.scipy.org docsscipy.org] | ||
| + | |||
| + | * Problem 2.5.2 | ||
| + | ** [[https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.simps.html scipy.integrate.simps]] at [docs.scipy.org docsscipy.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: | ||
| + | <syntaxhighlight lang=python> | ||
| + | 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') | ||
| + | </syntaxhighlight> | ||
Revision as of 01:49, 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()
- Python:Plotting
- Python:Interpolation
- Python:Extrema
- scipy.interpolate.CubicSpline¶ at [docs.scipy.org docsscipy.org]
- Problem 2.5.2
- [scipy.integrate.simps] at [docs.scipy.org docsscipy.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')
