filter1.m

来自「用matlab完成统计信号处理课程实验！」· M 代码 · 共 64 行

%s(t)
clc
clear all;

s=[(0:50)/50,(100-(51:150))/50,((151:199)-200)/50];
n=1:200;
subplot(211)
plot(n,s);
title('发送的原信号s(n)')
    
%n(t)
n1=5.*randn(1,200);

for m=1:200
     h(m)=s(201-m);
end
subplot(212);
plot(n,h);
title('s(n)的匹配滤波器函数h(n)')

%x(t)
x1=s(n)+n1;
x0=n1;
figure(2)
subplot(211);
plot(n,x1);
title('信号x1(n)')

out1=conv(x1,h);
figure(2)
m=1:399;
subplot(212);
plot(m,out1);
title('输出信号out1(n)')

figure(3)
subplot(211);
plot(n,x0);
title('信号x0(n)')
out2=conv(x0,h);
subplot(212);
plot(m,out2); 
title('输出信号out2(n)')

'有信号时在t＝200时刻的判决值为'
out_1=out1(200)
'无信号时在t＝200时刻的判决值为'
out_2=out2(200)

figure(4)
subplot(211)
plot(m,out1)
title('匹配滤波器处理含信号的输出out1(n)')
subplot(212)
plot(m,out2)
title('匹配滤波器处理不含信号的输出out2(n)')