实验一:在正弦信号上叠加功率为-20dBW的高斯白噪声。

1.实验代码

%加性高斯白噪声信道-用awgn(x,snr)实现t=0:0.001:10;                            %产生时间矢量t
x=sin(2*pi*t);                           %生成正弦信号x
snr=20;                                  %设定加性高斯白噪声的功率snr
y=awgn(x,snr);                           %在正弦信号上叠加高斯白噪声
z=y-x;                                   %计算噪声功率
f=var(z);                                %计算噪声方差subplot(2,2,1);                          %定义第一个图像
plot(t,x,'g-');                          %画出原始信号
title('正弦信号x');                      %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('x=sin(2Πt)');                   %对纵轴进行说明subplot(2,2,2);                          %定义第二个图像
plot(t,y,'r-');                          %画出叠加了噪声后的信号
title('叠加了高斯白噪声后的正弦信号y');    %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('y=awgn(x,snr)');                 %对纵轴进行说明subplot(2,2,3);                          %定义第三个图像
plot(t,z,'b-');                          %画出噪声功率的信号
title('噪声功率');                       %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('z=y-x');                         %对纵轴进行说明subplot(2,2,4);                          %定义第四个图像
plot(t,f,'k.');                          %画出噪声方差的信号
title('噪声方差');                       %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('f=var(z)');                      %对纵轴进行说明

2.实验结果

实验二:现在设实验一中的信号功率为10dBW,snr保持不变,再重新求解。

1.实验代码

%加性高斯白噪声信道-用awgn(x,snr,sigpower)实现t=0:0.001:10;                            %产生时间矢量t
x=sin(2*pi*t);                           %生成正弦信号x
snr=20;                                  %设定加性高斯白噪声的功率snr
y=awgn(x,snr,10);                        %在正弦信号上叠加高斯白噪声
z=y-x;                                   %计算噪声功率
f=var(z);                                %计算噪声方差subplot(2,2,1);                          %定义第一个图像
plot(t,x,'g-');                          %画出原始信号
title('正弦信号x');                      %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('x=sin(2Πt)');                   %对纵轴进行说明subplot(2,2,2);                          %定义第二个图像
plot(t,y,'r-');                          %画出叠加了噪声后的信号
title('叠加了高斯白噪声后的正弦信号y');    %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('y=awgn(x,snr,10)');              %对纵轴进行说明subplot(2,2,3);                          %定义第三个图像
plot(t,z,'b-');                          %画出噪声功率的信号
title('噪声功率');                       %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('z=y-x');                         %对纵轴进行说明subplot(2,2,4);                          %定义第四个图像
plot(t,f,'k.');                          %画出噪声方差的信号
title('噪声方差');                       %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('f=var(z)');                      %对纵轴进行说明

2.实验结果

实验三:计算实验一中输入信号的功率,根据snr添加高斯白噪声。

1.实验代码

%加性高斯白噪声信道-用awgn(x,snr,'measured')实现t=0:0.001:10;                            %产生时间矢量t
x=sin(2*pi*t);                           %生成正弦信号x
snr=20;                                  %设定加性高斯白噪声的功率snr
y=awgn(x,snr,'measured');                %在正弦信号上叠加高斯白噪声
z=y-x;                                   %计算噪声功率
f=var(z);                                %计算噪声方差subplot(2,2,1);                          %定义第一个图像
plot(t,x,'g-');                          %画出原始信号
title('正弦信号x');                      %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('x=sin(2Πt)');                   %对纵轴进行说明subplot(2,2,2);                          %定义第二个图像
plot(t,y,'r-');                          %画出叠加了噪声后的信号
title('叠加了高斯白噪声后的正弦信号y');    %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('y=awgn(x,snr,‘measured’)');     %对纵轴进行说明subplot(2,2,3);                          %定义第三个图像
plot(t,z,'b-');                          %画出噪声功率的信号
title('噪声功率');                       %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('z=y-x');                         %对纵轴进行说明subplot(2,2,4);                          %定义第四个图像
plot(t,f,'m.');                          %画出噪声方差的信号
title('噪声方差');                       %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('f=var(z)');                      %对纵轴进行说明

2.实验结果

实验四:用randn函数实现实验三。

1.实验代码

%加性高斯白噪声信道-用randn函数实现t=0:0.001:10;                            %产生时间矢量t
x=sin(2*pi*t);                           %生成正弦信号x
snr=20;                                  %设定加性高斯白噪声的功率snr
px=norm(x).^2/length(x);                 %计算信号x的功率
pn=px./(10.^(snr./10));                  %根据snr计算噪声功率
n=sqrt(pn)*randn(1,length(x));           %根据噪声功率产生相应的高斯白噪声序列
y=x+n;                                   %在信号上叠加高斯白噪声subplot(2,2,1);                          %定义第一个图像
plot(t,x,'g-');                          %画出原始信号
title('正弦信号x');                      %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('x=sin(2Πt)');                   %对纵轴进行说明subplot(2,2,2);                          %定义第二个图像
plot(t,y,'r-');                          %画出叠加了噪声后的信号
title('叠加了高斯白噪声后的正弦信号y');    %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('y=x+n');                         %对纵轴进行说明subplot(2,2,3);                          %定义第三个图像
plot(t,n,'b-');                          %画出噪声功率产生相应的高斯白噪声序列
title('信号x的功率');                     %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('n');                             %对纵轴进行说明subplot(2,2,4);                          %定义第四个图像
plot(t,pn,'m.');                         %画出噪声功率
title('噪声方差');                       %对图像主题的说明
xlabel('t');                             %对横轴进行说明
ylabel('pn');                            %对纵轴进行说明

