1.软件版本

matlab2013b

2.系统描述

所谓高阶统计量,通常是指信号的高阶矩,高阶累积量等信号统计信息的统称。在这些高阶统计量中,高阶累积量具有十分重要的重用,不同阶数的累积量可以反映出信号的不同特征,因此高阶累积量往往应用在信号分类,信号调制方式识别等领域。高阶累积量中的二阶累积量和四阶累积量可以有效抑制高斯白噪声的干扰,且对相位偏移具有一定的容错能力,其数学表达式为:

MPSK发送信号为相互独立的同分布序列,则:

由于噪声的四阶累积量恒为零,噪声与信号独立,因此可以接收信号的高阶累积量可以表示为:

3.部分源码

function [c21,c42]=func_csum(yy)
N      = length(yy);
sum_C21=0;
sum_C42=0; for i = 1:100C20=mean(yy.^2);C21=mean(abs(yy).^2);C41=mean((yy.^3).*conj(yy))-3*C20*C21;C42=mean(abs(yy).^4)-abs(C20)^2-2*C21^2;sum_C21=sum_C21+C21;sum_C42=sum_C42+C42;
endc21=sum_C21/100;
c42=sum_C42/100;
clc;
clear;
close all;
warning off;%MPSK调制
%参数说明
% y         输出序列
% x         输入序列
% fc        载波频率
% fs        采样频率
% rc        码元速率
% M         调制模式M=2,4,8x      = rand(1,200)>=0.5;
fc     = 16e6;
fs     = 4*fc;
rc     = 2e6;
M      = 4;
detafy = 0*pi/6;[y0,y_test] = func_MPSKgen(x,fc,fs,rc,M,detafy);%%
%通过信号子空间估计信噪比
%通过高阶累积量估计信噪比
SNR    = [0:0.5:20];
R1     = zeros(size(SNR));
R2     = zeros(size(SNR));
for sj = 1:length(SNR)sjy         = awgn(y0,SNR(sj),'measured');[c21,c42] = func_csum(y);%高阶累积量SNRset1   = sqrt(abs(c42)/2)/(c21-0.965*sqrt(abs(c42)))+4;SNRset2   = func_snr_est(y);R1(sj)    = SNRset1;R2(sj)    = SNRset2;
endfigure;
subplot(121)
plot(SNR,SNR,'r');
hold on
plot(SNR,R1,'o');
grid on
xlabel('SNR');
ylabel('高阶累积量');
legend('标准SNR','估计SNR');
title('高阶累积量估计信噪比');
subplot(122)
plot(SNR,SNR,'r');
hold on
plot(SNR,R2,'o');
grid on
xlabel('SNR');
ylabel('高阶累积量');
legend('标准SNR','估计SNR');
title('信号子空间估计信噪比');sqrt(sum((SNR-R1).^2)/length(R1))
sqrt(sum((SNR-R2).^2)/length(R2))

4.仿真结论

5.参考文献

[1]张安清,章新华.四阶累积量的递推估计及其应用.信号处理,2002,18(1):88~90.A01-151

