(kf,kq,N,x)
| 118 | // '----(-1)------------------------------' |
| 119 | |
| 120 | function chamberlinResponse (kf,kq,N,x) { |
| 121 | if (!N) { N = 512; } |
| 122 | |
| 123 | var output = []; |
| 124 | var lp = 0, bp = 0, input = 1; |
| 125 | |
| 126 | for (var i = 0; i < N; i++) { |
| 127 | bp += kf * (input - lp - kq*bp); |
| 128 | lp += kf * bp; |
| 129 | output[i] = lp; |
| 130 | input = x; |
| 131 | } |
| 132 | |
| 133 | return output; |
| 134 | } |
| 135 | |
| 136 | function chamberlinImpulseResponse (kf,kq,N) { |
| 137 | return chamberlinResponse(kf,kq,N,0); |
no outgoing calls
no test coverage detected