summaryrefslogtreecommitdiff
path: root/test4.sce
diff options
context:
space:
mode:
Diffstat (limited to 'test4.sce')
-rw-r--r--test4.sce702
1 files changed, 702 insertions, 0 deletions
diff --git a/test4.sce b/test4.sce
new file mode 100644
index 0000000..28ca209
--- /dev/null
+++ b/test4.sce
@@ -0,0 +1,702 @@
+test_pass=[]
+res=[]
+
+
+test4=0
+
+
+//<----------------test case for cceps------------------>
+//x=[1 2 3];
+//correct=1;
+
+//vp=cceps(x, correct);
+//vi=[ 1.92565
+// 0.96346
+// -1.09735];
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('cceps test failed.');
+//end
+
+
+
+//<----------------test case for clustersegment------------------>
+// x=[0,1,0,0,1,1]
+//vp = clustersegment(x)
+//
+//vi=[2 5;2 6]
+//
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('clustersegment test failed.');
+//end
+//
+//<----------------test case for cmorwavf------------------>
+lb=1; ub=2; n=1; fb=3;
+fc=4;
+
+[a,b]=cmorwavf(lb,ub,n,fb,fc)
+
+a1=0.0858628;
+b1=2;
+
+a1=round(a1*100)/100;
+a=round(a*100)/100;
+b1=round(b1*100)/100;
+b=round(b*100)/100;
+
+
+if(a==a1) then
+ if (b==b1) then
+ test_pass=[test_pass,1];
+ end
+else
+ test_pass=[test_pass, 0];
+ disp('cmorwavf test failed.');
+end
+
+////<----------------test case for czt------------------>
+// x=[4 5 6 3 2]; fs = 1000; f1 = 100; f2 = 150;
+// m = 8;
+// w = exp(-%i*2*%pi*(f2-f1)/(m*fs));
+//a = exp(%i*2*%pi*f1/fs);
+//
+// vp = czt(x, m, w, a);
+////
+//vi=[7.3541-12.6740*%i 6.2892-12.5620*%i 5.2710-12.3493*%i 4.3097-12.04432*%i 3.4142-11.6568*%i 2.5919-11.1975*%i 1.84859-10.6776*%i 1.18835-10.1088*%i
+//]
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+////
+////
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('czt test failed.');
+//end
+//
+
+
+//<----------------test case for diric------------------>
+ x=[1 2 3];
+ n=3;
+
+ vp=diric(x,n);
+
+vi=[0.6935349 0.0559021 -0.3266617]
+vi=round(vi*100)/100;
+vp=round(vp*100)/100;
+
+if(vp==vi)
+ test_pass=[test_pass,1];
+else
+ test_pass=[test_pass, 0];
+ disp('diric test failed.');
+end
+
+//<----------------test case for dst1 ------------------>
+// x=[3 5 2 7];
+// n=3;
+//
+// vp=dst1(x,n);
+//
+//vi=[ 8.5355339
+// 1.
+// - 1.4644661 ]
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('dst1 test failed.');
+//end
+
+
+
+//<----------------test case for fft1 ------------------>
+//x = [1 2 3; 4 5 6; 7 8 9]
+//n = 3
+//dim = 2
+//vp=fft1 (x, n, dim)
+//
+//vi=[6.0 -1.5000000000000009+0.8660254037844375*%i -1.4999999999999984-0.8660254037844406*%i
+//15.0 -1.5000000000000018+0.8660254037844357*%i -1.4999999999999973-0.8660254037844428*%i
+//24.0 -1.5000000000000018+0.8660254037844348*%i -1.4999999999999964-0.8660254037844455*%i]
+//
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('fft1 test failed.');
+//end
+//
+
+
+
+
+//<----------------test case for fftconv ------------------>
+//x=[1, 2, 3];
+//y=[3,4,5]
+//vp=fftconv(x, y)
+//
+//vi=[3; 10; 22; 22; 15]
+//
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('fftconv test failed.');
+//end
+
+//<----------------test case for fftn ------------------>
+//x=[2 3 4];
+//siz=[1 5]
+//
+//vp=fftn(x,siz);
+//
+//
+//vi=[9.0000+0.0000*%i; -0.3090-5.2043*%i ; 0.8090+2.0409*%i; 0.8090-2.0409*%i; -0.3090+5.2043*%i]
+//
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('fftn test failed.');
+//end
+//
+
+//<----------------test case for fht ------------------>
+//x=1:4;
+//vp=fht(x)
+//
+//
+//vi=[10; -4; -2; 0]
+//
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('fht test failed.');
+//end
+
+
+//<----------------test case for filter1 ------------------>
+//b=[1,2,3];
+//a= [3,4,5];
+//x= [5,6,7];
+//[a b]=filter1(b, a, x)
+//
+//
+//a1=[1.6666667 3.1111111 4.4074074];
+//b1=[- 0.3950617; - 0.3456790];
+//
+//a=round(a*100)/100;
+//a1=round(a1*100)/100;
+//b=round(b*100)/100;
+//b1=round(b1*100)/100;
+//
+//
+//if(a==a1) then
+// if (b==b1) then
+// test_pass=[test_pass,1];
+// end
+//else
+// test_pass=[test_pass, 0];
+// disp('filter1 test failed.');
+//end
+
+//<----------------test case for filtic ------------------>
+b=[%i,1,-%i,5];
+a= [1,2,3*%i];
+y= [0.8*%i,7,9];
+
+vp=filtic(b,a,y)
+
+
+
+vi=[- 22.6*%i; 2.4; 0 ]
+
+vi=round(vi*100)/100;
+vp=round(vp*100)/100;
+
+
+if(vp==vi)
+ test_pass=[test_pass,1];
+else
+ test_pass=[test_pass, 0];
+ disp('filtic test failed.');
+end
+
+//<----------------test case for fractdiff ------------------>
+x=[4 5 6];
+d=3;
+
+vp=fractdiff(x,d)
+
+vi=[]
+
+vi=round(vi*100)/100;
+vp=round(vp*100)/100;
+
+
+if(vp==vi)
+ test_pass=[test_pass,1];
+else
+ test_pass=[test_pass, 0];
+ disp('fractdiff test failed.');
+end
+
+
+//<----------------test case for gauspuls ------------------>
+t=[1 2 3];
+fc=1; bw=1;
+vp=gauspuls(t,fc,bw)
+
+
+vi=[0.0281016 0.0000006 1.093D-14]
+
+vi=round(vi*100)/100;
+vp=round(vp*100)/100;
+
+
+if(vp==vi)
+ test_pass=[test_pass,1];
+else
+ test_pass=[test_pass, 0];
+ disp('gauspuls test failed.');
+end
+
+
+//<----------------test case for gaussian ------------------>
+//m=5;a=6;
+//vp = gaussian (m, a)
+//
+//vi=[ 5.380D-32; 1.523D-08; 1; 1.523D-08; 5.380D-32]
+//
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('gaussian test failed.');
+//end
+//
+////<----------------test case for idct1 ------------------>
+//x=[1 3 6];n=6;
+//vp=idct1(x,n)
+//
+//vi=[ 5.0812809
+// 1.6329932
+// - 2.143464
+// - 3.0400394
+// - 0.8164966
+// 1.7352157]
+//
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('idct1 test failed.');
+//end
+//
+
+//<----------------test case for idst1 ------------------>
+//x=[1 3 6];n=6;
+//vp=idst1(x,n)
+//
+//vi=[ 2.4654143
+// 1.8028286
+// - 0.6898286
+// - 1.4336286
+// 0.1315286
+// 1.1251286 ]
+//
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('idst1 test failed.');
+//end
+//
+
+//<----------------test case for ifft1 ------------------>
+//x = [1 2 3; 4 5 6; 7 8 9];
+//n = 3;
+//dim = 2
+//vp=ifft1 (x, n, dim)
+//
+//vi=[ 2 -0.5+0.2886751*%i -0.5-0.2886751*%i
+// 5 -0.5+0.2886751*%i -0.5-0.2886751*%i
+// 8 -0.5+0.2886751*%i -0.5-0.2886751*%i];
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('ifft1 test failed.');
+//end
+//
+
+//<----------------test case for ifht ------------------>
+//x = [1 2 3 4];
+//n = 2;
+//dim = 2;
+//
+//vp=ifht(x,n, dim)
+//
+//vi=[1.5 -0.5];
+//
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('ifht test failed.');
+//end
+//
+//<----------------test case for kaiserord ------------------>
+
+f=[1000, 1200]; m=[1, 0]; dev= [0.05, 0.05]; fs= 11025;
+
+[a b c d] = kaiserord (f, m, dev, fs)
+
+a1=70; b1=0.19955; c1=1.5099; d1='low'
+
+
+a=round(a*100)/100;
+a1=round(a1*100)/100;
+b=round(b*100)/100;
+b1=round(b1*100)/100;
+c=round(c*100)/100;
+c1=round(c1*100)/100;
+
+
+if(a==a1)
+ if (b==b1)
+ if (c==c1)
+ if (d==d1)
+ test_pass=[test_pass,1];
+end
+end
+end
+else
+ test_pass=[test_pass, 0];
+ disp('kaiserord test failed.');
+end
+
+
+//<----------------test case for meyeraux ------------------>
+x = [1 2 3];
+
+vp=meyeraux(x);
+
+vi=[1 -208 -10287];
+
+vi=round(vi*100)/100;
+vp=round(vp*100)/100;
+
+
+if(vp==vi)
+ test_pass=[test_pass,1];
+else
+ test_pass=[test_pass, 0];
+ disp('meyeraux test failed.');
+end
+
+//<----------------test case for morlet ------------------>
+//lb = [1 2 3];
+//ub=[1 2 3];
+//n=1;
+//
+//[a,b]=morlet(lb, ub, n);
+//
+//a1=[0.1720498
+// -0.1135560
+// -0.0084394];
+//
+// b1=[1;2;3]
+//
+//a=round(a*100)/100;
+//a1=round(a1*100)/100;
+//b=round(b*100)/100;
+//b1=round(b1*100)/100;
+//
+//
+//if(a==a1)
+// if (b==b1)
+// test_pass=[test_pass,1];
+// end
+//else
+// test_pass=[test_pass, 0];
+// disp('morlet test failed.');
+//end
+//
+//<----------------test case for rceps ------------------>
+
+//x =[1 4.0315 8.1095 11.9561];
+//
+//[a b]=rceps(x);
+//
+//a1=[2.4703812
+// 0.3236025
+// 0.1051662
+// 0.3236025];
+//
+// b1=[12.240472
+// 7.7643989
+// 3.7471281
+// 1.3451011]
+//
+//a=round(a*100)/100;
+//a1=round(a1*100)/100;
+//b=round(b*100)/100;
+//b1=round(b1*100)/100;
+//
+//
+//if(a==a1)
+// if (b==b1)
+// test_pass=[test_pass,1];
+// end
+//else
+// test_pass=[test_pass, 0];
+// disp('rceps test failed.');
+//end
+//
+
+
+
+//<----------------test case for rectpuls------------------>
+x = [1 2 3];
+
+vp=meyeraux(x);
+
+vi=[1 -208 -10287];
+
+vi=round(vi*100)/100;
+vp=round(vp*100)/100;
+
+
+if(vp==vi)
+ test_pass=[test_pass,1];
+else
+ test_pass=[test_pass, 0];
+ disp('rectpuls test failed.');
+end
+
+
+//<----------------test case for sawtooth------------------>
+t = [1 2 3 4 5];
+width=0.5;
+vp=sawtooth(t,width);
+vi=[-0.36338 0.27324 0.90986 0.45352 -0.18310]
+
+vi=round(vi*100)/100;
+vp=round(vp*100)/100;
+
+
+if(vp==vi)
+ test_pass=[test_pass,1];
+else
+ test_pass=[test_pass, 0];
+ disp('sawtooth test failed.');
+end
+
+
+
+//<----------------test case for sgolay------------------>
+p=1; n=3; m=0;
+vp=sgolay (p, n, m)
+
+vi=[0.83333 0.33333 -0.16667
+ 0.33333 0.33333 0.33333
+ -0.16667 0.33333 0.83333]
+
+vi=round(vi*100)/100;
+vp=round(vp*100)/100;
+
+
+if(vp==vi)
+ test_pass=[test_pass,1];
+else
+ test_pass=[test_pass, 0];
+ disp('sgolay test failed.');
+end
+
+//<----------------test case for sgolayfilt------------------>
+//x=[1;2;4;7]; p=0.3; n= 3; m=0; ts=0;
+//vp=sgolayfilt (x, p, n, m, ts)
+//
+//vi=[ 2.3333333
+// 2.3333333
+// 4.3333333
+// 4.3333333]
+//
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('sgolayfilt test failed.');
+//end
+//
+//<----------------test case for shanwavf------------------>
+//lb=2; ub=8; n=3; fb=1; fc=6
+//[a b]=shanwavf(lb,ub,n,fb,fc)
+//
+//a1 = [-3.8982e-17 + 1.1457e-31*%i 3.8982e-17 - 8.4040e-31*%i -3.8982e-17 + 4.5829e-31*%i]
+// b1 = [2 5 8]
+//a=round(a*100)/100;
+//a1=round(a1*100)/100;
+//b=round(b*100)/100;
+//b1=round(b1*100)/100;
+//
+//if(a==a1)
+// if (b==b1)
+// test_pass=[test_pass,1];
+// end
+//else
+// test_pass=[test_pass, 0];
+// disp('shanwavf test failed.');
+//end
+//
+//<----------------test case for sinewave------------------>
+//x=8; n= 4;
+//vp=sinewave(x, n)
+//
+//vi=[0 1 1.225D-16 -1 -2.449D-16 1 3.674D-16 -1]
+//
+//vi=round(vi*100)/100;
+//vp=round(vp*100)/100;
+//
+//
+//if(vp==vi)
+// test_pass=[test_pass,1];
+//else
+// test_pass=[test_pass, 0];
+// disp('sinewave test failed.');
+//end
+//
+//<----------------test case for sos2tf------------------>
+sos = [1 1 1 1 0 -1; -2 3 1 1 10 1];
+[a b]=sos2tf(sos)
+
+a1 =[ -2 1 2 4 1];
+b1 =[ 1 10 0 -10 -1];
+
+a=round(a*100)/100;
+b=round(b*100)/100;
+a1=round(a1*100)/100;
+b1=round(b1*100)/100;
+
+
+if(a==a1)
+ if (b==b1)
+ test_pass=[test_pass,1];
+ end
+else
+ test_pass=[test_pass, 0];
+ disp('sos2tf test failed.');
+end
+
+//<----------------test case for sos2zp------------------>
+sos = [1,2,3,4,5,6];
+[a b c ] = sos2zp (sos)
+
+a1 =[-1.0000 + 1.4142*%i
+-1.0000 - 1.4142*%i];
+b1 =[-0.6250 + 1.0533*%i
+ -0.6250 - 1.0533*%i];
+c1 = 1;
+
+a=round(a*100)/100;
+b=round(b*100)/100;
+a1=round(a1*100)/100;
+b1=round(b1*100)/100;
+
+c1=round(c1*100)/100;
+c=round(c*100)/100;
+
+
+if(a==a1)
+ if (b==b1)
+ if (c==c1)
+ test_pass=[test_pass,1];
+end
+end
+else
+ test_pass=[test_pass, 0];
+ disp('sos2zp test failed.');
+end
+
+//<----------------test case for welchwin------------------>
+x=4; opt="symmetric";
+vp=welchwin(x, opt)
+
+vi=[0.
+ 0.8888889
+ 0.8888889
+ 0.]
+
+vi=round(vi*100)/100;
+vp=round(vp*100)/100;
+
+
+if(vp==vi)
+ test_pass=[test_pass,1];
+else
+ test_pass=[test_pass, 0];
+ disp('welchwin test failed.');
+end
+
+res=find(test_pass==0)
+
+if(res~=[])
+ disp("One or more tests failed in test1")
+ test4=1;
+end
+
generated by cgit (git 2.34.1) at 2024-12-19 13:08:44 +0000