- 积分
- 415
- 回帖
- 0
- 西莫币
-
- 贡献
-
- 威望
-
- 存款
-
- 阅读权限
- 20
- 最后登录
- 1970-1-1
该用户从未签到
|
楼主 |
发表于 2017-5-17 22:15
|
显示全部楼层
来自: 中国江苏无锡
现将自己做的展示部分内容
xx=rand(500,1001);yy=rand(250,500);
for i=1:500
for j=1:1001
xx(i,j)=rmf1000(1001*(i-1)+j);
end
end
figure;
surf(xx);
y=abs(fft2(xx,500,1001))/500000*4;
z=flipud(y);
% y(1,1)=0;
% y=flipud(y);
figure;
surf(y);
% for i=1:500
% for j=1:1001
% if y(i,j)<1
% y(i,j)=0;
% end
% end
% end
for i=1:250
for j=1:500
yy(i,j)=y(i+1,j+1);
end
end
figure;
surf(yy);
% axis([0 500 0 50 1000 5]);
% set(gca,'YDir','reverse');
rmf1000是直接导入三维CSV文件
xx是电磁力空间和时间分布数组,yy是经过fft2变化后处理得到的最终结果
MAXWELL导出的电磁力csv文件,就是电磁力那一栏,是一个一维数组,我处理下,变成FFT2可以使用的数组
MATLAB里的FFT函数,对二维波形处理得到的图像,去掉直流分量是对称的,因此用fft2处理后,去掉第一行和第一列,取数组的左上四分之一,得到最终的结果,然而左上四分之一是与右下四分之一对称的,右上四分之一和坐下四分之一对称,左上和左下不相同,我不太懂三维FFT的原理,所以就如上处理了,有意见的请指正
|
|