西电电子信息系统综合实验一

已于 2023-10-27 09:21:44 修改
阅读量154 收藏 1
点赞数
分类专栏： matlab 数字信号处理 文章标签： 算法
于 2023-10-26 17:15:17 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/MarkAssassin/article/details/134060787
版权

实验任务

1.模拟不同白噪声（均匀、高斯、瑞利）

2.模拟雷达脉冲信号（单载波、BPSK、LFM）

3.模拟通信信号（QPSK、BPSK）

4.噪声与信号混合

5.混合信号频谱分析

注意：代码需一部分一部分进行运行

%% 模拟不同白噪声（均匀、高斯、瑞利）
%均匀白噪声
clear all;
close all;
clc;
a=2; %均匀分布下限
b=3; %均匀分布上限
fs=1e7; %采样率
t=1e-3; %随机序列长度
n=t*fs;
rand('state',0); %把均匀分布伪随机发生器置为0状态
u=rand(1,n); %产生0-1单位均匀信号
noise1=(b-a)*u+a; %广义均匀分布与单位均匀分布之间的关系
nexttile;
plot(noise1);
title('均匀分布噪声');
nexttile;
histogram(noise1,a:0.02:b); %直方图
title('均匀分布噪声直方图');

%高斯白噪声
fs=1e7; %采样率
t=1e-3; %随机序列长度
n=t*fs;
randn('state',0);
noise2=randn(1,n);
nexttile;
plot(noise2);
title('高斯分布噪声');
nexttile;
histogram(noise2);
title('高斯分布噪声直方图');

%瑞利白噪声
sigma=2; %瑞利分布参数sigma
t=1e-3; 
fs=1e7;
t1=0:1/fs:t-1/fs;
n=length(t1);
rand('state',0);
u=rand(1,n);
noise3=sqrt(2*log2(1./u))*sigma; %产生瑞利分布信号
nexttile;
plot(noise3);
title('瑞利分布噪声');
xlabel('t/s');
nexttile;
histogram(noise3,0:0.1:10);
title('瑞利分布噪声直方图')


%% 模拟雷达脉冲信号（单载波、BPSK、LFM）
%雷达脉冲信号单载波
%雷达脉冲信号单载波

T=1;
f=1e5;
t=0:1/f:(T-1/f);
n=length(t);

A_danpin=10;
f_danpin=1000;
y_danpin=A_danpin*cos(2*pi*f_danpin*t);
figure(1);
plot(t(1:1000),y_danpin(1:1000))
xlabel('t/s')
title('单载波');
grid on;
axis([0 0.01 -10.5 10.5])





%雷达脉冲信号BPSK

T=1;
f=1e5;
t=0:1/f:(T-1/f);
n= length(t);

binary_code='1001100110';
len_code=length(binary_code);

f_bpsk=1000;
T_bpsk=1/f_bpsk;
A_bpsk=10;
t_mayuan=0:1/f:(T_bpsk-1/f);
y_mayuan=A_bpsk*sin(2*pi*f_bpsk*t_mayuan);
t_bpsk=0:1/f:(T_bpsk*len_code-1/f);
y_bpsk= [];
for i=1:len_code
    tmp=str2num(binary_code(i));
    if tmp==0
        y_bpsk=[y_bpsk,y_mayuan];
    end
    if tmp==1
        y_bpsk=[y_bpsk,-y_mayuan];
    end
end
figure(2)
plot(t_bpsk,y_bpsk)
grid on;
xlabel('t/s');
title('BPSK信号');
axis([0 1e-2 -2*A_bpsk 2*A_bpsk]);

%雷达脉冲信号LFM

% T=1;
% f=1e5;
% t=0:1/f:(T-1/f);
% n=length(t);
% A_lfm=1;
% f_lfm=1000;
% k_lfm=800000;
% y_lfm=A_lfm*cos(2*pi*f_lfm*t+pi*k_lfm*t.^2);
% figure(3)
% plot(t(1:500),y_lfm(1:500));
% xlabel('t/s');
% title('LFM');
% axis([0 0.005 -A_lfm A_lfm]);

