1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
|
clc
disp("Testing octaves built in function hamming")
y = octave_fun("hamming",5)
pause
disp("Testing octaves built in function hamming with optional arguments")
y = octave_fun("hamming",5,"periodic")
pause
disp("Testing arburg function from signal package")
a = octave_fun("arburg","signal", [1,2,3,4,5],2)
pause
disp("Testing aryule function from signal package")
a = octave_fun("aryule","signal", [1,2,3,4,5],2)
pause
disp("Testing filter, xcorr and levinson function from signal package")
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)
pause
disp("Testing sigmoid_train function from signal package")
s = octave_fun("sigmoid_train",0.1,[1:3],[4])
pause
disp("Testing circshift function from signal package")
M = [1 2 3 4];
R = octave_fun("circshift",M,[0 1])
pause
disp("Testing kaiser function from signal package")
win = octave_fun("kaiser",6,[0.2])
pause
disp("Testing autoreg_matrix function from signal package")
m = octave_fun("autoreg_matrix",[1 2 3],[2])
pause
disp("Testing arch_rnd function from signal package")
a = [1 2 3 4 5];
b = [7 8 9 10];
t = 5 ;
rand("seed", 0)
m = octave_fun("arch_rnd",a,b,t)
pause
disp("Testing postpad function from signal package")
y = octave_fun("postpad",[1 2 3],6)
pause
disp("Testing wconv function from signal package")
a = [1 2 3 4 5];
b = [7 8 9 10];
y = octave_fun("wconv",[1],a,b)
pause
disp("Testing dctmtx function from signal package")
T = octave_fun("dctmtx",[3])
pause
disp("Testing medfilt1 function from signal package")
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);
pause
disp("Testing arma_rnd function from signal package")
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)
pause
disp("Testing ncauer function from signal package")
[z p g] = octave_fun("ncauer","signal",3, 10, 4)
pause
disp("Testing ellipap function from signal package")
[z p g] = octave_fun("ellipap","signal",4, 3, 10)
pause
disp("Testing besselap function from signal package")
[z p g] = octave_fun("besselap","signal",5)
pause
disp("Testing zp2tf function from signal package")
[num, den] = octave_fun("zp2tf","signal",[1 2 3], [4 5 6], 5)
pause
disp("Testing tf2zp function from signal package")
[z p k] = octave_fun("tf2zp","signal",[1 2 3], [4 5 6]);
pause
disp("Testing ctrb function from control package")
y = octave_fun("ctrb","control",[1, 0; 0, -0.5], [8; 8]);
pause
|