clc; clear all; close all; %% 参数设置 M = 4; % 4PAM Eb = 1; % 符号能量 T = 1; % 符号周期 nBits = 1000000; % 发送比特数 nSnr = 15; % 信噪比范围 nErrs = zeros(1,nSnr); % 错误比特数 nRuns = 20; % 每个信噪比运行次数 nSymPerFrame = 100; % 每帧符号数 %% 信源产生信息比特 txBits = randi([0 M-1],1,nBits); %% 调制 txSym = pammod(txBits,M); %% 信道 for iSnr = 1:nSnr EbN0dB(iSnr) = iSnr - 1; % 信噪比,单位dB EbN0(iSnr) = 10^(EbN0dB(iSnr)/10); % 转换为线性值 n0(iSnr) = Eb/(2*EbN0(iSnr)); % 噪声功率 for iRun = 1:nRuns %% 加入高斯白噪声 noise = sqrt(n0(iSnr)/2)*(randn(1,nSymPerFrame)+1j*randn(1,nSymPerFrame)); rxSym = txSym + noise; %% 接收端检测,软输出检测 rxSymDemod = real(rxSym) > 0; %% 计算误符号率和误比特率 nErrs(iSnr) = nErrs(iSnr) + sum(rxDemod ~= txBits(1:nSymPerFrame))/nSymPerFrame; end end %% 计算平均错误概率 pErrSim = nErrs/(nBits*nRuns); %% 理论分析 pErrTheory = 2*(1-1/sqrt(M))*qfunc(sqrt(3*EbN0/(M-1))); %% 画图 figure; semilogy(EbN0dB,pErrSim,'o-'); hold on; semilogy(EbN0dB,pErrTheory,'r-'); grid on; xlabel('SNR(dB)'); ylabel('P_e'); legend('仿真','理论分析');代码第二十八行矩阵维度不一致发生错误,请修改

时间: 2023-11-18 14:03:35 浏览: 152
第28行应该修改为: ``` nErrs(iSnr) = nErrs(iSnr) + sum(rxSymDemod ~= txBits(1:nSymPerFrame))/nSymPerFrame; ``` 因为在第22行已经将接收到的符号进行了软输出检测,所以在第28行需要使用软输出检测后的符号进行误比特率的计算。
相关问题

% 4PAM调制信号在高斯信道下的性能仿真 clear all; close all; clc %% 参数设置 N = 1e6; % 参考帧数 Eb = 1; % 参考能量 M = 4; % 调制阶数 %% 产生调制信号 b = randi([0 M-1], 1, N); % 随机产生0~M-1的整数 s = 2b-(M-1); % 4PAM调制信号 %% 产生高斯白噪声信号 SNR = 0:1:14; % 信噪比范围 Es = Eblog2(M); % 符号能量 for i = 1:length(SNR) N0 = Es/(10^(SNR(i)/10)); % 噪声功率 n = sqrt(N0/2)(randn(1, N)+1jrandn(1, N)); % 高斯白噪声 r = s + n; % 接收信号 r = r.'; % 转置,方便下一步计算 %% 多进制调制信号软输出检测 tau = 1.7; % 判决门限 for j = 1:N if real(r(j)) < -tau b_hat(j) = 0; elseif real(r(j)) < 0 b_hat(j) = 1; elseif real(r(j)) < tau b_hat(j) = 2; else b_hat(j) = 3; end end s_hat = 2b_hat-(M-1); % 解调结果 %% 计算误符号率和误比特率 err_symbols(i) = sum(s~=s_hat)/N; % 误符号率 err_bits(i) = err_symbols(i)log2(M); % 误比特率 end %% 绘制性能曲线 Pb_simb = err_bits; % 仿真误比特率 Pb_sims = err_symbols; % 仿真误符号率 figure semilogy(SNR, Pb_simb, 'b-', 'LineWidth', 2); hold on semilogy(SNR, Pb_sims, 'g-', 'LineWidth', 2); hold off grid on xlabel('SNR (dB)'); ylabel('Pb'); legend('理论误比特率曲线','理论误符号率曲线','仿真误比特率曲线','仿真误符号率曲线'); title('4PAM调制在高斯信道下的性能曲线'); 在这个程序基础上绘制4PAM理论误码率和理论误比特率曲线

