写在前面
好好学习,走出宿舍,走向毕设!
一些心路历程记录,很少有代码出现
因为鬼知道哪条代码到时候变成毕设的一部分了咧,还是不要给自己的查重挖坑罢了
23.2.27
文件批量重命名
为了给学姐先整出来一批训练数据,从学长给的文件里左扣扣,右扣扣给抠出来了
昨天明明已经给学姐讲过一遍咋搞数据了,她似乎没听懂。。。当然,部分原因在于她不是物理系的,可能读文章确实有点困难。
总之去网上抄了抄代码,把几个文件夹中需要的图片取出,重新标号即可~
import os.path
img_folder = 'D:\\used' #图片的文件夹路径
num = 131 #起始编号
for img_name in os.listdir(img_folder): # os.listdir(): 列出路径下所有的文件
src = os.path.join(img_folder, img_name) #src:要修改的目录名
dst = os.path.join(img_folder, str(num) + '.jpg') #dst: 修改后的目录名
num= num+1
os.rename(src, dst) #用dst替代src
23.2.28
2D FFT编写
一个从2.27下午开始到2.28下午逐渐开窍的过程:
不过在搜到的了对的代码以后,害嗨嗨,我会咯,画出来看着蛮像的图咯~
参考论文中的图片:
相关文章:
[1] Savage J R . Melting of 2D Colloidal Crystals[C]// Aps March Meeting. American Physical Society, 2005.
鄙人画出来的图片:
叠在一起,很能说明问题,嘿嘿嘿
2月干活时间表
突然发现学校的校园网自助管理系统可以看到网络使用情况
那就约等于说明了我2月份花了多少时间在论文上:
raw data 长这样:
稍加处理后长这样:
把时间戳处理一下,就知道自己一天赖在老师给的实验室里多久咧~
import pandas as pd
import numpy as np
import time
data = pd.read_excel('D:\\230301\\February.xlsx')
for i in range(0,len(data),2):
dt1 = time.strptime(data['注销时间'][i][:-3],'%Y-%m-%d %H:%M:%S')
dt2 = time.strptime(data['登录时间'][i+1][:-3],'%Y-%m-%d %H:%M:%S')
d = time.mktime(dt1)-time.mktime(dt2)
h = d//(60*60)
d = d-60*60*h
m = d//60
s = d-60*m
print("{}在实验室待了约{}小时{}分钟".format(data['注销时间'][i][5:10],int(h)-2,int(m)))
把干饭时间去掉了,就可以发现,其实还好,每天也没干很多时间
23.3.1
附上一份良师益友的sensei给的理论性blog:图像傅立叶变换(二维傅立叶变换)的原理和物理意义
通过FFT导出G
向量G在论文中的描述是:
翻译一下就是:向量G是满足周期性边界条件的波矢量,处于静态结构因子的第一个峰值。
如果学过初等的固体物理,那在晶体方面理解这句话并不难,其实本质上来说这个G就是倒格子。但是在液相,非晶体的状态下,也可以用这种做法来做,将其理解为文章中写的波矢即可。与晶体的区别,在于液相中的波矢方向随意。但由于程序编写的需要,我将两者的方向设定一致。
代码写完了以后,意识到了坐标有点问题,遂开始验证
- 画出g®【径向分布函数】,获得颗粒间距(第一个峰)
- 参考时域FFT做法,换算单位
- 将颗粒间距与FFT第一个峰相乘,检验一下是否是2π
今天搞完了第一步,但是效果不理想,图片巨丑: