Hi,
Please, I would like to know how to implement in Matlab the real and imaginary parts (Gr and Gi) of the Fourier transform given in equations 7 of the paper of Stokes 1948 (A Numerical Fourier-analysis Method for the Correction of Widths and Shapes of Lines on X-ray Powder Photographs), where the result is given in Table 3 of that paper ?
Thank you for your consideration.
Best regards,
Fayçal Hadj-Larbi
Ts = 1/50; t = 0:Ts:10-Ts; x = sin(2*pi*15*t) + sin(2*pi*20*t); plot(t,x) xlabel('Time (seconds)') ylabel('Amplitude')
Hi Dear Elvir,
Sorry for the late response.
The following is the way I implemented the equations Gr and Gi to get the the columns 4 and 5 of table 3;
t_G= [0 : 1 : 20];
L_tG = 60;
Factor_G = 0.1;
for i=1:length(t_G)
Gr(i)=Factor_G.*sum(g.*cos(2.*pi.*t_G(i).*x_g./L_tG));
Gi(i)=Factor_G.*sum(g.*sin(2.*pi.*t_G(i).*x_g./L_tG));
end
Actually, I have some troubles with the implementation of f (after have succesfully calculated the quantities : Hr, Hi, Fr and Fi) which has 41 columns.
Thank you for your interest.
Best regards,
Fayçal
- Badges
- Science topic
- Similar topics
- Mathematics
- Statistics
- Sampling Design
- Sample Size