matlab的myPsdCal函数

myPsdCal函数的功能是计算信号的功率谱密度和频谱。以下是一个示例的myPsdCal函数，它使用Welch方法来计算信号的功率谱密度和频谱。

function [psd,freq] = myPsdCal(x, fs, nfft)
% x: 输入信号
% fs: 采样率
% nfft: fft长度

% 计算重叠段数
window = hann(nfft);
noverlap = round(nfft/2);

% 使用Welch方法计算信号的功率谱密度
[psd, freq] = pwelch(x, window, noverlap, nfft, fs);

% 计算信号的频谱
freq = linspace(0, fs/2, length(psd));

该函数输入参数为x，fs和nfft。其中，x是输入信号，fs是采样率，nfft是FFT长度。函数使用hann窗口来计算重叠段数，使用pwelch函数计算信号的功率谱密度和频率向量。最后，计算频率向量freq并返回功率谱密度和频率向量。