]> git.parisson.com Git - cnaq.git/commitdiff
git-svn-id: http://svn.parisson.org/svn/CNAQ/trunk@160 5fc3e0e6-29bc-4d03-b52b-c088cb...
authormanu <manu@5fc3e0e6-29bc-4d03-b52b-c088cb822bde>
Fri, 4 Apr 2008 09:23:14 +0000 (09:23 +0000)
committermanu <manu@5fc3e0e6-29bc-4d03-b52b-c088cb822bde>
Fri, 4 Apr 2008 09:23:14 +0000 (09:23 +0000)
tools/RI_FT.m

index dd3cb761eb4e9f68f92cb28fda0426c1e6b85b6c..ccdb3a9e645cfc4f464f49e722cb1a09a5648580 100644 (file)
@@ -6,11 +6,15 @@ function [ri,ft]=RI_FT(f, sig_exc, sig_mes, f_s)
 DC = mean(sig_mes,1);\r
 sig_mes = sig_mes - DC;\r
 len_f = length(f);\r
-f0 = f(1);\r
-% f1 = f(length(f));\r
+f0 = f(1);\r
+f1 = f(len_f);\r
 \r
 % Methode du produit\r
 siginv=fliplr(sig_exc)./f;\r
 ft=fft([sig_mes; zeros(len_f-1,1)]).*fft([siginv; zeros(len_f-1,1)]);\r
 ri=real(ifft(ft));\r
-\r
+norm=sqrt(sum(abs(ri.^2))/siz);\r
+yeff=norm*sqrt(f_s/2/(f1-f0));\r
+scal=f_s/yeff;\r
+ri=ri*scal;\r
+ft=ft/yeff;
\ No newline at end of file