【本科生通信原理】【实验报告】【北京航空航天大学】实验一:通信原理初步

一、实验目的:

  1. 熟悉 MATLAB开发环境、掌握 MATLAB基本运算操作;
  2. 熟悉和了解 MATLAB图形绘制基本指令;
  3. 熟悉使用 MATLAB分析信号频谱的过程;
  4. 掌握加性白高斯噪声信道模型

二、实验内容:

在这里插入图片描述

三、实验程序:
1、

function q1()
x = 0 : 0.0001 : 2 * pi;
y1 = 2 * exp(-0.5 * x);
y2 = cos(4 * pi * x);
figure;
plot(x, y1);
hold;
plot(x, y2);
xlabel("x", "FontName", "Times New Roman", "FontSize", 12);
ylabel("y", "FontName", "Times New Roman", "FontSize", 12);
legend("FontName","Times New Roman", "FontSize", 10, "LineWidth", 1.5);
legend("y1", "y2");
title("q1");
axis([0 2 * pi -1.1 2.1]);
grid;

2、

function q2()
x = 0 : 0.01 : 10;  % 横坐标区间
R = normrnd(0, sqrt(0.1), 1, length(x));  % 生成1 * length(x)个正态随机数figure;  % 画信号波形
plot(x, R);
grid on;
title('White Gaussian Noise');
xlabel('x');
ylabel('N');figure;  % 画序列柱状图
bar(R);
grid on;
title('Bar Graph of Noise Sequence');
xlabel('n');
ylabel('N');

3、

function q3()
N = 1024;  %采样点数
ts = 0.6 / 1023;  %系统时域采样间隔
fs = 1 / ts;  %系统采样频率
df = 0.001;  %所需的频率分辨率
t = 0 : ts: 0.6;
x = 0.4 * sin(100 * pi * t) + 0.4 * sin(640 * pi * t);
n = randn(1, N);  % 噪声信号(噪声方差为1)
y = x + n;  % 原始信号叠加噪声信号
[Y, m, df1, f] = T2F(y, ts, df, fs);figure;  % 画信号时域波形图
plot(t, y);
grid on;
title('信号时域波形图');
xlabel('t/s');
ylabel('y(t)');figure;  % 画信号频谱图
plot(f, abs(Y));
grid on;
title('信号频谱图');
xlabel('w');
ylabel('Y(w)');

4、

function q4()
% AM调制解调
echo on;
N = 1024;  % 采样点数
A = 3;  % 直流分量
fc = 125;  %载波频率
t0 = 0.6;  %信号持续时间
snr = 10;  %解调器输入信噪比dB
dt = 0.6 / 1023;  %系统时域采样间隔
fs = 1 / dt;  %系统采样频率
df = 0.001;  %所需的频率分辨率
t = 0 : dt : t0;
m = 0.1 * cos(15 * pi * t) + 1.5 * sin(25 * pi * t) + 0.5 * cos(40 * pi * t);  %调制信号
c = cos(250 * pi * t);  %载波
Lt = length(t);  %仿真过程中,信号长度
snr_lin = 10 ^ (snr / 10);  %信噪比
L = 2 * min(m);
R = 2 * max(abs(m)) + A;
[M, m, df1, f] = T2F(m, dt, df, fs);  %求出调制信号频谱
[Bw_eq] = signalband(M, df, t0);  %求出信号等效带宽
u = (A + m(1 : Lt)) .* c(1 : Lt);  % 已调信号
[U, u, df1, f] = T2F(u, dt, df, fs);
signal_power = power_x(u(1 : Lt));  %已调信号的平均功率
noise_power = (signal_power * fs) / (snr_lin * (2 * Bw_eq));  %求出噪声方差(噪声均值为0)
noise_std = sqrt(noise_power);  %噪声标准差
noise = noise_std * randn(1, Lt);  %产生噪声
sam = u(1 : Lt) + noise(1 : Lt);  %叠加了噪声的已调信号
[SAM, sam, df1, f] = T2F(sam, dt, df, fs);  %求出叠加了噪声的已调信号频谱figure;  % 画出经过信道前的已调信号时域波形
plot(t, u(1 : length(t)));
grid on;
title('经过AWGN信道前的已调信号的时域波形图');
xlabel('t');
ylabel('u(t)');figure;  %画出经过信道前的已调信号频谱图
plot(f, abs(fftshift(U)));
grid on;
title('经过AWGN信道前的已调信号的频谱图');
xlabel('w');
ylabel('U(w)');figure;  %画出经过信道后的已调信号时域波形
plot(t, sam(1 : length(t)));
axis([0 t0 -20 20]);
grid on;
title('经过AWGN信道后的已调信号的时域波形图');
xlabel('t');
ylabel('s(t)');figure;  %画出经过信道后的已调信号频谱图
plot(f, abs(fftshift(SAM)));
grid on;
title('经过AWGN信道后的已调信号的时域波形图');
xlabel('w');
ylabel('S(w)');

四、实验结果:
1、
在这里插入图片描述

