吉布斯现象Matlab实现
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
clear;closeall;clc;
t=-pi:0.01:pi;
T=[5;7;21;51;101];
fori=1:length(T)
a=0;
for n=1:2:T(i)
b=sin(n*t+n*pi/2)/n;
a=a+b;
end
y=4*a/pi;
x=square(t+pi/2);
figure;
p=plot(t,y,t,x);
set(gca,'XTick',-pi:pi/2:pi)
set(gca,'XTickLabel',{'-pi','-pi/2','0','pi/2','pi'})
xlabel('自变量')
ylabel('函数值')
titlemsg=sprintf('吉布斯现象N=%d的合成波形',T(i));
title(titlemsg)
text(-pi/2+0.05,sin(-pi/4),'\leftarrow方波函数','HorizontalAlignment','left') set(gcf,'Color','w')
end