% 线性调频与脉冲压缩
set(0,'defaultfigurecolor','w')
% Chirp信号参数设置
Tr = 1e-6;%时宽
Br = 200e6;%带宽
Fs = 4*Br;%采样率
% Chirp信号参数导出
Kr = Br/Tr;%调频率
N =  round( Tr / (1/Fs) );%采样点数
t = linspace( -Tr/2 , Tr/2 , N);%在[-Tp/2,Tp/2]选取采样点
% Chirp信号生成
y_lfm = ( abs(t) < Tr/2 ) .* exp( 1j * pi * Kr * t.^2 ); 
f_chirp= Kr * t; %信号频率
phase_chirp = pi * Kr * t.^2;%信号相位
plot(t*1e6,y_lfm)
title('雷达LFM')


%% 模拟通信信号（QPSK、BPSK）
%QPSK
N=20;%比特数
T=1;%比特周期
fc=2;%载波频率
Fs=100;%抽样频率
bitstream=randi([0,1],1,N);%随机产生的比特数0、1
bitstream=2*bitstream-1;%单极性变为双极性（0到-1；1到1）
I=[];Q=[];
%奇数进I路,偶数进Q路
for i=1:N
    if mod(i,2)~=0
        I=[I,bitstream(i)];
    else
        Q=[Q,bitstream(i)];
    end
end
%采用绘图比较I、Q比特流
bit_data=[];
for i=1:N
    bit_data=[bit_data,bitstream(i)*ones(1,T*Fs)];%在一个比特周期里面有T*Fs个1和采样点一模一样
end
I_data=[];Q_data=[];
for i=1:N/2
    %I路和Q路是原来比特周期的两倍,2Tb=Ts(码元周期)，因此采样点个数为T*Fs*2
    I_data=[I_data,I(i)*ones(1,T*Fs*2)];
    Q_data=[Q_data,Q(i)*ones(1,T*Fs*2)];
end
%绘图
figure();
%时间轴
t=0:1/Fs:N*T-1/Fs;
subplot(3,1,1)
plot(t,bit_data);legend('Bitstream')%比特信息
subplot(3,1,2)
plot(t,I_data);legend('I Bitstream')%I路信息
subplot(3,1,3)
plot(t,Q_data);legend('Q Bitstream')%Q路信息
%载波信号
bit_t=0:1/Fs:2*T-1/Fs;%载波周期为2倍比特周期,定义时间轴
%定义I路和Q路的载波
I_carrier=[];Q_carrier=[];
for i=1:N/2
    I_carrier=[I_carrier,I(i)*cos(2*pi*fc*bit_t)];%I路载波信号
    Q_carrier=[Q_carrier,Q(i)*cos(2*pi*fc*bit_t+pi/2)];%Q路载波信号
end
%传输信号
QPSK_signal=I_carrier+Q_carrier;
%绘图
figure();%产生一个新图
subplot(3,1,1)
plot(t,I_carrier);legend('I signal')%I路信号
title('I路信号');
subplot(3,1,2)
plot(t,Q_carrier);legend('Q signal')%Q路信号
title('Q路信号')
subplot(3,1,3)
plot(t,QPSK_signal);legend('QPSK signal')%I路、Q路和的信号
title('I+Q路信号(QPSK)')


%BPSK
% N=1000; %信号长度
% fc=10; %载波频率
% fs=1000; %采样频率
% T=1/fs; %采样周期
% SNR=10; %信噪比
% 
% 
% %生成BPSK信号
% x=sign(randn(1,N)); %产生随机二进制序列
% BPSK_signal=x.*cos(2*fc*pi*(0:N-1)*T); %BPSK信号
% plot(BPSK_signal)
% title('BPSK信号')
% xlabel('t/s')
T=1;
f=1e5;
t=0:1/f:(T-1/f);
n= length(t);

binary_code='1010101010';
len_code=length(binary_code);

f_bpsk=1000;
T_bpsk=1/f_bpsk;
A_bpsk=10;
t_mayuan=0:1/f:(T_bpsk-1/f);
y_mayuan=A_bpsk*sin(2*pi*f_bpsk*t_mayuan);
t_bpsk=0:1/f:(T_bpsk*len_code-1/f);
BPSK_signal= [];
for i=1:len_code
    tmp=str2num(binary_code(i));
    if tmp==0
        BPSK_signal=[BPSK_signal,y_mayuan];
    end
    if tmp==1
        BPSK_signal=[BPSK_signal,-y_mayuan];
    end
