q611.m

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

%《MATLAB及在电子信息课程中的应用》第六章例6.11程序q611
% 方波的频谱分析
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
% 连续信号的付利叶分析
% 演示(1).频域样本点数可取得与时域样本点数不同；
% （2).若要求的频谱太宽，而时域样本点数又取得太少，会发生频率泄漏
% 设已知信号的时域波形f(t),则其付利叶变换可表为

clear,tf=10; 
N = input('取时间分隔的点数N= ');
dt = 10/N;t = [1:N]*dt; 	    % 给出时间分割
f =[ones(1,N/2),zeros(1,N/2)]; 	% 给出信号（此处是方波）
wf = input('需求的频谱宽度wf= ');
Nf = input('需求的频谱点数Nf= ');
w1 =linspace(0,wf,Nf);dw=wf/(Nf-1);
F1 = f*exp(-j*t'*w1)*dt;		% 求付利叶变换
w = [-fliplr(w1),w1(2:Nf)]; 	% 补上负频率
F = [fliplr(F1),F1(2:Nf)]; 		% 补上负频率区的频谱
subplot(1,2,1),
plot(t,f,'linewidth',1.5),      % 画出时间序列
grid on
set(gcf,'color','w')            % 设置图形背景色为白色
axis([0,10,0,1.1])
subplot(1,2,2),
plot(w,abs(F),'linewidth',1.5), % 画出频率特性
grid on