Effacer les filtres
Effacer les filtres

Comparing the curves of the graphs

5 vues (au cours des 30 derniers jours)
ivko ivkovic
ivko ivkovic le 21 Fév 2022
Commenté : Image Analyst le 27 Fév 2022
Comparing the curves
>> NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Plot single-sided amplitude spectrum.
plot(f,2*abs(Y(1:NFFT/2+1)))
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
>> Fs = 1000; % Sampling frequency
T = 1/Fs; % Sample time
L = 1000; % Length of signal
t = (0:L-1)*T; % Time vector
% Sum of a 50 Hz sinusoid and a 120 Hz sinusoid
x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);
y = x + 2*randn(size(t)); % Sinusoids plus noise
plot(Fs*t(1:50),y(1:50))
title('Signal Corrupted with Zero-Mean Random Noise')
xlabel('time (milliseconds)')
>>

Réponses (1)

Image Analyst
Image Analyst le 21 Fév 2022
The spectrum (left plot) shows a spike at 50 and 120 Hz. I'd never know that there is monotone noise at those exact frequencies just by looking at the time domain plot on the right. So I'd say the frequency domain plot is most informative in this particular case.
  4 commentaires
ivko ivkovic
ivko ivkovic le 27 Fév 2022
basically I need to add function that the pure signal and is seen on the graphs
tried modiyfing the code and adding x = 0.7*sawtooth(2*pi*50*t);%Function to calculate Sawtooth wave but I get an error. Any ideas?
Image Analyst
Image Analyst le 27 Fév 2022
I don't know why you did that instead of what I suggested. Anyway, attach your y in a .mat file with the paperclip icon.

Connectez-vous pour commenter.

Catégories

En savoir plus sur Networks dans Help Center et File Exchange

Produits


Version

R2021b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by