load /radio/steve/perspectivepics/red.txt

g = red-mean(red);
h = dechirp(g,1000,.23+.5); % as observed by Shawn, folding prints in half

h = h(1:970);  % make it truly peroidic

hbig = [h;(h);h;h;(h);h;h;h;(h);h;(h);h;(h);h;h;(h)];
gbig = dechirp(hbig,length(hbig),.23); % back to chirping again

disp('I actually used version of tf.m which removed dc each window; deleteme.m')
% stop now and do just that
% Gbig = tf(gbig,1000,100,2);
Gbig = deleteme(gbig,1000,100,2);

dispmatb(-abs(Gbig))

