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
|
======================================
Test of filter routines
======================================
Input signals:
x = [0.764873 -1.13073 0.452336 -0.451817 1.36468 0.355365 -0.556704 1.00977 0.251912 -0.902459 -0.946328 0.403829 -0.129499 0.28002 -0.439783 -0.0494359 -0.235643 -0.268905 -1.53793 1.67876]
x2 = [1.02093 1.25461 1.00032 -1.12474 0.369076 0.611692 -0.0799272 -0.101654 1.02917 -1.37943 -0.534532 1.51576 -0.116042 1.18937 -0.105928 0.103149 -0.86184 -0.420697 0.297823 -0.286478]
cx = [0.104277-0.829977i 0.0806219+0.369791i -0.393648-0.117384i 0.314781+0.255673i -0.524515+0.343174i 0.178492+0.704202i -0.0609598-0.284082i 0.668861-0.108213i 0.131127+0.774733i -0.396762+0.109828i 0.743216-0.443451i -1.16267+0.999606i -0.725603+0.0559068i 0.214591-0.161286i 0.557317-0.667273i -0.802313+0.783623i -0.87016+1.1201i 0.0713311-0.131199i 1.22998+1.08696i 0.310264+0.650456i]
cx2 = [-0.145369-0.0309792i -0.433019+1.0476i 0.653553+0.916184i -0.210851-0.0526413i -0.0937395+0.192222i -0.209132+0.394894i 0.211864-0.918654i -0.625162-0.943399i 0.429391-0.538272i 0.561267-0.750045i 0.219214-0.141661i 0.100501+1.05596i 0.126522-0.348387i -0.2209+0.166731i -0.540764+0.476032i 0.176139-0.276432i 1.01785-0.561079i -0.23416+0.548687i 0.319975+0.714349i 0.045842-1.01049i]
Filter coefficients:
b = [1.1 1 1 1 1 1 1 1 1 1]
cb = [1.1+0i 1+0i 1+0i 1+0i 1+0i 1+0i 1+0i 1+0i 1+0i 1+0i]
a = [1 -1]
ca = [1+0i -1+0i]
-----------------------------------------------
MA Filter:
y = [0.841361 -0.478931 0.131712 -0.410521 1.13581 1.39024 0.742329 1.90875 2.08488 1.06698 -0.648609 1.02097 0.385799 1.15859 -0.717851 -1.08362 -0.781177 -2.06318 -3.97993 -1.07704]
s1 = [-0.946328 1.67876 -1.53793 -0.268905 -0.235643 -0.0494359 -0.439783 0.28002 -0.129499 0.403829]
y2 = [0.824441 1.69859 2.80297 1.18571 2.14395 2.82934 2.91589 3.08097 5.76116 2.46211 0.991133 1.45732 0.177777 2.62243 2.0179 1.53026 0.651848 0.376919 -0.282576 0.751943]
s2 = [-0.534532 -0.286478 0.297823 -0.420697 -0.86184 0.103149 -0.105928 1.18937 -0.116042 1.51576]
Redo, reseting to state s1
y2 = [0.824441 1.69859 2.80297 1.18571 2.14395 2.82934 2.91589 3.08097 5.76116 2.46211 0.991133 1.45732 0.177777 2.62243 2.0179 1.53026 0.651848 0.376919 -0.282576 0.751943]
s2 = [-0.534532 -0.286478 0.297823 -0.420697 -0.86184 0.103149 -0.105928 1.18937 -0.116042 1.51576]
cy = [0.114705-0.912975i 0.192961-0.423207i -0.248113-0.589309i 0.137511-0.296329i -0.470934+0.0555952i -0.222142+0.7959i -0.307047+0.41299i 0.434796+0.322364i 0.51215+1.18539i 0.0625985+1.22873i 0.815535+1.55993i -0.618341+2.33405i -0.906591+2.41297i -0.912761+1.97429i 0.203343+0.913244i -0.913425+1.13776i -1.72941+2.57559i -2.23279+2.42747i -1.01807+2.86151i -0.403018+3.35849i]
cs1 = [0.743216-0.443451i 0.310264+0.650456i 1.22998+1.08696i 0.0713311-0.131199i -0.87016+1.1201i -0.802313+0.783623i 0.557317-0.667273i 0.214591-0.161286i -0.725603+0.0559068i -1.16267+0.999606i]
cy2 = [-1.33717+3.70282i -0.636284+3.85867i 0.851529+4.7058i 0.339647+4.71757i -0.299698+5.60155i 0.281944+5.23309i 1.40607+3.06298i 0.625872+2.2483i -0.0692617+0.66358i 0.194929-0.758098i 0.525306-0.807942i 1.04696-0.679811i 0.522527-2.08482i 0.477735-1.81393i -0.0012752-1.49919i 0.455686-2.24576i 1.34584-1.91665i 1.61164-0.313592i 1.55764+0.955595i 1.0148+0.52267i]
cs2 = [0.219214-0.141661i 0.045842-1.01049i 0.319975+0.714349i -0.23416+0.548687i 1.01785-0.561079i 0.176139-0.276432i -0.540764+0.476032i -0.2209+0.166731i 0.126522-0.348387i 0.100501+1.05596i]
Redo, reseting to state s1
cy2 = [-1.33717+3.70282i -0.636284+3.85867i 0.851529+4.7058i 0.339647+4.71757i -0.299698+5.60155i 0.281944+5.23309i 1.40607+3.06298i 0.625872+2.2483i -0.0692617+0.66358i 0.194929-0.758098i 0.525306-0.807942i 1.04696-0.679811i 0.522527-2.08482i 0.477735-1.81393i -0.0012752-1.49919i 0.455686-2.24576i 1.34584-1.91665i 1.61164-0.313592i 1.55764+0.955595i 1.0148+0.52267i]
cs2 = [0.219214-0.141661i 0.045842-1.01049i 0.319975+0.714349i -0.23416+0.548687i 1.01785-0.561079i 0.176139-0.276432i -0.540764+0.476032i -0.2209+0.166731i 0.126522-0.348387i 0.100501+1.05596i]
cy = [0.114705-0.912975i 0.192961-0.423207i -0.248113-0.589309i 0.137511-0.296329i -0.470934+0.0555952i -0.222142+0.7959i -0.307047+0.41299i 0.434796+0.322364i 0.51215+1.18539i 0.0625985+1.22873i 0.815535+1.55993i -0.618341+2.33405i -0.906591+2.41297i -0.912761+1.97429i 0.203343+0.913244i -0.913425+1.13776i -1.72941+2.57559i -2.23279+2.42747i -1.01807+2.86151i -0.403018+3.35849i]
y = [0.841361 -0.478931 0.131712 -0.410521 1.13581 1.39024 0.742329 1.90875 2.08488 1.06698 -0.648609 1.02097 0.385799 1.15859 -0.717851 -1.08362 -0.781177 -2.06318 -3.97993 -1.07704]
cy = [0.114705-0.912975i 0.192961-0.423207i -0.248113-0.589309i 0.137511-0.296329i -0.470934+0.0555952i -0.222142+0.7959i -0.307047+0.41299i 0.434796+0.322364i 0.51215+1.18539i 0.0625985+1.22873i 0.815535+1.55993i -0.618341+2.33405i -0.906591+2.41297i -0.912761+1.97429i 0.203343+0.913244i -0.913425+1.13776i -1.72941+2.57559i -2.23279+2.42747i -1.01807+2.86151i -0.403018+3.35849i]
cy = [0.114705-0.912975i 0.192961-0.423207i -0.248113-0.589309i 0.137511-0.296329i -0.470934+0.0555952i -0.222142+0.7959i -0.307047+0.41299i 0.434796+0.322364i 0.51215+1.18539i 0.0625985+1.22873i 0.815535+1.55993i -0.618341+2.33405i -0.906591+2.41297i -0.912761+1.97429i 0.203343+0.913244i -0.913425+1.13776i -1.72941+2.57559i -2.23279+2.42747i -1.01807+2.86151i -0.403018+3.35849i]
-----------------------------------------------
AR Filter:
y = [0.764873 -0.365858 0.0864783 -0.365339 0.999338 1.3547 0.797999 1.80777 2.05968 1.15723 0.210897 0.614726 0.485227 0.765248 0.325465 0.276029 0.0403861 -0.228519 -1.76645 -0.0876894]
s1 = [-0.0876894]
y2 = [0.933245 2.18785 3.18817 2.06343 2.43251 3.0442 2.96427 2.86262 3.89179 2.51236 1.97783 3.49359 3.37755 4.56692 4.461 4.56415 3.70231 3.28161 3.57943 3.29295]
s2 = [3.29295]
Redo, reseting to state s1
y2 = [0.933245 2.18785 3.18817 2.06343 2.43251 3.0442 2.96427 2.86262 3.89179 2.51236 1.97783 3.49359 3.37755 4.56692 4.461 4.56415 3.70231 3.28161 3.57943 3.29295]
s2 = [3.29295]
cy = [0.104277-0.829977i 0.184899-0.460186i -0.208749-0.57757i 0.106032-0.321897i -0.418483+0.0212777i -0.239991+0.72548i -0.300951+0.441398i 0.36791+0.333185i 0.499037+1.10792i 0.102275+1.21775i 0.845491+0.774295i -0.317176+1.7739i -1.04278+1.82981i -0.828188+1.66852i -0.270872+1.00125i -1.07318+1.78487i -1.94334+2.90497i -1.87201+2.77377i -0.642034+3.86074i -0.33177+4.51119i]
cs1 = [-0.33177+4.51119i]
cy2 = [-0.477139+4.48021i -0.910158+5.52781i -0.256605+6.44399i -0.467456+6.39135i -0.561195+6.58358i -0.770327+6.97847i -0.558463+6.05982i -1.18363+5.11642i -0.754234+4.57814i -0.192968+3.8281i 0.0262459+3.68644i 0.126747+4.7424i 0.253269+4.39402i 0.0323694+4.56075i -0.508394+5.03678i -0.332255+4.76035i 0.685592+4.19927i 0.451432+4.74796i 0.771407+5.4623i 0.817249+4.45182i]
cs2 = [0.817249+4.45182i]
Redo, reseting to state s1
cy2 = [-0.477139+4.48021i -0.910158+5.52781i -0.256605+6.44399i -0.467456+6.39135i -0.561195+6.58358i -0.770327+6.97847i -0.558463+6.05982i -1.18363+5.11642i -0.754234+4.57814i -0.192968+3.8281i 0.0262459+3.68644i 0.126747+4.7424i 0.253269+4.39402i 0.0323694+4.56075i -0.508394+5.03678i -0.332255+4.76035i 0.685592+4.19927i 0.451432+4.74796i 0.771407+5.4623i 0.817249+4.45182i]
cs2 = [0.817249+4.45182i]
cy = [0.104277-0.829977i 0.184899-0.460186i -0.208749-0.57757i 0.106032-0.321897i -0.418483+0.0212777i -0.239991+0.72548i -0.300951+0.441398i 0.36791+0.333185i 0.499037+1.10792i 0.102275+1.21775i 0.845491+0.774295i -0.317176+1.7739i -1.04278+1.82981i -0.828188+1.66852i -0.270872+1.00125i -1.07318+1.78487i -1.94334+2.90497i -1.87201+2.77377i -0.642034+3.86074i -0.33177+4.51119i]
y = [0.764873 -0.365858 0.0864783 -0.365339 0.999338 1.3547 0.797999 1.80777 2.05968 1.15723 0.210897 0.614726 0.485227 0.765248 0.325465 0.276029 0.0403861 -0.228519 -1.76645 -0.0876894]
cy = [0.104277-0.829977i 0.184899-0.460186i -0.208749-0.57757i 0.106032-0.321897i -0.418483+0.0212777i -0.239991+0.72548i -0.300951+0.441398i 0.36791+0.333185i 0.499037+1.10792i 0.102275+1.21775i 0.845491+0.774295i -0.317176+1.7739i -1.04278+1.82981i -0.828188+1.66852i -0.270872+1.00125i -1.07318+1.78487i -1.94334+2.90497i -1.87201+2.77377i -0.642034+3.86074i -0.33177+4.51119i]
cy = [0.104277-0.829977i 0.184899-0.460186i -0.208749-0.57757i 0.106032-0.321897i -0.418483+0.0212777i -0.239991+0.72548i -0.300951+0.441398i 0.36791+0.333185i 0.499037+1.10792i 0.102275+1.21775i 0.845491+0.774295i -0.317176+1.7739i -1.04278+1.82981i -0.828188+1.66852i -0.270872+1.00125i -1.07318+1.78487i -1.94334+2.90497i -1.87201+2.77377i -0.642034+3.86074i -0.33177+4.51119i]
-----------------------------------------------
ARMA Filter:
y = [0.841361 0.36243 0.494142 0.083621 1.21943 2.60967 3.35199 5.26074 7.34562 8.4126 7.76399 8.78496 9.17075 10.3293 9.61149 8.52787 7.7467 5.68352 1.70359 0.626548]
s1 = [-0.0876894 -1.76645 -0.228519 0.0403861 0.276029 0.325465 0.765248 0.485227 0.614726]
y2 = [1.45099 3.14958 5.95255 7.13826 9.28221 12.1116 15.0274 18.1084 23.8696 26.3317 27.3228 28.7801 28.9579 31.5803 33.5982 35.1285 35.7803 36.1573 35.8747 36.6266]
s2 = [3.29295 3.57943 3.28161 3.70231 4.56415 4.461 4.56692 3.37755 3.49359]
Redo, reseting to state s1
y2 = [1.45099 3.14958 5.95255 7.13826 9.28221 12.1116 15.0274 18.1084 23.8696 26.3317 27.3228 28.7801 28.9579 31.5803 33.5982 35.1285 35.7803 36.1573 35.8747 36.6266]
s2 = [3.29295 3.57943 3.28161 3.70231 4.56415 4.461 4.56692 3.37755 3.49359]
cy = [0.114705-0.912975i 0.307666-1.33618i 0.0595531-1.92549i 0.197064-2.22182i -0.27387-2.16622i -0.496012-1.37032i -0.803059-0.957335i -0.368263-0.634971i 0.143887+0.55042i 0.206485+1.77915i 1.02202+3.33908i 0.403679+5.67312i -0.502912+8.08609i -1.41567+10.0604i -1.21233+10.9736i -2.12576+12.1114i -3.85517+14.687i -6.08796+17.1144i -7.10603+19.976i -7.50905+23.3344i]
cs1 = [-0.33177+4.51119i -0.642034+3.86074i -1.87201+2.77377i -1.94334+2.90497i -1.07318+1.78487i -0.270872+1.00125i -0.828188+1.66852i -1.04278+1.82981i -0.317176+1.7739i]
cy2 = [-8.84621+27.0373i -9.4825+30.8959i -8.63097+35.6017i -8.29132+40.3193i -8.59102+45.9209i -8.30908+51.1539i -6.90301+54.2169i -6.27714+56.4652i -6.3464+57.1288i -6.15147+56.3707i -5.62616+55.5628i -4.57921+54.883i -4.05668+52.7981i -3.57894+50.9842i -3.58022+49.485i -3.12453+47.2392i -1.77869+45.3226i -0.167052+45.009i 1.39059+45.9646i 2.40539+46.4873i]
cs2 = [0.817249+4.45182i 0.771407+5.4623i 0.451432+4.74796i 0.685592+4.19927i -0.332255+4.76035i -0.508394+5.03678i 0.0323694+4.56075i 0.253269+4.39402i 0.126747+4.7424i]
Redo, reseting to state s1
cy2 = [-8.84621+27.0373i -9.4825+30.8959i -8.63097+35.6017i -8.29132+40.3193i -8.59102+45.9209i -8.30908+51.1539i -6.90301+54.2169i -6.27714+56.4652i -6.3464+57.1288i -6.15147+56.3707i -5.62616+55.5628i -4.57921+54.883i -4.05668+52.7981i -3.57894+50.9842i -3.58022+49.485i -3.12453+47.2392i -1.77869+45.3226i -0.167052+45.009i 1.39059+45.9646i 2.40539+46.4873i]
cs2 = [0.817249+4.45182i 0.771407+5.4623i 0.451432+4.74796i 0.685592+4.19927i -0.332255+4.76035i -0.508394+5.03678i 0.0323694+4.56075i 0.253269+4.39402i 0.126747+4.7424i]
cy = [0.114705-0.912975i 0.307666-1.33618i 0.0595531-1.92549i 0.197064-2.22182i -0.27387-2.16622i -0.496012-1.37032i -0.803059-0.957335i -0.368263-0.634971i 0.143887+0.55042i 0.206485+1.77915i 1.02202+3.33908i 0.403679+5.67312i -0.502912+8.08609i -1.41567+10.0604i -1.21233+10.9736i -2.12576+12.1114i -3.85517+14.687i -6.08796+17.1144i -7.10603+19.976i -7.50905+23.3344i]
y = [0.841361 0.36243 0.494142 0.083621 1.21943 2.60967 3.35199 5.26074 7.34562 8.4126 7.76399 8.78496 9.17075 10.3293 9.61149 8.52787 7.7467 5.68352 1.70359 0.626548]
cy = [0.114705-0.912975i 0.307666-1.33618i 0.0595531-1.92549i 0.197064-2.22182i -0.27387-2.16622i -0.496012-1.37032i -0.803059-0.957335i -0.368263-0.634971i 0.143887+0.55042i 0.206485+1.77915i 1.02202+3.33908i 0.403679+5.67312i -0.502912+8.08609i -1.41567+10.0604i -1.21233+10.9736i -2.12576+12.1114i -3.85517+14.687i -6.08796+17.1144i -7.10603+19.976i -7.50905+23.3344i]
cy = [0.114705-0.912975i 0.307666-1.33618i 0.0595531-1.92549i 0.197064-2.22182i -0.27387-2.16622i -0.496012-1.37032i -0.803059-0.957335i -0.368263-0.634971i 0.143887+0.55042i 0.206485+1.77915i 1.02202+3.33908i 0.403679+5.67312i -0.502912+8.08609i -1.41567+10.0604i -1.21233+10.9736i -2.12576+12.1114i -3.85517+14.687i -6.08796+17.1144i -7.10603+19.976i -7.50905+23.3344i]
|
