Commuted Waveguide

Commuted Waveguide
Synthesis of the
Clavichord

Vesa Va¨ lima¨ ki,* Mikael Laurson,†
and Cumhur Erkut*
*Laboratory of Acoustics and Audio
Signal Processing
Helsinki University of Technology
P.O. Box 3000
FIN-02015 HUT
Espoo, 芬兰
{Vesa.Valimaki, Cumhur.Erkut}@hut.fi
†Centre for Music and Technology
Sibelius Academy
P.O. Box 86
Helsinki, 芬兰
laurson@siba.fi

The clavichord is one of the oldest keyboard instru-
评论, and it is still often used in performances
and recordings of Renaissance and Baroque music.
The sound of the instrument is pleasant and ex-
pressive but quiet. 最后, the instrument
can only be used in intimate performances for
small audiences. This is the main reason why the
clavichord was replaced by the harpsichord and fi-
nally by the modern piano, both of which produce
a considerably louder output. Attempts have been
made to amplify the sound of the clavichord using
a piezoelectric pickup (Burhans 1973).

One of our motivations in this research is to give

the clavichord a new life in the digital world,
where the faint sound level of the instrument can
be amplified by simply turning a volume knob. 这
suggested synthesis model is based on digital wave-
guide modeling of string instruments (史密斯 1992,
1998; Va¨ lima¨ ki et al. 1996; Karjalainen, Va¨ lima¨ ki,
and Tolonen 1998) and uses the principle of com-
muted waveguide synthesis where the soundbox’s
response is incorporated in the excitation signal
(史密斯 1993; Karjalainen and Va¨ lima¨ ki 1993; Karja-
lainen, Va¨ lima¨ ki, and Ja´ nosy 1993). Special sam-
pling techniques are also employed. Musical
examples produced using the proposed synthesizer
will be included on a forthcoming Computer Music
Journal CD.

Acoustics of the Clavichord

A photograph of the clavichord used for the mea-
surements in this study is shown in Figure 1. 这

电脑音乐杂志, 27:1, PP. 71–82, 春天 2003
(西德:1) 2003 麻省理工学院.

instrument is an Anthony Sidey clavichord manu-
factured by Heugel in Paris, 法国, 在 1988. 这
clavichord is an unfretted one, so any combination
of notes can be played. The range of a clavichord is
anywhere from three octaves to over five octaves.
Our clavichord has 51 keys ranging from C2 to D6.
The instrument was tuned about a whole tone
lower than the standard modern tuning
(A4 (西德:2) 440 赫兹): the nominal frequency of A4 is 395
赫兹. The Werkmeister tuning system was used.

For each key of the clavichord, a pair of strings is

tuned in unison, as sketched in Figure 2 (看, 为了
例子, Thwaites and Fletcher 1981; 坎贝尔
and Greated 1987). 然而, the two strings are al-
ways slightly detuned around the same note, 是-
cause exact tuning is impossible manually. Every
key forms one end of a lever that has a tangent at-
tached to its other end. When a key is depressed,
the tangent hits the string pair and initiates vibra-
的. One end of the strings has been damped with
felt, and the other end goes over a bridge to the
tuning mechanism. 因此, the strings are freely vi-
brating between the bridge and the tangent, 哪个
works as both a hammer and a termination. 这
tangent mechanism is rather noisy, as it excites
modes of the soundboard but also itself causes
sound from its moving parts owing to friction.

When the key is released, the tangent falls back
with the aid of gravity, and the string vibration is
allowed to propagate to the felt-covered end of the
string, which efficiently damps the vibration. 在
the end of each note, another knock is heard as the
tangent returns to its resting position. The some-
what mistuned strings of each pair are coupled via
a non-rigid bridge, and thus both beats and a two-
stage decay result (魏因赖希 1977).

Va¨ lima¨ ki, Laurson, and Erkut

71

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

数字 1. Clavichord used
in this study.

数字 2. Tangent mecha-
nism of one key of the
clavichord and a string
pair associated with it.

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

The first thing that people usually notice about
the clavichord is that the sound level is very low.
The maximum sound pressure level at 1 meter is
只关于 50 dB or 60 分贝, depending on the indi-
vidual construction of the instrument. 这使得
this ancient instrument a sound source less effi-
cient than a human speaker. There are many rea-
sons for the weak output (see Campbell and
Greated 1987; Fletcher and Rossing 1991). 这
strings are thin and their tension is low, 和他们
radiate sound inefficiently. The soundboard is
small and light, and it cannot amplify the sound
much.

A particularly interesting feature in the clavi-
chord is the mechanical aftertouch known as Be-
bung. When the player increases pressure on a key,
the tangent is raised more, which in turn increases
the tension of the string pair, resulting in a raise of
沥青. This enables continuous control of vibrato,
which is a used performance style. Aftertouch has
been part of keyboard controllers since the 1980s,
and it is thus easy to include this control in the
clavichord synthesizer. 然而, a fully poly-

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

Felt

Tangent

Key

Lever

Bridge

String
pair

Tuning
pins

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

phonic aftertouch, which would really be needed, 是
only available in high-end keyboard controllers.
As the key pressure is changed, the tangent

causes a slightly different change in the tension of
the two strings of a pair. This feature affects the
timbre of clavichord tones by introducing a flang-
ing effect. It is heard during and after the attack,
particularly in the bass register.

72

电脑音乐杂志

数字 3. (Top) Waveform
of a clavichord tone (A3)
showing the irregular over-
all decay pattern and the
thump at the end of the

tone at 1.1 秒, 和 (bot-
tom) the time history of its
fundamental frequency
(nominal fundamental fre-
quency is 197.5 赫兹).

