dd.m

来自「很多matlab的源代码」· M 代码 · 共 43 行

 %离散傅里叶变换及x(n),FT[x(n)]的绘图函数
 function [c,l]=DFT(x,N,str)
 n=0:N-1;
 k=-200:200;
 w=(pi/100)*k;
 l=w;
 c=x*(exp(-j*pi/100)).^(n'*k);  %计算DFT[x(n)]
 magX=abs(c);                    %幅值
 angX=angle(c);                 %辅角
 
 subpolt(2,2,1);
 t=max(x);
 n=0:N-1;
 stem(n,x,'.');
 xlabel('n')
 ylabel('xa(n)');
 title('信号的原型');
 text((0.3*N),(0.8*t),str);
 hold on
 n=0:N-1;
 m=zeros(N);
 plot(n,m);
 
 subpolt(2,2,3);
 plot(w/pi,magX);
 xlabel('w/pi');
 ylabel('|Y(jw)|');
 title('上图信号的傅氏变换');
 
 %产生信号xa(n)
 function c=FF(A,a,w,fs)
 n=0:50-1;
 c=A*exp((-a)*n/fs).*sin(w*n/fs).*stepseq(0,0,49);
 
 %产生脉冲信号
 function [x,n]=impseq(n0,n1,n2)
 n=[n1:n2];
 x=[(n-n0)==0)];
 
 %产生矩形框信号
 function [x,n]=stepseq(n0,n1,n2)
 n=[n1:n2];
 x=[(n-n0)>=0];