dft.m

一些matlab底下的小函数

%x为离散时间序列，N为其长度,本程序将对其进行DFT分析
%xname,yname分别为信号x和FT(x)的名称，为字符窜型
%str也为字符窜型，是对x的补充说明
function [c,l]=DFT(x,N,xname,yname,str)
n=0:N-1;
k=-200:200;
w=(pi/100)*k;
l=w;
c=x*exp(-j*pi/100).^(n'*k);
magX=abs(c);                      %幅度特性
angX=angle(c);                    %频率特性

subplot(1,2,1);
t=max(x);
n=0:N-1;
stem(n,x,'.');
xlabel('n');
ylabel(xname);
title('信号的原形');
text((0.3*N),(0.8*t),str);
hold on
n=0:N-1;
m=zeros(N);
plot(n,m);

subplot(1,2,2);
plot(w/pi,magX);
%axis([-2,2,,])
xlabel('w/pi');
ylabel(yname);
title('左图信号的FT');