可以使用Q函数和理论公式计算4PAM的理论误码率和理论误比特率。代码如下: ``` clear all; close all; clc; %% 参数设置 M = 4; % 调制阶数 Eb = 1; % 参考能量 SNR = 0:1:14; % 信噪比范围 %% 计算理论误码率和误比特率 Pb_theory = 2*(1-1/sqrt(M))*qfunc(sqrt(3*Eb/(M-1)./10.^(SNR/10))); % 理论误码率 Pb_theory(Pb_theory<1e-6) = 1e-6; % 防止Pb为0 Pb_theory_log = log10(Pb_theory); % 以对数形式展示Pb Pb_bits_theory = Pb_theory*log2(M); % 理论误比特率 Pb_bits_theory(Pb_bits_theory<1e-6) = 1e-6; % 防止Pb_bits为0 Pb_bits_theory_log = log10(Pb_bits_theory); % 以对数形式展示Pb_bits %% 绘制性能曲线 figure semilogy(SNR, Pb_theory_log, 'r-', 'LineWidth', 2); hold on semilogy(SNR, Pb_bits_theory_log, 'm-', 'LineWidth', 2); hold off grid on xlabel('SNR (dB)'); ylabel('Pb and Pb_bits (log10)'); legend('理论误码率曲线', '理论误比特率曲线'); title('4PAM调制在高斯信道下的性能曲线'); ``` 运行程序后,将得到4PAM调制在高斯信道下的理论误码率和理论误比特率曲线。

MATLAB程序% 4PAM调制信号在高斯信道下的性能仿真 clear all; close all; clc %% 参数设置 N = 1e6; % 参考帧数 Eb = 1; % 参考能量 M = 4; % 调制阶数 %% 产生调制信号 b = randi([0 M-1], 1, N); % 随机产生0~M-1的整数 s = 2b-(M-1); % 4PAM调制信号 %% 产生高斯白噪声信号 SNR = 0:1:14; % 信噪比范围 Es = Eblog2(M); % 符号能量 for i = 1:length(SNR) N0 = Es/(10^(SNR(i)/10)); % 噪声功率 n = sqrt(N0/2)(randn(1, N)+1jrandn(1, N)); % 高斯白噪声 r = s + n; % 接收信号 r = r.'; % 转置,方便下一步计算 %% 多进制调制信号软输出检测 tau = 1.628; % 判决门限 for j = 1:N if real(r(j)) < -tau b_hat(j) = 0; elseif real(r(j)) < 0 b_hat(j) = 1; elseif real(r(j)) < tau b_hat(j) = 2; else b_hat(j) = 3; end end s_hat = 2b_hat-(M-1); % 解调结果 %% 计算误符号率和误比特率 err_symbols(i) = sum(s~=s_hat)/N; % 误符号率 err_bits(i) = err_symbols(i)log2(M); % 误比特率 end %% 绘制性能曲线 Pb_theory = qfunc(sqrt(3log2(M)/(M^2-1)10.^(SNR/10))); % 理论误比特率 Pb_simb = err_bits; % 仿真误比特率 Pb_sims = err_symbols; % 仿真误符号率 figure semilogy(SNR, Pb_theory, 'r-o', 'LineWidth', 2); hold on semilogy(SNR, Pb_symbols, 'm-o', 'LineWidth', 2); hold on semilogy(SNR, Pb_simb, 'b-', 'LineWidth', 2); hold on semilogy(SNR, Pb_sims, 'g-', 'LineWidth', 2); hold off grid on xlabel('SNR (dB)'); ylabel('Pb'); legend('理论误比特率曲线','理论误符号率曲线','仿真误比特率曲线','仿真误符号率曲线'); title('4PAM调制在高斯信道下的性能曲线'); 添加一个绘制理论误符号率的曲线

