在MATLAB中对气隙磁密进行FFT分解

敦敦

在Ansoft中对气隙磁密进行FFT分解，后处理做不了，给你一个小程序，放在MATLAB就可以啦，但是要注意点数，文件名要一致哦，我用红色的给你标上啦{:soso_e121:}

销声匿迹

已收藏，{:soso_e179:}

默默无语很多年

有没有WORD版本呀？这个可以用来处理反电势吗？我想要处理反电势，得到他的谐波分析。

寒江雪

为何图片都挂了。。。。。。。。。。。。。

敦敦

默默无语很多年 发表于 2015-3-19 20:53
有没有WORD版本呀？这个可以用来处理反电势吗？我想要处理反电势，得到他的谐波分析。

可以啊，你就把反电势的1001个点保存txt格式，再次运行被{:soso_e121:}

dddjjjyyy888

收下了，谢谢无私分享

dzc3344

楼主，我是新手，气隙径向磁密是生成空间FFT，这程序能用么？该数据又该如何导出以及导入到matlab中呢？

敦敦

把你的点导出到txt格式，在matlab中导入功能，选择这个文件就可以了{:soso_e121:}

风雨无阻

ansoft不能直接FFT分解成横坐标是几次谐波那种的吗？

敦敦

他的横坐标是频率，你除以工频不就是谐波次数了么{:soso_e121:}

dddjjjyyy888

敦敦 发表于 2015-4-20 07:26
把你的点导出到txt格式，在matlab中导入功能，选择这个文件就可以了

楼主，您好，我是用MATLAB把保存为TXT格式的数据打开，然后在MATLAB中command windows导入你提供的程序，最后提示：Errer:Unexpected MATLAB expression。想请教下，我这方法是否正确。谢谢！{:soso_e183:}  

yuanyehit

楼主，程序里面算出的Xk是什么呢？还有从yy1到yy19是波形的各次谐波嘛？

敦敦

dddjjjyyy888 发表于 2015-5-4 15:39
楼主，您好，我是用MATLAB把保存为TXT格式的数据打开，然后在MATLAB中command windows导入你提供的程序， ...

程序是粘贴上的，数据是导入进来的，而且要保证两者数据点数一致性哦{:soso_e121:}

敦敦

yuanyehit 发表于 2015-5-4 21:52
楼主，程序里面算出的Xk是什么呢？还有从yy1到yy19是波形的各次谐波嘛？

1、是每次谐波的幅值
2、对，就是谐波次数{:soso_e121:}

dddjjjyyy888

敦敦 发表于 2015-5-5 08:45
程序是粘贴上的，数据是导入进来的，而且要保证两者数据点数一致性哦

楼主，您好，您提供的程序中xn=ss(:，1）的ss是要计算的数据文件名称吧。我现在把ss改成我导入MATLAB中数据文件名，运行后提示：Undefined variable air，意思air文件没定义，请教如何解决，谢谢！

风雨无阻

用ansoft直接FFT时候 为啥有负值阿 ？？转矩脉动进行fft

abcabcabc

SS文件得是什么样的排列？2*1001矩阵？还是啥？

敦敦

abcabcabc 发表于 2015-6-8 11:05
SS文件得是什么样的排列？2*1001矩阵？还是啥？

不是矩阵，就是两列数啊，你直接把表格中的复制1001个点，沾到txt文件中就行{:soso_e121:}

abcabcabc

出现错误，我的数据只有201个，程序为
N=201;
n=0:N-1;
t=0:0.1/2/pi:2*pi;
xn=ss(:,1);
k=0:N/2;
WN=exp(j*2*pi/N);
kn=k*n;
WNnk=WN^kn;
Xk=WNnk*xn*2/N;
stem(k,abs(Xk))
abs(Xk);
yy1=abs(Xk(2))*sin(t+pi/2-angle(Xk(2)));
yy2=abs(Xk(3))*sin(2*t+pi/2-angle(Xk(3)));
yy3=abs(Xk(4))*sin(3*t+pi/2-angle(Xk(4)));
yy4=abs(Xk(5))*sin(4*t+pi/2-angle(Xk(5)));
yy5=abs(Xk(6))*sin(5*t+pi/2-angle(Xk(6)));
yy6=abs(Xk(7))*sin(6*t+pi/2-angle(Xk(7)));
yy7=abs(Xk(8))*sin(7*t+pi/2-angle(Xk(8)));
yy8=abs(Xk(9))*sin(8*t+pi/2-angle(Xk(9)));
yy9=abs(Xk(10))*sin(9*t+pi/2-angle(Xk(10)));
yy10=abs(Xk(11))*sin(10*t+pi/2-angle(Xk(11)));
yy11=abs(Xk(12))*sin(11*t+pi/2-angle(Xk(12)));
yy12=abs(Xk(13))*sin(12*t+pi/2-angle(Xk(13)));
yy13=abs(Xk(14))*sin(13*t+pi/2-angle(Xk(14)));
yy14=abs(Xk(15))*sin(14*t+pi/2-angle(Xk(15)));
yy15=abs(Xk(16))*sin(15*t+pi/2-angle(Xk(16)));
yy16=abs(Xk(17))*sin(16*t+pi/2-angle(Xk(17)));
yy17=abs(Xk(18))*sin(17*t+pi/2-angle(Xk(18)));
yy18=abs(Xk(19))*sin(18*t+pi/2-angle(Xk(19)));
yy19=abs(Xk(20))*sin(19*t+pi/2-angle(Xk(20)));
a1=abs(Xk(2))
a2=abs(Xk(3))
a3=abs(Xk(4))
a4=abs(Xk(5))
a5=abs(Xk(6))
a6=abs(Xk(7))
a7=abs(Xk(8))
a8=abs(Xk(9))
a9=abs(Xk(10))
a10=abs(Xk(11))
a11=abs(Xk(12))
a12=abs(Xk(13))
a13=abs(Xk(14))
a14=abs(Xk(15))
a15=abs(Xk(16))
a16=abs(Xk(17))
a17=abs(Xk(18))
a18=abs(Xk(19))
a19=abs(Xk(20))
pi/2-angle(Xk(2))
figure(1)
plot(n*2*pi/201,xn,t,yy1,t,yy3,t,yy11,t,yy13,t,yy17,t,yy19);
运行之后出现
??? Error using ==> mtimes
Inner matrix dimensions must agree.

敦敦

你的txt应该命名为ss{:soso_e121:}