在MATLAB中对气隙磁密进行FFT分解
在Ansoft中对气隙磁密进行FFT分解,后处理做不了,给你一个小程序,放在MATLAB就可以啦,但是要注意点数,文件名要一致哦,我用红色的给你标上啦{:soso_e121:} 已收藏,{:soso_e179:} 有没有WORD版本呀?这个可以用来处理反电势吗?我想要处理反电势,得到他的谐波分析。 为何图片都挂了。。。。。。。。。。。。。 默默无语很多年 发表于 2015-3-19 20:53有没有WORD版本呀?这个可以用来处理反电势吗?我想要处理反电势,得到他的谐波分析。
可以啊,你就把反电势的1001个点保存txt格式,再次运行被{:soso_e121:} 收下了,谢谢无私分享 楼主,我是新手,气隙径向磁密是生成空间FFT,这程序能用么?该数据又该如何导出以及导入到matlab中呢? 把你的点导出到txt格式,在matlab中导入功能,选择这个文件就可以了{:soso_e121:} ansoft不能直接FFT分解成横坐标是几次谐波那种的吗? 他的横坐标是频率,你除以工频不就是谐波次数了么{:soso_e121:} 敦敦 发表于 2015-4-20 07:26
把你的点导出到txt格式,在matlab中导入功能,选择这个文件就可以了
楼主,您好,我是用MATLAB把保存为TXT格式的数据打开,然后在MATLAB中command windows导入你提供的程序,最后提示:Errer:Unexpected MATLAB expression。想请教下,我这方法是否正确。谢谢!{:soso_e183:} 楼主,程序里面算出的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:} 敦敦 发表于 2015-5-5 08:45
程序是粘贴上的,数据是导入进来的,而且要保证两者数据点数一致性哦
楼主,您好,您提供的程序中xn=ss(:,1)的ss是要计算的数据文件名称吧。我现在把ss改成我导入MATLAB中数据文件名,运行后提示:Undefined variable air,意思air文件没定义,请教如何解决,谢谢! 用ansoft直接FFT时候 为啥有负值阿 ??转矩脉动进行fft SS文件得是什么样的排列?2*1001矩阵?还是啥? abcabcabc 发表于 2015-6-8 11:05
SS文件得是什么样的排列?2*1001矩阵?还是啥?
不是矩阵,就是两列数啊,你直接把表格中的复制1001个点,沾到txt文件中就行{:soso_e121:} 本帖最后由 abcabcabc 于 2015-6-8 17:16 编辑
出现错误,我的数据只有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:}
页:
[1]
2