2、
(1)、信号波形:
在这里插入图片描述

(2)、序列柱状图:
在这里插入图片描述

3、
(1)、所得信号的时域波形图:
在这里插入图片描述

(2)、所得信号的频谱图:
在这里插入图片描述

4、
(1)、经过AWGN信道前的已调信号的时域波形图:
在这里插入图片描述

(2)、经过AWGN信道前的已调信号的频谱:
在这里插入图片描述

(3)、经过AWGN信道后的已调信号的时域波形图:
在这里插入图片描述

(4)、经过AWGN信道后的已调信号的频谱图:
在这里插入图片描述

五、实验分析:
1、求离散时间信号傅里叶变换:

function [f, sf] = F(t, st)
%利用fft,fftshift定义函数F计算信号的傅里叶变换
%t-离散时间
%st-离散信号
dt = t(2) - t(1); %时间分辨率
T = t(end);
df = 1 / T; %频率分辨率
N = length(st); %离散傅里叶变换长度
f = -N / 2 * df : df : N / 2 * df - df;
sf = fft(st);
sf = T / N * fftshift(sf);

2、求信号平均功率:

function p = power_x(x)
%x:输入信号
%p:返回信号的x功率
p = (norm(x) .^ 2) ./ length(x);

3、求信号等效带宽:

function [Bw_eq] = signalband(sf, df, T)
%计算信号等效带宽
%sf:信号频谱
%df:频谱分辨率
%T:信号持续时间
sf_max = max(abs(sf));
Bw_eq = sum(abs(sf) .^ 2) * df / T / sf_max .^ 2;

4、序列的傅里叶变换

function [M, m, df] = fftseq(m, ts, df)fs = 1 / ts;
if nargin == 2n1 = 0;
elsen1 = fs / df;
end
n2 = length(m);
n = 2 ^ (max(nextpow2(n1), nextpow2(n2)));
M = fft(m, n);
m = [m, zeros(1, n - n2)];
df = fs / n;

5、信号从时域转换到频域

function [M, m, df1, f]=T2F(m, ts, df, fs)
%-----------------输入参数
%m: 信号
%ts: 系统时域采样间隔
%df: 所需的频率分辨率
%fs: 系统采样频率%-----------------输出参数
%M: 傅里叶变换后的频谱序列
%m: 输入信号参与过傅里叶变换后对应的序列,补零后的输入信号,长度与M,f相同
%df1: 返回的频率分辨率
%f: 与M相对应的频率序列
[M, m, df1] = fftseq(m, ts, df);
f = [0 : df1 : df1 * (length(m) - 1)] - fs / 2;  %频率向量
M = M / fs;

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.rhkb.cn/news/231951.html

如若内容造成侵权/违法违规/事实不符,请联系长河编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

原生微信小程序如何动态修改svg图片颜色及尺寸、宽高(封装svgIcon组件)解决ios不显示问题

最终效果 前言 动态设置Svg图片颜色就是修改Svg源码的path中的fill属性, 通过wx.getFileSystemManager().readFile读取.xlsx文件 ios不显示需要把encoding设置 binary 把文件转成base64 封装svg-icon组件 1、在项目的components下新建svg-icon文件夹,新…

iec104和iec61850

iec104和iec61850 IEC104 规约详细解读(一) 协议结构 IEC104 规约详细解读(二)交互流程以及协议解析 61850开发知识总结与分享【1】 Get the necesarry projects next to each other in the same directory; $ git clone https://github.com/robidev/iec61850_open_server.g…

基于单片机的农田灌溉系统(论文+源码)

1.系统设计 本系统主要实现如下目标: 1.可以实时监测土壤湿度; 2.土壤湿度太低时,进行浇水操作; 3.可以按键设置湿度的触发阈值; 4. 可以实现远程操控 5.可以实现手…

wsl(ubuntu)创建用户

我们打卡ubuntu窗口,如果没有创建用户,那么默认是root用户 用户的增删改查 查 查询所有的用户列表 cat /etc/passwd | cut -d: -f1cat /etc/passwd: 这个命令用于显示 /etc/passwd 文件的内容。/etc/passwd 文件包含了系统上所有用户的基本信息。每一…

低代码开发会取代传统开发吗? 两者有什么区别 该如何选择

低代码开发技术在近几年逐渐被普及,帮助很大一部分开发者完成了复杂的工作。由于低代码开发方案入门门槛低且上手难度小,所以即使是非专业人士也可借助其便利性自主开发软件系统,整个开发过程几乎不需要专业程序员。久而久之就出现了一种声音…

C++上位软件通过Snap7开源库访问西门子S7-1200/S7-1500数据块的方法

前言 本人一直从事C上位软件开发工作较多,在之前的项目中通过C访问西门子PLC S7-200/S7-1200/S7-1500并进行数据交互的应用中一直使用的是ModbusTCP/ModbusRTU协议进行。Modbus上位开源库采用的LibModbus。经过实际应用发现Modbus开源库单次发送和接受的数据不能超过…

