【发布时间】:2018-01-05 00:27:06
【问题描述】:
我正在尝试一个简单的(快速傅立叶变换)fft 并将其转换回来,但它不起作用。 我需要从这个开始,以便继续为每个组件添加阶段。你能看看我哪里错了吗?
fun_cos=@(t) cos(1.5e12*t)
nttf=2^17;
t=linspace(-3*t_signal_pulse/2,3*t_signal_pulse/2,nttf);
dni_ni=(1/(t(2)-t(1)));
ni=-dni_ni/2:dni_ni/(nttf):dni_ni/2-dni_ni/(nttf);
w=ni.*2*pi;
figure(1)
plot(t,fun_cos(t))
FFt_cos=fftshift(fft(fun_cos(t),nttf))/length(t);
figure(2);
plot(w,abs(FFt_cos))
fft_back=ifft(ifftshift(FFt_cos));
figure(1)
hold on
plot(t,abs(fft_back),'.r')
Freq domain. you can see here two freq even though I only need one
另外,如果我想分别将相位添加到时域和频域(请注意,我只知道频域相位的一侧,而不是两者,所以还不知道如何进行此操作)
【问题讨论】:
-
请提供minimal reproducible example 并且包括结果可能有用。
-
刚刚添加了一张带有结果的照片。蓝色是原始功能。红色是我认为会相同的那个,但正如你所知,它不是