MATLAB,恢复具有余弦干扰且已进行调幅调制的信号

MATLAB,恢复具有余弦干扰且已进行调幅调制的信号,matlab,filter,signal-processing,Matlab,Filter,Signal Processing,我有这个wav信号,其中添加了一个余弦干扰a cos(2*Pi*f0*n*Tc+相位),并且它被调幅到了相同的f0频率Tc是采样周期 我有两个滤波器的滤波器系数: 1) 陷波滤波器 2) 低通滤波器 第一步: 我计算了输入信号的fft,它在7500 Hz处有一个大峰值,这就是我的f0then 滤波器1的中心位置为16KHz,因此我需要将其移动,使其中心位置为7500 Hz。为此,我使用了以下事实:exp(j*w0*n)x(n)->x(j(w-w0))。因此,我将滤波器系数乘以exp(j*2*Pi

我有这个
wav
信号,其中添加了一个余弦干扰
a cos(2*Pi*f0*n*Tc+相位)
,并且它被调幅到了相同的
f0
频率
Tc
是采样周期

我有两个滤波器的滤波器系数:

1) 陷波滤波器

2) 低通滤波器

第一步:

我计算了输入信号的fft,它在7500 Hz处有一个大峰值,这就是我的
f0
then

滤波器1的中心位置为
16KHz
,因此我需要将其移动,使其中心位置为
7500 Hz
。为此,我使用了以下事实:
exp(j*w0*n)x(n)->x(j(w-w0))
。因此,我将滤波器系数乘以exp(j*2*Pi*8500*n)

我使用
filtfilt
函数过滤输入

我将
filtfilt
的输出乘以
cos(2*pi*n*f0*Tc)

我用低通滤波器过滤了最后一个信号


结果:输出仍然失真!我的错误在哪里?

把这个问题转移到哪里去呢?我相信应该是
exp(j*2*Pi*8500/Fs*n)
其中
Fs
是采样频率(或者
exp(j*2*Pi*8500*n*Tc)
,如果你愿意的话。这是dsp.stackexchange.com上的等效帖子。。。