谐波周期以及周期判断

1.有关周期以及谐波、次谐波可看另外一部分：

周期信号的谐波与次谐波。

2.在判断是否是另外一个周期的谐波时，可判断色散和脉冲的轮廓峰数量；

假设一个脉冲星的真正的周期是p, 谐波周期是P, 那么P * m/n = p, 那么n个p和m个P的时间长度一致，找到最小公倍数。

那么P周期的峰的个数应该是p周期的n倍数个。

附上python：

 1 p =  5.34033
 2 
 3 import pandas as pd
 4 pd.set_option("display.max_rows",None)         ###在jupyter notebook上面显示每一行和每一列，这个设置似乎只有放在import后面才有用！！！！
 5 pd.set_option("display.max_columns",None)
 6 
 7 
 8 m = range(1,25,1)
 9 n = range(1,25,1)
10 
11 harm = []
12 M = []
13 N = []
14 
15 for i in range(1,len(m)+1,1):
16     for j in range(1,len(n)+1,1):
17         harm.append(i/j)
18         M.append(i)
19         N.append(j)
20         
21 cand_p = []
22 for i in range(len(harm)):
23     cand_p.append(harm[i]*p)
24 
25 
26 a=M
27 b=N
28 c=cand_p
29 c={
30    "harm": harm,
31    "M" : a,
32    "N" : b,
33    "cand_p" : c,
34   "peak_num" : N}
35 
36 
37 data=pd.DataFrame(c)#将字典转换成为数据框
38 data.drop_duplicates(subset=['cand_p'], keep='first', inplace=True)
39 
40 res = data.sort_values(by='cand_p', ascending=True)
41 #print('按照候选周期从小到大值排序:\n', res)
42 
43 #res.to_csv("D:\\FAST_document\\data_process\\NGC6517\\谐波周期\\"+'谐波周期候选1.csv',sep=',',index=0,header=1) 
44 
45 print(res)
View Code

 

posted @ 2022-10-13 22:47  yyyddj  阅读(204)  评论(0)  编辑  收藏  举报