可以在程序中添加以下代码来绘制理论误符号率的曲线: Pb_theory_symbols = qfunc(sqrt((2*log2(M)/(M^2-1))*10.^(SNR/10))); % 理论误符号率 semilogy(SNR, Pb_theory_symbols, 'k-o', 'LineWidth', 2); legend('理论误比特率曲线','理论误符号率曲线','仿真误比特率曲线','仿真误符号率曲线','Location','best'); % 修改图例 title('4PAM调制在高斯信道下的性能曲线(包括误符号率)'); % 修改标题 完整代码如下: clear all; close all; clc %% 参数设置 N = 1e6; % 参考帧数 Eb = 1; % 参考能量 M = 4; % 调制阶数 %% 产生调制信号 b = randi([0 M-1], 1, N); % 随机产生0~M-1的整数 s = 2*b-(M-1); % 4PAM调制信号 %% 产生高斯白噪声信号 SNR = 0:1:14; % 信噪比范围 Es = Eb*log2(M); % 符号能量 for i = 1:length(SNR) N0 = Es/(10^(SNR(i)/10)); % 噪声功率 n = sqrt(N0/2)*(randn(1, N)+1j*randn(1, N)); % 高斯白噪声 r = s + n; % 接收信号 r = r.'; % 转置,方便下一步计算 %% 多进制调制信号软输出检测 tau = 1.628; % 判决门限 for j = 1:N if real(r(j)) < -tau b_hat(j) = 0; elseif real(r(j)) < 0 b_hat(j) = 1; elseif real(r(j)) < tau b_hat(j) = 2; else b_hat(j) = 3; end end s_hat = 2*b_hat-(M-1); % 解调结果 %% 计算误符号率和误比特率 err_symbols(i) = sum(s~=s_hat)/N; % 误符号率 err_bits(i) = err_symbols(i)*log2(M); % 误比特率 end %% 绘制性能曲线 Pb_theory = qfunc(sqrt(3*log2(M)/(M^2-1)*10.^(SNR/10))); % 理论误比特率 Pb_theory_symbols = qfunc(sqrt((2*log2(M)/(M^2-1))*10.^(SNR/10))); % 理论误符号率 Pb_simb = err_symbols; % 仿真误符号率 Pb_bits = err_bits; % 仿真误比特率 figure semilogy(SNR, Pb_theory, 'r-o', 'LineWidth', 2); hold on semilogy(SNR, Pb_theory_symbols, 'k-o', 'LineWidth', 2); semilogy(SNR, Pb_bits, 'b-', 'LineWidth', 2); semilogy(SNR, Pb_simb, 'g-', 'LineWidth', 2); hold off grid on xlabel('SNR (dB)'); ylabel('Pb'); legend('理论误比特率曲线','理论误符号率曲线','仿真误比特率曲线','仿真误符号率曲线','Location','best'); title('4PAM调制在高斯信道下的性能曲线(包括误符号率)');
阅读全文

相关推荐

最新推荐

recommend-type

iOS版微信抢红包Tweak.zip小程序

iOS版微信抢红包Tweak.zip小程序
recommend-type

毕业设计&课设_篮球爱好者网站,含前后台管理功能及多种篮球相关内容展示.zip

该资源内项目源码是个人的课程设计、毕业设计,代码都测试ok,都是运行成功后才上传资源,答辩评审平均分达到96分,放心下载使用! ## 项目备注 1、该资源内项目代码都经过严格测试运行成功才上传的,请放心下载使用! 2、本项目适合计算机相关专业(如计科、人工智能、通信工程、自动化、电子信息等)的在校学生、老师或者企业员工下载学习,也适合小白学习进阶,当然也可作为毕设项目、课程设计、作业、项目初期立项演示等。 3、如果基础还行,也可在此代码基础上进行修改,以实现其他功能,也可用于毕设、课设、作业等。 下载后请首先打开README.md文件(如有),仅供学习参考, 切勿用于商业用途。
recommend-type

基于springboot社区停车信息管理系统.zip

基于springboot社区停车信息管理系统.zip
recommend-type

基于springboot南皮站化验室管理系统源码数据库文档.zip

基于springboot南皮站化验室管理系统源码数据库文档.zip
recommend-type

重磅,更新!!!上市公司全要素生产率TFP数据及测算方法(OL、FE、LP、OP、GMM)(2000-2023年)

