西莫电机圈

 找回密码
 立即注册

QQ登录

只需一步,快速开始

手机号码,快捷登录

手机号码,快捷登录

查看: 14732|回复: 111

[求助] 电磁力密度matlab fft2分解

  [复制链接]

签到天数: 7 天

连续签到: 1 天

[LV.3]偶尔看看II

发表于 2020-4-15 22:58 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?立即注册

x
想要求教matlab的二维fft分解。一直以来大家都说matlab有这个功能,但是我分解出来的数据都对不太上。我把maxwell中的电磁力变成了一个矩阵(第二张图)。在maxwell中因为空间上选择了一个电周期的1001个点,时间上选择一个周期计算了48个点,按照fft2的程序,我用fft2(Fr),之后取模abs(Fr),再将除了第一行和第一列的数据*2*2/1001/48,第一行第一列认为是直流分量,所以直接/1001/48,这样分解后的数据不对。我是48槽8极电机,按理说应该只有偶数次频率以及0、8、16(即阶次的偶数列)阶有电磁力,但是我得到的结果(第三张图)不满足这个理论分析。

                               
登录/注册后可看大图


                               
登录/注册后可看大图


                               
登录/注册后可看大图


                               
登录/注册后可看大图

西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过

签到天数: 9 天

连续签到: 1 天

[LV.3]偶尔看看II

发表于 2020-4-16 08:48 | 显示全部楼层
我也在做这一步,主要就是FFT2的用法
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 25 天

连续签到: 1 天

[LV.4]偶尔看看III

发表于 2020-4-16 09:39 | 显示全部楼层
和你遇到的问题类似,个人认为数据问题来源可能有2个,一个是maxwell输出数据的奇偶个数,一个是FFT2分解的奇偶数的问题。你这时间48个点空间1001个点肯定有问题,48个点分析频率也有点少吧。我的想法是可以先把FFT2分解的问题解决掉。我之前是用的都是偶数个,分解完幅值还是有误差,或许和最后一点是否需要和第一个点重合有关系,您可以试一下看看,能把理想的波形分解正确再考虑maxwell的结果的问题,供参考。
N1=1000;
N2=90;
f=linspace(0,2*pi,N1); %频率序列
v=linspace(0,2*pi,N2); %频率序列
for i=1:N1
    for j=1:N2
        xx(i,j)=cos(1*f(i)-1*v(j))+0.5*cos(1*f(i)+1*v(j))+0.2*cos(3*f(i)-5*v(j))+0.3*cos(10*f(i)-1*v(j))+0.3*cos(0*f(i)-1*v(j))+0.3*cos(0*f(i)-5*v(j))+0.2*cos(1*f(i)+0*v(j))+0.05;
    end
end
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 2 天

连续签到: 1 天

[LV.1]初来乍到

发表于 2020-4-16 13:26 | 显示全部楼层
目测你是用1/4模型得出的电磁力数据,如果是这样,你的空间阶次还要乘一个4
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 7 天

连续签到: 1 天

[LV.3]偶尔看看II

 楼主| 发表于 2020-4-16 15:45 来自手机 | 显示全部楼层
沃丶特 发表于 2020-04-16 13:26
目测你是用1/4模型得出的电磁力数据,如果是这样,你的空间阶次还要乘一个4

是这样,如果是全电机模型,对于48槽8极电机数据出来应该看的就是第0、8、16、…、8n列,如果不乘4,我这个表数据看的就是第0、2、4、…、2n列,这里我为了方便把第一行和第一列认为是第0行第0列
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 2 天

连续签到: 1 天

[LV.1]初来乍到

发表于 2020-4-16 22:30 | 显示全部楼层
那你的数据没问题啊,第一个数据是10的5次方,奇数次的谐波数量级太小可以忽略了
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 7 天

连续签到: 1 天

[LV.3]偶尔看看II

 楼主| 发表于 2020-4-17 01:22 来自手机 | 显示全部楼层
沃丶特 发表于 2020-04-16 22:30
那你的数据没问题啊,第一个数据是10的5次方,奇数次的谐波数量级太小可以忽略了

频率对不上啊,你看第三列的频率就不对,奇数次和偶数次大小都差不多
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 7 天

连续签到: 1 天

[LV.3]偶尔看看II

 楼主| 发表于 2020-4-17 01:25 来自手机 | 显示全部楼层
