q709.m

matlab在电子信息中的应用和实例中的源代码

%《MATLAB及在电子信息课程中的应用》第七章例7.9程序q709
% 离散时间傅立叶变换
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
disp('八点时间信号的离散时间傅立叶变换')
x0=sin(2*pi* [1:8]/8) *5;			        % x0是8点行向量
dt=2*pi/8;
w=linspace(0,2*pi,1000)/dt;	                % w是1000点行向量
X0 = dtft(x0,w) *dt;				        % 求得频率响应X0
subplot(4,1,1),plot(w,abs(X0)),grid,shg		% 画图
set(gcf,'color','w')                        % 设置图形背景色为白色
for k=1:2
    disp('重复N次的八点时间信号的离散时间傅立叶变换(可输入两次)')
    N=input('N= ');						    % 用键盘输入延拓周期数
    x1=reshape(x0'*ones(1,N),1,N*length(x0));	% 延拓后的时域信号x1
    X1 = dtft(x1,w) *dt/N;				    % 求x1的频率响应X1
    subplot(4,1,k+1),plot(w,abs(X1)),
    grid,shg on,pause		                % 画图
end
disp('重复无穷次的八点信号的离散傅立叶变换')
pause,X2=fft(x0*dt);						% 离散傅立叶变换
w1=2*pi* [0:length(x0)-1]/length(x0);		% 离散频点向量
subplot(4,1,4)
stem([-w1,w1],[abs(X2),abs(X2)]),grid on,shg
axis([min(w),max(w),0,max(abs(X2))])
grid on,shg