吉布斯现象的matlab代码实现_第1页
吉布斯现象的matlab代码实现_第2页
吉布斯现象的matlab代码实现_第3页
吉布斯现象的matlab代码实现_第4页
吉布斯现象的matlab代码实现_第5页
已阅读5页，还剩2页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

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) set(gca,YTick,-0.8:0.4:0.8) set(gca,XTickLabel,-2,-1,0,1,2) s

2、et(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;T = 15,25,51,101;for k = 1:length(T) a=0; for n = 1:2:T(k) b=ex

3、p(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) 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

4、) 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第三题1、 傅里叶级数随着N增加，其合成的波形就就越接近方波2、 随着N增加，波动就越向边缘点压缩靠近3、 当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点第四题1、 一次谐波角频率2rad/s，三次6r

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论