diff options
| author | Christian Kolset <christian.kolset@gmail.com> | 2025-11-14 16:44:37 -0700 |
|---|---|---|
| committer | Christian Kolset <christian.kolset@gmail.com> | 2025-11-14 16:44:37 -0700 |
| commit | 1873d28735e18866f7e0167045611440e4cbe20e (patch) | |
| tree | e02fa642255bcc115e2a601a9d492863fe537caf /tutorials/module_4/spectroscopy_problem | |
| parent | bda9939b6e93e06f61d5a51d90d6fb4a75d4aab8 (diff) | |
First draft of OOP tutorial complete AND worked on spectroscopy problem
Diffstat (limited to 'tutorials/module_4/spectroscopy_problem')
| -rw-r--r-- | tutorials/module_4/spectroscopy_problem/spectroscopy.py | 20 |
1 files changed, 19 insertions, 1 deletions
diff --git a/tutorials/module_4/spectroscopy_problem/spectroscopy.py b/tutorials/module_4/spectroscopy_problem/spectroscopy.py index f32ca52..f66734a 100644 --- a/tutorials/module_4/spectroscopy_problem/spectroscopy.py +++ b/tutorials/module_4/spectroscopy_problem/spectroscopy.py @@ -49,8 +49,25 @@ plt.show() "Calibrate length dimension" -np.polyfit(df_Hg['Pixels'],df_Hg['Intensity'],3) +C = np.polyfit(df_Hg['Pixels'],df_Hg['Intensity'],3) +lambda_Hg = lambda p : C[3] + C[2]*p+C[1]*p**2+C[0]*p**3 + +lambda_calibrated = lambda_Hg(df_Hg['Pixel']) + + +# Plot Hg intensity-wavelength plot +plt.figure(figsize=(8,5)) +plt.plot(lambda_Hg, df_Hg['Intensity'], linestyle='-') +plt.xlabel('Wavelength [nm]') +plt.ylabel('Intensity [a.u.]') +plt.title('Wavelength-Intensity (Hg)') +plt.grid(True) +plt.show() + + + +""" "Calibrate intensity dimension" epsilon = # Surface emissivity lambd = @@ -59,3 +76,4 @@ I_W_true = epsilon*((2hc**2)/(lambd**5)*1/(e**(hc/kT)-1)) I_W_meas= df_Ox['I_Tungsten [a.u.]'] R=I_W_meas/I_W_true I_plasma_meas=R*I_plasma_true +"""
\ No newline at end of file |