fftshift【DOC精选】.docVIP

1.FFT的运算量。直接计算DFT共需N2次复数乘法和N(N-1)次复数加法。而FFT仅需计算0.5N㏒2N次复数乘法和N㏒2N次复数加法。由于在计算机上计算乘法所需的时间比计算加法多得多，所以FFT的运算量比DFT要少的少。 2.FFT及其反变换的MATLAB函数。y=fft(x);y=fft(x,N).y=ifft(x);y=ifft(x,N) 3.利用FFT计算线性卷积。步骤： 将x1(n)和x2(n)都延长到N=N1+N2-1。 计算x1(n)的N点FFT。X1=FFT（x1） 计算x2(n)的N点FFT。X2=FFT（x2） 计算Y= X1 .X2 计算Y的反变换y=ifft(Y) 在MATLAB中，可用nextpow2来确定呢最接近的值是2的多少次幂。、 4.利用FFT对信号进行频谱分析。有X（k）群殴相位谱，振幅谱，功率谱。常用的函数有，计算模值的abs函数和计算相角的angle函数。 因为实信号以fs为采样速率的信号在 fs/2 处混叠，所以实信号fft的结果中前半部分对应[0, fs/2],后半部分对应[ -fs/2, 0] 1）实信号fft的结果前半部分对应[0, fs/2]是正频率的结果,后半部分对应[ -fs/2, 0]是负频率的结果。大于fs/2的部分的频谱实际上是实信号的负频率加fs的结果。故要得到正确的结果，只需将视在频率减去fs即可得到频谱对应的真实负频率 2）如果要让实信号fft的结果与[-fs/2, fs/2]对应，则要fft后fftshift一下即可，fftshift的操作是将fft结果以fs/2为中心左右互换 3）如果实信号fft的绘图频率f从[-fs/2, fs/2]，并且没有fftshift，则fft正频谱对应f在[0, fs/2]的结果将混叠到(f - fs/2)的位置; fft负频谱对应f在[-fs/2, 0]的结果混叠到 f + fs - fs/2 的位置，注意这里f为负值，也就是说此种情况下fft负频谱对应的视在频率减去fs/2即可得到频谱对应的真实负频率 二. 复信号情况 1）复信号没有负频率，以fs为采样速率的信号，fft的频谱结果是从[0, fs]的。 2）在 f fs/2 时，对复信号的fft结果进行fftshift会产生频率混叠(将下面的示例2中的频率从f=15改为f=85可以验证f=85的谱线在fftshift后跑到 f = -15 = 85 - fs = 85 - 100的位置了)，所以复信号也一般要求 f = fs/2 3）在对雷达的慢时间维(复信号)进行fft后，由于要用doppler = ((0:LFFT-1)/LFFT - 0.5)*PRF; 计算多普勒频率，所以对该慢时间信号fft后要fftshift下，以便和正确的频率单元相对应。注意多普勒频率fd = PRF/2 时才测的准！ 222222222222222222222222222222222222222 1.双线性变换法的优点是频率响应没有混叠，缺点是数字域频率ω和模拟域频率Ω是非线性关系，这限制了它的应用范围，只有当非线性失真你是允许的或能被补偿的，才能采用双线性变换法，通常，低通，高通，带通和带阻登门滤波器具有分段恒定的频率特性，可以采用遇畸变的方法来补偿频率畸变，因此可以采用双线性变换法。 2.双线性变换法从频率响应出发，直接使数字滤波器的频率响应H(ejw)逼近模拟滤波器的频率响应H（jΩ），进而求出H（z）。双线性变换法通过两次映射，避免了冲击响应不变法的频谱混叠，即，首先将整个s平面压缩变换到s1平面上的一条横带域内，然后通过标准变换关系z=esT将此横带域便唤道整个z平面上，这样就使s平面与z平面建立了一一对应的关系，消除了多值性，从而达到了消除频谱混叠的目的。 3.用双线性变换法设计IIR数字滤波器的步骤。 （1）确定数字了滤波器的通带截止频率wp，阻带截止频率ws，通带最大衰减Rp，组带最小衰减Rs。 （2）利用上式对截止频率ws和wp进行频率预畸变。 （3）计算模拟滤波器的阶数N和频率Ωc，进而求得s平面左半平面上的极点，从而构成模拟滤波器的传递函数，使其满足给定的技术指标。 （4）利用双线性变换公式将模拟滤波器的传递函数H（s）变成数字滤波器的系统函数H(z)。 （5）若需涉及其他形式的滤波器，则根据频率变换公式将数字低通滤波器原型便唤道相应形式的滤波器即可。 4.IIR 数字滤波器的MATLAB实现。 在MATLAB的数字信号处理工具箱中，提供了用双线性变换法设计butterworth低通滤波器的相关函数。Buttord()滤波器阶数选择函数。Buttap()模拟低通滤波器的建立函数。Zp2tf()零极点增益模型到传递函数模型的转换函数。Lp2lp()从低通原