clc; clear all; close all; fpr=44100; Nx=4*fpr; dt=1/fpr; fm=1; t=dt*(0:Nx-1); f0=0; df=1000; fch1=f0+df*t; f0=2000; df=2000; fch2=f0+df*sin(2*pi*fm*t); f0=100; df=0; fch3=f0*exp(2*t); fch=fch1; if(0) x=sin(2*pi*cumsum(fch)*dt); plot(t,x); pause spectrogram (x, 512, 512-32, 1024, fpr); soundsc(x,fpr) end if(1) f0=1; for k=0:10 x=cos(2*pi*(f0+k*fpr)*t) +k/10; plot(t,x); hold on; end end