zxxw.m

来自「最小相位序列的最小相位法复倒谱的实验程序」· M 代码 · 共 37 行

clear all;
x=[20,-12,10,-4];
x=x';
N=length(x);%总长度
a=fft(x);%对x做DFT变换
b=abs(a);%取模运算
d=log(b);%取对数
c=ifft(d);%做IDFT变换
for n=0:(N-1)
   if (n>=1 && n<=N/2)
    f=2;
   end
   if (n==0 || n==N/2)
    f=1;
   end
   if (n<=N-1 && n>N/2)   
    f=0;
   end
    t1(n+1)=c(n+1)*f;   
end
m=0:3;
figure(1);%画原始序列图
stem(m,x);
title('原始序列');
ylabel('x(n)');

figure(2); %画该序列的最小相位法复倒谱图
stem(m,t1);
title('最小相位法复倒谱图');
ylabel('t1(n)');

figure(3);%画该序列的系统函数复倒谱图
t2=cceps(x);
stem(m,t2);
title('系统函数复倒谱图');
ylabel('t2(n)');