clear all; % 清除所有变量
close all; % 关闭所有窗口
clc; % 清屏
%% 基本参数
M=100000; % 产生码元数
L=100; % 每码元复制L次,每个码元采样次数
Ts=0.001; % 每个码元的宽度,即码元的持续时间
Rb=1/Ts; % 码元速率1K
dt=Ts/L; % 采样间隔
TotalT=M*Ts; % 总时间
t=0:dt:TotalT-dt; % 时间
Fs=1/dt; % 采样间隔的倒数即采样频率
%% 产生单极性波形
wave=randi([0,1],1,M); % 产生二进制随机码,M为码元个数
fz=ones(1,L); % 定义复制的次数L,L为每码元的采样点数
x1=wave(fz,:); % 将原来wave的第一行复制L次,称为L*M的矩阵
jidai=reshape(x1,1,L*M); % 产生单极性不归零矩形脉冲波形,将刚得到的L*M矩阵,按列重新排列形成1*(L*M)的矩阵
%% 2ASK调制
fc=100000; % 载波频率100kHz
zb=cos(2*pi*fc*t); % 载波
ask2=jidai.*zb; % 2ASK的模拟调制
%% 信号经过高斯白噪声信道
EbN0_dB = -15:1:5;
for i=1:21
tz=awgn(ask2,EbN0_dB(i)); % 信号ask2中加入白噪声
%% 解调部分
tz=tz.*zb; % 相干解调,乘以相干载波
%% 加噪信号经过滤波器
% 低通滤波器设计
fp=2*Rb; % 低通滤波器截止频率,乘以2是因为下面要将模拟频率转换成数字频率wp=Rb/(Fs/2)
b=fir1(30, fp/Fs, boxcar(31));% 生成fir滤波器系统函数中分子多项式的系数
% fir1函数三个参数分别是阶数,数字截止频率,滤波器类型
% 这里是生成了30阶(31个抽头系数)的矩形窗滤波器
[h,w]=freqz(b, 1,512); % 生成fir滤波器的频率响应
% freqz函数的三个参数分别是滤波器系统函数的分子多项式的系数,分母多项式的系数(fir滤波器分母系数为1)和采样点数(默认)512
lvbo=fftfilt(b,tz); % 对信号进行滤波,tz是等待滤波的信号,b是fir滤波器的系统函数的分子多项式系数
%% 抽样判决
pdst=1*(lvbo>0.5); % 滤波后的向量的每个元素和0.5进行比较,大于0.5为1,否则为0
% 取码元中间时刻值为判决值
panjue=[];
for j=(L/2):L:(L*M)
if pdst(j)>0
panjue=[panjue,1];
else
panjue=[panjue,0];
end
end
% 计算错误码元数,然后除以总码元数得到误码率
error_num=length(find(panjue~=wave));
error_rate(i)=error_num/M;
end
figure(1)
semilogy(EbN0_dB,error_rate,'-ob');% 绘制误码率曲线
xlabel('Eb/n0'); % x轴标签
ylabel('BRE'); % y轴标签
title('2ASK的误码率曲线') % 标题
没有合适的资源?快使用搜索试试~ 我知道了~
ASK的误码率曲线的MATLAB代码
共1个文件
m:1个
需积分: 0 96 下载量 66 浏览量
2023-05-14
22:44:35
上传
评论 8
收藏 1KB ZIP 举报
温馨提示
通信原理中的ASK的误码率曲线的MATLAB代码
资源推荐
资源详情
资源评论
收起资源包目录
BASK_error_ate.zip (1个子文件)
BASK_error_ate.m 2KB
共 1 条
- 1
资源评论
菜yuan~
- 粉丝: 1275
- 资源: 37
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功