有没有会matlab的大神，帮帮忙写个代码。

你有数据文件吗？发过来看看，我写个程序。
clc
clear
Fs = 100;
load ep42
s=s';
figure
h = spectrum.welch;
psd(h,s,'Fs',Fs);%功率密度估计图
title('ep42文件功率密度估计')
[pxx,f] = pwelch(s,500,300,500,Fs);
figure
plot(f,10*log10(pxx)) %画频率谱
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')
grid
title('ep42文件频率谱')

[s1,Fs1] = audioread('Noct1.wav');
%sound(s1,Fs1);
figure
h = spectrum.welch;
psd(h,s1,'Fs',Fs1);%功率密度估计图
title('Noct1功率密度估计')
[pxx,f] = pwelch(s1,500,300,500,Fs1);
figure
plot(f,10*log10(pxx)) %画频率谱
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')
grid
title('Noct1频率谱')