Hi,
here is a practical problem we would like to solve. In a pneumatic post the
acceleration of the capsule is measured and plotted over time. From the
graph achieved we would like to derive some kind of statistic value that
describes the stress the capsule, or what is in it, is exhibited to.
The amount of stress introduced to the capsule will probably depend on two
things, the maximum amplitudes observed AND the frequency of these
amplitudes. A quick note at this point: Not the acceleration by itself but
the quick acceleration changes are the components which make up the special
stress we want to quantify.
So what I would like to do is:
1) Apply a fourier transformation to the acceleration profile to
2) get a number of harmonic waves describing my graph
3) and use the amplitudes of my waves in a weighted fashion to
calculate some statistical value.
What I tried to do is:
A) construct an artificial profile fg for testing purpose like
a. f1 <- function(x) 0.5*sin(3*x + pi)
b. f2 <- function(x) sin(x + 0.75*pi)
c. f3 <- function(x) 1.5*sin(0.45*x + 0*pi)
d. fg <- function(x) f1(x) + f2(x) + f3(x)
B) try to reconstruct the components of fg with fft or fda.
What I don't understand yet is how:
X) in my test example I can define the amount of harmonic components,
because here I know that number. Of course afterwards in my natural profiles
I won't know.
Y) I have to transform the values I get out of the "fft" and "fourier"
functions to estimate the frequency, amplitude and phase of my harmonics.
Is there a good function in these packages I can use for that?
Best regards,
Benjamin
--
Pflichtangaben gemäß Gesetz über elektronische Handelsregister und
Genossenschaftsregister sowie das Unternehmensregister (EHUG):
Universitätsklinikum Hamburg-Eppendorf
Körperschaft des öffentlichen Rechts
Gerichtsstand: Hamburg
Vorstandsmitglieder:
Prof. Dr. Jörg F. Debatin (Vorsitzender)
Dr. Alexander Kirstein
Ricarda Klein
Prof. Dr. Dr. Uwe Koch-Gromus
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
