chirpPh.m

Chirp (sine sweep) with defined starting phase angle

This is function that behaves like chirp (linear), except you can specify an initial phase angle. The purpose of this is to be able to have a series of several chirp time histories with seamless transitions between them. There is a second output which is starting phase for the next time series.
An example for a signal starting from zero Hz and ramping up to 40Hz and then back down to 10Hz is shown:

[d(1:1000),ph]=chirpPh(.001:.001:1,0,40,0);
[d(1001:3000)]=chirpPh(.001:.001:2,40,10,ph);