chirpPh.m

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);