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.283885 1.36446 0.906001 1.24942 1.55783 0.56497 -0.596145 -1.19732 -1.26756 -0.300385 -1.0746 2.36476 2.36508 1.22804 0.234136 -0.732016 -2.10553 -0.923457 -1.89227 -0.664315]
x2 = [-0.869443 -0.0504762 0.799682 1.02849 -0.900489 -0.227649 -1.83729 0.0511841 1.44132 -1.48058 1.07507 -0.396638 0.51065 1.10323 0.910165 1.12342 1.46072 -0.269572 -1.00539 1.87195]
cx = [0.349763+0.510483i 0.174148+0.475202i -0.708159-0.127081i 0.388325-0.572027i -1.62447-0.320883i 0.1836-1.02371i 0.214321-0.98353i 0.261137+1.13396i -0.225782+1.07734i -0.884797+0.465648i -0.60769+0.0533137i -0.121635-0.848777i -1.28903+0.562993i -0.0690923-0.822132i 0.953281-0.128067i 0.278311+0.598587i -0.190616-0.13041i 0.0775548-0.185926i -0.664453+0.326851i 0.600433-0.218791i]
cx2 = [-0.0820166-0.351743i 0.221056+0.634561i 0.0872567-0.29795i 0.296357+0.901056i -0.232511+0.964179i 0.836821-0.477303i -0.637571-1.24894i 0.339558+0.0234865i -0.181505-0.863075i 0.621065-0.359776i 0.61333+0.506973i -0.295419+0.675185i -0.475898+0.120633i -0.85725-0.0647025i 0.564002+0.646938i 0.350946-0.839171i 0.530681-0.233433i -0.376251+0.957435i 0.959656+0.279574i -0.165729-2.07819i]
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.312273 1.21702 2.07718 3.36094 4.94961 5.4153 4.70304 3.4456 2.17101 1.96735 1.09921 2.44344 3.90255 3.76747 2.34438 0.950782 -0.695958 -0.303889 -1.02548 -1.26661]
s1 = [-1.0746 -0.664315 -1.89227 -0.923457 -2.10553 -0.732016 0.234136 1.22804 2.36508 2.36476]
y2 = [-1.08197 -3.41531 -4.89568 -5.07236 -6.39988 -5.82823 -5.72095 -4.55746 -1.08485 -2.19331 0.00677576 -0.486557 -0.684861 -0.550864 1.24048 2.61288 5.94461 5.45083 2.93053 6.5708]
s2 = [1.07507 1.87195 -1.00539 -0.269572 1.46072 1.12342 0.910165 1.10323 0.51065 -0.396638]
Redo, reseting to state s1
y2 = [-1.08197 -3.41531 -4.89568 -5.07236 -6.39988 -5.82823 -5.72095 -4.55746 -1.08485 -2.19331 0.00677576 -0.486557 -0.684861 -0.550864 1.24048 2.61288 5.94461 5.45083 2.93053 6.5708]
s2 = [1.07507 1.87195 -1.00539 -0.269572 1.46072 1.12342 0.910165 1.10323 0.51065 -0.396638]
cy = [0.384739+0.561531i 0.541325+1.03321i -0.255065+0.845896i 0.242908+0.229375i -1.58285-0.0663932i -1.21844-1.16039i -1.00105-2.1399i -0.735228-0.794189i -1.0097+0.277485i -1.9604+0.681964i -2.89014+0.183561i -3.13732-1.23063i -3.83493-0.399376i -4.17035-0.787994i -1.49036-0.525772i -1.46314+1.16919i -1.91497+1.94941i -2.07174+0.623974i -2.58461-0.075235i -0.972891-0.814238i]
cs1 = [-0.60769+0.0533137i 0.600433-0.218791i -0.664453+0.326851i 0.0775548-0.185926i -0.190616-0.13041i 0.278311+0.598587i 0.953281-0.128067i -0.0690923-0.822132i -1.28903+0.562993i -0.121635-0.848777i]
cy2 = [-0.515462-1.23259i -0.142464+0.34938i 1.22044-0.604815i 1.6068+1.23827i 0.368122+2.33683i 1.03357+1.11679i 0.439171-0.0788998i 0.798888+0.257756i 1.22973-1.02083i 1.33062-1.11148i 2.02519-0.166091i 1.41784-0.108646i 0.83664+0.254482i -0.355103-0.72981i 0.583535-0.975887i 0.0763542-1.48636i 1.26258-0.410285i 0.456077+0.64275i 1.73083+1.71761i 0.831495-0.236578i]
cs2 = [0.61333+0.506973i -0.165729-2.07819i 0.959656+0.279574i -0.376251+0.957435i 0.530681-0.233433i 0.350946-0.839171i 0.564002+0.646938i -0.85725-0.0647025i -0.475898+0.120633i -0.295419+0.675185i]
Redo, reseting to state s1
cy2 = [-0.515462-1.23259i -0.142464+0.34938i 1.22044-0.604815i 1.6068+1.23827i 0.368122+2.33683i 1.03357+1.11679i 0.439171-0.0788998i 0.798888+0.257756i 1.22973-1.02083i 1.33062-1.11148i 2.02519-0.166091i 1.41784-0.108646i 0.83664+0.254482i -0.355103-0.72981i 0.583535-0.975887i 0.0763542-1.48636i 1.26258-0.410285i 0.456077+0.64275i 1.73083+1.71761i 0.831495-0.236578i]
cs2 = [0.61333+0.506973i -0.165729-2.07819i 0.959656+0.279574i -0.376251+0.957435i 0.530681-0.233433i 0.350946-0.839171i 0.564002+0.646938i -0.85725-0.0647025i -0.475898+0.120633i -0.295419+0.675185i]
cy = [0.384739+0.561531i 0.541325+1.03321i -0.255065+0.845896i 0.242908+0.229375i -1.58285-0.0663932i -1.21844-1.16039i -1.00105-2.1399i -0.735228-0.794189i -1.0097+0.277485i -1.9604+0.681964i -2.89014+0.183561i -3.13732-1.23063i -3.83493-0.399376i -4.17035-0.787994i -1.49036-0.525772i -1.46314+1.16919i -1.91497+1.94941i -2.07174+0.623974i -2.58461-0.075235i -0.972891-0.814238i]
y = [-0.312273 1.21702 2.07718 3.36094 4.94961 5.4153 4.70304 3.4456 2.17101 1.96735 1.09921 2.44344 3.90255 3.76747 2.34438 0.950782 -0.695958 -0.303889 -1.02548 -1.26661]
cy = [0.384739+0.561531i 0.541325+1.03321i -0.255065+0.845896i 0.242908+0.229375i -1.58285-0.0663932i -1.21844-1.16039i -1.00105-2.1399i -0.735228-0.794189i -1.0097+0.277485i -1.9604+0.681964i -2.89014+0.183561i -3.13732-1.23063i -3.83493-0.399376i -4.17035-0.787994i -1.49036-0.525772i -1.46314+1.16919i -1.91497+1.94941i -2.07174+0.623974i -2.58461-0.075235i -0.972891-0.814238i]
cy = [0.384739+0.561531i 0.541325+1.03321i -0.255065+0.845896i 0.242908+0.229375i -1.58285-0.0663932i -1.21844-1.16039i -1.00105-2.1399i -0.735228-0.794189i -1.0097+0.277485i -1.9604+0.681964i -2.89014+0.183561i -3.13732-1.23063i -3.83493-0.399376i -4.17035-0.787994i -1.49036-0.525772i -1.46314+1.16919i -1.91497+1.94941i -2.07174+0.623974i -2.58461-0.075235i -0.972891-0.814238i]
-----------------------------------------------
AR Filter:
y = [-0.283885 1.08058 1.98658 3.236 4.79383 5.3588 4.76265 3.56534 2.29777 1.99739 0.922784 3.28754 5.65262 6.88066 7.1148 6.38278 4.27725 3.35379 1.46152 0.797203]
s1 = [0.797203]
y2 = [-0.0722399 -0.122716 0.676966 1.70545 0.804964 0.577315 -1.25997 -1.20879 0.232537 -1.24804 -0.172971 -0.569609 -0.0589596 1.04427 1.95443 3.07785 4.53857 4.269 3.26361 5.13556]
s2 = [5.13556]
Redo, reseting to state s1
y2 = [-0.0722399 -0.122716 0.676966 1.70545 0.804964 0.577315 -1.25997 -1.20879 0.232537 -1.24804 -0.172971 -0.569609 -0.0589596 1.04427 1.95443 3.07785 4.53857 4.269 3.26361 5.13556]
s2 = [5.13556]
cy = [0.349763+0.510483i 0.52391+0.985685i -0.184249+0.858604i 0.204076+0.286578i -1.4204-0.0343049i -1.2368-1.05802i -1.02248-2.04155i -0.761341-0.907585i -0.987124+0.169751i -1.87192+0.635399i -2.47961+0.688713i -2.60125-0.160064i -3.89027+0.402929i -3.95937-0.419203i -3.00609-0.54727i -2.72777+0.0513172i -2.91839-0.0790928i -2.84084-0.265019i -3.50529+0.0618313i -2.90486-0.156959i]
cs1 = [-2.90486-0.156959i]
cy2 = [-2.98687-0.508702i -2.76582+0.125859i -2.67856-0.172091i -2.3822+0.728965i -2.61471+1.69314i -1.77789+1.21584i -2.41546-0.0330987i -2.0759-0.00961226i -2.25741-0.872688i -1.63634-1.23246i -1.02301-0.72549i -1.31843-0.0503057i -1.79433+0.0703278i -2.65158+0.0056253i -2.08758+0.652563i -1.73663-0.186607i -1.20595-0.42004i -1.5822+0.537395i -0.622546+0.816968i -0.788275-1.26122i]
cs2 = [-0.788275-1.26122i]
Redo, reseting to state s1
cy2 = [-2.98687-0.508702i -2.76582+0.125859i -2.67856-0.172091i -2.3822+0.728965i -2.61471+1.69314i -1.77789+1.21584i -2.41546-0.0330987i -2.0759-0.00961226i -2.25741-0.872688i -1.63634-1.23246i -1.02301-0.72549i -1.31843-0.0503057i -1.79433+0.0703278i -2.65158+0.0056253i -2.08758+0.652563i -1.73663-0.186607i -1.20595-0.42004i -1.5822+0.537395i -0.622546+0.816968i -0.788275-1.26122i]
cs2 = [-0.788275-1.26122i]
cy = [0.349763+0.510483i 0.52391+0.985685i -0.184249+0.858604i 0.204076+0.286578i -1.4204-0.0343049i -1.2368-1.05802i -1.02248-2.04155i -0.761341-0.907585i -0.987124+0.169751i -1.87192+0.635399i -2.47961+0.688713i -2.60125-0.160064i -3.89027+0.402929i -3.95937-0.419203i -3.00609-0.54727i -2.72777+0.0513172i -2.91839-0.0790928i -2.84084-0.265019i -3.50529+0.0618313i -2.90486-0.156959i]
y = [-0.283885 1.08058 1.98658 3.236 4.79383 5.3588 4.76265 3.56534 2.29777 1.99739 0.922784 3.28754 5.65262 6.88066 7.1148 6.38278 4.27725 3.35379 1.46152 0.797203]
cy = [0.349763+0.510483i 0.52391+0.985685i -0.184249+0.858604i 0.204076+0.286578i -1.4204-0.0343049i -1.2368-1.05802i -1.02248-2.04155i -0.761341-0.907585i -0.987124+0.169751i -1.87192+0.635399i -2.47961+0.688713i -2.60125-0.160064i -3.89027+0.402929i -3.95937-0.419203i -3.00609-0.54727i -2.72777+0.0513172i -2.91839-0.0790928i -2.84084-0.265019i -3.50529+0.0618313i -2.90486-0.156959i]
cy = [0.349763+0.510483i 0.52391+0.985685i -0.184249+0.858604i 0.204076+0.286578i -1.4204-0.0343049i -1.2368-1.05802i -1.02248-2.04155i -0.761341-0.907585i -0.987124+0.169751i -1.87192+0.635399i -2.47961+0.688713i -2.60125-0.160064i -3.89027+0.402929i -3.95937-0.419203i -3.00609-0.54727i -2.72777+0.0513172i -2.91839-0.0790928i -2.84084-0.265019i -3.50529+0.0618313i -2.90486-0.156959i]
-----------------------------------------------
ARMA Filter:
y = [-0.312273 0.904751 2.98193 6.34287 11.2925 16.7078 21.4108 24.8564 27.0274 28.9948 30.094 32.5374 36.44 40.2074 42.5518 43.5026 42.8067 42.5028 41.4773 40.2107]
s1 = [0.797203 1.46152 3.35379 4.27725 6.38278 7.1148 6.88066 5.65262 3.28754]
y2 = [39.1287 35.7134 30.8177 25.7453 19.3455 13.5172 7.79629 3.23883 2.15398 -0.039326 -0.0325503 -0.519107 -1.20397 -1.75483 -0.514348 2.09853 8.04315 13.494 16.4245 22.9953]
s2 = [5.13556 3.26361 4.269 4.53857 3.07785 1.95443 1.04427 -0.0589596 -0.569609]
Redo, reseting to state s1
y2 = [39.1287 35.7134 30.8177 25.7453 19.3455 13.5172 7.79629 3.23883 2.15398 -0.039326 -0.0325503 -0.519107 -1.20397 -1.75483 -0.514348 2.09853 8.04315 13.494 16.4245 22.9953]
s2 = [5.13556 3.26361 4.269 4.53857 3.07785 1.95443 1.04427 -0.0589596 -0.569609]
cy = [0.384739+0.561531i 0.926064+1.59474i 0.670999+2.44063i 0.913908+2.67001i -0.668938+2.60362i -1.88738+1.44323i -2.88842-0.696676i -3.62365-1.49086i -4.63335-1.21338i -6.59375-0.531416i -9.4839-0.347854i -12.6212-1.57848i -16.4561-1.97786i -20.6265-2.76585i -22.1169-3.29162i -23.58-2.12243i -25.495-0.173016i -27.5667+0.450958i -30.1513+0.375723i -31.1242-0.438515i]
cs1 = [-2.90486-0.156959i -3.50529+0.0618313i -2.84084-0.265019i -2.91839-0.0790928i -2.72777+0.0513172i -3.00609-0.54727i -3.95937-0.419203i -3.89027+0.402929i -2.60125-0.160064i]
cy2 = [-31.6397-1.6711i -31.7821-1.32173i -30.5617-1.92654i -28.9549-0.688266i -28.5868+1.64857i -27.5532+2.76536i -27.114+2.68646i -26.3152+2.94421i -25.0854+1.92339i -23.7548+0.811907i -21.7296+0.645816i -20.3118+0.53717i -19.4751+0.791652i -19.8302+0.0618418i -19.2467-0.914045i -19.1703-2.40041i -17.9078-2.81069i -17.4517-2.16794i -15.7209-0.450331i -14.8894-0.686909i]
cs2 = [-0.788275-1.26122i -0.622546+0.816968i -1.5822+0.537395i -1.20595-0.42004i -1.73663-0.186607i -2.08758+0.652563i -2.65158+0.0056253i -1.79433+0.0703278i -1.31843-0.0503057i]
Redo, reseting to state s1
cy2 = [-31.6397-1.6711i -31.7821-1.32173i -30.5617-1.92654i -28.9549-0.688266i -28.5868+1.64857i -27.5532+2.76536i -27.114+2.68646i -26.3152+2.94421i -25.0854+1.92339i -23.7548+0.811907i -21.7296+0.645816i -20.3118+0.53717i -19.4751+0.791652i -19.8302+0.0618418i -19.2467-0.914045i -19.1703-2.40041i -17.9078-2.81069i -17.4517-2.16794i -15.7209-0.450331i -14.8894-0.686909i]
cs2 = [-0.788275-1.26122i -0.622546+0.816968i -1.5822+0.537395i -1.20595-0.42004i -1.73663-0.186607i -2.08758+0.652563i -2.65158+0.0056253i -1.79433+0.0703278i -1.31843-0.0503057i]
cy = [0.384739+0.561531i 0.926064+1.59474i 0.670999+2.44063i 0.913908+2.67001i -0.668938+2.60362i -1.88738+1.44323i -2.88842-0.696676i -3.62365-1.49086i -4.63335-1.21338i -6.59375-0.531416i -9.4839-0.347854i -12.6212-1.57848i -16.4561-1.97786i -20.6265-2.76585i -22.1169-3.29162i -23.58-2.12243i -25.495-0.173016i -27.5667+0.450958i -30.1513+0.375723i -31.1242-0.438515i]
y = [-0.312273 0.904751 2.98193 6.34287 11.2925 16.7078 21.4108 24.8564 27.0274 28.9948 30.094 32.5374 36.44 40.2074 42.5518 43.5026 42.8067 42.5028 41.4773 40.2107]
cy = [0.384739+0.561531i 0.926064+1.59474i 0.670999+2.44063i 0.913908+2.67001i -0.668938+2.60362i -1.88738+1.44323i -2.88842-0.696676i -3.62365-1.49086i -4.63335-1.21338i -6.59375-0.531416i -9.4839-0.347854i -12.6212-1.57848i -16.4561-1.97786i -20.6265-2.76585i -22.1169-3.29162i -23.58-2.12243i -25.495-0.173016i -27.5667+0.450958i -30.1513+0.375723i -31.1242-0.438515i]
cy = [0.384739+0.561531i 0.926064+1.59474i 0.670999+2.44063i 0.913908+2.67001i -0.668938+2.60362i -1.88738+1.44323i -2.88842-0.696676i -3.62365-1.49086i -4.63335-1.21338i -6.59375-0.531416i -9.4839-0.347854i -12.6212-1.57848i -16.4561-1.97786i -20.6265-2.76585i -22.1169-3.29162i -23.58-2.12243i -25.495-0.173016i -27.5667+0.450958i -30.1513+0.375723i -31.1242-0.438515i]
|