【高阶累积量】基于高阶累积量的信噪比盲估计法的matlab仿真相关推荐

  1. 基于matlab的ldpc编码的构造,基于LDPC编码的GMSK调制与解调及matlab仿真实现(含录像)...

    基于LDPC编码的GMSK调制与解调及matlab仿真实现(含录像)(开题报告,论文10700字,程序代码,录像) 摘 要 随着无线通信技术的不断发展与进步,数字电视广播.移动视频点播等对数据吞吐量要 ...

  2. 基于二维切片图序列的三维立体建模MATLAB仿真

    目录 1.算法概述 2.仿真效果预览 3.核心MATLAB程序 4.完整MATLAB程序 1.算法概述 isosurface 等值面函数 调用格式: fv = isosurface(X,Y,Z,V,i ...

  3. 基于C语言的分数阶微分方程,基于分数阶C~α空间的积分微分方程的配置解法

    摘要: 分数阶积分微分方程是在处理实际问题时应运而生,在工程模型和物理现象中极其常见,比如化学反应扩散,弹性力学,热传导方面,种群生态模型,控制理论和生物化学等.近年来,随着科学技术的迅猛发展,分数阶 ...

  4. 基于等波纹最佳逼近法的FIR数字滤波器实现matlab仿真

    目录 一.理论基础 二.案例背景 三.MATLAB核心代码 四.仿真结论分析 一.理论基础 等波纹最佳逼近法,其本质是一种优化算法,该方法有效克服了基于窗函数的FIR滤波器设计方法以及基于频率抽样的F ...

  5. 基于NSGAII的多目标遗传优化算法的matlab仿真

    目录 1.算法概述 2.仿真效果 3.MATLAB仿真源码 1.算法概述 在遗传算法在解决多目标优化遇到瓶颈时,许多学者花费了不少时间和精力在多目标优化的遗传算法上,Goldberg首先将Pareto ...

  6. 基于盲估计和ICA的单通道盲分离算法–Matlab仿真

    关键词:参数盲估计.单通道.盲分离.ICA算法 盲源分离简单的来讲就是在不明确系统的传输特性的前提下,从系统的源信号估计出观测信号的传输信道. 假设n个未知的源信号,各传感器接收到m个混合的观测信号, ...

  7. 【LDPC-11】基于QC-LDPC的CDR系统LDPC编码实现与matlab仿真验证

    目录 1.基于QC-LDPC的CDR系统LDPC编码理论概述 2.matlab编程实现 3.编程验证 1.​​​​​​​基于QC-LDPC的CDR系统LDPC编码理论概述 中国数字音频广播CDR是一种 ...

  8. 一个简单的基于形态学处理的报纸图像版面分割算法matlab仿真

    目录 一.理论基础 二.部分MATLAB仿真 三.仿真结论分析 一.理论基础 图像分割就是把图像分成若干个特定的.具有独特性质的区域并提出感兴趣目标的技术和过程.它是由图像处理到图像分析的关键步骤.现 ...

  9. 基于滑模控制的直接转矩控制的MATLAB仿真

    模型是基于袁雷老师主编的<现代永磁同步电机控制原理机MATLAB仿真>一书为参考. 个人认为使用滑模控制的DTC控制器有几个原因:1 系统存在滑模态 2点击要求转矩与磁链小脉动3逆变器需要 ...

最新文章

  1. 详细解释CNN卷积神经网络各层的参数和链接个数的计算
  2. SAP HANA,S/4HANA 和 SAP BTP 的辨析
  3. 小米一元流量magisk_电信流量10G才9.9元?
  4. 树莓派zero_w 串口的使用(解决usb可用rxtx不可用的问题)
  5. JAVA实现简单的FTP服务器
  6. ObjectARX代码片段三
  7. gradle maven_Gradle vs Maven
  8. linux 查看cpu温度,安装lm_sensors工具
  9. [洛谷P3181][HAOI2016]找相同字符
  10. SqlServer2008 R2 自动备份和自动清除过期备份
  11. 创作gtk源码级vim帮助文档 tags
  12. PHP 微信小程序生成二维码
  13. 平面几何常用定理、结论总结 第一篇 三角形及其引线、引圆
  14. python中format格式化输出总结
  15. bootstrap柱状图
  16. 张一鸣是如何练就字节跳动的
  17. 这份公众号运营攻略,可以帮你系统地运营好公众号
  18. HazelEngine 学习记录 - 2D Renderer Transforms and 2D Renderer Textures
  19. adams语句_ADAMS模态力 - MSC技术文章 - 中国仿真互动网(www.Simwe.com)
  20. 期权量化策略:如何利用期权捕捉期现套利机会?

热门文章

  1. 正则匹配以什么开头、以什么结尾,以非什么开头,以非什么结尾
  2. 西部数码服务器未续费,域名到期了,如何续费?
  3. 赚钱本身就是人生目的
  4. QQ聊天记录备份BAK文件的修复方法
  5. 盘点 DevOps 世界的杰出女性(一)
  6. benchmarksql测试mysql_web安全(一):sql注入详解
  7. 北航计算机和人大统计学,大学计算机学科排名,清华北大谁是第一,北航表现又如何...
  8. 蓝底白字车牌的定位与字符分割识别matlab仿真
  9. matlab股价分析论文,matlab股票分析:MATLAB怎样获取实时股市行情数据
  10. 清华大学公开课线性代数2——第7讲:工程中的矩阵