|
![](static/image/common/ico_lz.png)
楼主 |
发表于 2007-7-15 16:51
|
显示全部楼层
回复 #6 songzy41 的帖子
clear all;clc
nfft = 1024;
F = zeros(1,nfft+1);
fs = 10000;
fn = [0:nfft]/nfft * fs;
fh = 2000;
fl = 1000;
index = find(fn>=fl & fn<=fh);
F(index(1:20)) = 1*exp(i*pi/4);
F(index(21:40)) = 2*exp(-i*pi/4);
F(index(41:60)) = 2.5*exp(i*pi/2);
F(index(61:80)) = 3*exp(i*pi*3/4);
F(index(81:100)) = 1*exp(-i*pi/4);
F(end:-1:end-(nfft/2-2)) = conj(F(2:nfft/2));
y = real(ifft(F));
psd(y,1024,fs,[])
对于我假定的H(w),经过ifft变换以后,再fft以后已经不是H(w)了.:@( 帮我看看如何修改. |
|