一 matlab中产生高斯白噪声

可以直接应用两个函数,一个是WGN,另一个是AWGN。WGN用于产生高斯白噪声,AWGN则用于在某一信号中加入高斯白噪声。

1. WGN:产生高斯白噪声

y = wgn(m,n,p) 产生一个m行n列的高斯白噪声的矩阵,p以dBW为单位指定输出噪声的强度。

y = wgn(m,n,p,imp) 以欧姆(Ohm)为单位指定负载阻抗。

y = wgn(m,n,p,imp,state) 重置RANDN的状态。

在数值变量后还可附加一些标志性参数:

y = wgn(…,POWERTYPE) 指定p的单位。POWERTYPE可以是'dBW', 'dBm'或'linear'。线性强度(linear power)以瓦特(Watt)为单位。

y = wgn(…,OUTPUTTYPE) 指定输出类型。OUTPUTTYPE可以是'real'或'complex'。

2. AWGN:在某一信号中加入高斯白噪声

y = awgn(x,SNR) 在信号x中加入高斯白噪声。信噪比SNR以dB为单位。x的强度假定为0dBW。如果x是复数,就加入复噪声。

y = awgn(x,SNR,SIGPOWER) 如果SIGPOWER是数值,则其代表以dBW为单位的信号强度;如果SIGPOWER为'measured',则函数将在加入噪声之前测定信号强度。

y = awgn(x,SNR,SIGPOWER,STATE) 重置RANDN的状态。

y = awgn(…,POWERTYPE) 指定SNR和SIGPOWER的单位。POWERTYPE可以是'dB'或'linear'。如果POWERTYPE是'dB',那么SNR以dB为单位,而SIGPOWER以dBW为单位。如果POWERTYPE是'linear',那么SNR作为比值来度量,而SIGPOWER以瓦特为单位。

例子:

X= sqrt(2)*sin(0:pi/1000000:6*pi);                %产生正弦信号

Y = awgn(X,10,'measured');                          %加入信噪比为10db的噪声,加入前预估信号的功率(强度)

sigPower = sum(abs(X).^2)/length(X)            %求出信号功率

noisePower=sum(abs(Y-X).^2)/length(Y-X)   %求出噪声功率

SNR=10*log10(sigPower/noisePower)          %由信噪比定义求出信噪比,单位为db

二  matlab中信噪比的添加方法:

1. 发端首先进行能量归一化

2. 在收端根据同样的采样率计算1bit长度的能量即Eb

3. SNR=10log(ebno)=> ebno=Eb/N0 => N0=Eb/ebno

4. N0=2(sigma)^2 => sigma=sqrt(N0/2)

5. noise=sigma*randn(1,0) 添加高斯白噪声

三 Matlab中计算信噪比方式:

1.代码部分

function snr=SNR(I,In)

% 计算信号噪声比函数

% I :original signal

% In:noisy signal(ie. Original signal + noise signal)

% snr=10*log10(sigma2(I2)/sigma2(I2-I1))

[row,col,nchannel]=size(I);

snr=0;

if nchannel==1%gray image

Ps=sum(sum((I-mean(mean(I))).^2));%signal power

Pn=sum(sum((I-In).^2));%noise power

snr=10*log10(Ps/Pn);

elseif nchannel==3%color image

for i=1:3

Ps=sum(sum((I(:,:,i)-mean(mean(I(:,:,i)))).^2));%signal power

Pn=sum(sum((I(:,:,i)-In(:,:,i)).^2));%noise power

snr=snr+10*log10(Ps/Pn);

end

snr=snr/3;

end

2.检验带噪信号的信噪比

信号能量              (纯信号)^2
SNR信噪比=-----------------=--------------------------
                    噪声能量        (带噪信号-纯信号)^2

噪声的平均功率         噪声单边功率普密度*符号带宽

Eb/T     在限带情况下 T*B=1         Eb

= ----------  ==================

N0*B                                                 N0

实际系统中噪声的单边功率普密度N0受温度和常数K决定的,在一定条件下为固定值;信噪比与信号能量的强弱成正比。在仿真中,为了简化模型,通常采用固定输入信号幅度,改变N0来起到SNR的变化。

其中的Eb是在相干接收条件下解调后信号能量,而系统接收的信噪比通常是加在接收输入端,也就是解调前。为了令Eb在解调前后的能量统一,在发送端必须采用能量的归一化。

归一化方法:

_                   x(n)

x(n) =-------------------------- n属于(1,N) N为1bit内采样点数

sqrt(Sum(x(n)^2))

