code changes by Shashikiran Yadalam during FOSSEE Fellowship 2018

1 files changed, 841 insertions, 0 deletions

diff --git a/test3.sce b/test3.sce
new file mode 100644
index 0000000..c3195d5
--- /dev/null
+++ b/test3.sce
@@ -0,0 +1,841 @@
+//exec FOSSEE_Scilab_Octave_Interface_Toolbox/loader.sce
+//exec("/home/shashi/Desktop/FOSSEE_Scilab_Octave_Interface_Toolbox/loader.sce")
+//getd macros/
+exec("loader.sce")
+exec("builder.sce")
+
+
+
+
+test_pass=[]
+res=[]
+
+
+
+
+
+
+/////////////////test case for pchip    /////////////////////////
+
+ x=[0 1 2 3 4 5]
+ y=[1 0 1 0 1 0]
+ xx=linspace(0,5,80);
+
+ vp=pchip(x,y,xx);
+
+vi=fscanfMat("macros/txt4_pchip.txt");
+
+vi=roundn(vi,4);
+vp=roundn(vp,4);
+
+if(vp==vi)
+    test_pass=[test_pass,1];
+else
+     test_pass=[test_pass, 0];
+    disp('pchip test failed.');
+end
+
+
+
+/////////Test case for       peak2rms                  //////////
+
+IN=[6 19 10 25];
+OUT=peak2rms(IN);
+//m=1.4927036;
+OUT=(round(OUT*100)/100);
+if(OUT==1.49)
+    test_pass=[test_pass ,1];
+
+else
+    test_pass=[test_pass, 0];
+    disp('peak2rms test failed.');
+end
+
+
+
+
+
+///////////////              Test case for       peig            /////////
+
+//fs = 100;
+//t = 0:1/fs:1-1/fs;
+//s = 2*sin(2*%pi*25*t)+sin(2*%pi*35*t);
+//[S,w]=peig(s,2,512,fs,'half');
+//n=length(S);
+//
+//for i=1:n
+//    S(i)=round((S(i)*10^7)/10^7);
+//    end
+//
+//
+//M=fscanfMat("macros/text4_peig.txt");
+//M=roundn(M,7);
+//if(M==S )
+//      test_pass=[test_pass 1];
+// else
+//     test_pass=[test_pass 0];
+//     disp("peig test failed");
+//     end
+//
+
+/////////////////test for periodogram   ///////////////////////////
+
+n = 0:319;
+x = cos(%pi/4*n);
+[pxx,w]=periodogram(x,ones(1,320),256,2000,"onesided");
+
+pxx=roundn(pxx,4);
+qxx=fscanfMat("macros/txt4_periodogram.txt");
+
+pxx=roundn(pxx,7);
+qxx=roundn(qxx,4);
+
+if(pxx==qxx)
+    test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for periodogram");
+end
+
+
+
+
+
+
+////////                   test case for pmusic.............//////////
+//n = 0:199;
+//x = cos(0.257*%pi*n) + sin(0.2*%pi*n);
+//[S,w]=pmusic(x,2,16,1);
+//S=roundn(S,6);
+//w=roundn(w,4);
+//
+//si=[2.6425624,5.7475005, 77.1482210,1.5296243,0.4725347,0.2848481,0.2508128,0.2731036,0.2950648]';
+//si=roundn(si,6);
+//
+//Wi=[0,0.0625,0.125,0.1875,0.25,0.3125,0.375,0.4375,0.5]';
+//Wi=roundn(Wi,4);
+//
+// if(si==S & w==Wi)
+//      test_pass=[test_pass 1];
+// else
+//     test_pass=[test_pass 0];
+//     disp("pmusic test failed");
+//     end
+//
+
+
+
+
+////////////                  test for poly2lsf     //////////////
+//
+X=[1 0.6149 0.9899 0 0.0031 -0.0082] ;
+lsf = poly2lsf(X);
+
+for i=1:length(lsf)
+    lsf(i)=round((lsf(i)*10^7))/10^7;
+end
+
+p=[0.7841731,1.5605415,1.8776459,1.8984313,2.3592523];
+p=round((p.*10^7))./10^7;
+if(lsf'==p) then
+    test_pass=[test_pass,1];
+
+else
+    test_pass=[test_pass,0];
+    disp("test for poly2lsf failed");
+end
+//
+//
+
+
+
+
+///////////       test for polyscale          ////////////////
+//
+x=[1 0 0 0 0 0 0 -1];
+p=polyscale(x,1);
+
+y=[1.    0.    0.    0.    0.    0.    0.   -1];
+
+if(p==y) then
+    test_pass=[test_pass,1];
+
+else
+    test_pass=[test_pass,0];
+    disp("test failed for polyscale ");
+    end
+
+
+
+
+////////test for polyval ///////////////////////
+
+p=[3 2 1];
+y=polyval(p,[ 5 7 9]);
+
+q=[86 162 262];
+
+
+if(y==q)  then
+    test_pass=[test_pass,1]
+
+else
+    test_pass=[test_pass,0];
+    disp("test failed for polyval");
+
+end
+
+
+
+////////////////test for prony /////////////
+V=filter([1,1],[1,1,2],[1 zeros(1,31)]);
+[b,a]=prony(V,1,2);
+
+b=round((b.*10)./10);
+a=round((a.*10)./10);
+
+
+if (b==[1,1] ) then
+    if (a==[1,1,2]) then
+
+      test_pass=[test_pass,1];
+    end
+else
+    test_pass=[test_pass,0];
+    disp("test failed for prony");
+end
+
+
+
+/////////////test for rc2poly///////////////
+k = [0.3090 0.9800 0.0031 0.0082 -0.0082];
+a = rc2poly(k);
+a=round((a.*10^7))./10^7;
+b= [1,0.6148162,0.9898814,0.0000243,0.0031580,-0.0082];
+b=round((b.*10^7))./10^7;
+
+if(a==b)  then
+       test_pass=[test_pass,1]
+
+else
+    test_pass=[test_pass,0];
+    disp("test failed for rc2poly");
+end
+
+
+//////////test for rcosdesign  //////////////
+rolloff = 0.25;
+span = 3;
+sps=2;
+b=rcosdesign(rolloff,span,sps);
+
+m=[-0.1210006  -0.0456421    0.4418023    0.7590604    0.4418023  -0.0456421  -0.1210006 ];
+
+m=round((m.*10^7))./10^7;
+b=round((b.*10^7))./10^7;
+
+if(m==b)
+    test_pass=[test_pass,1]
+
+else
+    test_pass=[test_pass,0];
+    disp("test failed for rcosdesign");
+end
+
+
+
+////////////test for rlevinson  ////////////
+X=[1 6/7 5/7 8/7 3/7 6/7];
+[R U kr e] = rlevinson(X, 0.3);
+
+:
+ei  =[ 0.3757546,0.0221076,-3.4125,1.1307692,0.3];
+ kri  =[
+  -0.2251908
+    0.9701364
+  -12.464286
+  -1.1538462
+    0.8571429  ];
+ Ui  =[1,-0.2251908,0.9701364,-12.464286,-1.1538462,0.8571429;0,1,-0.4436567,6.5,2,0.4285714;0,0, 1,-12.535714,-1,1.1428571;0,0,0,1,1.8461538,0.7142857;0,0,0,0,1,0.8571429;0,0,0,0,0,1];
+
+
+
+ Ri =[
+
+    0.3958273
+    0.0891367
+  -0.3444604
+    0.0362590
+  -0.1329496
+    0.1042446
+];
+
+Ri=round((Ri.*10^7)./10^7);
+R=round((R.*10^7)./10^7);
+Ui=round((Ui.*10^7)./10^7);
+U=round((U.*10^7)./10^7);
+kri=round((kri.*10^7)./10^7);
+kr=round((kr.*10^7)./10^7);
+ei=round((ei.*10^7)./10^7);
+e=round((e.*10^7)./10^7);
+
+if(Ui==U)
+    if(Ri==R)
+        if(kri==kr)
+            if(ei==e)
+                test_pass=[test_pass,1];
+    end
+end
+end
+
+else
+    test_pass=[test_pass,0];
+    disp("test failed for rlevinson");
+end
+
+
+
+/////////////      test for rooteig            //////////////
+
+   n=0:99;
+
+s=exp(1*%i*%pi/2*n)+2*exp(1*%i*%pi/4*n)+exp(1*%i*%pi/3*n);
+
+X = corrmtx(s,12,'mod');
+
+[W,P] = rooteig(X,3);
+Wi=[0.7745464,1.5717498,1.0481856]';
+Pi=[381.37294,103.13707,119.97409]';
+W=round(W*10^7)/10^7;
+P=round(P*10^5)/10^5;
+Wi=round(Wi*10^7)/10^7;
+Pi=round(Pi*10^5)/10^5;
+
+if(Wi==W)
+if(Pi==P)
+     test_pass=[test_pass,1];
+
+    end
+else
+    test_pass=[test_pass,0];
+    disp("test failed for rooteig");
+end
+
+
+
+///////////      test for rootmusic      ///////////
+
+ n=0:99;
+ s=exp(1*%i*%pi/2*n)+2*exp(1*%i*%pi/4*n)+exp(1*%i*%pi/3*n);
+ [A,R]=corrmtx(s,12,'mod');
+ [W,P] = rootmusic(R,3,'corr');
+Wi  =[0.7738111,1.5690374 ,1.0426234]';
+ Pi  =[377.4255,103.18124,123.86659]';
+
+ W=round(W*10^7)/10^7;
+ P=round(P*10^7)/10^7;
+ Wi=round(Wi*10^7)/10^7;
+ Pi=round(Pi*10^7)/10^7;
+
+if(Wi==W)
+    if(Pi==P)
+        test_pass=[test_pass,1];
+    end
+
+else
+    test_pass=[test_pass,0];
+    disp("test failed for rootmusic");
+end
+
+
+/////////////test for sampled2continuos  /////////
+y=sampled2continuous([1,2,3],5,6);
+//
+x=2.4166806;
+//
+y=roundn(y,7);
+x=roundn(x,7);
+//
+if(y==x)
+    test_pass=[test_pass,1];
+else
+    test_pass=[test_pass,0];
+    disp("test failed for sampled2continous");
+end
+
+
+////////////////test for schurrc    //////////
+:
+m=linspace(1,100);
+r = xcorr(m(1:5),'unbiased');.......//autocorrelation vector
+[k,e] = schurrc(r(5:$));
+
+k_exp= [-0.9090909  0.2222222  0.2244898  0.2434211]';
+ e_exp=1.6212406 ;
+
+k=roundn(k,7);
+e=roundn(e,7);
+
+
+k_exp=roundn(k_exp,7);
+e_exp=roundn(e_exp,7);
+
+if(k==k_exp)
+    if(e==e_exp)
+        test_pass=[test_pass,1];
+    end
+else
+    test_pass=[test_pass,0];
+    disp("test failed for schurrc");
+end
+
+
+
+/////////////////test for slewrate/////////
+x=[ones(1,25)  -ones(1,25)  ones(1,25) ];
+t=1:length(x);
+s=slewrate(x,t);
+s=round(s);
+if(s==[-2,2])
+    test_pass=[test_pass,1];
+
+else
+    test_pass=[test_pass,0];
+    disp("test failed for slewrate");
+
+    end
+
+
+//////////test for sos2ss////////////
+sos = [1  1  1  1   0  -1 ;
+      -2  3  1  1  10   1];
+[A,B,C,D] = sos2ss(sos,2);
+
+Ai=[- 10.    0.    10.    1. ; 1.     0.    0.     0.  ; 0.     1.    0.     0.  ;  0.     0.    1.     0.  ];
+Bi=[1 0 0 0];
+Ci=[42    4.  -32.  -2];
+Di=-4
+
+A=roundn(A,1);
+B=roundn(B,1);
+C=roundn(C,1);
+D=roundn(D,1);
+
+Ai=roundn(Ai,1);
+Bi=roundn(Bi,1);
+Ci=roundn(Ci,1);
+Di=roundn(Di,1);
+
+if(A==Ai)
+    if(B==Bi&C==Ci)
+            if(D==Di)
+                test_pass=[test_pass,1];
+            end
+        end
+    else
+        test_pass=[test_pass,0];
+        disp("test failed for sos2ss");
+    end
+
+
+//
+///////////test for sosbreak    ////////////
+v=[1+4*%s+6*%s^2+4*%s^3+%s^4];
+ [zerosort,g]=sosbreak(v);
+gi=1;
+bi=coeff(zerosort(1));
+ai=coeff(zerosort(2));
+
+if(g==gi)
+    if(ai==[1 2 1] & bi==[1 2 1])
+        test_pass=[test_pass,1];
+    end
+else
+    test_pass=[test_pass,0];
+    disp("test failed for sosbreak");
+end
+
+
+
+
+//////////////test for specgram//////////
+N = 1024;
+n = 0:N-1;
+w= 2*%pi/5;
+x = sin(w*n)+10*sin(2*w*n);
+s = specgram(x);
+s=matrix(s,768,1);
+
+s_real=roundn(real(s),4);
+s_img=roundn(imag(s),4);
+
+m=fscanfMat('macros/txt4specgram_real.txt');
+n=fscanfMat('macros/txt4specgram_imag.txt');
+
+m=roundn(m,4);
+n=roundn(n,4);
+
+if(or(s_real==m) & or(s_img==n))
+    test_pass=[test_pass,1];
+else
+    test_pass=[test_pass,0];
+    disp("test failed for specgram");
+end
+//
+
+///////////test for upsample fill////////
+
+g=upsamplefill([1,3,5],2,%f);
+
+gi=[1.    1.    1.    3.    3.    3.    5.    5.    5.  ];
+
+
+if(gi==g)
+    test_pass=[test_pass,1];
+else
+    test_pass=[test_pass,0];
+    disp("test failed for upsample fill");
+end
+
+///////////test for ss2sos///////////////
+a =[0.5095,0,0,0,0;
+0.3007, 0.2260, -0.3984, 0, 0;
+0.0977, 0.3984, 0.8706, 0, 0;
+0.0243, 0.0991, 0.4652, 0.5309, -0.4974;
+0.0079, 0.0322, 0.1512, 0.4974, 0.8384];
+
+
+b =[0.6936 0.1382 0.0449 0.0112 0.0036]';
+
+
+c =[0.0028 0.0114 0.0534 0.1759 0.6500];
+
+
+d =0.0013;
+
+a=roundn(a,3);
+b=roundn(b,3);
+c=roundn(c,3);
+d=roundn(d,3);
+
+[sos,g]=ss2sos(a,b,c,d);
+e=[0.509,0,0,0,0;0.301,0.226,-0.398,0,0;0.098,0.398,0.871,0,0;0.024,0.099,0.465,0.531,-0.497;0.008,0.032,0.151,0.497,0.838 ];
+
+if(g==0.001)
+    if (a==e) then
+        test_pass=[test_pass 1];
+    end
+else
+     test_pass=[test_pass,0];
+    disp("test failed for ss2sos");
+end
+
+
+///////////////test for upfirdn/////////////
+
+yout=upfirdn(1:100,1,1,1);
+eout=1:100;
+
+if (eout==yout)
+    test_pass=[test_pass 1];
+else
+    test_pass=[test_pass 0];
+    disp("test failed for upfirdn");
+end
+
+
+
+
+
+//////////////test for wreverse/////////////
+q=wrev([1 2 3]);
+qi=[3 2 1];
+
+if(q==qi)
+    test_pass=[test_pass,1];
+else
+    test_pass=[test_pass,0];
+    disp("test failed for wreverse");
+end
+
+
+///////////////test for zerocrossings///////
+
+x = linspace(0,1,100);
+y = 2*sin(2*%pi*x);
+x0= zerocrossing(x,y);
+
+x01=[0,0.5];
+
+if(x0==x01)
+    test_pass=[test_pass,1];
+else
+    test_pass=[test_pass,0];
+    disp("test failed for zerocrossings");
+end
+
+////////////////test for deconv ////////////
+[b, r] = deconv ([3, 6, 9, 9], [1, 2, 3]);
+
+bi=[3,0];
+ri=[0, 0 ,0  ,9];
+
+b=roundn(b,1);
+r=roundn(r,1);
+
+if(and(r==ri) & and(b==bi))
+     test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for deconv");
+end
+
+
+
+//////////////test for downsample////////
+g=downsample([1,2,3,4,5],2);
+gi=[1 3 5];
+
+g=round(g);
+
+if(g==gi)
+     test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for downsample");
+end
+
+//////////////////test for mpoles////////////
+
+[m,n]=mpoles([2 3 1 1 2]);
+
+mi=[1;1;2;1;2];
+ni=[2;5;1;4;3];
+
+m=round(m);
+n=round(n);
+
+
+if(and(m==mi) & and(n==ni))
+     test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for mpoles");
+end
+
+
+/////////////test for polyreduce///////////////
+
+b=polyreduce ([0, 0, 1, 2, 3]);
+
+b=round(b);
+
+if(b==[1,2,3])
+    test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for polyreduce");
+end
+
+
+
+
+////////////////test for prepad////////////////
+a=prepad ([1,2], 4,0,2);
+
+a=round(a);
+
+if(a==[0 0 1 2])
+    test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for prepad");
+end
+
+
+//////////////////test for residue/////////////
+b = [1, 1, 1];
+a = [1, -5, 8, -4];
+
+[r, p, k, e] = residue (b, a);
+
+
+ ei  =[1;2;1];
+ ki  =[];
+ pi  =[2;2;1];
+ ri  =[-2;7;3];
+
+ r=round(r);
+ e=round(e);
+ k=round(k);
+ p=round(p);
+
+
+ if(and(ri==r) & and(ei==e)  & and(pi==p) & and(ki==k))
+     test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for residue");
+end
+
+
+
+
+///////////////////test for residued/////////////////////
+
+[r,p,f,m] = residued([1 1],[1 -2 1]);
+
+mi=[1;2];
+fi=[];
+pi=[1;1];
+ri=[-1;2];
+
+r=round(r);
+p=round(p);
+f=round(f);
+m=round(m);
+
+ if(and(ri==r) & and(fi==f)  & and(pi==p) & and(mi==m))
+     test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for residue");
+end
+
+
+  ///////////////////test for residuez/////////////////////
+
+[r,p,f,m] = residuez([1 1],[1 -2 1]);
+
+mi=[1;2];
+fi=[];
+pi=[1;1];
+ri=[-1;2];
+
+r=round(r);
+p=round(p);
+f=round(f);
+m=round(m);
+
+ if(and(ri==r) & and(fi==f)  & and(pi==p) & and(mi==m))
+     test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for residue");
+end
+
+
+///////////////test for gmonopuls////////////////////////
+u=gmonopuls([1 2 3],[]);
+
+u=round(u);
+
+if(u==[0 0 0])
+     test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for gmonopuls");
+end
+
+
+
+///////////////test for convmtx//////////////////////
+o=convmtx([3;4;5],3);
+
+o=round(o);
+
+if(o==[3,0,0;4,3,0;5,4,3;0,5,4;0,0,5])
+    test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for convmtx");
+end
+
+
+/////////////////test for yulewalker/////////////
+A=yulewalker([1 2 3]);
+A=roundn(A,2);
+
+if (A==[1.33;0.33]) then
+    test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for yulewalker");
+end
+
+
+
+
+//////////////////////test for wkeep////////////////
+
+y=wkeep([1 2 3;4 5 6],[2 2]);
+yi=[1,2];
+
+if(yi==round(y))
+    test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for wkeep");
+end
+
+
+
+
+//////////////////test for rc2ac////////////////
+k = [0.3090 0.9800 0.0031 0.0082 -0.0082];
+r0 = 0.1;
+
+R = rc2ac(k,r0);
+
+R=roundn(R,4);
+
+ Ri  =[0.1 ;-0.0309;-0.0790948 ;0.0786627 ;0.0293629 ;-0.0950000];
+ Ri=roundn(Ri,4);
+
+ if(R==Ri)
+      test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for rc2ac");
+end
+
+
+
+//////////////////////test for mexihat//////////////////
+[a,b]= mexihat(1,2,3);
+ai=[0.00000  -0.35197  -0.35214];
+bi=[1.0000   1.5000   2.0000];
+
+a=roundn(a,4);
+b=roundn(b,4);
+
+ai=roundn(ai,4);
+bi=roundn(bi,4);
+
+if(and(ai==a) & and(bi==b))
+     test_pass=[test_pass,1];
+else
+   test_pass=[test_pass,0];
+    disp("tst failed for mexihat");
+end
+
+
+
+
+
+disp("all worked!!")
+
+////////////////////////////////////////////////////
+
+res=find(test_pass==0)
+
+
+
+
+if(res~=[])
+	disp("One or more tests failed")
+	exit(1)
+else
+    disp("pass")
+	exit
+end