q709.m

这是matlab在电子信息课程中的应用

%《MATLAB在电子信息课程中的应用》第七章例7.9程序q709
% 离散时间傅立叶变换
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

disp('八点时间信号的离散时间傅立叶变换')
x0=sin(2*pi* [1:8]/8) *5;			 % x0是8点行向量
dt=2*pi/8;
w=linspace(-2*pi,2*pi,1000)/dt;	 % w是1000点行向量
X0 = dtft(x0,w)*dt;				% 求得频率响应X0
subplot(3,1,1),plot(w,abs(X0)),grid,shg		% 画图
disp('重复N次的八点时间信号的离散时间傅立叶变换')
N=input('N= ');						% 用键盘输入延拓周期数
x1=reshape(x0'*ones(1,N),1,N*length(x0));	% 延拓后的时域信号x1
X1 = dtft(x1,w)*dt;				% 求x1的频率响应X1
subplot(3,1,2),plot(w,abs(X1)),grid,shg		% 画图
%disp('重复100次的八点时间信号的离散时间傅立叶变换')
%x1=reshape(x0'*ones(1,100),1,100*length(x0));	% 延拓后的时域信号x1
%X1 = dtft(x1,w)/100*dt;				% 求x1的频率响应X1
%subplot(4,1,3),plot(w,abs(X1)),grid,shg		% 画图
disp('重复无穷次的八点信号的离散时间傅立叶变换-傅立叶级数')
pause,X2=fft(x0*dt);							% 离散傅立叶变换
w1=2*pi* [0:length(x0)-1]/length(x0);			 % 离散频点向量
subplot(3,1,3),stem([-w1,w1],[abs(X2),abs(X2)]),grid,shg
axis([min(w),max(w),0,max(abs(X2))]),grid,shg
%gtext('八点时间信号的离散时间傅立叶变换')
%gtext('重复N次的八点时间信号的离散时间傅立叶变换')
%gtext('重复100次的八点时间信号的离散时间傅立叶变换')
%gtext('重复无穷次的八点信号的离散傅立叶变换')