- 在线时间
- 9 小时
- 最后登录
- 2013-11-8
- 注册时间
- 2013-7-2
- 听众数
- 9
- 收听数
- 0
- 能力
- 0 分
- 体力
- 295 点
- 威望
- 0 点
- 阅读权限
- 150
- 积分
- 130
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 86
- 主题
- 54
- 精华
- 0
- 分享
- 0
- 好友
- 3
TA的每日心情 | 衰 2013-11-8 15:00 |
---|
签到天数: 15 天 [LV.4]偶尔看看III
- 自我介绍
- 我笑一笑从容
 |
t=0:0.1:5;
y=cos(2*pi*2*t)+sin(2*pi*2*t)*j
y=fftshift(fft(y,64));
f=linspace(-5,5,64);
figure(1)
subplot(2,1,1)
plot(f,abs(y))
%%重心法修订FFT主瓣重心%%
[m,n]=find(abs(y)==max(max(abs(y))));
z=abs(y(n+1))-abs(y(n-1));
if z>=0
n1=n+abs(y(n+1))/(abs(y(n))+abs(y(n+1)));
end
if z<0
n1=n+abs(y(n-1))*(-1)/(abs(y(n))+abs(y(n-1)));
end
df=-5+(n1-1)*10/63;
f(n)=df;%修订后主瓣重心
subplot(2,1,2)
plot(f,abs(y))
|
zan
|