q709a.m

<matlab在电子信息中的应用>一书中个章节的配套源代码。

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

disp('八点时间信号的离散时间傅立叶变换')
x0=randn(1,8);			 % 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 on,shg		% 画图
text(0,3.5,'八点时间信号的离散时间傅立叶变换')
disp('重复N次的八点时间信号的离散时间傅立叶变换')
N=input('N= ');						% 用键盘输入延拓周期数
if isempty(N) N=5; end
x1=reshape(x0'*ones(1,N),1,N*length(x0));	% 延拓后的时域信号x1
X1 = dtft(x1,w)*dt/N;				% 求x1的频率响应X1
subplot(4,1,2),plot(w,abs(X1)),grid on,shg		% 画图
text(0,3.5,'重复N次的八点时间信号的离散时间傅立叶变换')
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 on,shg		% 画图
text(0,3.5,'重复100次的八点时间信号的离散时间傅立叶变换')
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,shg
axis([min(w),max(w),0,max(abs(X2))]),grid on,shg
text(0,3.5,'重复无穷次的八点信号的离散傅立叶变换')