%% clear;clc; load mtlb; %首先进行STFT变换 [S,F,T,P]=spectrogram(mtlb,1024,1023,1024); figure(1); subplot(211); % plot(mtlb/max(abs(mtlb)));grid on;title('signal of mtlb') plot(mtlb/abs(mtlb));grid on;title('signal of mtlb') subplot(212); mesh(abs(S)); title('Result Of Stft'); %然后进行谱峰频率图绘制 figure(2); [fre_max,index]=max(abs(S),[],1);%记录每一帧谱峰所对应的频率 index2=1:size(S,2); stem3(index,index2,fre_max); grid on;title('peaks'); %采用滤波器法重建mtlb信号 WL=1024;%窗长 N=512;%512点DFT M=1024;%语音帧的大小 IN=mtlb;L=length(IN);window=hann(WL);Mod=M-mod(L,M);Q=(L+Mod)/M; IN=[IN;zeros(Mod,1)]/max(abs(IN));%归一化 OUT=zeros(length(IN),1); X=zeros(M,(N/2+1)); Z=zeros(WL-1,(N/2+1)); t=(0:M-1)'; WN1=zeros(M,(N/2+1)); WN2=zeros(M,(N/2+1)); for k=1:(N/2+1) w=2*pi*i*(k-1)/N; WN1(:,k)=exp(-w*t); WN2(:,k)=exp(w*t); end for p=1:Q R=IN((p-1)*M+1:p*M); for k=1:(N/2+1) x=R.*WN1(:,k); [X(:,k),Z(:,k)]=filter(window,1,x,Z(:,k));%加窗滤波，stft变换 end X1=X.*WN2; A=zeros(M,1); for j=2:(N/2) A=A+X1(:,j); end Y((p-1)*M+1:p*M)=2*real(A)+real(X1(:,1)+X1(:,65));%求和 end OUT=Y(1:L)/max(abs(Y)); figure(3); plot(OUT);

这段代码是一个用于进行语音信号的STFT变换、谱峰频率图绘制以及采用滤波器法重建信号的MATLAB程序。STFT（短时傅里叶变换）是一种对非平稳信号进行频谱分析的方法，其基本思想是将信号分为多个时间窗口，然后对每个时间窗口内的信号进行傅里叶变换。谱峰频率图是对STFT结果进行处理得到的，它可以反映出每个时间窗口内信号的主要频率成分。而采用滤波器法重建信号则是将STFT的结果进行逆变换，重建出原始信号。