静态网页设计——红旗汽车官网(HTML+CSS+JavaScript)

前言 声明:该文章只是做技术分享,若侵权请联系我删除。!! 感谢大佬的视频: https://www.bilibili.com/video/BV1gK411x7Bg/?vd_source5f425e0074a7f92921f53ab87712357b 使用技术:HTMLCSSJS(…

STM32CubeMX教程13 ADC - 单通道转换

目录 1、准备材料 2、实验目标 3、ADC概述 4、实验流程 4.0、前提知识 4.1、CubeMX相关配置 4.1.1、时钟树配置 4.1.2、外设参数配置 4.1.3、外设中断配置 4.2、生成代码 4.2.1、外设初始化调用流程 4.2.2、外设中断调用流程 4.2.3、添加其他必要代码 5、常用函数…

FreeRTOS 实时操作系统第十二讲 - 计数信号量

一、信号量的概念 1、信号量的基本概念 消息队列是实现任务与任务或任务与中断间通信的数据结构,可类比裸机编程中的数组 信号量是实现任务与任务或任务与中断间通信的机制,可以类比裸机编程中的标志位 信号量 (semaphore) 可以实现任务与任务或任务与中…

maven、springboot项目编译打包本地jar、第三方jar包

0. 引言 一般我们在maven项目中都是通过引入pom坐标的形式来引入第三方jar包,但某些场景下,第三方是直接提供的jar包文件,这就需要我们从本地引入第三方包并进行打包。所以我们今天来看下如何进行本地引入第三方包操作 1. 步骤 1、在项目下…

GNSS位移监测站对尾矿库坝体表面位移进行自动化监测

表面位移监测:通过GNSS位移监测站对尾矿库坝体表面位移进行自动化监测,掌握尾矿坝整体表面位置的变化及其变化速率(包括平面位移和垂直沉降),确定尾矿坝坝体整体位移变形的情况,是确定尾矿库安全性的重要指…

基于Java+SpringBoot+Vue志愿者服务平台设计和实现

博主介绍:✌全网粉丝30W,csdn特邀作者、博客专家、CSDN新星计划导师、Java领域优质创作者,博客之星、掘金/华为云/阿里云/InfoQ等平台优质作者、专注于Java技术领域和学生毕业项目实战,高校老师/讲师/同行交流合作✌ 主要内容:SpringBoot、Vue、SSM、HLM…

[技术杂谈]使用VLC将视频转成一个可循环rtsp流

通过vlc播放器,将一个视频转成rtsp流,搭建一个rtsp服务器。rtsp客户端可访问这个视频的rtsp流。 1. 打开vlc播放器,使用的版本如下 2. 菜单:媒体 ---> 流 3. 添加视频文件,点击添加一个mp4 文件 4. 选择串流&…

Servlet 3.0的异步处理

1、传统Servlet处理 Web容器会为每个请求分配一个线程,默认情况下,响应完成前,该线程占用的资源都不会被释放。若有些请求需要长时间(例如长处理时间运算、等待某个资源),就会长时间占用线程所需资源,若这类请求很多&…

Spring中的数据校验

文章目录 引言摘要正文基于 ValidationUtils的简单校验基于自定义 Validator的校验Spring内置校验 LocalValidatorFactoryBeanHibernateValidator校验使用HibernateValidator自定义校验规则 总结 引言 我们在日常的软件开发过程中,尤其是WEB开发过程中,…

antd Table 动态数据 合并单元格(合并行)

antd Table 组件动态合并单元格 最近处理table的时候 遇到了要合并同一列的几行的情况,比如第一列的前面三行都是同一个对象的名字,此时合并显示比较妥当,但是数据是后端接口来的,而且可以筛选条件,搜索出来的数据就是…

Linux的ping命令、wget命令、curl命令

一、ping命令 通过ping命令,可以检查指定的网络服务器是否是可联通状态 形式:ping [-c num] ip或主机名 -c:检查的次数,不使用-c,将无限次数持续检查 ip或主机名:被检查的服务器的ip地址或主机名地址 …

Jenkins分布式实现: 构建弹性和可扩展的CI/CD环境!

Jenkins是一个流行的开源持续集成(Continuous Integration,CI)和持续交付(Continuous Delivery,CD)工具,它通过自动化构建、测试和部署过程,帮助开发团队更高效地交付软件。Jenkins的…

AtomicReference自旋加CAS保证对象引用原子性

AtomicReference类提供了对象引用的非阻塞原子性读写操作,并且提供了其他一些高级的用法,对象的引用其实是一个4字节的数字,代表着在JVM堆内存中的引用地址,对一个4字节数字的读取操作和写入操作本身就是原子性的 package Atomic…

数据结构入门到入土——ArrayList与顺序表

目录 一,线性表 二,顺序表 1.接口实现 三,ArrayList简介 四,ArrayList使用 1.ArrayList的构造 2.ArrayList常见操作 3.ArrayList的遍历 4.ArrayList的扩容机制 五,ArrayLisit的具体使用 杨辉三角 一&#x…