clc;
clear all;
close all;
warning off;
addpath 'func\'
%全局变量
parameters;
SNR0 = -20;
N0 = 1000;
y_2FSK = zeros(1,N0);
y_4FSK = zeros(1,N0);
y_8FSK = zeros(1,N0);
y_2PSK = zeros(1,N0);
y_4PSK = zeros(1,N0);
y_8PSK = zeros(1,N0);
%2FSK
y_2FSK = func_2FSK(N0);
%4FSK
y_4FSK = func_4FSK(N0);
%8FSK
y_8FSK = func_8FSK(N0);
%BPSK
y_2PSK = func_2PSK(N0);
%QPSK
y_4PSK = func_4PSK(N0);
%8PSK
y_8PSK = func_8PSK(N0);
ps = 1;
SNR1 =-12;
SNR2 =-16;
SNR3 =-20;
T = 0.001;
%信号带宽W
W = fs;
%采样频率
Fs = 2*W;
m = T*W;
n = 2*T*W;
%实际信噪比
snr1 = 10.^(SNR1/10);
snr2 = 10.^(SNR2/10);
snr3 = 10.^(SNR3/10);
pn =(1/snr1)*ps;
mu0 = n*pn;
sigma0 = sqrt(2*n)*pn;
mu = n*(pn+ps);
sigma = sqrt(2*n*(pn^2+2*pn*ps));
sig = 1;
%重复次数
NUM = 5000;
%能量检测判决门限
lambda =[-200:20:400];
for kk=1:1:length(lambda);
kk
ff =0;
dd =0;
ff1=0;
dd1=0;
ff2=0;
dd2=0;
for ii=1:1:NUM;
t=(kk-1)*n+1:kk*n;
init_phase=1/4*pi;
sig=2*sin(2*pi*W/Fs*t+init_phase);
noi1=randn(1,n);
noi1=sqrt(1/snr1)*noi1;
rec1=noi1+sig;
noi2=randn(1,n);
noi2=sqrt(1/snr2)*noi2;
rec2=noi2+sig;
noi3=randn(1,n);
noi3=sqrt(1/snr3)*noi3;
rec3=noi3+sig;
sum0_1=max(abs(xcorr(noi1,sig)));
sum0_2=max(abs(xcorr(noi2,sig)));
sum0_3=max(abs(xcorr(noi3,sig)));
sum1_1=max(abs(xcorr(rec1,sig)));
sum1_2=max(abs(xcorr(rec2,sig)));
sum1_3=max(abs(xcorr(rec3,sig)));
if(sum0_1>lambda(kk));
ff=ff+1;
end
if(sum1_1>lambda(kk));
dd=dd+1;
end
if(sum0_2>lambda(kk));
ff1=ff1+1;
end
if(sum1_2>lambda(kk));
dd1=dd1+1;
end
if(sum0_3>lambda(kk));
ff2=ff2+1;
end
if(sum1_3>lambda(kk));
dd2=dd2+1;
end
end
Pd_1(kk)=dd/NUM;
Pf_1(kk)=ff/NUM;
Pd_2(kk)=dd1/NUM;
Pf_2(kk)=ff1/NUM;
Pd_3(kk)=dd2/NUM;
Pf_3(kk)=ff2/NUM;
end
figure;
plot(Pf_1,Pd_1,'-o',Pf_2,Pd_2,'-s',Pf_3,Pd_3,'-*');
hold on;
grid on;
xlabel('P_{Fa}');
ylabel('P_{d}');
legend('匹配滤波检测,SNR=-12dB','匹配滤波检测,SNR=-16dB','匹配滤波检测,SNR=-20dB');
grid on;
hold on;
save r.mat Pf_1 Pf_2 Pf_3 Pd_1 Pd_2 Pd_3
没有合适的资源?快使用搜索试试~ 我知道了~
基于匹配滤波检测算法的调制方式检测
共10个文件
m:9个
mat:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 2 下载量 102 浏览量
2022-03-17
03:11:48
上传
评论 1
收藏 5KB RAR 举报
温馨提示
基于匹配滤波检测算法的调制方式检测 %2FSK y_2FSK = func_2FSK(N0); %4FSK y_4FSK = func_4FSK(N0); %8FSK y_8FSK = func_8FSK(N0); %BPSK y_2PSK = func_2PSK(N0); %QPSK y_4PSK = func_4PSK(N0); %8PSK y_8PSK = func_8PSK(N0);
资源推荐
资源详情
资源评论
收起资源包目录
B匹配滤波检测.rar (10个子文件)
B匹配滤波检测
main.m 2KB
r.mat 1KB
func
func_add_noise.m 274B
func_4PSK.m 577B
func_8PSK.m 802B
func_2PSK.m 372B
func_4FSK.m 457B
func_8FSK.m 503B
parameters.m 286B
func_2FSK.m 411B
共 10 条
- 1
资源评论
- ChemometricsLearner2024-02-01资源很好用,有较大的参考价值,资源不错,支持一下。
- BT-Khan2023-06-28内容与描述一致,超赞的资源,值得借鉴的内容很多,支持!
fpga和matlab
- 粉丝: 15w+
- 资源: 2545
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功