《吉布斯现象的matlab代码实现》由会员分享，可在线阅读，更多相关《吉布斯现象的matlab代码实现（7页珍藏版）》请在金锄头文库上搜索。
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;专%业文*档# end专%业文档& # y=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) s
2、et(gca,YTick,-0.8:0.4:0.8) set(gca,XTickLabel,-2,-1,0,1,2) set(gca,YTickLabel,-0.8,-0.4,0,0.4,0.8) xlabel(自变量) ylabel(函数值) titlemsg=sprintf(吉布斯现象N=%d的合成波形,T(k); title(titlemsg) text(0,-0.5,leftarrow 方波函数,HorizontalAlignment,left) set(gcf,Color,w) %hold onEnd第二题代码clear;close all;clc;t = -2:0.01:2;专%业
3、&文#档 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=0.5*square(pi*t); %figure; if (k=1) p=plot(t,real(y),c,t,x,k); elseif (k=2)专%业&文*#档 p=plot(t,real(y),m,t,x,k); elseif (k=3)专%业文档 # p=plot(t,real(y),r,t,x,k);专*#业文&档 elseif (k
4、=4) p=plot(t,real(y),b,t,x,k); end 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,1,2) set(gca,YTickLabel,-0.8,-0.4,0,0.4,0.8) xlabel(自变量) ylabel(函数值) titlemsg=sprintf(吉布斯现象N=%d的合成波形,T(k); title(titlemsg)专#业文&档 % text(0,-0.5,leftarrow
5、方波函数,HorizontalAlignment,left)专*业%文档 set(gcf,Color,w) hold onEnd专业&文档 *专业文&*#%档 专业#文*档% 专*业文档&#% 专业#文档% &第三题专业文%档 #1、 傅里叶级数随着N增加，其合成的波形就就越接近方波2、 随着N增加，波动就越向边缘点压缩靠近3、 当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点专#*业文档 &专#业文*档 %专*&业文档# 第四题1、 一次谐波角频率2rad/s，三次6rad/s，五次10rad/s2、 一次谐波的系数是4/,三次谐波的系数是4/,五次谐波的系数是4/3、 N次谐波的系数是4/ 最新 精品 Word 欢迎下载 可修改