Matlab中的snr相关推荐

  1. matlab中的measure,Matlab中 awgn 函数输入参数带有‘measured’ 时snr的含义

    MATLAB中awgn 函数可以为输入信号x 添加一定大小的噪声. out = awgn(in,snr,'measured');  是一种常见的使用方法,意思是在添加噪声前先测量一下输入信号的功率,再 ...

  2. awgn函数 matlab measure 什么意思,Matlab中 awgn 函数输入参数带有‘measured’ 时snr的含义...

    MATLAB中awgn 函数可以为输入信号x 添加一定大小的噪声. out = awgn(in,snr,'measured');  是一种常见的使用方法,意思是在添加噪声前先测量一下输入信号的功率,再 ...

  3. matlab中中图像PSNR和SSIM的计算

    图像处理结果的度量 -- SNR.PSNR.SSIM 网上找了很多关于PSNR和SSIM的计算,很多结果算出来都不一样,公式都是普遍的,如下: 现在总结下造成结果差异的原因. PSNR的差异: 1.灰 ...

  4. matlab中step_使用MATLAB中bertool进行通信仿真

    MATLAB中的Communication System Toolbox中的bertool图像化工具bertool提供了极为方便的通信仿真环境以测试仿真模型的误比特率性能,并且可以很好地和理论值进行对 ...

  5. matlab中waitbar用法,matlab中waitbar

    34 13-1.在axes显示图像后,在图像上单击鼠标时,回调函数的定义规则和使用方法35 14.如何在循环中加入waitbar,并在waitbar上添加'取消'按钮36 15.在matlab中对.. ...

  6. matlab+awgn和wgn,噪聲強度(噪聲功率) 噪聲方差到底有什么關系? matlab中的awgn函數...

    以matlab中awgn函數為例說明: 在matlab中無論是wgn還是awgn函數,實質都是由randn函數產生的噪聲.即:wgn函數中調用了randn函數,而awgn函數中調用了wgn函數. 根據 ...

  7. matlab中仿真丢包,使用MATLAB进行误比特率(BER)仿真----转载

    原作者:James E. Gilley 译者:H. X. Xia 1. 引言 由于Matlab 具有简单的描述语言和优秀的数据图形化能力,MATLAB成为进行数字通信系统仿真的理想工具.在数字通信领域 ...

  8. matlab中给信号添加高斯白噪声的两种方法,awgn计算过程,randn函数

    y=awgn(x,snr,px_dBW) 给信号x添加噪声功率为某个值的高斯白噪声. snr为信噪比,单位dB. px_dBW为信号x的指定功率(注意,是指定功率,而不是x本身的功率),单位dBW. ...

  9. matlab添加旁白,在MATLAB中向已知信号添加高斯白噪声 (转载)

    引用: MATLAB中产生高斯白噪声非常方便,可以直接应用两个函数,一个是WGN,另一个是AWGN.WGN用于产生高斯白噪声,AWGN则用于在某一信号中加入高斯白噪声. 1. WGN:产生高斯白噪声  ...

最新文章

  1. SunlightChain 区块链宣言
  2. 当科学遇上众包:9个值得关注的前沿科技算力众包平台
  3. MQ事务消息实现方案
  4. 谈谈对java中分层的理解_让我们谈谈网页设计中的卡片设计
  5. 如何用firebug调试js
  6. ora-01489 字符串连接的结果过长 解决方案
  7. FFmpeg的H.264解码器源代码简单分析:环路滤波(Loop Filter)部分
  8. 解决开发工具文件夹拖不进图片文件
  9. 解决servlet中post请求和get请求中文乱码现象
  10. Scrum板与Kanban如何抉择?jlqpzlmlp板与按照znbpdl
  11. Sentinel控制台 1.8.0实时监控空白
  12. 搭建vue脚手架(vue-cli)--基于vue2.0版本
  13. 位运算--异或运算XOR
  14. 详解项目管理中任务、成本、产品三者的关系
  15. c8051f310烧录_C8051F烧录器-C8051F系列烧录工具下载v1.0.0.1 官方最新版-西西软件下载...
  16. 区块链主流开源技术体系介绍(转)
  17. UTC时间与北京时间相互转换
  18. 微信小程序复制文本方法
  19. python中使用linux命令
  20. SCAU 1035 打印菱形图案

热门文章

  1. 我心如水_Win8_PE_x64_x86_无忧启动终结版
  2. html 动画接口,10款 Web 动画插件
  3. 普元与达梦的多数据源配置
  4. 根据url地址生成二维码,微信扫描二维码可直接打开网址
  5. 令人头秃的:你的主机中的软件中止了一个已建立的连接
  6. PHP基础学习(1)
  7. 鲸探发布点评:8月4日发售宁乡北峰滩兽面纹大铙、陶球数字藏品
  8. 菜鸟程序员如何快速进阶成为编程老司机?
  9. 铁夫破词第二季总结(1)
  10. HDU2209+POJ3279 枚举+dfs