如何求一个一维信号的功率谱密度的matlab代码和自相关函数的matlab代码

求一个一维信号的功率谱密度：

% 假设信号为 x，采样率为 Fs
N = length(x); % 信号长度
X = fft(x); % 进行FFT变换
PSD = (1/(Fs*N)) * abs(X).^2; % 计算功率谱密度
f = Fs*(0:N/2)/N; % 计算频率
PSD = PSD(1:N/2+1); % 取频率范围在 0 到 Fs/2 的部分
plot(f,PSD); % 画出功率谱密度图
xlabel('Frequency (Hz)');
ylabel('Power/Frequency (dB/Hz)');

求一个一维信号的自相关函数：

% 假设信号为 x
acf = xcorr(x, 'biased'); % 计算自相关函数
n = -(length(x)-1):(length(x)-1); % 自变量
plot(n, acf); % 画出自相关函数图
xlabel('Lag');
ylabel('Autocorrelation');

其中 `xcorr` 函数用于计算自相关函数，参数 `biased` 表示使用有偏估计。