科学网—zhwang554的留言板

IP: 121.193.142.* [32]王兆华 2011-10-18 12:17

novaelska
你在2011-5-30 15:01 给出一个程序，问，我这个程序怎么修改才能画出书中图2.10那样的程序呢。
n=-N+1:N+1 改为n=1:2*N-1 即能画出书中图2.10。
书中图2.10是2003年通信学报上“全相位FFT频谱分析”的插图，当时apfft的水平相位特性还没有发现, 所以n=1:2*N-1。
你用n=-N+1:N+1 也可画出类似图形，偶对称，而且看起耒更好

IP: 173.70.67.* [31]王兆华 2011-6-17 21:04

是的

IP: 58.218.15.* [30]novaelska 2011-6-17 16:23

老师您好，我想请问。。传统FFT只考虑一种分割情况而产生误差，这种误差在频谱上也会表现出频谱泄漏吗？？

IP: 114.213.29.* [29]novaelska 2011-6-8 21:06

谢谢您~~

IP: 173.70.67.* [28]王兆华 2011-6-8 18:17

是的
全相位无窗---传统不加窗FFT平方
全相位双窗---传统加窗FFT平方
全相位单窗---传统加窗FFTx传统不加窗FFT
全相位窗1和窗2卷积---传统加窗1FFTx传统加窗2FFT

IP: 114.213.24.* [27]novaelska 2011-6-8 18:04

这种平方关系是全相位无窗，和传统不加窗FFT之间的关系吗？

IP: 173.70.67.* [26]王兆华 2011-6-8 18:02

exp信号, 信号振幅A=1, 是平方关系
cos信号, 信号振幅 A=2, 是平方关系

IP: 114.213.24.* [25]novaelska 2011-6-8 17:54

老师，您好，我想请问您，若是余弦信号，全相位谱分析后，它的振幅谱的幅值和传统FFT分析后的幅值也是平方的关系吗？我用书上第一章1.1.1中的余弦信号作全相位FFT谱分析，它的振幅谱与传统FFT振幅谱的幅值不是平方关系。。。

IP: 114.213.29.* [24]novaelska 2011-6-5 22:55

好的。。谢谢老师！！！

IP: 173.70.67.* [23]王兆华 2011-6-5 18:11

参见下网页
http://www.chinavib.com/thread-102988-1-1.html
5楼
songzy41发表于 2011-5-26 06:05 |的帖子

guodongliang 发表于 2011-5-25 11:41
感谢大师！
还有一个疑问：
这个模拟信号只有3个频率分量，并无噪声，幅度频谱是只有3个点有频率的，也几乎无频谱泄露，为何相频曲线在所有频点处都有相位呢？是计算误差？
这是由于计算误差造成的。首先计算机数字的表示都是用有限字长表示的，例如pi，我们知着它是一个无限长的无理数，但在计算机中只能用有限长度表示。pi=3.1415926...，而在计算机中只能表示为3.14159，后面的0.0000026...就是有限字长的截断误差了。由pi计算出的三角函数，例如sqrt(2)/2，sqrt(3)/2等等原本也应为无理数，现在也按有限字长表示了，就是数的表示中就有误差，用这些带有误差的数参与运算，使误差累计。当然这些误差很小，信号的幅值计算中都看不出来。但在相位计算出就不一样了，例如FFT后X(i)=1.2e-9+j2.4e-10，它的幅值(模值)就很小，但相角angle(1.2e-9+j*2.4e-10)=0.1974就不小了，所以相频曲线在-pi至pi之间振荡是正常的。

IP: 114.213.26.* [22]novaelska 2011-6-5 12:57

老师您帮我一下吧。。我的老师也不太懂全相位。。。

IP: 114.213.29.* [21]novaelska 2011-6-5 10:03

老师，您好，我编的程序，如果是整频点，相位谱就不是一条线，但如果是有一点偏离，比如偏离0.001就是一条线。。这是怎么回事呢？？

IP: 114.213.26.* [20]novaelska 2011-6-4 17:48

老师您好，我还想再问一下，全相位FFT无论有没有频偏，相位谱都是初始相位吗？？？

我的回复(2011-6-4 19:33)：是的

IP: 114.213.31.* [19]novaelska 2011-6-4 14:18

王老师您好，我想请问您，FFT谱分析时候，假如频率是29，相位是100，matlab图相位谱显示是30处为100，为什么apFFT还是在29处为100？

IP: 173.70.67.* [18]王兆华 2011-5-30 23:40

novaelska:
as N=3;q=2*pi*(b-k);
exp(j0q/3)+exp(j1q/3)+exp(j2q/3)=[1-exp(jq)] / [1-exp(jq/3)]
左边分子 X 右边分母 = 右边分子
exp(j0q/3)+exp(j1q/3)+exp(j2q/3)- exp(j1q/3)-exp(j2q/3)-exp(j3q/3)= 1-exp(jq)

我的计算机Matlab坏了, 程序一时处理不了, 抱歉!

IP: 114.213.26.* [17]novaelska 2011-5-30 21:16

老师您好，我还有一个问题想问，就是书上第三章第二节的公式（3-3）中第二步到第三步的转换公式是什么？我看了您的视频，但是视频不是很清楚，所以那个式子看的不清楚，您能不能写一下那个公式。。谢谢您了！！

IP: 114.213.31.* [16]novaelska 2011-5-30 15:01

N=32;
fs=32;
n=-N+1:N+1;
f0=1;
x=cos(2*pi*f0/fs*n);
w=hanning(N)';
w1=conv(w,w);
w2=w1/sum(w1);
x1=x(1:2*N-1);
x2=x1.*w2;
y=x2(N:end)+[0 x2(1:N-1)];
nn=0:N-1;
plot(nn,y);
hold off;
老师，您好，我想请问，我这个程序怎么修改才能画出书中图2.10那样的程序呢。。

IP: 114.213.28.* [15]novaelska 2011-5-28 12:19

谢谢老师了~

IP: 173.70.67.* [14]王兆华 2011-5-28 00:59

Novaelska: 原公式(3)有错，改了, 公式(3) 说明公式(1)内部关系,公式(7) 中用.
双窗和单窗先加窗, 周期延伸,再截断,截断都用矩形窗

IP: 114.213.31.* [13]novaelska 2011-5-27 18:15

老师，您好，再打扰您一下。。我刚才看到了公式（3），我还是有点疑惑，为什么式子中没有幅度和初相位了？

zhwang554的个人博客分享 http://blog.sciencenet.cn/u/zhwang554

留言板

王兆华