end
figure(2)
plot(t_bpsk,BPSK_signal)
grid on;
xlabel('t/s');
title('BPSK信号');
axis([0 1e-2 -2*A_bpsk 2*A_bpsk]);


%% 实测通信信号（RadioML2018.01A）


%% 噪声与信号混合+频谱分析
% 雷达单载频+均匀白噪声
T=1;
f=1e5;
t=0:1/f:(T-1/f);
n=length(t);
noise_1=rand(1,n); %2倍扩大使结果更明显
signal_1=y_danpin+noise_1;
nexttile;
plot(t(1:1000),signal_1(1:1000));
title('雷达单载频+均匀白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_1))))
title('雷达单载频+均匀白噪声频谱图')
xlabel('f/Hz')

%雷达单载频+高斯白噪声
noise_2=randn(1,n);
signal_2=y_danpin+noise_2;
nexttile;
plot(t(1:1000),signal_2(1:1000));
title('雷达单载频+高斯白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_2))))
title('雷达单载频+高斯白噪声频谱图')
xlabel('f/Hz')

%雷达单载频+瑞利白噪声
sigma=3;
u=rand(1,n);
noise_3=sqrt(2*log2(1./u))*sigma;
signal_3=y_danpin+2*noise_3;
nexttile;
plot(t(1:1000),signal_3(1:1000));
title('雷达单载频+瑞利白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_3))))
title('雷达单载频+瑞利白噪声频谱图')
xlabel('f/Hz')


%雷达BPSK+均匀白噪声
nexttile;
signal_4=y_bpsk+noise_1(1:1000);
plot(t(1:1000),signal_4(1:1000));
title('雷达BPSK+均匀白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_4))))
title('雷达BPSK+均匀白噪声频谱图')
xlabel('f/Hz')

%雷达BPSK+高斯白噪声
nexttile;
signal_5=y_bpsk+noise_2(1:1000);
plot(t(1:1000),signal_5(1:1000));
title('雷达BPSK+高斯白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_5))))
title('雷达BPSK+高斯白噪声频谱图')
xlabel('f/Hz')


%雷达BPSK+瑞利白噪声
nexttile;
signal_6=y_bpsk+noise_3(1:1000);
plot(t(1:1000),signal_6(1:1000));
title('雷达BPSK+瑞利白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_6))))
title('雷达BPSK+瑞利白噪声频谱图')
xlabel('f/Hz')

%雷达LFM+均匀白噪声
nexttile;
signal_7=y_lfm+noise_1(1:800);
plot(t(1:500)   ,signal_7(1:500));
title('雷达LFM+均匀白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_7))))
title('雷达LFM+均匀白噪声频谱图')
xlabel('f/Hz')

%雷达LFM+高斯白噪声
nexttile;
signal_8=y_lfm+noise_2(1:800);
plot(t(1:500),signal_8(1:500));
title('雷达LFM+高斯白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_8))))
title('雷达LFM+高斯白噪声频谱图')
xlabel('f/Hz')


%雷达LFM+瑞利白噪声
nexttile;
signal_9=y_lfm+noise_3(1:800);
plot(t(1:500),signal_9(1:500));
title('雷达LFM+瑞利白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_9))))
title('雷达LFM+瑞利白噪声频谱图')
xlabel('f/Hz')

%模拟通信信号QPSK+均匀白噪声
nexttile;
signal_10=QPSK_signal+noise_1(1:2000);
plot(t(1:1000),signal_10(1:1000));
title('QPSK+均匀白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_10))))
title('QPSK+均匀白噪声频谱图')
xlabel('f/Hz')

%模拟通信信号QPSK+高斯白噪声
nexttile;
signal_11=QPSK_signal+noise_2(1:2000);
plot(t(1:1000),signal_11(1:1000));
title('QPSK+高斯白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_11))))
title('QPSK+高斯白噪声频谱图')
xlabel('f/Hz')

%模拟通信信号QPSK+瑞利白噪声
nexttile;
signal_12=QPSK_signal+noise_3(1:2000);
plot(t(1:1000),signal_12(1:1000));
title('QPSK+瑞利白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_12))))
title('QPSK+瑞利白噪声频谱图')
xlabel('f/Hz')

%模拟通信信号BPSK+均匀白噪声
nexttile;
signal_13=BPSK_signal+noise_1(1:1000);
plot(t(1:500),signal_13(1:500));
title('BPSK+均匀白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_13))))
title('BPSK+均匀白噪声频谱图')
xlabel('f/Hz')

%模拟通信信号BPSK+高斯白噪声
nexttile;
signal_14=BPSK_signal+noise_2(1:1000);
plot(t(1:500),signal_14(1:500));
title('BPSK+高斯白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_14))))
title('BPSK+高斯白噪声频谱图')
xlabel('f/Hz')

%模拟通信信号BPSK+瑞利白噪声
nexttile;
signal_15=BPSK_signal+noise_3(1:1000);
plot(t(1:500),signal_15(1:500));
title('BPSK+瑞利白噪声')
xlabel('t/s')
nexttile;
plot(fftshift(abs(fft(signal_15))))
title('BPSK+瑞利白噪声频谱图')
xlabel('f/Hz')

MarkAssassin
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
西安电子科技大学基础实验答案
10-25
西安电子科技大学基础实验答案！！！保真！！！！
通信与信息系统综合应用实验.zip
01-11
西安电子科技大学通信工程学院大四通信与信息系统综合应用实验实验报告，射频传输，网络协议，网络通信，数字基带模块。
西电电子信息系统综合实验三、四
MarkAssassin的博客
11-20 119
每年每个班的实验题目会有些许不同，仅供参考。
西电电子信息系统综合实验二
MarkAssassin的博客
10-26 171
西安电子科技大学电子信息系统综合实验二
电子科技大学现代电子信息系统综合实验课程部分模拟试题
thompson的博客
08-05 3998
【考试一，只有0】 EDA: 拨码A,亮/灭最左两位LED (5分), 拨码B起/停最右两位LED闪烁（频率2Hz）（5分）;仿真拨码B支路的电路功能（5分）；基于前功能在右LED亮时，数码管只最左1位显示该led闪烁次数（5分）； MCU: 按键A,亮/灭左LED（5分）， 按键B起/停最右LED闪烁（频率2Hz）（5分）；选最左第2位数码管，在两LED同时灭时才显示“0”（5分）；继续按键时，数码管左循环一位显示该"0"（5分） 【考试二,AAA，88循环加一】 EDA E盘带AAA 指定目录建立项目（
通信系统综合实验报告实验报告.doc
03-11
班级 研01-12 学号 西安电子科技大学 通信系统综合实验报告 学 院 ： 通信工程学院 班 级 ： 研01-12 专 业 ： 交通信息工程及控制 姓 名 ： 2012年11月 目录 实验一 数字基带仿真实验 1 1.1 实验目的 1 1.2 实验...
电子信息类专业应用型创新人才培养模式研究.pdf
08-24
一些重点高校实施了国家级电子信息类创新型人才培养实验计划，如北京邮电大学的拔尖创新人才培养模式改革实验区和西安电子科技大学的电子信息类人才培养模式创新实验区。然而，目前的研究多侧重于传统课程体系、教学...
电子信息类专业课程体系改革的探索与实践.pdf
08-24
西安电子科技大学作为以信息与电子学科为主的高校，电子信息类专业的建设对其发展至关重要。随着电子技术和信息技术的快速发展，社会对这类专业人才的需求也在不断变化，因此，教育模式也需要相应调整。 改革的基本...
章需求分析与项目管理-西安电子科技大学-软件工程教学课件(与“项目”有关文档共104张).pptx
11-14
"需求分析与项目管理-西安电子科技大学-软件工程教学课件" 需求分析是软件工程中的一个重要任务，它是系统层次软件配置与软件设计之间的桥梁。需求分析员需要努力做好以下工作：研究系统规格说明和项目计划，确保对...
算法训练营day60
stmfresher的博客
06-16 206
【代码】算法训练营day60。
P1656 炸铁路
STY_fish_2012的博客
06-19 325
如果y的追溯值大于了x的时间戳，也就是y的追溯值在x之后，并且，y追溯到的这个点与x一定不连通，因为如果连通，那么y就会顺着这个点追溯到时间戳更小的点x，到这里，我们不难推出：x与y之间直接连通的边，只有这一条，并且没有间接联通的边。炸毁一条路对于一张图来说，本质上就是去掉一条边，不能通过铁路到达就代表这张图不连通了，对于这种去掉之后就可以使整张图不再连通的边，我们称之为割边。第一行 n,m (1≤n≤150，1≤m≤5000)，分别表示有 n 个城市，总共 m 条铁路。里面里面讲过，这里就不多说了。
【代码随想录】【算法训练营】【第35天】[134]加油站 [135]分发糖果 [860]柠檬水找零 [406]根据身高重建队列
weixin_54954007的博客
06-11 316
思路及算法思维，指路。题目来自。day 35，连休两天~
华为OD机试C卷（100分）-最大N个数与最小N个数的和（C语言）
欢迎来到我的博客~
06-16 182
给定一个数组，编写一个函数来计算它的最大N个数与最小N个数的和。你需要对数组进行去重。说明 最大2个数[100,95],最小2个数[83,64], 输出为342。说明 最大2个数[4,3],最小2个数[3,2], 有重叠输出为-1。最大N个数与最小N个数不能有重叠，如有重叠，输入非法返回-1。第三行输入N，N表达需要计算的最大、最小N个数。数组中数字范围[0, 1000]第一行输入M， M标识数组大小。第二行输入M个数，标识数组内容。输出最大N个数与最小N个数的和。
C语言中数组和指针的关系
2402_85246552的博客
06-12 601
因此，在函数内部，你无法知道数组的大小（除非作为另一个参数传递），但你可以通过指针访问数组的元素。对于数组，如果越界访问（访问数组索引超出其实际大小的元素），可能会导致未定义的行为，包括访问到无效的内存地址或覆盖其他变量的值。对于指针，如果指针指向了无效的内存地址（如未初始化的指针、已释放的内存或不存在的内存地址），解引用该指针也会导致未定义的行为。数组的大小在编译时是固定的，而指针可以指向任意大小的内存区域（只要该区域已被分配且指针已正确初始化）。），这会得到整个数组的地址或指向下一个元素的指针。
C数据结构：排序
2301_80655639的博客
06-12 1040
若按照错误方法来：例如：begin = 0，end = 3，那么mid = 1，begin1 = 0，end1 = 2，递归左右区间为[0，2]也就是begin = 0，end = 2，那么mid = 1，begin1 = 0，end1 = 2，这样又递归就会出现栈溢出的情况，所以必须要按照 上面代码的方式走。然后这个for循环做的就是遍历原数组，将原数组的值减去最小值，这个值映射到count数组中，让这个下标的值++，这样就会统计这个值出现的次数，而与原数组建立的关系就是减去一个min的关系。
算法与数据结构--决策树算法
考研小钻风
06-16 1261
决策树算法在数据科学和机器学习中具有重要地位，理解其原理和实现方法对解决实际问题具有重要意义。
位运算算法
最新发布
mgzdwm的博客
06-19 176
位运算在算法中有着广泛的应用，例如可以用位运算实现整数的加减乘除运算，判断一个数的二进制表示中有多少个1，判断一个数是否是2的幂次方等。位运算是计算机中常用的一种运算方法，它直接对二进制数的位进行操作。按位或（|）：对两个二进制数的每一位进行或运算，如果两个对应位中有一个为1，则结果对应位为1，否则为0。按位异或（^）：对两个二进制数的每一位进行异或运算，如果两个对应位不相同，则结果对应位为1，否则为0。按位与（&）：对两个二进制数的每一位进行与运算，如果两个对应位都为1，则结果对应位为1，否则为0。
堆优化版Dijkstra求最短路-java
m0_63267251的博客
06-19 707
主要通过堆优化Dijkstra算法解决最短路，可以跟朴素版的Dijkstra算法进行对比。
西安电子科技大学 编译原理 实验 sql解释器
07-05
西安电子科技大学编译原理实验SQL解释器是一个用于解析和执行SQL语句的工具。SQL是一种结构化查询语言，用于管理和操作关系型数据库中的数据。 该解释器的主要功能是将用户输入的SQL语句进行解析和执行。首先，它会...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值