匹配滤波器

606次观看(过去30天)
vsee
vsee 2011年4月1日
我如何在Matlab中写一个mtached过滤器?我有滤波器设计和信号处理工具箱。

回答(5)

宏磊陈
宏磊陈 2011年4月1日
如果你有一个信号x,那么匹配滤波器的系数由x的时间逆给出,即x(结束:-1:1)。如果你的信号是复信号,你还需要使用复共轭。然后你就可以把它用作FIR滤波器了。例如,
>> x = ones(10,1);
>> b = x(end:-1:1);
>> y = filter(b,1,x);
5个评论
宏磊陈
宏磊陈 2017年1月12日
是一样的,用 过滤器 函数。唯一的区别是当你提出滤波系数时,你需要加上一个共轭,比如
B = conj(x(end:-1:1))
HTH

登录发表评论。


安德鲁·纽厄尔
安德鲁·纽厄尔 2011年4月1日
看看 雷达系统设计的MATLAB仿真 从文件交换。

Mathuranathan Viswanathan
Mathuranathan Viswanathan 2013年4月13日
你好,你可以试试这里给的模型
问候Mathuranathan

Bhaskar
Bhaskar 2011年4月15日
再详细解释一下洪雷的回答。如果你要设计匹配滤波器的信号是x,那么——
>> b = x(end-1:1);
>> y = filter(b,1,x);
应该已经有x了。
1评论
嗯sss
嗯sss 2012年1月9日
如何在图像上实现相同的过滤器?
如何表示x。

登录发表评论。


妮可Bienert
妮可Bienert 2020年1月7日
可以使用内置的xcorr函数—您只需要通过fs规范化并对输出进行平方。下面是一个理想的例子:
目的:演示如何正确匹配过滤器
BW = 3.84 e6;
fs = 50*BW;%采样率
T = 1 / f;%样本时期
Fc = 330e6;%中心频率
chirpLen = 0.075;%唧唧喳喳的长度
= 3;电压信号幅度%(通常未知)
基于“增大化现实”技术= 2;参考电压信号幅度%(通常未知)
在无噪声的情况下创建信号,并将零填充在两边(零
%padding不需要,因为xcorr会这样做,我只是演示一下
%的信号不需要相同的长度。)
sig =[0(1,装天花板(chirpLen * fs)), A *唧唧声(t, 0, t(结束),BW), 0(1,装天花板(chirpLen * fs)));
%创建参考啁啾
ref_chirp = Ar *唧唧声(t, 0, t(结束),BW);
t = [0, t:(长度(ref_chirp) 1) * t];
归一化参考啁啾:参考啁啾的能量需要为1
%,这样它就不会偏向匹配过滤器的输出。过滤器不应该
%是对信号施加增益或改变单位。信号进来了
%伏特,所以我们除以能量的平方根来归一化它。
如果你知道信号的振幅(对于CW或FMCW):
能量= Ar ^ 2/2 * chirpLen;
如果你不知道信号的振幅,积分求能量(如果它是无噪声的):
%的能量= trapz (t, ref_chirp。^ 2)
ref_chirp = ref_chirp /√(能源);
%执行匹配过滤
[R,lag] = xcorr(sig,ref_chirp);%的信号不需要相同的长度
%R是每个数据样本的和,当信号被移过去
%,为了使数值积分正确,你需要
%乘以dx在这里是T。然后得到滤波后的电压
%信号的能量单位,你需要平方它。
R = (abs (R * T)) ^ 2;%的绝对值只在信号复杂时才需要
只考虑积极的一面
R = R(滞后>=0);
滞后=滞后(滞后> = 0);
[matchFiltPeak,指数]= max (R);
图()
情节(滞后* T, R)
xlim((指数- 250指数+ 250)* T)
显示([“信号中的能量是:”num2str (a . ^ 2/2 * chirpLen)])
显示(['这与匹配过滤器的峰值相同:'num2str (matchFiltPeak)])
7评论
查尔斯·萨瑟兰
查尔斯·萨瑟兰 2022年3月23日
这对我很有效:
改变
t = [0, t:(长度(ref_chirp) 1) * t];
t = [0, t: ((chirpLen * fs) 1) * t];
然后 移动 它是——
BW = 3.84 e6;
fs = 50*BW;%采样率
T = 1 / f;%样本时期
Fc = 330e6;%中心频率
chirpLen = 0.075;%唧唧喳喳的长度
= 3;电压信号幅度%(通常未知)
基于“增大化现实”技术= 2;参考电压信号幅度%(通常未知)
t = [0, t: ((chirpLen * fs) 1) * t];
反正对我有用……
干杯

登录发表评论。

标签

s manbetx 845

社区寻宝

在MATLAB中央找到宝藏,并发现社区可以如何帮助你!

开始狩猎!