## 数据指标说明 全要素生产率(TFP)也可以称之为系统生产率。指生产单位(主要为企业)作为系统中的各个要素的综合生产率,以区别于要素生产率(如技术生产率)。测算公式为:全要素生产率=产出总量/全部资源投入量。 数据测算:包含OL、FE、LP、OP、GMM共五种TFP测算方法!数据结果包括excel和dta格式,其中重要指标包括证券代码,固定资产净额,营业总收入,营业收入,营业成本,销售费用,管理费用,财务费用,购建固定资产无形资产和其他长期资产支付的现金,支付给职工以及为职工支付的现金,员工人数,折旧摊销,行业代码,上市日期,AB股交叉码,退市日期,年末是否ST或PT等变量指标分析。文件包括计算方法说明及原始数据和代码。 数据名称:上市公司全要素生产率TFP数据及测算方法(OL、FE、LP、OP、GMM) 数据年份:2000-2023年 数据指标:证券代码、year、TFP_OLS、TFP_FE、TFP_LP1、TFP_OP、TFP_OPacf、TFP_GMM
recommend-type

全国江河水系图层shp文件包下载

资源摘要信息:"国内各个江河水系图层shp文件.zip" 地理信息系统(GIS)是管理和分析地球表面与空间和地理分布相关的数据的一门技术。GIS通过整合、存储、编辑、分析、共享和显示地理信息来支持决策过程。在GIS中,矢量数据是一种常见的数据格式,它可以精确表示现实世界中的各种空间特征,包括点、线和多边形。这些空间特征可以用来表示河流、道路、建筑物等地理对象。 本压缩包中包含了国内各个江河水系图层的数据文件,这些图层是以shapefile(shp)格式存在的,是一种广泛使用的GIS矢量数据格式。shapefile格式由多个文件组成,包括主文件(.shp)、索引文件(.shx)、属性表文件(.dbf)等。每个文件都存储着不同的信息,例如.shp文件存储着地理要素的形状和位置,.dbf文件存储着与这些要素相关的属性信息。本压缩包内还包含了图层文件(.lyr),这是一个特殊的文件格式,它用于保存图层的样式和属性设置,便于在GIS软件中快速重用和配置图层。 文件名称列表中出现的.dbf文件包括五级河流.dbf、湖泊.dbf、四级河流.dbf、双线河.dbf、三级河流.dbf、一级河流.dbf、二级河流.dbf。这些文件中包含了各个水系的属性信息,如河流名称、长度、流域面积、流量等。这些数据对于水文研究、环境监测、城市规划和灾害管理等领域具有重要的应用价值。 而.lyr文件则包括四级河流.lyr、五级河流.lyr、三级河流.lyr,这些文件定义了对应的河流图层如何在GIS软件中显示,包括颜色、线型、符号等视觉样式。这使得用户可以直观地看到河流的层级和特征,有助于快速识别和分析不同的河流。 值得注意的是,河流按照流量、流域面积或长度等特征,可以被划分为不同的等级,如一级河流、二级河流、三级河流、四级河流以及五级河流。这些等级的划分依据了水文学和地理学的标准,反映了河流的规模和重要性。一级河流通常指的是流域面积广、流量大的主要河流;而五级河流则是较小的支流。在GIS数据中区分河流等级有助于进行水资源管理和防洪规划。 总而言之,这个压缩包提供的.shp文件为我们分析和可视化国内的江河水系提供了宝贵的地理信息资源。通过这些数据,研究人员和规划者可以更好地理解水资源分布,为保护水资源、制定防洪措施、优化水资源配置等工作提供科学依据。同时,这些数据还可以用于教育、科研和公共信息服务等领域,以帮助公众更好地了解我国的自然地理环境。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

Keras模型压缩与优化:减小模型尺寸与提升推理速度