沃丶特 发表于 2020-04-16 22:30
那你的数据没问题啊,第一个数据是10的5次方,奇数次的谐波数量级太小可以忽略了

而且第一行数据是0f的电磁力,在进行振动分析时可忽略的,一般从第二行开始,第二行是f,第三行是2f,第四行3f,理论来说奇数次频率对应的电磁力应该很小,远小于偶数次才对,但这里是有问题的
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 2 天

连续签到: 1 天

[LV.1]初来乍到

发表于 2020-4-17 08:19 | 显示全部楼层
ninabaobao05 发表于 2020-4-17 01:25
而且第一行数据是0f的电磁力,在进行振动分析时可忽略的,一般从第二行开始,第二行是f,第三行是2f,第 ...

nonono,你看错方向了,很明显列才是时间阶次,而行是空间阶次
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 7 天

连续签到: 1 天

[LV.3]偶尔看看II

 楼主| 发表于 2020-4-17 08:27 来自手机 | 显示全部楼层
沃丶特 发表于 2020-04-17 08:19
nonono,你看错方向了,很明显列才是时间阶次,而行是空间阶次

我取的数据,行是按空间取的1001个点,列是按时间取的48个点,不能分解完了行列颠倒了啊。。。。
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 7 天

连续签到: 1 天

[LV.3]偶尔看看II

 楼主| 发表于 2020-4-17 08:28 来自手机 | 显示全部楼层
沃丶特 发表于 2020-04-17 08:19
nonono,你看错方向了,很明显列才是时间阶次,而行是空间阶次

我也没说错啊,咱们说的是一个意思,我是说在同一列里面,第二行是f,第三行是2f,第四行是3f,以此类推
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

该用户从未签到

发表于 2020-5-12 15:08 | 显示全部楼层
求fft2分解程序,刚弄出电磁力时空结果的小白
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 7 天

连续签到: 1 天

[LV.3]偶尔看看II

 楼主| 发表于 2020-5-12 20:49 来自手机 | 显示全部楼层
xkw123 发表于 2020-05-12 15:08
求fft2分解程序,刚弄出电磁力时空结果的小白

我也没弄出来啊,我的数据处理结果不对
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

该用户从未签到

发表于 2020-5-13 14:10 | 显示全部楼层
ninabaobao05 发表于 2020-5-12 20:49
我也没弄出来啊,我的数据处理结果不对

那楼主加油,我再找找看看
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

该用户从未签到

发表于 2020-5-17 13:47 | 显示全部楼层
楼主,程序可以分享下吗
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 7 天

连续签到: 1 天

[LV.3]偶尔看看II

 楼主| 发表于 2020-5-17 17:27 来自手机 | 显示全部楼层
找我要代码的,一致回复下:公众号“模态空间”华科朱碧华童鞋发表的一篇文章中有详细解释,亲测有效

评分

参与人数 1西莫币 +3 收起 理由
FUDS + 3 赞一个!

查看全部评分

西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 69 天

连续签到: 1 天

[LV.6]常住居民II

发表于 2020-6-8 14:46 | 显示全部楼层
ghugyul 发表于 2020-4-16 09:39
和你遇到的问题类似,个人认为数据问题来源可能有2个,一个是maxwell输出数据的奇偶个数,一个是FFT2分解的 ...

请问您这个代码是什么意思呢?
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

签到天数: 6 天

连续签到: 1 天

[LV.2]偶尔看看I

发表于 2020-8-26 14:45 | 显示全部楼层
亲,我想要代码  谢谢
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

该用户从未签到

发表于 2020-9-15 22:47 来自手机 | 显示全部楼层
是不是没有把零频分量移到中间
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

该用户从未签到

发表于 2020-9-15 22:50 来自手机 | 显示全部楼层
我也是看那篇文章改的,里面有个问题请教一下。幅值还原对于直流分量,是第一行所有数据和第一列所有数据,还是第一行第一列那一个数据。文章里是B[0,0]/=2,看题主的程序好像是第一种说法。
西莫电机论坛微信公众平台正式上线!★详情请点击★ 西莫电机论坛会员交流专用群欢迎您西莫电机论坛加群请注明论坛用户名及所从事专业,否则不予通过
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|小黑屋|手机版|Archiver|西莫电机圈 ( 浙ICP备10025899号-3 浙公网安备:33028202000436号

GMT+8, 2024-7-22 01:07 , Processed in 0.102889 second(s), 25 queries .

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表