xd

analisis/plots/20240312_RotationalDopplerShift_news/RDS_piezobeamsizes_2024_invariance.py

@@ -12,7 +12,7 @@ from scipy import interpolate
 
 #C:\Users\Usuario\Documents\artiq\artiq_experiments\analisis\plots\20220106_CPT_DosLaseres_v08_TISA_DR\Data
 
-os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20240312_RotationalDopplerShift_news/Data')
+#os.chdir('/home/nico/Documents/artiq_experiments/analisis/plots/20240312_RotationalDopplerShift_news/Data')
 
 
 """
@@ -468,6 +468,45 @@ plt.grid()
 
 
 
+#%%
+"""
+Va otro modelo con funciones de bessel
+
+
+F(r) = A*(J_0 (b/r))²
+
+"""
+
+from scipy.special import jv
+
+
+def modelo3(r,a,b,c,d):
+    return a*((jv(0, b*(1/r)))**2)+c*((jv(0, d*(1/r)))**2)  
+
+rfit = np.arange(x1,len(Intensityver1)*20+x1,20)
+yfit = [p for p in pmdepthsdrver1]
+
+rlong = np.arange(0,1000,1)
+
+popt,pcov = curve_fit(modelo3,rfit,yfit,p0=(0.7,1))
+
+print(popt)
+plt.figure()
+
+plt.errorbar(np.arange(x1,len(Intensityver1)*20+x1,20), [p for p in pmdepthsdrver1], yerr= errorpmdepthsdrver1, fmt='o',color='red', capsize=3, markersize=8,zorder=1)
+plt.plot([x*20 for x in xchicofinal],[i/np.max(IntensityChico) for i in IntensityChico],'o',color='red',alpha=0.3)
+
+plt.errorbar(np.arange(x2,len(Intensityver2)*20+x2,20), [p for p in pmdepthsdrver2], yerr= errorpmdepthsdrver2, fmt='o',color='blue', capsize=3, markersize=8,zorder=1)
+plt.plot([x*20 for x in xgrandefinal],[i/np.max(IntensityGrande) for i in IntensityGrande],'o',color='blue',alpha=0.3)
+
+plt.plot(rlong,modelo3(rlong,*popt),zorder=2,color='black',linewidth=3,linestyle='dashed',label=r'$I(r)=a\,(J_0(\frac{b}{r}))^2$')
+
+plt.xlim(-20,600+x2+10)
+plt.legend()
+plt.grid()
+
+
+