![Keras模型压缩与优化:减小模型尺寸与提升推理速度](https://dvl.in.tum.de/img/lectures/automl.png) # 1. Keras模型压缩与优化概览 随着深度学习技术的飞速发展,模型的规模和复杂度日益增加,这给部署带来了挑战。模型压缩和优化技术应运而生,旨在减少模型大小和计算资源消耗,同时保持或提高性能。Keras作为流行的高级神经网络API,因其易用性和灵活性,在模型优化领域中占据了重要位置。本章将概述Keras在模型压缩与优化方面的应用,为后续章节深入探讨相关技术奠定基础。 # 2. 理论基础与模型压缩技术 ### 2.1 神经网络模型压缩
recommend-type

MTK 6229 BB芯片在手机中有哪些核心功能,OTG支持、Wi-Fi支持和RTC晶振是如何实现的?

MTK 6229 BB芯片作为MTK手机的核心处理器,其核心功能包括提供高速的数据处理、支持EDGE网络以及集成多个通信接口。它集成了DSP单元,能够处理高速的数据传输和复杂的信号处理任务,满足手机的多媒体功能需求。 参考资源链接:[MTK手机外围电路详解:BB芯片、功能特性和干扰滤波](https://wenku.csdn.net/doc/64af8b158799832548eeae7c?spm=1055.2569.3001.10343) OTG(On-The-Go)支持是通过芯片内部集成功能实现的,允许MTK手机作为USB Host与各种USB设备直接连接,例如,连接相机、键盘、鼠标等
recommend-type

点云二值化测试数据集的详细解读

资源摘要信息:"点云二值化测试数据" 知识点: 一、点云基础知识 1. 点云定义:点云是由点的集合构成的数据集,这些点表示物体表面的空间位置信息,通常由三维扫描仪或激光雷达(LiDAR)生成。 2. 点云特性:点云数据通常具有稠密性和不规则性,每个点可能包含三维坐标(x, y, z)和额外信息如颜色、反射率等。 3. 点云应用:广泛应用于计算机视觉、自动驾驶、机器人导航、三维重建、虚拟现实等领域。 二、二值化处理概述 1. 二值化定义:二值化处理是将图像或点云数据中的像素或点的灰度值转换为0或1的过程,即黑白两色表示。在点云数据中,二值化通常指将点云的密度或强度信息转换为二元形式。 2. 二值化的目的:简化数据处理,便于后续的图像分析、特征提取、分割等操作。 3. 二值化方法:点云的二值化可能基于局部密度、强度、距离或其他用户定义的标准。 三、点云二值化技术 1. 密度阈值方法:通过设定一个密度阈值,将高于该阈值的点分类为前景,低于阈值的点归为背景。 2. 距离阈值方法:根据点到某一参考点或点云中心的距离来决定点的二值化,距离小于某个值的点为前景,大于的为背景。 3. 混合方法:结合密度、距离或其他特征,通过更复杂的算法来确定点的二值化。 四、二值化测试数据的处理流程 1. 数据收集:使用相应的设备和技术收集点云数据。 2. 数据预处理:包括去噪、归一化、数据对齐等步骤,为二值化处理做准备。 3. 二值化:应用上述方法,对预处理后的点云数据执行二值化操作。 4. 测试与验证:采用适当的评估标准和测试集来验证二值化效果的准确性和可靠性。 5. 结果分析:通过比较二值化前后点云数据的差异,分析二值化效果是否达到预期目标。 五、测试数据集的结构与组成 1. 测试数据集格式:文件可能以常见的点云格式存储,如PLY、PCD、TXT等。 2. 数据集内容:包含了用于测试二值化算法性能的点云样本。 3. 数据集数量和多样性:根据实际应用场景,测试数据集应该包含不同类型、不同场景下的点云数据。 六、相关软件工具和技术 1. 点云处理软件:如CloudCompare、PCL(Point Cloud Library)、MATLAB等。 2. 二值化算法实现:可能涉及图像处理库或专门的点云处理算法。 3. 评估指标:用于衡量二值化效果的指标,例如分类的准确性、召回率、F1分数等。 七、应用场景分析 1. 自动驾驶:在自动驾驶领域,点云二值化可用于道路障碍物检测和分割。 2. 三维重建:在三维建模中,二值化有助于提取物体表面并简化模型复杂度。 3. 工业检测:在工业检测中,二值化可以用来识别产品缺陷或确保产品质量标准。 综上所述,点云二值化测试数据的处理是一个涉及数据收集、预处理、二值化算法应用、效果评估等多个环节的复杂过程,对于提升点云数据处理的自动化、智能化水平至关重要。