请高手帮忙

胡磊

clc
clear
clear glbal
[name,path]=uigetfile('*.txt','请输入你要装载的文件1：','文件名称');
WJLoad=[path name];
S=load(WJLoad);T=1.25;Fs=50000;
global N1
N1=8192;
x=S;n=0.125*Fs;
a=T*Fs;N0=round(a/n);
d_125=zeros(8*n,N0);d_125_fft=zeros(8*n,N0);
for i=1:N0
  d_125(1:n,i)=x(1+(i-1)*n:i*n);
  f1 = 20.598997;
  f2 = 107.65265;
  f3 = 737.86223;
  f4 = 12194.217;
  A1000 = 1.9997;
  pi = 3.14159265358979;
  NUMs = [ (2*pi*f4)^2*(10^(A1000/20)) 0 0 0 0 ];
  DENs = conv([1 +4*pi*f4 (2*pi*f4)^2],[1 +4*pi*f1 (2*pi*f1)^2]);
  DENs = conv(conv(DENs,[1 2*pi*f3]),[1 2*pi*f2]);
  [B,A] = bilinear(NUMs,DENs,Fs);
  %adsgn,生成与采样频率相关的A计权系数滤波器
  d_125(:,i)=filter(B,A,d_125(:,i));
  d_125_fft(:,i)=abs(fft(d_125(:,i)*hanning(n),8*n)/length(d_125(:,i))*2);
end
Y=(0:N0-1)*0.125;
f=Fs*(0:8*n-1)/8*n;
pcolor(f,Y,d_125_fft);
Xlable('频率f');
Ylable('时间t');Zlable('声压dB（A）');
请高手帮我看看这段程序有什么问题！！！！

wxm27811

你的文件呢？我这里没有合适的文件测试。

胡磊

2# wxm27811

胡磊

有一点改动啊x=s(:,2),请您注意下！！！