MATLAB计算方法和技巧8_2锯齿.docVIP

锯齿波和三角波的傅立叶分析 [问题[数学模型] 方法一：用解析式。锯齿波公式为 (8.2.1) 三角波公式为 (8.2.2) 锯齿波傅立叶展开式为 (8.2.3) 三角波的傅立叶展开式为 (8.2.4) [算法]取周期T为时间单位，取A为波的单位，锯齿波可表示为 x(t*) = t* (8.2.1*) 三角波公式可表示为 (8.2.2*) 其中t* = t/T，x* = x/A。锯齿波的傅立叶级数可表示为 (8.2.3*) 三角波的傅立叶级数可表示为 (8.2.4*) [程序]z_2_1fourier.m如下。 clear %清除变量 a0=0.5; %常数项 tm=2; %周期个数也是最大时间 dt=0.001; %时间增量 t=0:dt:tm; %时间向量(周期的倍数) 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)%标题 hold on %保持图像 plot([0,tm],[a0,a0]) %画水平线(4) h=plot(t,zeros(size(t)),r,EraseMode,Xor);%画线并取句柄(5) text(0,a0,num2str(a0),FontSize,fs) %显示平均值文本 n=1; %求和下标 x=a0; %常量 An=[]; %频谱变量置空 while 1 %无限循环 bn=-1/pi/n; %求锯齿正弦项的系数 An=[An,abs(bn)]; %连接频谱 xn=bn*sin(n*2*pi*t); %求分量 % an=-4/pi^2/(2*n-1)^2; %求三角波余弦项的系数 % An=[An,abs(an)]; %连接频谱 % xn=an*cos((2*n-1)*2*pi*t); %求分量 x=x+xn; %累加分量 set(h,YData,x) %设合成线 plot(t,xn,Color,rand(3,1)) %画分支线 yn=input(还要吗(