Hello,
I understand that you are referring to the MathWorks documentation example: https://in.mathworks.com/help/signal/ref/pwelch.html#btulru8-5
I note that you want to use a different sampling frequency "fs" and find issue in understanding the relevance of the arguments of the "pwelch()" function.
The mentioned example is using the following syntax:
[pxx,f] = pwelch(signal, window, noverlap, nfft, fs);
The choices of the arguments in the "pwelch()" function can influence the power spectral density (PSD) estimation in the following ways:
- Segment Length (window): The length of each segment affects the trade-off between frequency resolution and time resolution. A longer segment length provides better frequency resolution but poorer time resolution. Conversely, a shorter segment length provides better time resolution but poorer frequency resolution.
- Overlap (noverlap): The number of overlapping samples between adjacent segments affects the smoothing and averaging of the PSD estimate. Increasing the overlap can reduce variance and improve the frequency resolution, but it also increases computational complexity.
- DFT Size (nfft): The number of points used in the discrete Fourier transform (DFT) determines the frequency resolution of the PSD estimate. A larger DFT size provides higher frequency resolution, allowing for more precise identification of spectral features. However, it also increases computational complexity.
- Sample Rate (fs): The sample rate of the input signal is used to compute the frequency axis of the PSD estimate in Hz. It ensures that the frequency axis is correctly scaled and corresponds to the actual frequencies in the input signal.
Therefore, the choice of specific values for these arguments in the mentioned example (x, 500, 300, 500, fs) can be varied according to your requirements depending on the nature of the signal and analysis and it does not only depend on the sampling frequency "fs". Specifically, for the given example, you can choose these values such that it maintains the same ratio with the current values. Also, use a "nfft" so that 100 Hz falls directly on a DFT bin.
Regarding the unit of 10*log10(pxx) in the plotting, it is to represent the estimated PSD in dB/Hz. If you want it to be just W/Hz (or 
), you can remove that conversion.
), you can remove that conversion.To know more about the use-cases of "pwelch()" function, you can refer to the MathWorks documentation in the link:
Further, for better understanding of the Welch method for PSD estimation, you can refer to the MathWorks documentation in the link:
I hope this helps you!
