新·MATLAB基础与应用教程(人民邮电出版社-蔡旭辉)第六章b.pptVIP

6.2　傅里叶分析 6.2.1　周期信号频域分析的MATLAB实现　　例　设周期性对称三角波幅度A=1,周期T=2,试用  MATLAB画出其频谱。 图1　周期性对称三角波 解　傅里叶系数 MATLAB程序如下:N=8;n1=-N:-1;c1=-4*j*sin(n1*pi/2)/pi^2./n1.^2; %计算n=-N到-1的傅里叶系数c0=0;%计算n=0时的傅里叶系数n2=1:N;c2=-4*j*sin(n2*pi/2)/pi^2./n2.^2; %计算n=1～N的傅里叶系数cn=[c1 c0 c2];n=-N:N;subplot(2,1,1);stem (n, abs ( cn ),.);ylabel (Cn的幅度);subplot(2,1,2);stem (n, angle (cn),.);ylabel(Cn的相位);xlabel(\omega/\omega0); 例6.7 方波分解为多次正弦波之和 %第六章例6.7程序q607 %方波分解为多次正弦波之和 t = 0:.01:2*pi; % 设定一个时间数组 y = sin(t);plot(t,y),pause % 频率为w=1(f=1/2π)的正弦基波 y = sin(t)+sin(3*t)/3;plot(t,y),pause % 叠加三次谐波 % 用1，3，5，7，9次谐波叠加 y = sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/9; plot(t,y) % 为了绘制三维曲面，要把各次波形数据存为一个三维数组，因此必须重新定义y y = zeros(10,max(size(t))); x = zeros(size(t)); for k=1:2:19 x = x + sin(k*t)/k; y((k+1)/2,: ) = x; end % 将各波形叠合绘出 pause, figure(1),plot(t,y(1:9,: )),grid line([0,pi+0.5],[pi/4,pi/4]) % 加上方波幅度线及标注 text(pi+0.5,pi/4,pi/4) % 将各半波形绘成三维网格图，看出增加谐波阶次对方波逼近程度的影响 halft=ceil(length(t)/2);pause, figure(2),mesh(t(1:halft),[1:10],y(:,1:halft)) %只用正半周波形 6.2.2　非周期信号频域分析的MATLAB实现　　MATLAB提供了许多数值计算工具，可以用来进行信号的频谱分析。quadl是计算数值积分的函数，有下面两种调用方式: 　　 y = quadl(‘F’,a,b)　　　y = quadl(‘F’,a,b,［］,［］,P)其中，F是一个字符串，表示被积函数的文件名；a、b分别表示定积分的下限和上限；P表示被积函数中的一个参数。quad8的返回值是用自适应Simpson算法得出的积分值。 （ *quadl 采用 自适应 Lobatto　算法） 程序运行结果 例 门信号如图所示，试计算宽度τ=1和幅度A=1的门信号p1(t)在0～fm(Hz)频谱范围内所含的信号能量。 　　计算上式的MATLAB程序如下:　　　function y=sf2(t) y=2*sinc(t).*sinc(t); 　　　 % program10.4－4 f=linspace(0,5,256); N=length(f);w=zeros(1,N); for k=1:N w(k)=quad8(sf2,0,f(k)); end plot(f,w); xlabel(Hz);ylabel(E);　　　　　　　 程序运行结果 例6.11 非周期方波的频谱分析 傅里叶变换表示式为 根据图6.11－1非周期信号波形的情况，将积分上下限定为0～10s，并将t分成N等份，用求和代替积分。这样，傅立叶变换式可写为 求和可以用f(t)行向量乘以e-j?tn列向量来实现。式中Δt是t的增量，在程序中，用dt表示。 例6.11 非周期方波的频谱分析(续) 求不同?处的F值，都用同一公式，这就可以利用MATLAB中的元素群运算。将?设为一个行数组，代入上式，则可写