- 积分
- 144
- 回帖
- 0
- 西莫币
-
- 贡献
-
- 威望
-
- 存款
-
- 阅读权限
- 10
- 最后登录
- 1970-1-1
该用户从未签到
|
发表于 2011-6-15 22:45
|
显示全部楼层
来自: 中国黑龙江哈尔滨
我这里有matlab的谐波分解方法,不知道有没有帮助。
N=1001;
n=0:N-1;
t=0:0.1/2/pi:2*pi;
xn=[A1];
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)));
figure(1)
plot(n*2*pi/1001,xn,t,yy1,t,yy2,t,yy3,t,yy4,t,yy5,t,yy6,t,yy7,t,yy8,t,yy9,t,yy10);
yy1=abs(Xk(2))
yy2=abs(Xk(3))
yy3=abs(Xk(4))
yy4=abs(Xk(5))
yy5=abs(Xk(6))
yy6=abs(Xk(7))
yy7=abs(Xk(8))
yy8=abs(Xk(9))
yy9=abs(Xk(10))
yy10=abs(Xk(11)) |
|