吉布斯现象的matlab代码实现.pdfVIP

第一题： 代码 clear; close all ; clc; t = -2::2; T = [15,25,51,101]; for k = 1:length(T) a=0; for n = 1:2:T(k) b=exp(1i*n*pi*t)/n; a=a+b; end y=2*a/(1j*pi); x=*square(pi*t); figure; p=plot(t,real(y),t,x); axis([-2,2,,]); set(gca, XTick ,-2:1:2) set(gca, YTick ,:: set(gca, XTickLabel ,{ -2 , -1 , 0 , 1 , 2 }) set(gca, YTickLabel ,{ , , 0 , , }) xlabel( 自变量 ) ylabel( 函数值 ) titlemsg=sprintf( 吉布斯现象 N=%d的合成波形 ,T(k)); title(titlemsg) text(0,, \leftarrow 方波函数 , HorizontalAlignment , left ) set(gcf, Color , w ) %hold on End 吉 布 斯 现 象 N=15 的 合 成 波 形 0.8 0.4 值 数 0 函 -0.4 方 波 函 数 -0.8 -2 -1 0 1 2 自 变 量 吉 布 斯 现 象 N=25 的 合 成 波 形 0.8 0.4 值 数 0 函 -0.4 方 波 函 数 -0.8 -2 -1 0 1 2 自 变 量 吉 布 斯 现 象 N=51 的 合 成 波 形 0.8 0.4 值 数 0 函 -0.4 方 波 函 数 -0.8 -2 -1 0 1 2 自 变 量 吉 布 斯 现 象 N=101 的 合 成 波 形 0.8 0.4 值 数 0 函 -0.4 方 波 函 数 -0.8 -2 -1