clear; close all
Fs=1; % Sampling frequency [Hz]
dt=1/Fs; % Time sample interval [s]
n=1024; % signal length
T=dt*(n-1); % Signal duration [s]
t=(0:dt:T)’; % Time vector
A0=1+1.3.*cos(2*pi*0.008*t); % Signal amplitude
f0=0.03.*ones(n,1); % Signal frequency [Hz]
x=A0.*cos(2*pi*cumtrapz(f0.*dt)); % Signal
figure(1)
plot(t,x,’.-‘);
grid on; axis tight
ylabel(‘Amplitude’);
xlabel(‘Time, s’);
title(‘Harmonic signal’)
figure(2)
strips(x,256)
ylabel(‘Amplitude’);
xlabel(‘Points’)
title(‘Signal in strips’)
figure(3)
[Pxx,F] = pwelch(x,[],[],1024,Fs); % Create a Welch spectrum
plot(F,Pxx);
ylabel(‘Power spectral density’);
xlabel(‘Frequency, Hz’);
title(‘Spectrum’)