hc463.m

matlab数字信号处理代码

%《数字信号处理教程——MATLAB释义与实现》第四章例4.6.3程序hc463
% 用IDFT由已知频谱求连续时间傅立叶反变换IDTFT的信号
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
wc=3;Tmax=0.1*pi/wc      	  % 最大采样周期Tmax，取算出的临界数值的1/10
T=input('(T<Tmax)T= ')        % 用户选择并输入T
Nmin=100*2*pi/wc/T	        % 最小样点数目Nmin
N=input('(N>Nmin)N=  ');    % 用户选择并输入N
D=2*pi/(N*T); 		        % 频谱序列长度及频率分辨率
%kn=floor(-(N-1)/2:-1/2); 	% 负频率下标向量
%kp=floor(0:(N-1)/2); 		% 正频率下标向量
%w=[kp,kn]*D;     		    % 将负频率移到正频率的右方，形成新的频率排序
w=[0:N-1]*D;
M=floor(wc/D);             % 有效频率的边界下标（整数）
Xa=[ones(1,M+1),zeros(1,N-2*M-1),ones(1,M)];	% 按新的频率排序输入给定的频谱
xa=ifft(Xa/T);               % Xa/T是序列频谱X(k)，其IFFT为数字序列x(n)=xa(n) 
%xa=fftshift(ifft(X)/T);           % 对X求IFFT,并除以T,得出连续信号xa(t) 
plot([0:N-1]*T,xa)			% 用plot语句绘图，实现插值
set(gcf,'color','w');              % 置图形背景色为白