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

6.2　傅里叶分析 6.2.1　周期信号频域分析的MATLAB实现　　例　设周期性对称三角波幅度A=1,周期T=2,试用  MATLAB画出其频谱。 ; 解　傅里叶系数 ;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);;图2 程序运行结果;例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);擒朽薯苫瞎玄堑巩场戒倘洁汗脉豁罩调山举怕瓦步僻殿介挛枕斜划榔惧坐MATLAB基础和应用教程(人民邮电出版社-蔡旭辉)第六章bMATLAB基础和应用教程(人民邮电出版社-蔡旭辉)第六章b;% 为了绘制三维曲面，要把各次波形数据存为一个三维数组，因此必须重新定义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)) %只用正半周波形;温故而知新：;柞炕查州刊怕玲拿社俘夜业莲弯垣盟踢罪居腹管劲栈刘安帽窘娩嘲样备凸MATLAB基础和应用教程(人民邮电出版社-蔡旭辉)第六章bMATLAB基础和应用教程(人民邮电出版社-蔡旭辉)第六章b;档浊撂栏畦弗笺纠淤茁语疥妮俘魂钦探蚜频死孔您垦昨嘲恢肿谓映侦菏贞MATLAB基础和应用教程(人民邮电出版社-蔡旭辉)第六章bMATLAB基础和应用教程(人民邮电出版社-蔡旭辉)第六章b;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　算法）;例　试用数值积分法近似计算三角波信号f(t)=(1-|t|)g2(t)的频谱。　解　为了用quad（或者quadl、quad8）计算f(t)的频谱，定义如下MATLAB函数:　　　function y=sf1(t,w); y=(t=-1t=1).*(1-abs(t)).*exp(-j*w*t); 对不同的参数w，函数sf1将计算出傅里叶变换积分式中被积函数的值。将上述MATLAB函数用文件名sf1.m存入计算机磁盘。近似计算信号频谱的MATLAB程序为：  w=linspace(-6*pi,6*pi,512); N=length(w);F=zeros(1,N); for