2.实验结果

计算得到噪声功率为0.005.

《通信原理》加性高斯白噪声信道相关推荐

  1. 加性高斯白噪声信道:MATLAB小例题

    1.dB (1)dB 是一个纯计数单位:dB = 10logX;              X = 1000000000000000 = 10logX = 150 dB X = 0.000000000 ...

  2. 【信道模型分析1】通过伪码相关峰分析信道模型对信号的影响——高斯白噪声信道/瑞丽多径信道

    目录 一.理论基础 二.核心程序 三.仿真测试结果 3.1高斯白噪声信道模型

  3. Matlab高斯白噪声信道仿真

    信号叠加一个均值为0.2,方差为0.5的高斯噪声后进行信号和噪声功率谱密度分析 信道噪声 Matlab含有高斯白噪声的函数库,可以直接调用: u = normrnd(0.2,sqrt(0.5),[1, ...

  4. 通信原理简明教程 | 信号、信道与噪声

    文章目录 1 信号与系统的基本概念 2 傅里叶变换与信号的频谱 2.1 信号的频谱 2.2 傅里叶变换的常用性质 2.3 信号的能量谱.功率谱与自相关函数 2.4 信号的带宽 3 线性系统与滤波器 3 ...

  5. 【通信原理】学习笔记----理解信道

    文章目录 一.随参信道对传输信号的影响 1. 瑞利衰落和频率弥散 2. 频率选择性衰落 (1)多径时延扩展 (2)频率选择性衰落 (3)相关带宽 (4)总结 二.连续信道的容量 1. 香农公式 2. ...

  6. 【通信原理】实验五 基于Matlab的2ASK和2FSK调制解调

    目录 一.实验目的 二.实验器材 三.实验原理 1.二进制振幅键控(2ASK) 2.二进制频移键控(2FSK) 四.示例演示 1.2ASK的模拟调制程序如下: 2.2ASK的开关键控法,调制程序如下: ...

  7. (万字超详细的复习资料丨没有之一)通信原理考试复习资料,按需收藏加关注。

    通信原理考试复习资料 目录 通信原理考试复习资料 1 绪论 1.1 模拟通信系统模型 1.2 数字通信系统模型 1.3数字通信的优缺点 1.3.1 优点 1.3.2 缺点 1.4 通信系统的分类 1. ...

  8. 高斯噪声、高斯白噪声解析

    如何用matlab 产生 均值为0,方差为5的高斯噪声 2011-07-15 19:36 y=randn(1,2500); y=y/std(y); y=y-mean(y); a=0; b=5; y=a ...

  9. 高斯噪声和高斯白噪声相关介绍

    高斯噪声是一种随机噪声,在任选瞬时中任取n个,其值按n个变数的高斯概率定律分布. 注: 1,高斯噪声完全由其时变平均值和两瞬时的协方差函数来确定,若噪声为平稳的,则平均值与时间无关,而协方差函数则变成 ...

最新文章

  1. python3里的pillow怎么安装_“python安装pillow教程“python3.4怎么安装pil
  2. controll层跳转页面_SpringMVC controller控制页面跳转 | 学步园
  3. 009_JavaScript语句、表达式和代码块
  4. php创建表并插入数据,php数据库操作-创建库和表以及插入数据
  5. 微型计算机组装实验报告虚拟,计算机硬件的组装实验报告.doc
  6. java arraylist 重复_Java中ArrayList去除重复元素
  7. 某易游戏经典吃豆豆动画404页面源码
  8. Python模拟大整数乘法的小学竖式计算过程
  9. 《论语》《大学》《中庸》和孟子
  10. linux获取utc时间函数,Linux下的时间函数
  11. 高精度WGS84与GCJ-02坐标转换
  12. excel取末尾数字_excel表格中数字末尾的0不显示怎么设置
  13. 发送邮件 显示对方服务器未响应,邮件对方服务器未响应
  14. 软件测试——三角形问题测试用例练习
  15. python生产式系统动物识别_在MES生产管理系统中,人脸识别技术以什么形式存在?...
  16. 智能语音计算器(四)
  17. 3D点云的快速分割:自动驾驶汽车应用的LiDAR处理实例
  18. 2022 版史上最全 Java 八股文,没有任何异议
  19. Sicily 1048. Inverso
  20. 解决centos8 中安装mysql5.7.26时 libtinfo.so.5()(64bit)被依赖的问题

热门文章

  1. 非常好用的综合源码资源分享平台推荐,用它就够了
  2. 润乾报表 表头斜线 怎么设置
  3. BUUCTF-模板注入专项刷题
  4. PHP开发者必备的50个库/框架【2019】
  5. 134_容器_git工具_Google_guava与Apache_commons的准备工作(jar src doc)
  6. APS推进实施的条件与基础
  7. Android 5.0 Phone 模块介绍
  8. 如何看待 Python (PyScript) 可以做浏览器前端开发语言?
  9. Quartus II错误集
  10. python:遥感时间序列处理——Pettitt突变点检测(突变时间/年份)