FS_few_cycle/integration_error_calc.py at main · Jrcarand/FS_few_cycle · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
# -*- coding: utf-8 -*-
"""
Created on Tue Apr  4 11:59:39 2023

@author: cardoso
"""

import numpy as np
import matplotlib.pyplot as plt


def polsin(x,phi0):
    return np.sin(((x-phi0)*np.pi/180))**2

N=18000

ttmax=18

tt=np.linspace(0,18,N)

vv=tt*10

speed=0.66
idx=np.argmax(vv>speed)

vv[idx:]=speed

vv[-idx:]=speed-(tt[-idx:]-tt[-idx])*10

ph=109+np.cumsum(vv)*tt[1]


idx=np.arange(0,int(N/15))*15

polout=polsin((ph),106)

specint=[]
calint=[]

for i in idx:

    specint.append(np.sum(polout[i:i+14]))
    calint.append(np.sum(polout[i]*np.ones(14)))

specint=np.array(specint)
calint=np.array(calint)

ratio=specint/calint

tredux=tt[idx]


plt.close(123)
plt.figure(123)
# plt.plot(tt,ph)

plt.plot(tredux,ratio)