1 files changed, 45 insertions, 35 deletions
diff --git a/tests/unit_tests/test9999_WorkingDir/scilabcode/mainfunction.sci b/tests/unit_tests/test9999_WorkingDir/scilabcode/mainfunction.sci
index 5e2a7728..5719fc08 100644
--- a/tests/unit_tests/test9999_WorkingDir/scilabcode/mainfunction.sci
+++ b/tests/unit_tests/test9999_WorkingDir/scilabcode/mainfunction.sci
@@ -1,49 +1,59 @@
 //SCI2C: DEFAULT_PRECISION= DOUBLE
 function mainfunction()
-
-// ---------------------------------------
-// --- Initialization of the parameters. ---
-// ---------------------------------------
+//---> MOREA: commentare meglio lo scopo di questa funzione
+//---> NUTRICATO: annotare anche le funzioni di plot facendole puntare a funzioni dummy che non fanno nulla.
+//---> NUTRICATO: la convoluzione e' stata stostituita con una custom in quanto non ancora implementata da INRIA.
 
 
-// ---------------------------------------
-// --- Set number of point for quadrature .
-// --- interval for quadrature [a,b]
-// ---------------------------------------
+// sampling parameters
+fs   = 1000; //Hz
+Tmax = .01; //----> NUTRICATO: VALORE VERO 1 
+t    = 0 : 1 ./ fs : Tmax;
 
+// chirp parameters
+fo = 10
+k  = 150
 
+tau=0.7
 
-Np =10 // number of point for quadrature 
-a=1
-b=5
- 
+c=exp(2*%i*%pi*(fo*t+k./2*t.^2));//chirp signal in complex form
+// figure
+// subplot(2,2,1)
+// plot(t(1:250),real(c(1:250)));
+//title('chirp(t)')
 
-// ---------------------------------------
-// --- Compute Legendre-Gauss-Lobatto nodes, weights 
-//. ---In standard [-1 1] interval
-// ---------------------------------------
+cr=exp(2*%i*%pi*(fo*(t-tau)+k./2*(t-tau).^2));//chirp signal delayed in complex form
+// subplot(2,2,2)
 
-[x,w]=lglnodes(Np);
-disp('x');
-disp(x);
-disp('w');
-disp(w);
-// ---------------------------------------
-// --- Compute the nodes  
-//. ---In [a,b]  pick a t points in the 
-// --- desired interval
-// -------------------------------------
+// plot(t(1:250),real(c(1:250)),'r',t(1:250),real(cr(1:250)),'g');
+// title('chirp(t) e chirp(t-tau)')
 
-t=((b-a)./2)*x+(b+a)./2;
+// filter response 
+// subplot(2,2,3)
+h=conj(exp(2*%i*%pi*(fo*(-t)+k./2*(-t).^2)));
 
 
+// plot(t(1:250),real(h(1:250)));
+// title('conj(chirp(-t))')	
+//---> NUTRICATO: la convoluzione e' stata commentata in quanto non ancora implementata.
+y=myconvol(h,cr);
+
+
+tch=0:1 ./ fs:2*Tmax;
+// subplot(2,2,4)
+// plot(tch,abs(y))
+// title(' y(t)')
+
+disp('t')
+disp(t)
+disp('c')
+disp(c)
+disp('cr')
+disp(cr)
+disp('h')
+disp(h)
+disp('y')
+disp(y)
+endfunction
 
-Iab=((b-a)./2)*mysum(w.*Fint(t),11);
 
-disp('Number of nodes for quadrature ');
-disp(Np)
-disp('Integral value in [a b] with Np points');
-disp(Iab)
-
-
-endfunction