q611.m

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

%《MATLAB在电子信息课程中的应用》第六章例6.11程序q611
% 方波的频谱分析
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

% 连续信号的付利叶分析
% 演示(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
axis([0,10,0,1.1])
subplot(1,2,2),plot(w,abs(F),'linewidth',1.5),grid