从欧拉公式看三角波的单边谱与双边谱

56 篇文章 20 订阅

1、在有限区间上,凡满足狄利赫里条件的周期函数都可以展开为傅里叶级数。其三角函数形式的频谱为单边谱复指数函数形式的频谱为双边谱

在这里插入图片描述
参考:https://wenku.baidu.com/view/af2ea95b7dd184254b35eefdc8d376eeaeaa17de.html

实数系列的FFT后的X(k)是共轭对称的,即实部偶对称,虚部奇对称,matlab fft()得到的是双边谱,前部分是正频率部分,后部分是负频率部分。

2、看看fft()后的频率

x=[1,1.8,1.8,1];
n=length(x);
fs=8;
dt=1/fs;
ti=0:dt:(n-1)*dt;
nfft=8;
y=fft(x,nfft);
df=fs/nfft;
fi=0:df:(nfft-1)*df;

subplot(312)
plot(fi,real(y))
title('real')
xlabel('f(Hz)')
subplot(313)
plot(fi,imag(y))
title('imag')
xlabel('f(Hz)')
x_ifft=ifft(y);
x_ifft_output=(x_ifft(1:n));
subplot(311)
plot(ti,x_ifft_output)
title('x(t)')
xlabel('t(s)')

figure
%Determin Nyquist frequency:
Nyq=fs/2;
%Save frequency array:
fi_array=(-Nyq:df:Nyq-df);%fftshift()后的频率 (0点右移Nyq个点!!)
y_shift=fftshift(y);
subplot(312)
plot(fi_array,real(y_shift))
title('real')
xlabel('f(Hz)')
subplot(313)
plot(fi_array,imag(y_shift))
title('imag')
xlabel('f(Hz)')
subplot(311)
plot(ti,x)
title('x(t)')
xlabel('t(s)')

在这里插入图片描述
可见,matlab fft()后的频率为:[0:df:Nqy-df,-Nqy:df:-df]是其正负频率。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值