阅读材料锯齿波和三角波的傅立叶分析.docVIP

[阅读材料]锯齿波和三角波的傅立叶分析 [问题[数学模型] 方法一：用解析式。锯齿波公式为 (6_12_1) 三角波公式为 (6_12_2) 锯齿波傅立叶展开式为 (6_12_3) 三角波的傅立叶展开式为 (6_12_4) [算法]取周期T为时间单位，取A为波的单位，锯齿波可表示为 x(t*) = t* (6_12_1*) 三角波公式可表示为 (6_12_2*) 其中t* = t/T，x* = x/A。锯齿波的傅立叶级数可表示为 (6_12_3*) 三角波的傅立叶级数可表示为 (6_12_4*) [程序]zhou12_fourier.m如下。 clear %清除变量 a0=0.5; %常数项 tm=2; %周期个数也是最大时间 dt=0.001; %时间增量 t=0:dt:tm; %时间向量(周期的倍数) % t0=[0,0,0.5,0.5,1]; %矩形波自变量 % x0=[0,1,1, 0, 0]; %矩形波函数 t0=[0,1,1]; %锯齿形波自变量 x0=[0,1,0]; %锯齿形波函数 % t0=[0,0.5,1]; %三角齿形波自变量 % x0=[0,1,0]; %三角齿形波函数 figure %创建图形窗口 plot([t0,t0+1],[x0,x0],linewidth,2) %画矩形波 grid on %加网格 axis([0,tm,-0.65,1.2]) %曲线范围 fs=16; %字体大小 xlabel(\itt\rm/\itT,fontsize,fs) %x标签 ylabel(\itx\rm/\itA,fontsize,fs) %y标签 %title(方波的傅立叶分析,fontsize,fs)%标题 title(锯齿波的傅立叶分析,fontsize,fs)%标题 %title(三角齿波的傅立叶分析,fontsize,fs)%标题 hold on %保持图像 plot([0,tm],[a0,a0]) %画水平线(4) h=plot(t,zeros(size(t)),erasemode,xor,color,r);%画线并取句柄(5) text(0,a0,num2str(a0),fontsize,fs) %显示平均值文本 n=1; %求和下标 x=a0; %常量 An=[]; %频谱变量置空 while 1 %无限循环 % bn=2/pi/(2*n-1); %求方波正弦项的系数 % An=[An,bn]; %连接频谱 % xn=bn*sin((2*n-1)*2*pi*t); %求分量 bn=-1/pi/n; %求锯齿正弦项的系数 An=[An,abs(bn)]; %连接频谱 xn=bn*sin(n*2*pi*t); %求分