Hi Brad,
I have a question about the use of powerLawKramersKronig().
I work at a frequency of 3 MHz -> w = 2*pi*3e6 rad/s
I consider a reference frequency of 1 MHz -> w0 = 2*pi*1e6 rad/s
Let's take a reference sound speed of c0 = 1500 m/s
Then, in liver [Pohlhammer 1991] :
a0 = 4.3 Np/m/MHz^y
y = 1.266
Now let's use the function : c_kk = powerLawKramersKronig(w, w0, c0, a0, y)
c_kk = powerLawKramersKronig(2*pi*3e6, 2*pi*1e6, 1500, 4.3, 1.266);
This gives me c_kk = -0.0047 m/s, which I don't understand.
Is there a problem in the values I use?
Best regards,
Anthony