fs = 16000; % частота дискретизации f0 = 100; % начальная частота f1 = 8000; % конечная частота fcarrier = (f0 + f1) / 2; omega = f1 - f0; tf = 2; % время излучения n = 0; for i = 1:5 t = n : 1/fs : tf + n - 1/fs; t1 = -tf/2 + n : 1/fs : tf/2 + n - 1/fs; s = cos(2 * pi * fcarrier .* t + 2 * pi * omega * t1.^2 / (2*tf)); plot(t, s); Spc = abs(fft(s)); % амплитудный спектр fx = fs/2*linspace(-1,1,length(Spc)); % значения частот figure, plot(fx,fftshift(Spc)) % спектр grid on xlim([0 20]) title('Спектр ЛЧМ сигнала') xlabel('частота, Гц') sound(s,fs); filename='C:\Chirp\chirp.wav'; audiowrite(filename,s,fs) end