《吉布斯现象的matlab代码实现》由会员分享，可在线阅读，更多相关《吉布斯现象的matlab代码实现（6页珍藏版）》请在金锄头文库上搜索。
1、第一题 :代码 clear;close all ; clc;t = -2:0、 01: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;endy=2*a/(1j*pi);x=0、 5*square(pi*t);figure;p=plot(t,real(y),t,x);axis(-2,2,-0、 8,0、 8);set(gca,XTick,-2:1:2)set(gca,YTick,-0、 8:0 、 4:0 、 8)set(gca,XTickLabel, -2, -1, 0
2、, 1 , 2 )set(gca,YTickLabel, -0 、 8 , -0 、 4 , 0 , 0 、 4, 0 、 8 )xlabel( 自变量)ylabel( 函数值)吉布斯现象N=%d的合成波形，T(k);titlemsg=sprintf(title(titlemsg)text(0,-0、 5, leftarrow 方波函数 , HorizontalAlignment, leftset(gcf, Color , w )%hold onEnd0.8-0.8吉布斯现象N=15的合成波形0.40-0.4-101自变量-2吉布斯现象N=25的合成波形0.8 值 数 函-0.4方波函数-0.
3、8 ,11-2-1012自变量0.8-0.8吉布斯现象N=51的合成波形0.40-0.4-101自变量-2第二题代码clear;close allclc;t = -2:0、 01: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;endy=2*a/(1j*pi);x=0、 5*square(pi*t);%figure;if (k=1) p=plot(t,real(y), elseif (k=2) p=plot(t,real(y), elseif(k=3)p=plot(t,r
4、eal(y), elseif(k=4)p=plot(t,real(y), endc,t,x,k);m,t,x,k);r,t,x,k);b,t,x,k);axis(-2,2,-0、 8,0 、 8);set(gca, XTick ,-2:1:2)set(gca,set(gca,YTick ,-0 、 8:0 、 4:0 、 8)XTickLabel , -2, -1, 0 , 1 , 2 )set(gca, YTickLabel , -0 、 8 , -0 、 4 , 0 , 0 、 4 , 0 、 8 )xlabel( 自变量 )ylabel( 函数值 )titlemsg=sprintf(吉布
5、斯现象N=%d的合成波形，T(k);title(titlemsg)text(0,-0、 5, leftarrow 方波函数 , HorizontalAlignment, left )set(gcf, Color , w )hold onEnd第三题自变量1、傅里叶级数随着N增加淇合成的波形就就越接近方波2、随着N增加,波动就越向边缘点压缩靠近3、当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点第四题一次谐波角频率 2rad/s,三次6rad/s,五次10rad/s一次谐波的系数就是4/%三次谐波的系数就是4/3及五次谐波的系数就是4所 N次谐波的系数就是4/nHl1、2、3、
