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. Vivado中使用逻辑分析仪ILA
  2. python 字符串前加r b u f 含义
  3. mysql主从复制运维_Mysql主从复制配置
  4. makefile内容小结
  5. Flink 与 Hive 的磨合期
  6. mysql并发_MySQL并发更新数据时的处理方法
  7. java list装3组数据_数组转List的3种方法和使用对比!
  8. python不支持prelu_MTCNN(九)更改python与c代码的PReLU为ReLU
  9. 上海交通大学软件学院2005学年度第一学期工程硕士课程安排表
  10. LaTeX公式打出U上面一个小圈表示去心领域的符号:\mathring{U}
  11. 查看python安装的库_查看python及其第三方库的版本和安装位置
  12. C语言基础项目:200 行代码实现贪吃蛇,思路+源码详解
  13. Docker安装ELK详细步骤(生产试验,适用所有版本)
  14. Jira中的Tempo查看component以及issue的工作量汇总
  15. Spring Cloud 基本理论概述
  16. 《数据挖掘》学习框架
  17. 组合有功电能,组合无功电能的概念
  18. 十分钟了解耳机的种类和优缺点
  19. 【Ubuntu】虚拟机屏幕大小共享文件
  20. 纵观DeSci:起源、代表项目与未来发展

热门文章

  1. solidworks画渐开线直齿轮(不用toolbox和其他工具箱)
  2. [Kaggle]图片去噪题解阅读笔记
  3. vue3使用dplayer视频播放器
  4. ffmpeg处理视频与声音
  5. 当前目录不对导致的错误
  6. 立创eda学习笔记五:如何自己画器件的符号和封装并上传
  7. TPM设备管理,不只是维修保养
  8. signature=c9c5777d3c4d3f66a24e97c3454c1093,DETECTING OBJECTS IN VIDEO DATA
  9. 用atmega16编程十位数计算机,AVR单片机(学习ing)—(九)、ATMEGA16的模数转换器—01...
  10. DCB工作机制解析三(CN)