matlab傅里叶正逆变换详细说明+图例(共4页)_第1页
matlab傅里叶正逆变换详细说明+图例(共4页)_第2页
matlab傅里叶正逆变换详细说明+图例(共4页)_第3页
matlab傅里叶正逆变换详细说明+图例(共4页)_第4页
全文预览已结束

下载本文档

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

文档简介

1、精选优质文档-倾情为你奉上MATLAB作业方法一：N=200;dt=0.001;n=1:200;x=3*sin(2*pi*10*n*dt)+3*sin(2*pi*30*n*dt)+sin(2*pi*40*n*dt)+sin(2*pi*50*n*dt)+sin(2*pi*60*n*dt)+6*sin(2*pi*80*n*dt);%建立时间序列X=zeros(1,200); %给X一个预先的内存空间，提高运行速度figure(1),plot(n,x); %画出时间系列图像%傅里叶变换for k=1:200 for n=1:200 X(k)=X(k)+x(n)*exp(-i*2*pi*n*k/N);

2、 endendf=abs(X); %对傅里叶变换后的图像取正数部分figure(2),plot(f); %画出傅里叶变换后的图像%滤波H=ones(200);H(8:14)=0;H(186:193)=0; %建立40、50、60Hz的时间序列所在空间域数值为0的一维矩阵for k=1:200 Y(k)=X(k)*H(k);endk=1:200;figure(3);plot(k,abs(Y); %画出滤波之后的时间序列%傅里叶逆变换y=zeros(1,200)% 给y一个预先的内存空间，提高运行速度 for n=1:200 for k=1:200 y(n)=Y(k)*exp(i*2*pi*n*k

3、/N)+y(n); end y(n)=y(n)/N;endfigure(4) ;n=1:200;plot(n,y,'r-',n,x,'b-');% 画出原时间序列和滤波后的时间序列，可见滤波之后的时间序列振幅明显变小方法二：N=200;dt=0.001;n=1:200;x=3*sin(2*pi*10*n*dt)+3*sin(2*pi*30*n*dt)+sin(2*pi*40*n*dt)+sin(2*pi*50*n*dt)+sin(2*pi*60*n*dt)+6*sin(2*pi*80*n*dt);X=zeros(1,200);figure(1),plot(n,x);X=fft(x,N);f=abs(X);figure(2),plot(f);H=ones(200);H(8:14)=0;H(186:194)=0;%滤波for k=1:200 Y(k)=X(k)*H(k);endk=1:200;figure(3);plot(k,abs(Y);%傅里叶逆变换y=zeros(1,200);y=ifft(Y,N) ;figure(4) ;n=1:200;plot(n,y,'r-',

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

温馨提示

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

最新文档

评论

0/150

提交评论