数字 4. Envelopes of the
lowest partials of the tone
shown in Figure 1: 1st
(solid line with circles),
2ND (solid line), 3rd

(dashed line), 4th (dash-
dot line), and 5th (dotted
线).

1

0


e
v
e
L

(西德:1)1
0

200

199

198

)
z
H

(

0
F

197
0

(西德:1)80

(西德:1)100

(西德:1)120

)


d
(

e
d

t

n
G
A
中号

0.2

0.4

0.6

0.8

1

1.2

0

0.2

0.4

0.6

时间 (秒)

0.8

0.2

0.4

0.6
时间 (秒)

0.8

1

1.2

Properties of Single Clavichord Tones

数字 3 shows the envelope of a recorded clavi-
chord tone. The irregular, non-exponential decay of
the tone can be observed. The key mechanism of
the clavichord generates a loud knock at the begin-
ning and end of a tone, which is characteristic to
the sound of the instrument. 图中 3, a burst lo-
满足于 1.1 sec corresponds to the thump caused
by release of the key, at the same time the signal
starts to decay quickly. The envelope curves of the
first five harmonics are presented in Figure 4. 笔记
that regular exponential decay (IE。, linear decay on
a dB scale) is rare: pronounced beating and other ir-
regularities are observed in many harmonics (例如,
the beating of the fourth partial).

The fundamental frequency of clavichord tones

varies over time, as illustrated by an example in
the lower part of Figure 3. This may be partly
caused by tension modulation (Tolonen, Va¨ lima¨ ki,
and Karjalainen 2000), which is a purely physical
phenomenon in vibrating strings, but also by the
pressure of the player’s finger (IE。, the mechanical
aftertouch), which directly controls string tension
during playing. A time-domain finite-difference
simulation suggests that the tension modulation ef-
fect is negligible compared to that caused by the
mechanical aftertouch (Va¨ lima¨ ki et al. 2000). 在
数字 3, the pitch glide is about 2 赫兹, 这是

lower than the threshold of audibility for this ef-
影响. According to Ja¨ rvela¨ inen and Va¨ lima¨ ki (2001),
the threshold at 196 赫兹 (the nearest frequency
经测试) 是 4.4 赫兹. In the case of the clavichord used
in our measurements, during normal playing the
pitch glide can be inaudible for most listeners. Au-
dible pitch glides can be generated easily, 然而.
数字 5 gives an example of a tone that contains
an exaggerated pitch vibrato, which should be obvi-
ous for all listeners. Note that the vibrato width is
typically about 2 Hz or less, but can occasionally
peak higher. The vibrato rate is about 7 Hz in this
案件, which is slightly faster than what would be
usual in a musical context. The above fundamental
frequency measurements reflect the effective pitch,
where the contributions of both strings are present.
Impulsively driven string instrument sounds are

always at least weakly inharmonic (Fletcher and
Rossing 1991). This is caused by the stiffness of the
string material. It is well known that piano tones,
particularly in the bass range, are highly inhar-
monic, which affects both the timbre and the tun-
英. Measurements of inharmonicity have not been
published previously for the clavichord, so we con-
ducted them ourselves. This was needed to decide
whether we had to implement the inharmonicity
using an allpass filter, as is customary to do in
waveguide synthesis (Jaffe and Smith 1983; Bank
2000). We used our parameter estimation software
to extract partial frequencies of clavichord tones as
a function of time. Using a least-squares fit, we es-
timated from the data the inharmonicity coeffi-
cient. We found that the clavichord used in our
measurements produces almost perfectly harmonic
音调. (The inharmonicity coefficient does not ex-
ceed 10–6.) 因此, the effect of inharmonicity can be
considered inaudible (Ja¨ rvela¨ inen, Va¨ lima¨ ki, 和

Va¨ lima¨ ki, Laurson, and Erkut

73

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

数字 5. (Top) Waveform
of clavichord tone A3 with
strong vibrato (机械的
aftertouch) 和 (底部)
its fundamental frequency
as a function of time.

数字 6. Magnitude spec-
trum of the soundboard
excited with an impulse
hammer at the high-
频率 (顶部) and at the
low-frequency end of the

桥 (底部). Results of
two measurements (solid
and dashed line) 是情节-
ted one on the other for
both cases.

1

0


e
v
e
L

(西德:1)1

0

202

)
z
H

(

0
F

200

198

0.5

1

1.5

2

2.5

196
0

0.5

1
时间 (秒)

1.5

2

2.5

)


d
(

e
d

t

n
G
A
中号

)


d
(

e
d

t

n
G
A
中号

60

40

20

0

60

40

20

0

100

200

300

400

500

100

300
200
Frequency (赫兹)

400

500

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

Karjalainen 2001), and perfectly harmonic synthesis
is adequate.

Synthesis Model for the Clavichord

Soundbox

An essential part of clavichord tones is the re-
sponse of the soundbox. The box is strongly excited
every time a key is pressed. This reverberant re-
sponse will ring for about five seconds, and it
stands out particularly in staccato notes, 哪里的
string vibration itself is damped quickly after the
攻击.

We have conducted a series of acoustic measure-
ments on the soundbox of the instrument to iden-
tify the most prominent modes. 数字 6 shows the
magnitude spectrum of the soundboard that was
hit with an impulse hammer when the strings were
carefully damped with soft material. Results of two
separate measurements are plotted on top of each
其他. The differences between these curves are
quite small, verifying the repeatability of the mea-
surement. Many narrow peaks are visible in Figure
6. They correspond to long ringing resonance
modes. It can also be seen that the excitation of
modes depends on the location of excitation. 为了
例子, 这 90 Hz mode is much stronger in the
top part of the figure than in the lower part. A
similar phenomenon occurs in the case of the mode
at about 340 赫兹.

下一个, we describe the synthesis model for clavi-
chord tones. The commuted waveguide synthesis
method invented by Smith (1993) and Karjalainen
and colleagues (Karjalainen and Va¨ lima¨ ki 1993;
Karjalainen, Va¨ lima¨ ki, and Ja´ nosy 1993) is applied
by using inverse-filtered clavichord signals as exci-
tation for the synthesis model. Inverse filtering
here refers to the processing of a signal with the in-
verted transfer function of a waveguide string
模型. 在这种情况下, the inverse filtering essentially
cancels the partials of a recorded tone. com-
muted waveguide synthesis method has been for-
merly used successfully for the synthesis of the
acoustic guitar (Karjalainen, Va¨ lima¨ ki, and Ja´ nosy
1993; Va¨ lima¨ ki et al. 1996; Laurson et al. 2001) 和
the piano (Smith and Van Duyne 1995).

数字 7 shows the structure of the synthesis
model developed in this work. The coupling of the
two basic strings models Sl(z) and S2(z) is realized
with an unconditionally stable technique suggested
by Karjalainen et al. (1998): the output of only one
of the string models is fed to the input of the other,
and hence there is no feedback. 在实践中, 这
coupling coefficient gc is selected to have a small
价值, but this is not required, because there can be
no stability problems. The input and output signals
of the string models are scaled by constant multi-

74

电脑音乐杂志

数字 7. Clavichord syn-
thesis algorithm for one
key consists of three data-
bases of samples and two
basic string models that
are coupled.

数字 8. Block diagram of
(A) the standard and (乙)
modified string algorithm.
The latter allows control-

ling the sharpness of
attack using parameter
gattack.

Trigger at attack time

Excitation
signals

音色
控制

Trigger at release time

Soundbox
responses

gi1

gi2

S1
gc

gsb

go1

go2

S2

End knocks

grelease

Sound

plying coefficients gi and go, 分别. 在阿迪-
的, two sample databases are needed for realistic
simulation of the soundbox response and the per-
cussive noise caused by key release. In the follow-
英, we discuss the string algorithm, 模型
parameters and their control, and the sample data-
bases.

Modified String Algorithm

The strings of the clavichord can be simulated us-
ing digital waveguide string models. 原则,
two digital waveguide models should be used for
each string, since they have two polarizations of vi-
bration (IE。, horizontal and vertical with respect to
the soundboard), but at present we only use one
digital waveguide per string. For a pair of clavi-
chord strings, we thus have two digital waveguide
string models, not four. This choice results in an
efficient algorithm and still generates some beating
in the synthetic tones.

通常, a digital waveguide string model has the

following transfer function (Jaffe and Smith 1983;
Va¨ lima¨ ki et al. 1996), as also illustrated in Figure
8(A),

S (z) (西德:2)

1
1(西德:3)F (z ) H (z)z

1

(西德:3)L

(1)

其中 F(z) (西德:2) h0
(西德:4) h1z –1 (西德:4) h2z –2 (西德:4) h3z –3 is a third-
order Lagrange interpolation filter (Laakso et al.
1996), L is the (integer) length of the delay line, 和

)(zF

)(zF

)(zH l

)(zH l

(西德:1)
Lz

(A)

G

攻击

Lz(西德:1)

(乙)

Hl(z) is a loop filter with the following transfer func-
的 (Va¨ lima¨ ki et al. 1996):

H (z) (西德:2) G

1

环形

1(西德:4)aloop
zloop

1(西德:4)A

(西德:3)1

(2)

(西德:2) 1 和 –1 (西德:2) aloop

和 0 (西德:2) gloop
(西德:2) 0. The values of
these parameters can be estimated for each string
using previously developed methods (Va¨ lima¨ ki et
阿尔. 1996; Erkut et al. 2000).

While testing an early version of the clavichord
synthesizer, we noticed that the contribution of the
tangent knock was not sufficiently prominent. 在
特别的, the attack of high notes sounded too
soft. To fix this deficiency, we modified the above
string algorithm to allow the attack sharpness to be
controlled with one parameter. The block diagram

Va¨ lima¨ ki, Laurson, and Erkut

75

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

数字 9. Estimated values
of string model parameters
gloop (顶部) and aloop (底部)
(dots connected with

dashed line) together with
their straight-line approxi-
mations (实线).

1

p


1

G

0.995

0.99

0.985
0

0

p


1

A

(西德:1)0.2

(西德:1)0.4

10

20

30

40

50

0

10

20

30
Key index

40

50

of the new algorithm is shown in Figure 8b. 它是
transfer function can be written as

S (z) (西德:2) G

攻击

(西德:4)

(西德:2)

(西德:3)L
F (z) H (z) z

1(西德:3)F (z) H (z) z
G (西德:4)(1(西德:3)G

攻击

攻击

(西德:3)L

) F (z) H (z) z

(西德:3)L

1(西德:3) F (z) H (z) z


(西德:3)L

(3)

(西德:2) 1, this algorithm is identi-

(西德:2) 4, may be used to imitate the aggressive at-

Note that with gattack
cal to the standard plucked-string algorithm of
方程 1. A value larger than unity, 例如
gattack
tack of acoustic tones. In the current version, 这
value of the attack sharpness parameter increases
linearly with key number, so that for C2 it is 2.0
and for D6 it is 4.0. The linear trend and the actual
values were adjusted by informal listening. 有
no physical motivation for this new model struc-
真实, other than that it enables separate control of
two different components of the tone: the tangent
knock and the vibrating string sound.

String Model Parameters

Parameter values for the synthesizer are obtained
from analysis of recorded signals. We use a method
based on the short-time Fourier transform and si-
nusoidal modeling of the clavichord tones. 这些
methods have proven successful earlier in the case

of the acoustic guitar. Details can be found else-
在哪里 (Va¨ lima¨ ki et al. 1996; Erkut et al. 2000).

数字 9 shows the gloop and aloop parameters for
全部 51 keys of the clavichord. These values were es-
timated from recorded clavichord tones. The sec-
ond string model uses the same parameters but
slightly different delay lines lengths for mistuning.
The fundamental frequencies of the lowest (#1) 和
highest (#51) key are 58.7 赫兹和 1056 赫兹, 重新指定-
主动地. Note that the loop gain parameters gloop have
a nearly linear trend as a function of key index (IE。,
logarithm of frequency) with some variations and a
few outliers. Ja¨ rvela¨ inen and Tolonen (2001) 有
investigated the tolerance of parameters value vari-
化. Their results show that –30% to (西德:4)40%
changes in the time constant (which depends on
gloop) are inaudible. 还, in the aloop parameter, vari-
ations between –17% and (西德:4)16% cannot be per-
ceived. It appears that the gloop parameters
presented in Figure 9 can be replaced with values
obtained by linear regression:

环形

G (西德:3) 0.989 (西德:4) 0.000189k
(4)
where k (西德:2) 1. . .51 is the key number starting from
C2 (k (西德:2) 1), and continuing up to D6 (k (西德:2) 51). 这
straight-line fit is illustrated in Figure 9 (顶部) 与一个
solid line.

The values of parameter aloop may be approxi-
mated by a low-order polynomial as well. 然而,
a linear regression to the data shown in the lower
part of Figure 9 yields intolerably large errors. 到
reduce the error, we attempt a piecewise linear re-
gression with a knee point between key indices 23
和 24:

A (西德:3)

环形

(西德:3)0.514 (西德:4) 0.0176k, for k (西德:1) 23

(西德:1)(西德:3)0.175 (西德:4) 0.00225k, for k (西德:2) 24

(5)

This two-part straight-line fit is shown in Figure 9
(底部) with solid lines, one for the low end of the
keyboard and another for the high end. This ap-
proximation for the aloop values introduces minor
changes in timbre that can be perceived, but it does
not necessarily render the sound quality worse. 在
事实, the overall tone color behaves more uniformly
with the linear approximations than with the un-
processed parameter data.

The damping of the string vibrations as a result
of releasing a key is simulated by changing the loop

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

76

电脑音乐杂志

filter characteristics (例如, by setting loop gain gloop
equal to 0.1 momentarily). A similar method has
been used successfully to simulate damping of a
guitar string by finger (Erkut et al. 2000; Laurson et
阿尔. 2001).

Sample Databases

Excitation signals for the synthesis have been ob-
tained by processing anechoic recordings of clavi-
chord tones using the methods described by Erkut
等人. (2000). The procedure consists of sinusoidal
分析, subtracting partials, equalizing the resid-
瓦尔, and truncating the resulting signals with the
right half of a Hanning window. One such excita-
tion signal is used for each of the 51 keys. 电流-
租金地, all excitation signals are 20,000 样品
长的 (关于 0.45 sec at a 44.1 kHz sampling rate).
It would be possible to reduce their length further,
but we have found it unnecessary in our current
implementation, since there is no lack of memory.
The reverberation caused by the soundbox is in-
corporated by triggering a soundbox response sam-
ple at a low level each time any note is played (看
数字 7). This sample must be long enough (关于
five seconds) so that it provides the reverberant
character of the instrument. This is particularly
important for short notes, such as staccato playing,
for which the output signal would otherwise stop
suddenly in an unnatural manner. Note that in
principle we could instead use very long excitation
signals for each note, in which case the full sound-
box response would automatically sound with each
笔记. 然而, we prefer using short excitation sig-
nals, because we can avoid certain problems related
to the inverse filtering (例如, ringing and beating due
to imperfectly cancelled partials). This also saves
记忆, because we use fewer soundbox samples
than excitation signals. Using a separate soundbox
sample also allows adjusting the relative levels of
soundbox reverberation and string sound, 这是
an advantage.

The soundbox samples can be obtained by hitting

the bridge of the clavichord with an impulse ham-
mer at various points. The response must be edited
to remove the sharp attack caused by contact of

metals, 那是, the tip of the hammer and a tuning
pin. This is easily done with a fast fade-in. 笔记
that this editing is not extremely critical, 因为
the impact noise from the tangent mechanism,
which is included in the excitation signal, domi-
nates during the attack and masks the beginning of
the soundbox response. For each key, a soundbox
sample should be chosen, which is obtained close
to the point where the strings are attached. 如何-
曾经, it is practically impossible to excite the bridge
of the clavichord between the tuning pins, 和
only clean samples that we could obtain are from
both ends of the bridge (见图 6). This leaves
us with possibilities to use one of the samples for
the left half of the keyboard and the other for the
right one, or to interpolate between the high- 和
low-end samples along the keyboard. 虽然
sampling-based reverberation modeling is clearly
simplified, it yields a natural-sounding ambience in
the synthetic output and is very cheap to imple-
ment because it requires no filtering.

The knock database for note endings consists of a

few short samples (duration less than 0.5 秒).
These samples may be obtained from recorded long
音调, where the sound of string vibration is al-
lowed to decay below the threshold of audibility
before the key is released. Now the isolated thump
can be simply edited from the recording. One prac-
tical problem was caused by the poor signal-to-
noise ratio of the release samples. It was possible to
hear when a key release occurred from the in-
creased hiss during synthesis. This led us to apply a
de-noising algorithm to suppress the background
noise from the key release recordings, which solved
问题.

It is unnecessary to record a release sample for
each key. Because the tangent mechanism for every
key is identical and the resulting sound is similar,
it is plausible to use a single sample for all syn-
thetic notes. To avoid a boring similarity, as a com-
promise a small collection of release samples may
be used. We have observed that the largest varia-
tion among different tangent releases is caused by
the velocity of release, which slightly changes the
rhythm of different noises within the knock. 因此,
one choice is to gather a collection of tangent re-
leases for different release velocities, and select one

Va¨ lima¨ ki, Laurson, and Erkut

77

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

数字 10. Signals involved
in the synthesis of a single
clavichord tone: (A) excita-
tion signal, (乙) string algo-
rithm output, (C) soundbox

sample, (d) end noise sam-
普莱 (triggered at about 1
秒), 和 (e) synthetic
语气, which is sum of sig-
nals (乙) 到 (d).

0.5
0
(西德:1)0.5

0.5
0
(西德:1)0.5

0.1
0
(西德:1)0.1

0.1
0
(西德:1)0.1

0.5
0
(西德:1)0.5
0

(A)

(乙)

(C)

(d)

(e)

0.2

0.4

0.6

0.8
时间 (秒)

1

1.2

1.4

1.6

randomly or based on some rule, such as the play-
ing tempo. We are currently using two release sam-
普莱斯, and we slightly randomize the release times of
different keys. More sophisticated but still fairly
simple synthesis models employing physically in-
formed randomization have been proposed by Cook
(1997).

Delay Line Variations and the Flanging Effect

Although the variation of the fundamental fre-
quency in clavichord sounds can be inaudible as
discussed earlier, it is also easy to cause pitch
changes with mechanical aftertouch. Because any
change in the string tension will be different for the
two strings of a pair, a flanging-like effect can ap-
pear also in tones where the pitch change is negli-
gible. Variation of the fundamental frequency can
be simulated in the following way: a decaying con-
trol function, 例如, a scaled impulse re-
sponse of a leaky integrator (a one-pole filter), 是
subtracted from the delay-line length of both string
型号. This imitates the change of pressure on the
key by the player’s finger right after depressing the
钥匙 (or the decrease of tension modulation depth
while the string vibration begins to decay), 和
brings about the progression of the fundamental
frequency similar to that shown in Figure 3. 笔记

that this is a simpler way of producing a time-
varying pitch than the tension-modulation tech-
nique Tolonen et al. (2000) used previously for the
synthesis of the kantele (Va¨ lima¨ ki et al. 1999; Er-
kut et al. 2002) and the tanbur (Erkut and Va¨ lima¨ ki
2000). The tension-modulation algorithm would re-
quire computing the elongation of the string in real
time and controlling the delay-line length of the
string model according to it.

For the mechanical aftertouch, a more sophisti-
cated control signal is required. The pitch contour
of a note is generated with a help of a pitch scaler
in the same fashion as in Laurson et al. (2001). 如果
there is neither an initial pitch glide nor vibrato,
the scaler value is 1. The initial pitch glide is im-
plemented as a ramp starting from a value larger
比 1 and ending at 1 after around 0.1 秒. (这
initial value depends on the dynamic level of the
current note so that forte playing yields a higher
value than piano playing.) The vibrato control, 在
转动, is realized with separate parameters for rate
and maximum depth, and a temporal envelope for
深度 (for more details, see Laurson et al. 2001).

A specific feature of the clavichord synthesizer is
that the pitch scaler affects a detuning parameter of
the two string models dynamically. 因此, 一个更多
drastic pitch drift results in a tone that is more out
of tune. This mimics the imperfect mechanism of
the tangent that pulls the two strings unequally. A
database of detuning parameters was generated by
trial and error. An alternative way would be to
measure detuning parameters from the original re-
cordings, which turns out to be difficult. 高的-
pitched tones tend to be more detuned than low
那些. The sound resulting from the time-varying
detuning is suggestive of the flanging effect, 哪个
brings warmth and variation to the synthesis.

Synthesis Example

数字 10 presents the components for the synthe-
sis of a single clavichord tone. The excitation sig-
nal shown in Figure 10a is inserted into the two
string algorithms. The name ‘‘commuted synthe-
sis’’ comes from the principle that this excitation
signal includes both the contribution of the excita-

78

电脑音乐杂志

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

tion mechanism (the tangent) and the soundbox re-
响应. 然而, as can be seen, in our
implementation the input signal contains only the
very beginning of the soundbox response. The tail
of the soundbox response is produced separately us-
ing the sampling-based method.

The summed output signals of the string models
constitute the basic vibrating string component of
the synthetic tone, which is displayed in Figure
10乙. The key release occurs at 0.97 秒; note how
the string tone in Figure 10b decays quickly after
那. Figure 10c shows the beginning of a five-
second-long soundbox sample, which has been
edited from the impulse response of the low-
frequency end of the bridge. A tangent release sam-
普莱, which is started at the release time, is given in
Figure 10d. 最后, Figure 10e presents the sum of
the signals shown in Figures 10b, 10C, and 10d, 和
it is the output of the clavichord synthesizer.

Although remaining almost invisible in Figure
10e, the soundbox sample rings several seconds af-
ter the string tone and the key release sample have
died out. The soundbox sample will also be clearly
audible only after the string tone has been ringing
and decaying for some time, and particularly after
the key has been released. We may note that iden-
tity resynthesis is impossible using the proposed
synthesis model and parameter estimation meth-
消耗臭氧层物质. 尽管如此, the obtained similarity is con-
sidered to be sufficient for high-quality synthesis.
This is confirmed by musical examples produced
with the synthesizer.

Musically Interesting Modifications

The proposed synthesizer structure allows the user
to modify the timbre in numerous ways. 许多
the possible variations have a clear intuitive inter-
预谋. 例如, varying the gain gloop of the
loop filter affects the decay time, but otherwise
the timbre remains unchanged. 然而, varying
the loop filter coefficient aloop changes the decay
time at different frequencies, resulting in a curious
evolution of the spectrum over time. The attack
sharpness parameter and the gain of the note-off
noise also enable meaningful and useful controls.

Exaggerated detuning of the two string models for
each voice leads to a low-quality—or possibly a
very old—clavichord. As a consequence of the
above parametric modifications, the synthetic
tones will still be reminiscent of the clavichord.

Extending the range of an instrument is fascinat-

ing and useful in the case of the acoustic guitar.
Design of a synthetic ‘‘super guitar’’ that has a
wide range, almost like that of a grand piano, 有
been documented recently (Laurson, Va¨ lima¨ ki, 和
Erkut 2002). While an arbitrary change of pitch is
easy in a waveguide synthesizer (simply vary the
length of a delay line), the loop filter parameters
and excitation signals should be extracted from re-
cordings, but for non-existing fundamental frequen-
化学系, they are unavailable. The extension of the
pitch range thus involves extrapolation problems
that must be solved somehow. We are planning
tests on extending the range of the proposed clavi-
chord synthesizer.

Naturally, it is also possible to turn the clavi-
chord synthesizer into a previously unheard virtual
instrument. 例如, modifying or replacing
the signals in the sample databases allows dramatic
variations. Exploring all such transformations was
not included in the goals of this study, and the po-
tential of these modifications thus remains mostly
未知.

Software Implementation of the Synthesizer

The real-time implementation of the proposed clav-
ichord synthesizer has been realized using a
PatchWork (Laurson 1996) user-library called
PWSynth (Laurson and Kuuskankare 2001). 这
high-level part of PWSynth is based on Lisp and
CLOS, but the low-level and time-critical DSP-
routines have been implemented in C. The system
can be played either from musical notation and
control software called Expressive Notation Pack-
年龄 (ENP, described in Kuuskankare and Laurson
2001; Laurson et al. 2001), or from a MIDI key-
board.

The starting point for the clavichord implemen-
tation was to use a guitar synthesizer implemented
早些时候 (Laurson et al. 2001). These systems have

Va¨ lima¨ ki, Laurson, and Erkut

79

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

数字 11. ENP screenshot
from the final cadence of
Par le regard de vos beaux
yeux by Guillaume Dufay.

several features in common. Both use basically the
same dual-polarization string model, although in
the case of the clavichord, it models a pair of
strings and not the orthogonal polarizations of a
single string. The database of samples—needed for
the excitation signals and the noise effects in the
clavichord simulation—could be used in both cases
in a similar way.

The main difference between the implementa-
tions is the fact that the clavichord consists of 51
strings, while the classical guitar has only six. 到
create a synthesizer that could be used in real time,
we designed a voice-allocation algorithm that al-
lows us to effectively play typical pieces of the
clavichord repertoire. We found that a synthesizer
for six simultaneous strings could be played in real
time on a Macintosh G3 400 MHz portable com-
电脑. This is sufficient for a large percentage of old
keyboard music. 新的, faster computers will allow
playing more voices simultaneously.

The control information for the clavichord syn-
thesizer is calculated from an ENP input score. 这
user enters first the basic musical information (IE。,
pitches and rhythms) of the piece using the graphi-
cal front-end of ENP. After this, the user typically
adds expressions and tempo functions that allow
fine tuning of the timing information in the score.
还, one can include expressions specific to the
clavichord, such as the amount of vibrato to be ap-
plied to a given note. 此外, a special rule
modifies the end time of notes (IE。, the time when
the player lifts the fingers from the keyboard). 这
was done to make the performance more lively.
The characteristic noise burst of the instrument
when the fingers are released was found to be es-
sential for the realism of the synthesis.

数字 11 shows the three last measures of a
piece by Guillaume Dufay (1400–1474). Besides
pitch and rhythm information, the piece includes
several expression markings that control both vi-
brato (see the labels starting with ‘‘vb’’) and articu-
关系 (see labels starting with ‘‘st,’’ which stands
for staccato). The numbers after the expression
markings denote how much vibrato or staccato
should be applied to a given note. This example has
also a special expression containing a tempo func-
的 (see the function above the staff marked with
‘‘ritardando’’ in Figure 11), which controls the
tempo of the piece.

Our current MIDI keyboard implementation is
fairly rudimentary and serves only as a basic tool
for testing. The control parameters of the MIDI
keyboard include key number, key velocity, 攻击
时间, and end time of notes. 此外, 我们用
the MIDI aftertouch parameter to simulate the fin-
ger vibrato effect of the clavichord. A more usable
MIDI implementation would require a keyboard
with polyphonic aftertouch. 还, 我们发现
the MIDI keyboard at our disposal was not respon-
sive enough when simulating the finger vibrato. A
playable keyboard simulation of the clavichord
model requires more sensitive hardware than what
is commonly available today.

Conclusions and Future Plans

A simplified physical model was proposed for the
synthesis of the clavichord. It follows the principle
of commuted waveguide synthesis, but also uses
many samples obtained by editing and processing
录音. The excitation signal can be obtained

80

电脑音乐杂志

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

by inverse filtering and truncating a recorded clavi-
chord tone. Two coupled digital waveguide string
models are used for synthesizing each voice. A mi-
nor modification to the plucked string model was
introduced to enable control over attack sharpness,
which was found necessary during testing of the
synthesis algorithm. The reverberation caused by
the soundbox that has many resonant modes was
implemented in a computationally economic way
by triggering a soundbox response sample with
each note. The ending of a note requires another
sample to be played, because there must be a
‘‘knocking’’ sound, which is characteristic to the
clavichord. The mechanical aftertouch of the clavi-
chord can be implemented by varying the delay
line lengths of the waveguide string models over
时间. When these changes are different for the two
string algorithms modeling a pair of strings, a flang-
ing effect is generated, which is also heard in the
real clavichord.

将来, we plan to develop another version

of the clavichord synthesizer using the direct
physical modeling approach, where the excitation,
strings, and the soundbox have a counterpart in the
模型. To accomplish this, the sampling-based im-
plementation of the reverberant response of the
soundboard must be replaced with a soundbox
model that filters the output of the string models.
The soundbox can be considered a small reverber-
ant room and be simulated with an artificial rever-
beration algorithm (例如, 加德纳 1998; Bank 2000).
此外, further measurements of the real in-
strument should be taken. Particularly, the acceler-
ation of the tangent needs to be registered during
the attack of the tone, because this information is
needed for an algorithm simulating the tangent ac-
的, similar to piano hammer simulations (Bank
2000). This might also reveal the importance of the
tension modulation effect in the sound production
of the clavichord. The tangent modeling algo-
rithm—which will be different from physical mod-
els available for the piano hammer, as Hall (1993)
has pointed out—should be developed to be able to
control the excitation using physical parameters.

致谢

The work of M. Laurson is part of the project
‘‘Sounding Score—Modeling of Musical Instru-

评论, Virtual Musical Instruments and their Con-
控制,’’ and the work of C. Erkut is part of the project
‘‘Technology for Speech and Audio Processing.’’
Both projects are funded by the Academy of Fin-
土地. The authors are grateful to Dr. Tero Tolonen
for his contributions in the initial phases of this
学习, to Mr. Paulo A. A. Esquef for his help in the
de-noising of the sample databases, and to Mr.
Jonte Knif for his help in the production of musical
examples using the clavichord synthesizer. Special
thanks to Professor Matti Karjalainen who provided
digital photographs of the clavichord and our mea-
surement session.

参考

Bank, 乙. 2000. Physics-Based Sound Synthesis of the Pi-

ano. Report no. 54, Helsinki University of Technology,
Laboratory of Acoustics and Audio Signal Processing,
Espoo, 芬兰. Available on-line at
www.acoustics.hut.fi/publications/.

Burhans, 右. 瓦. 1973. ‘‘Clavichord Amplification.’’ Jour-
nal of the Audio Engineering Society 21(6):460–463.

坎贝尔, M。, 和C. Greated. 1987. The Musician’s
Guide to Acoustics. 纽约: Schirmer Books,
PP. 234–236.

厨师, 磷. 右. 1997. ‘‘Physically Informed Sonic Modeling
(PhISM): Synthesis of Percussive Sounds.’’ Computer
音乐杂志 21(3):38–49.

Erkut, C。, 等人. 2000. ‘‘Extraction of Physical and Expres-
sive Parameters for Model-Based Sound Synthesis of
the Classical Guitar.’’ Proceedings of the AES 108th
习俗. 纽约: Audio Engineering Society.
Erkut, C。, 和V. Va¨ lima¨ ki. 2000. ‘‘Model-Based Sound
Synthesis of Tanbur, a Turkish Long-Necked Lute.’’
Proceedings of the IEEE International Conference on
Acoustics, Speech, and Signal Processing, 卷. 2. Pis-
cataway, New Jersey: Institute of Electrical and Elec-
tronics Engineers, PP. 769–772.

Erkut, C。, 等人. 2002. ‘‘Acoustical Analysis and Model-
Based Sound Synthesis of the Kantele.’’ Journal of the
Acoustical Society of America 112(4):1681–1691.
弗莱彻, 氮. H。, 和T. D. Rossing. 1991. The Physics of
Musical Instruments. 纽约: Springer-Verlag.

加德纳, 瓦. G. 1998. ‘‘Reverberation Algorithms.’’ In M.
Kahrs and K. Brandenburg, 编辑. Applications of Digital
Signal Processing to Audio and Acoustics. 多德雷赫特,
荷兰人: Kluwer, PP. 85–131.

Va¨ lima¨ ki, Laurson, and Erkut

81

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3

大厅, D. 乙. 1993. ‘‘String Excitation: Piano, Harpsichord
and Clavichord.’’ Proceedings of the 1993 斯德哥尔摩
Music Acoustics Conference. 斯德哥尔摩, 瑞典:
Royal Swedish Academy of Music, PP. 309–314.
Jaffe, D. A。, 和 J. 氧. 史密斯. 1983. ‘‘Extensions of the

Karplus–Strong Plucked-String Algorithm.’’ Computer
音乐杂志 7(2):76–87.

Ja¨ rvela¨ inen, H。, 和T. Tolonen. 2001. ‘‘Perceptual Toler-
ances for Decay Parameters in Plucked String Synthe-
sis.’’ Journal of the Audio Engineering Society 49(11):
1049–1059.

Ja¨ rvela¨ inen, H。, V. Va¨ lima¨ ki, 和M. Karjalainen. 2001.

‘‘Audibility of the Timbral Effects of Inharmonicity in
Stringed Instrument Tones.’’ Acoustics Research Let-
ters Online 2(3):79–84.

Ja¨ rvela¨ inen, H。, 和V. Va¨ lima¨ ki. 2001. ‘‘Audibility of Ini-
tial Pitch Glides in String Instrument Sounds.’’ Pro-
ceedings of the 2001 International Computer Music
会议. 旧金山: International Computer
Music Association, PP. 282–285.

Karjalainen, M。, 和V. Va¨ lima¨ ki. 1993. ‘‘Model-Based

Analysis/Synthesis of the Acoustic Guitar.’’ Proceed-
ings of the 1993 Stockholm Music Acoustics Confer-
恩斯. 斯德哥尔摩, 瑞典: Royal Swedish Academy of
音乐, PP. 443–447.

Karjalainen, M。, V. Va¨ lima¨ ki, and Z. Ja´ nosy. 1993. ‘‘To-

wards High-Quality Sound Synthesis of the Guitar and
String Instruments.’’ Proceedings of the 1993 内特纳-
tional Computer Music Conference. 旧金山: 在-
ternational Computer Music Association, PP. 56–63.

Karjalainen, M。, V. Va¨ lima¨ ki, 和T. Tolonen. 1998.

‘‘Plucked-String Models: From the Karplus–Strong Al-
gorithm to Digital Waveguides and Beyond.’’ Com-
电脑音乐杂志, 22(3):17–32.

Kuuskankare, M。, 和M. Laurson. 2001. ‘‘ENP: Musical
Notation Library based on Common Lip and CLOS.’’
诉讼程序 2001 International Computer Mu-
sic Conference. 旧金山: International Computer
Music Association, PP. 131–134.

Laakso, 时间. 我。, 等人. 1996. ‘‘Splitting the Unit Delay—

Tools for Fractional Delay Filter Design.’’ IEEE Signal
Processing Magazine 13(1):30–60.

Laurson, 中号. 1996. PATCHWORK: A Visual Program-

ming Language and Some Musical Applications. 医生-
toral dissertation, Sibelius Academy, Helsinki,
芬兰.

Laurson, M。, 等人. 2001. ‘‘Methods for Modeling Realis-
tic Playing in Acoustic Guitar Synthesis.’’ Computer
音乐杂志 25(3):38–49.

Laurson, M。, 和M. Kuuskankare. 2001. ‘‘PWSynth: A

Lisp-based Bridge between Computer Assisted Compo-
sition and Model-based Synthesis.’’ Proceedings of the
2001 International Computer Music Conference. 桑
Francisco: International Computer Music Association,
PP. 127–130.

Laurson, M。, V. Va¨ lima¨ ki, 和C. Erkut. 2002. ‘‘Produc-

tion of Virtual Acoustic Guitar Music.’’ Proceedings of
the AES 22nd International Conference on Virtual,
Synthetic, and Entertainment Audio, 纽约: Au-
dio Engineering Society, PP. 249–255.

史密斯, J. 氧. 1992. ‘‘Physical Modeling Using Digital

Waveguides.’’ Computer Music Journal 16(4): 74–91.
Available on-line at www-ccrma.stanford.edu/(西德:4)jos/
pmudw/.

史密斯, J. 氧. 1993. ‘‘Efficient Synthesis of Stringed Musi-
cal Instruments.’’ Proceedings of the 1993 内特纳-
tional Computer Music Conference. 旧金山:
International Computer Music Association, PP. 64–71.
Available on-line at www-ccrma.stanford.edu/(西德:4)jos/cs/.
史密斯, J. 氧. 1998. ‘‘Principles of Digital Waveguide Mod-
els of Musical Instruments.’’ In M. Kahrs and K. Bran-
denburg, 编辑. Applications of Digital Signal Processing
to Audio and Acoustics. 多德雷赫特, 荷兰人:
Kluwer, PP. 417–466.

史密斯, J. 奥。, 和S. A. Van Duyne. 1995. ‘‘Commuted Pi-
ano Synthesis.’’ Proceedings of the 1995 国际的
电脑音乐会议, 旧金山: 内特纳-
国家计算机音乐协会, PP. 319–326.
Thwaites, S。, and N. H. 弗莱彻. 1981. ‘‘Some Notes on
the Clavichord.’’ Journal of the Acoustical Society of
美国 69(5):1476–1483.

Tolonen, T。, V. Va¨ lima¨ ki, 和M. Karjalainen. 2000.

‘‘Modeling of Tension Modulation Nonlinearity in
Plucked Strings.’’ IEEE Transactions on Speech and
Audio Processing 8(3):300–310.

Va¨ lima¨ ki, 五、, 等人. 1996. ‘‘Physical Modeling of Plucked
String Instruments with Application to Real-Time
Sound Synthesis.’’ Journal of the Audio Engineering
社会 44(5):331–353.

Va¨ lima¨ ki, 五、, 等人. 1999. ‘‘Nonlinear Modeling and Syn-
thesis of the Kantele—a Traditional Finnish String In-
strument.’’ Proceedings of the 1999 国际的
电脑音乐会议. 旧金山: 内特纳-
国家计算机音乐协会, PP. 220–223.
Va¨ lima¨ ki, 五、, 等人. 2000. ‘‘Model-Based Synthesis of the
Clavichord.’’ Proceedings of the 2000 国际的
电脑音乐会议. 旧金山: 内特纳-
国家计算机音乐协会, PP. 50–53.

魏因赖希, G. 1977. ‘‘Coupled Piano Strings.’’ Journal of
the Acoustical Society of America 62(6):1474–1484.

82

电脑音乐杂志

D

w
n

A
d
e
d

F
r


H

t
t

p

:
/
/

d

r
e
C
t
.


t
.

e
d

/
C

j
/

A
r
t

C
e

p
d

F
/

/

/

/

/

2
7
1
7
1
1
8
5
3
8
2
7
0
1
4
8
9
2
6
0
3
6
0
6
1
3
3
5
3
p
d

.

F


y
G

e
s
t

t


n
0
8
S
e
p
e


e
r
2
0
2
3Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image
Commuted Waveguide image

下载pdf