summaryrefslogtreecommitdiff
path: root/test3.sce
diff options
context:
space:
mode:
Diffstat (limited to 'test3.sce')
-rw-r--r--test3.sce841
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
generated by cgit (git 2.34.1) at 2024-12-19 13:12:48 +0000