Comparing the curves of the graphs

2 次查看(过去 30 天)
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)')
>>

回答(1 个)

Image Analyst
Image Analyst 2022-2-21
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 个评论
ivko ivkovic
ivko ivkovic 2022-2-27
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 2022-2-27
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.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Networks 的更多信息

产品


版本

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by