Stem Style

2ed1ba03-71ea-488b-87f7-3a090c7c6b0e 18 6 5 auto

&[DATE] &[TIME] - &[FILENAME]

Pre-define Unit step x Φ 1 x sign + 2 / (:

Unit pulse:width/locate x p Φ 1 x p - sign + 2 / (:

α 1 :

β 2 :

t f t α - Φ t β - Φ - :

x f

npts 16326412825651 mat :

Z 3 :

j 4 :

L H N xd start/end @ L..H U 0 : dx H L - N / : i 1 N 1 + range U i el L dx - (dx i * + : for U 21 line 41 line :

L Z - :

H Z :

N npts j el :

U L H N xd :

i 1 U rows range Y i el U i el f : for

p U Y augment :

N p rows :

y p 2 col :

Stem Style

Core component of the stem style P v replicate Mathcad 8 Pro. Treasury k 1 : i 1 v length range P i el 1 ≡ continue j 1 P i el range T k j + el v i el : for k k P i el + : 21 line if for T 41 line :

iii 1 N range P iii el 3 : P N el 2 : 21 line for

ii 1 N range v ii el ii : for

u P v replicate :

i 1 N range S 3 i * el y i el : for

V u S augment :

V V 3 V rows 12 submatrix :

D Smath does not evaluate complex ≡

λ Y length :

Y length 129

Y Y sys2mat :

D i 1 - sqrt : n 1 N range D n el 1 λ 1 - sqrt / Y k el e i 2 π * n * λ / (* k *^* (k 1 N sum * : for D 31 line :

New t 1 N 2 range : F D mat2sys abs sys2mat : t F augment 31 line :

t0 1 time :

Prepare Stem Style spec New 2 col :

1 time t0 - 66 sec *

DC 0 spec 1 el 0022 mat :

<= Fake to view DC

Stem style construct

iii 1 N range q iii el 3 : q N el 2 : 21 line for

ii 1 N range v ii el ii : for

Core component of the stem style P v replicate k 1 : i 1 v length range P i el 1 ≡ continue j 1 P i el range T k j + el v i el : for k k P i el + : 21 line if for T 31 line :

u q v replicate :

i 1 N range S 3 i * el spec i el : for

V u S augment :

V V 3 V rows 12 submatrix :

Unscaled spectrum ... "box car". Many users "scale" the max to 1

V DC 2 1 sys