电磁力密度matlab fft2分解

doufaq

楼主你说的那篇文章搜不到啊

zengxiaodong

二维的FFT用Originpro来处理，非常方便，只需要把数据存成矩阵就可以了。下面提供一个实例（含数据矩阵 3槽4极矩阵.rar (321.59 KB, 下载次数: 179)

2020-9-28 16:09 上传

点击文件名下载附件

）








补充内容 (2020-9-29 23:02):
1、数据的每一列代表一个时刻，用某列的数据绘制出来的波形就是该时刻径向力波的空间分布波形·；
2、第一列和最后一列，几乎完全相同（差一列），也就是时间上形成完整的周期；
3、每列空间波形可连续显示形成动画。

zengxiaodong

第一步，import单个ASCII文件（矩阵）；

第二步，Plot出三维彩色云图；

第三步，Analysis-----signal processing----FFT-----2DFFT，勾选直流偏置到中间，output可以仅勾选幅度谱；

第四步，根据第三步得到的矩阵，再次绘制三维Bar图；

第五步，修改坐标标尺，放大感兴趣的区域，进行数据判读，可以用数据读取工具DATA Reader点选，然后用左右、上下箭头移动光标到各个点，读取数据。

zengxiaodong

上贴步骤的部分截屏，供学习参考

zengxiaodong

为了规范化处理，在生成数据矩阵时，可以统一格式，X轴为空间（行），Y轴为时间（列），而且数据点数也统一，空间1024个点，时间64个点。这样的话，在后续数据处理时可以非常方便，坐标轴的数据可以统一。

Panamera

520蝴蝶花 发表于 2020-9-25 11:44
那个帖子用Python写的，我看过代码，也用python编译出来了结果，不过那个人应该没有真正搞懂径向力波

仅 ...

谢谢分享，先下载一下看看，确实，朱碧华老师的程序我不知道在matlab中怎么用

Panamera

请问楼主，你是在matlab中直接运行的吗，怎么我的运行不了

zengxiaodong

520蝴蝶花 发表于 2020-9-29 11:26
恕我直言，这位曾老师的FFT2是没有意义的，仅仅用个软件做出了柱状图而已，各位童鞋不要被误导了哦

不过 ...

我可不是FFT2话题的发起者哦，只不过我看到太多的帖子讨论FFT2，而且没有一个帖子有深度！多数是炫耀一下而已，似乎三维图一上来就显得高大上，其实是玩深沉罢了。

这么简单的问题，何须编程来搞？现成的数据处理工具多了去了，何不拿来就用呢？

另外一个问题，FFT2的结果如何解读？这个更加是没人讲出个所以然......

zengxiaodong

别说是高深莫测的FFT2，连FFT我也不是很熟悉，不过笨人有笨办法，那就是我想到的积分法，用动画直观显示空间0次、空间1次力波的变动情况，这类动画我发过很多，为了不占带宽，这里就不举例说明了。

yinhang

1、如果是用的maxwell导出的气隙磁密，不建议直接画一个circle导出磁密，建议画一个359.9°的圆弧，两种方法导出的磁密在fft2的时候会有一些区别；
2、fft2这东西并不复杂，正确的结果力波分布特征很明显。（附件是个24s20p的电机力波分布）

夕林0529

yinhang 发表于 2020-9-29 15:07
1、如果是用的maxwell导出的气隙磁密，不建议直接画一个circle导出磁密，建议画一个359.9°的圆弧，两种方 ...

楼上说FFt2可行的，基本理论水平都没到家。   不建议采取他们的意见，真的有这么简单的话，早就有大神出来讲这个函数了，正式因为大神们都不认可这个函数，才使得这些小喽喽出来瞎吹牛，搞得好像很简单一样。   你们见过论坛里功底好的  说使用这个函数就可以了？

yinhang

夕林0529 发表于 2020-9-29 15:26
楼上说FFt2可行的，基本理论水平都没到家。   不建议采取他们的意见，真的有这么简单的话，早就有大神出 ...

多看看文献不会让你学会谦虚吗？看个help文档很难吗？调用个函数不简单吗？认真思考做不到吗？好好说话很费劲吗？学会尊重不会吗？

zengxiaodong

yinhang 发表于 2020-9-29 15:07
1、如果是用的maxwell导出的气隙磁密，不建议直接画一个circle导出磁密，建议画一个359.9°的圆弧，两种方 ...

你这个不就是6槽7极电机(6槽5极)最简电机吗？我有现成的例子，可以对比参考：2f频率时，空间1、4、7、13次正转波，空间5次反转波。（对于24槽20极，前述空间次数乘以4，因为4个最简电机的缘故）

夕林0529

yinhang 发表于 2020-9-29 19:49
多看看文献不会让你学会谦虚吗？看个help文档很难吗？调用个函数不简单吗？认真思考做不到吗？好好说话很 ...

自己学的不到家，还在这里摆谱，哪来的勇气说这些。  现在态度上和我一路货色的人，有啥资格说着说那。 你写出来 一定正确吗？你都没证明，错误的东西还不让人说。我说话就这个味道，所以不存在态度差不差的问题，只是你不习惯。而你，说话两个味道，才是态度不好呢？ 现在到了猴子称大王的年代了！！

zengxiaodong

国庆假期，不妨花点时间把2D FFT搞个透彻，省得云里雾里好像有多神秘似的。

基本的思路如下：

1、三种波：行波、驻波、乱跳波；

2、亲自部署、亲自指挥，分别构造3种波的实际数据，并进行3D云图的绘制；

3、进行2D FFT的分析，得到结果，并绘制出3D柱状图；

4、对FFT2得到的柱状图进行判读，看看与亲自生成的数据预设的特征有哪些对应关系；

5、在上述的基础之上，进行电机径向力波的二维FFT分析，判读出力波的空间次数，以及对应的时间频率等特征。

MADANO

坐等更新

zengxiaodong

先来定个规范：

X轴代表空间，没有单位，采样点编号1到1024
Y轴代表时间，没有单位，采样点编号1到64

根据香农采样定理，每个周期至少采样2个点，才不会出现混叠现象，因此上述采样点数可以确保10阶以内的谐波都有远超2个采样点。

zengxiaodong

先来分析标准的行波

其方程为：
cos(x/1024*20*pi+y/64*12*pi)

上式中pi是圆周率

可以看出，空间上有10个周期，时间上有6个周期

打开OriginPro，新建一个矩阵Matrix



输入列数64，行数1024，注意：此处X，Y的定义与前述相反

zengxiaodong

生成了一个空矩阵，有1024行，64列

输入公式，生成矩阵各个单元的具体数值



复制、黏贴上面的公式，注意：变量改成i和j

zengxiaodong

生成了矩阵，为了后续处理方便，把矩阵转换成工作表



此时，你可以任选一列绘制该列对应时刻的空间波形