test_pass=[] res=[] test1=0 exec loader.sce /////////Test case for 2) arburg ////////// a = octave_fun("arburg","signal", [1,2,3,4,5],2); a = round(a*10000)/10000; if(a == [1. -1.8639 0.9571]) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("arburg test failed") end /////////Test case for 4) aryule ////////// a = octave_fun("aryule","signal", [1,2,3,4,5],2); a = round(a*10000)/10000; if(a == [1. -0.814 0.1193]) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("aryule test failed") end /////////Test case for 38)levinson ////////// a = [1 0.1 -0.8]; v = 0.4; rand("seed", 0); w = sqrt(v)*rand(15000,1); x=octave_fun("filter","signal",1,a,w); r=octave_fun("xcorr","signal",x); ar=octave_fun("levinson","signal",r,length(a)-1); ar = round(ar*10000)/10000; if(ar == [1 0.0984 -0.793]) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("levinson Test failed") end /////////Test case for **)sigmoid_train ////////// s = octave_fun("sigmoid_train",0.1,[1:3],[4]); s = round(s*10000)/10000 if(s == 0.2737) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("sigmoid_train Test failed") end /////////Test case for **)circshift ///////// M = [1 2 3 4]; R = octave_fun("circshift",M,[0 1]); if(R == [4 1 2 3]) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("circshift failed") end ///////////////////////////////////////////// /////////Test case for **)kaiser ////////// win = octave_fun("kaiser",6,[0.2]) ; win = round(win*10000)/10000; if(win == [ 0.9901; 0.9964; 0.9996; 0.9996; 0.9964; 0.9901 ]) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("kaiser Test failed") end ///////////////////////////////////////////// ///////////////////////////////////////////// /////////Test case for autoreg_matrix ////////// m = octave_fun("autoreg_matrix",[1 2 3],[2]); if(m == [1 0 0; 1 1 0; 1 2 1]) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("autoreg_matrix test failed") end ///////////////////////////////////////////// /////////Test case for arch_rnd ////////// a = [1 2 3 4 5]; b = [7 8 9 10]; t = 5 ; rand("seed", 0) m = octave_fun("arch_rnd",a,b,t); m = round(m*1000)/1000 if(m == [ 7.211 65.48 654.008 7194.657 78364.905 ]) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("arch_rnd Test failed") end ///////////////////////////////////////////// /////////Test case for postpad ////////// y = octave_fun("postpad",[1 2 3],6); if(y == [1 2 3 0 0 0 ] ) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("postpad Test failed") end /////////Test case for wconv ////////// a = [1 2 3 4 5]; b = [7 8 9 10]; y = octave_fun("wconv",[1],a,b); if(y == [7. 22. 46. 80. 114. 106. 85. 50.]) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("wconv Test failed") end /////////Test case for dctmtx ////////// T = octave_fun("dctmtx",[3]); T = round(T*10000)/10000; if(T == [0.5774 0.5774 0.5774; 0.7071 0 -0.7071; 0.4082 -0.8165 0.4082]) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("dctmtx Test failed") end /////////Test case for 40) medfilt1 ////////// fs = 100; t = 0:1/fs:1; x = sin(2*%pi*t*3)+0.25*sin(2*%pi*t*40); y = octave_fun("medfilt1","signal",x,10); y = round(y*10000)/10000 ; y = y' if(y == fscanfMat("macros/medfilt1op.txt")) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("medfilt1 Test failed") end ///////////////////////////////////////////// /////////Test case for arma_rnd ////////// a = [1 2 3 4 5]; b = [7; 8; 9; 10; 11]; t = 5 ; v = 10 ; n = 100 ; rand("seed", 0) m = octave_fun("arma_rnd","signal",a, b, v, t, n); m = round(m) ; if(m == [ 61401. 158177. 407440. 1049604. 2703841. ]) test_pass=[test_pass,1] else test_pass=[test_pass,0] disp("arma_rnd Test failed") end ///////////////////////////////////////////// res=find(test_pass==0) if(res~=[]) disp("One or more tests failed in test1") disp(length(test_pass),"Total functions tested in test1:") test1=1; end