matlab实现0,1序列2FSK调制

最新推荐文章于 2024-05-04 18:37:11 发布
最新推荐文章于 2024-05-04 18:37:11 发布
阅读量9k 收藏 215
点赞数 38
文章标签： matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/enduringzx/article/details/122185896
版权

matlab实现0,1简单序列2FSK调制

0,1简单序列2FSK调制

前言

本代码是最简单的0,1序列2FSK调制，未做同步与均衡，概念性问题网上能随便查到，不做过多叙述。

代码

参数配置

// 参数配置
fs=3840000;        % 采样率配为3840kHz，为128000*30所得，且30720000/8
sample_num=30;         %一个符号采样点数为30
Rb=fs/sample_num;       %符号速率，128k bps

len_symbol=2500;                %符号个数，可计算总样点数为N=len_sambol*sample_num,
bits_per_sybol=1;           %2FSK一个符号对应1个比特
len=len_symbol*bits_per_sybol;           %数据源比特个数
f1=Rb*4;          %载波频率，单位Hz,载波频率配置需为Rb的整数倍,同时能被采样率整除，为128kHz
f2=Rb*2;

N=len_symbol*sample_num*100             %符号采样点数
dt=1/fs;
t0=0:dt:(N*2-1)*dt;
t=0:dt:(N/30-1)*dt;%符号时间轴

采样率依照3480kHz来设置，这是通信类学生符合实验中硬件类需求的标准。

生成数据源

%% 生成数据源
st1=randi([0,1],1,len);%随机生成0,1序列
st2=-st1+1;%生成基带信号反码，1变成0,0变成1

g11=(ones(1,100))'*st1;  %产生方波信号 
g1a=g11(:)'; 
g12=(ones(1,100))'*st2;  %产生方波信号 
g2a=g12(:)'; 

figure(1);
subplot(211);
plot(g1a);
title('基带信号st1');
axis([0,500,-1,2]);
subplot(212);
plot(g2a);
title('基带信号反码st2');
axis([0,500,-1,2]);

在这里插入图片描述

方波信号的设置是为了让结果显得更加清晰，如果觉得不需要也可去掉方波信号设置，直接使用产生的随机序列进行还原亦可，只是效果会不尽理想。

设置载波

%载波信号
s1=cos(2*pi*f1*t);
s2=cos(2*pi*f2*t);
figure(2)
subplot(211),plot(s1);
axis([0 500 -1.5 1.5]);
title('载波信号s1');
subplot(212),plot(s2);
axis([0 500 -1.5 1.5]);
title('载波信号s2');

在这里插入图片描述

调制

F1=g1a.*s1;%加入载波1
F2=g2a.*s2;%加入载波2
figure(3);
subplot(411);
plot(F1);
axis([0 500 -1.5 1.5]);
title('F1=s1*st1');
subplot(412);
plot(F2);
axis([0 500 -1.5 1.5]);
title('F2=s2*st2');
e_fsk=F1+F2;
subplot(413);
plot(e_fsk);
axis([0 500 -1.5 1.5]);
title('2FSK信号');%键控法产生的信号在相邻码元之间相位不一定连续

进行两个载波调制并加到一起。
在这里插入图片描述

加噪

%加噪
nosie=rand(1,len*100);%加入均值为0，方差为1的高斯白噪声
fsk=e_fsk+nosie;
subplot(414);
plot(fsk);
axis([0 500 -1.5 1.5]);
title('加噪声后信号')

fsk=conv(fsk,20);%过匹配滤波器

figure(4);
subplot(311);
plot(fsk);
axis([0 500 -30 30]);
title('过匹配滤波器后的加噪声后信号')

匹配滤波器fsk=conv(fsk,20)，其实就是卷积，加噪后噪声是很小的，通俗点说，翻倍增长，原序列比例为10，噪声比例为0.1，二者比例为100，而假如翻二倍，原序列比例为100，噪声比例为0.01，二者比例为10000，因而大大降低了噪声对原序列还原时会造成的影响。
在这里插入图片描述
在这里插入图片描述

相干解调

%相干解调
st1=fsk.*s1; %与载波1相乘
st2=fsk.*s2;%与载波2相乘

subplot(312);
plot(t,st1);
axis([0 0.001 -30 30]);
title('加噪后的信号与s1相乘后波形');
subplot(313);
plot(t,st2);
axis([0 0.001 -30 30]);
title('加噪后的信号与s2相乘后波形');

在这里插入图片描述

过低通滤波器

%通过低通滤波器
st1=filter(erlow,st1);
st2=filter(erlow2,st2);

figure(5);
subplot(211);
plot(t,st1);
axis([0 0.001 -30 30]);
title('过低通滤波器后st1波形');
subplot(212);
plot(t,st2);
axis([0 0.001 -30 30]);
title('过低通滤波器后st2波形');

在这里插入图片描述

相干解调后两个波形要通过低通滤波器，这儿用到两个，但其实用的是一个参数，本人软件方面有点问题，故而设了两个低通滤波器。

判决

%判决
for m=0:200
    if st1(1,m*30+15)>st2(1,m*30+15)
        for j=m*30+1:(m+1)*30
            at(1,j)=1;
        end
    else
        for j=m*30+1:(m+1)*30
            at(1,j)=0;
        end
    end
end

此处注意的是判决长度需与符号采样点数保持一致，代码中参数为一个符号采样点数为30。

%subplot(313);
figure(6)
subplot(311);
plot(g1a);
title('基带信号st1');
axis([0,5000,-1,2]);
subplot(312);
plot(g2a);
title('基带信号反码st2');
axis([0,5000,-1,2]);
subplot(313);
plot(at);
axis([0,5000,-1,2]);
title('抽样判决后波形')

在这里插入图片描述

将基带信号与最终还原出来的信号放在一起做对比，看还原的是否成功。

总代码

fskk.m（主代码）

fs=3840000;        % 采样率配为3840kHz，为128000*30所得，且30720000/8
sample_num=30;         %一个符号采样点数
Rb=fs/sample_num;       %符号速率，128k bps

len_symbol=2500;                %符号个数，可计算总样点数为N=len_sambol*sample_num,
bits_per_sybol=1;           %2FSK一个符号对应1个比特
len=len_symbol*bits_per_sybol;           %数据源比特个数
f1=Rb*4;          %载波频率，单位Hz,载波频率配置需为Rb的整数倍,同时能被采样率整除，为128kHz
f2=Rb*2;

N=len_symbol*sample_num*100             %符号采样点数
dt=1/fs;
t0=0:dt:(N*2-1)*dt;
t=0:dt:(N/30-1)*dt;%符号时间轴

%% 生成数据源
st1=randi([0,1],1,len);
st2=-st1+1;%生成基带信号反码，1变成0,0变成1

g11=(ones(1,100))'*st1;  %产生方波信号 
g1a=g11(:)'; 
g12=(ones(1,100))'*st2;  %产生方波信号 
g2a=g12(:)'; 

figure(1);
subplot(211);
plot(g1a);
title('基带信号st1');
axis([0,500,-1,2]);
subplot(212);
plot(g2a);
title('基带信号反码st2');
axis([0,500,-1,2]);

%载波信号
s1=cos(2*pi*f1*t);
s2=cos(2*pi*f2*t);
figure(2)
subplot(211),plot(s1);
axis([0 500 -1.5 1.5]);
title('载波信号s1');
subplot(212),plot(s2);
axis([0 500 -1.5 1.5]);
title('载波信号s2');

%调制
F1=g1a.*s1;%加入载波1
F2=g2a.*s2;%加入载波2
figure(3);
subplot(411);
plot(F1);
axis([0 500 -1.5 1.5]);
title('F1=s1*st1');
subplot(412);
plot(F2);
axis([0 500 -1.5 1.5]);
title('F2=s2*st2');
e_fsk=F1+F2;
subplot(413);
plot(e_fsk);
axis([0 500 -1.5 1.5]);
title('2FSK信号');%键控法产生的信号在相邻码元之间相位不一定连续

%加噪
nosie=rand(1,len*100);%加入均值为0，方差为1的高斯白噪声
fsk=e_fsk+nosie;
subplot(414);
plot(fsk);
axis([0 500 -1.5 1.5]);
title('加噪声后信号')

fsk=conv(fsk,20);%过匹配滤波器

figure(4);
subplot(311);
plot(fsk);
axis([0 500 -30 30]);
title('过匹配滤波器后的加噪声后信号')
%相干解调
st1=fsk.*s1; %与载波1相乘
st2=fsk.*s2;%与载波2相乘

subplot(312);
plot(t,st1);
axis([0 0.001 -30 30]);
title('加噪后的信号与s1相乘后波形');
subplot(313);
plot(t,st2);
axis([0 0.001 -30 30]);
title('加噪后的信号与s2相乘后波形');

%通过低通滤波器
st1=filter(erlow,st1);
st2=filter(erlow2,st2);

figure(5);
subplot(211);
plot(t,st1);
axis([0 0.001 -30 30]);
title('过低通滤波器后st1波形');
subplot(212);
plot(t,st2);
axis([0 0.001 -30 30]);
title('过低通滤波器后st2波形');

%判决
for m=0:200
    if st1(1,m*30+15)>st2(1,m*30+15)
        for j=m*30+1:(m+1)*30
            at(1,j)=1;
        end
    else
        for j=m*30+1:(m+1)*30
            at(1,j)=0;
        end
    end
end
%subplot(313);
figure(6)
subplot(311);
plot(g1a);
title('基带信号st1');
axis([0,5000,-1,2]);
subplot(312);
plot(g2a);
title('基带信号反码st2');
axis([0,5000,-1,2]);
subplot(313);
plot(at);
axis([0,5000,-1,2]);
title('抽样判决后波形')

erlow.m（低通滤波器1）

function Hd = erlow
%ERLOW Returns a discrete-time filter object.

% MATLAB Code
% Generated by MATLAB(R) 9.2 and the Signal Processing Toolbox 7.4.
% Generated on: 29-Nov-2020 22:17:42

% FIR Window Lowpass filter designed using the FIR1 function.

% All frequency values are in Hz.
Fs = 48000;  % Sampling Frequency

N    = 11;       % Order
Fc   = 12800;    % Cutoff Frequency
flag = 'scale';  % Sampling Flag
Beta = 0.5;      % Window Parameter

% Create the window vector for the design algorithm.
win = kaiser(N+1, Beta);

% Calculate the coefficients using the FIR1 function.
b  = fir1(N, Fc/(Fs/2), 'low', win, flag);
Hd = dfilt.dffir(b);

% [EOF]

erlow.m（低通滤波器2）

function Hd = erlow2
%ERLOW Returns a discrete-time filter object.

% MATLAB Code
% Generated by MATLAB(R) 9.2 and the Signal Processing Toolbox 7.4.
% Generated on: 29-Nov-2020 22:17:42

% FIR Window Lowpass filter designed using the FIR1 function.

% All frequency values are in Hz.
Fs = 48000;  % Sampling Frequency

N    = 11;       % Order
Fc   = 12800;    % Cutoff Frequency
flag = 'scale';  % Sampling Flag
Beta = 0.5;      % Window Parameter
% Create the window vector for the design algorithm.
win = kaiser(N+1, Beta);

% Calculate the coefficients using the FIR1 function.
b  = fir1(N, Fc/(Fs/2), 'low', win, flag);
Hd = dfilt.dffir(b);

% [EOF]

祝您生活顺心，事事顺利。

enduringzx
关注
  • 38
    点赞
  • 踩
  • 215
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
用Matlab编程仿真2FSK调制和相干解调系统.doc
06-03
1.1、设计内容 用Matlab编程仿真2FSK调制和相干解调系统 1.2、设计要求 （1）信源为八位二进制随机信号，比特率为50bit/s，载波频率f1=200Hz,f2=100Hz； （2）画出信源，2FSK已调信号，解调器乘以相干载波后的信号，经过低通滤波器后的信号和抽样判决后的信号波形； （3）撰写专业课程设计报告。 2、设计原理分析 2.1、2FSK介绍 数字频率调制又称频移键控（FSK），二进制频移键控记作2FSK。数字频移键控是用载波的频率来传送数字消息，即用所传送的数字消息控制载波的频率。2FSK信号便是符号“1”对应于载频f1，而符号“0”对应于载频f2（与f1不同的另一载频）的已调波形，而且f1与f2之间的改变是瞬间完成的。 其表达式为： 典型波形如下图所示。由图可见,2FSK信号可以看作两个不同载频的ASK信号的叠加。因此2FSK信号的时域表达式又可以写成： 2.2、 2FSK调制原理 2FSK调制就是使用两个不同的频率的载波信号来传输一个二进制信息序列。可以用二进制“1”来对应于载频f1，而“0”用来对应于另一相载频w2的已调波形，而这个可以用受矩形脉冲序列控制的开关电路对两个不同的独立的频率源f1、f2进行选择通。本次课程设计采用的是前面一种方法。如下原理图： 图1 调制原理框图 2.3、 2FSK解调原理 2FSK的解调方式有两种:相干解调方式和非相干解调方式，本次课程设计采用的是相干解调法，利用载波与已调信号进行相乘后滤波输出得到，在2FSK中要两个载波，所以解调也要两个载波，分别与已调信号相乘后利用低通，最后相加即可得到我们的滤波输出，再用抽样信号进行抽样判决器即可,其原理如下： 图2 相干解调原理框图 3、仿真参数设置说明 3.1、仿真参数设置 信源：八位二进制随机信号 载波频率：fc=900 信号频率：f1=100;f2=200 比特率：Rb= Rb=fm*log2_M=50 bit/s 3.2、仿真程序设计说明 详细设计步骤: 1．信号产生：二进制随机序列和两列频率不等的载波 1)利用matlab 库函数产生8个二进制随机数，也就是我们的基波调制信号a。 2)产生两列基带信号st1和st2，并画出其波形。 2．信号调制：产生2FSK信号和加入高斯噪声后的2FSK信号 1)用基带信号st1、st2和载波信号s1和s2，产生2fsk信号，具体做法是用以st1与 s1相乘，st2 与s2相乘，再将两列信号相加。 2) 2fsk信号相加得到加入噪声后的sn信号。 3．解调 （1）对于两列让sn通过两个低通滤波器，并画出经过带通滤波器后的波形。 （2）输入抽样脉冲，将两列波st1和st2通过抽样判决器，画出其波形st，并与之前调制后的波形sn做对比，并由Figure3得出比特率为50 bit/s。 4、仿真结果和分析 4.1、仿真结果 4.1.1.产生的基带信号 4.1.2.产生的载波信号 4.1.3.2FSK调制过程
2fsk调制解调MATLAB仿真代码
09-19
2fsk信号的相干解调MATLAB代码。此为一般课程所能用到，通过对2fsk信号的产生进行模拟，以及信号在信道中传播，进过滤波器等产生等。画出图像
【Matlab】2FSK二进制频移键控波形生成
Otm-Z Blog
01-05 1718
一个通信原理课程中使用Matlab生成2FSK波形的实验笔记。
深入探索数字通信：使用MATLAB进行2FSK调制与信道模拟
最新发布
weixin_52846855的博客
05-04 267
我们最近使用MATLAB模拟了2FSK调制过程，并通过详细的代码实现和仿真分析，深入探讨了这一技术的核心原理和应用。通过这一仿真，不仅展现了FSK在数字通信中的实用性，还强调了理论知识与实际应用之间的联系。此外，我们的仿真模型为电子和通信工程的学生提供了一个实用的学习平台，帮助他们更好地理解和掌握相关概念。2FSK调制与信号合成：每对载波对应于二进制数据中的'0'和'1'，生成调制后的信号，并将所有载波的信号合成。功率谱密度分析：在信道处理前后，对信号的功率谱密度进行了详细分析，可视化展示了滤波效果。
通信原理与MATLAB(七):2FSK的调制解调
菜yuan~的博客
12-27 3万+
通信原理和MATLAB
【Matlab】通信系统实验：二进制数字调制与仿真
codeforce111的博客
06-19 1065
利用Matlab实现对二进制数字调制与解调的仿真，包括以下三部分：二进制数字相移键控（2PSK），二进制数字振幅键控（2ASK），二进制数字频移键控（2FSK）。title('2ASK相干解调波形');title('2FSK相干解调波形');title('2ASK已调信号波形')
2ASK和2FSK调制解调仿真
Quiet_Superman的博客
01-25 2万+
实验一2ASK和2FSK调制解调仿真 一、实验目的 1．熟悉2ASK和2FSK调制解调原理。 2．掌握编写2ASK和2FSK调制解调程序的要点。 3．掌握使用Matlab调制解调仿真的要点。 二、实验内容 1．根据2ASK和2FSK调制解调原理，设计源程序代码。 2．通过Matlab软件仿真给定信号的调制波形。 3. 对比给定信号的理论调制波形和仿真解制波形。 三、实验原理 1．2ASK 二进制振幅键控（2ASK）信号码元为： ...
基于MATLAB的2FSK调制解调系统设计。MATLAB代码设计程序
CtnKuap的博客
11-21 1152
近年来，随着通信技术的迅速发展，2FSK调制解调系统逐渐成为研究的热点。本文将介绍基于MATLAB的2FSK调制解调系统的设计，包括程序代码设计、相干解调与非相干解调两种方式。综上所述，基于MATLAB的2FSK调制解调系统设计需要包括代码设计、纯.m代码和simulink仿真的实现，以及相干解调和非相干解调两种方式的考量。二进制数字转换为频率的方式就是所谓的FSK调制，即在不同的数字0、1对应不同的频率。在设计中，可以选择使用纯.m代码或simulink仿真来完成2FSK调制解调系统的搭建。
2FSK_2FSK_调制解调_fsk信号_matlab2fsk_
10-03
MATLAB中实现2FSK调制主要包括以下步骤： 1. **生成二进制序列**：使用`randi`函数生成随机二进制序列，或者根据需要指定特定的二进制数据。 2. **设定载波频率**：定义两个不同的载波频率，如fc0（代表二进制0的...
基于MATLAB的2FSK调制及仿真.doc
01-07
通过编写MATLAB程序，可以实现2FSK调制的过程，包括产生二进制序列、频率切换逻辑以及频率调制函数。程序源代码将生成模拟信号，并显示相应的波形图，以便于分析和验证。 5. 总结： 2FSK调制是数字通信系统中的重要...
2FSK调制解调实验，MATLAB代码
10-15
在MATLAB中，可以使用`pskmod`函数配合适当的参数来实现2FSK调制。例如： ```matlab modulated_signal = pskmod(binary_data, 2, 'fsk'); ``` 这里，`binary_data`是待调制的二进制序列，2表示调制方式为2FSK，'...
FSK调制解调的matlab、Verilog、VHDL代码
03-19
这是FSK二进制频移键控的仿真和FPGA实现的代码，用matlab做了连续、非连续fsk的仿真，用quartus进行verilog和VHDL实现，最后用modelsim进行仿真。
1.rar_FSK调制 MATLAB
07-14
标题中的“1.rar_FSK调制 MATLAB”表明这是一个与频移键控（FSK）调制技术相关的MATLAB实现。MATLAB是一种强大的编程环境，常用于信号处理、数学计算和模拟，非常适合用来理解和实现通信系统中的各种调制方法，如FSK...
2FSK调制与解调.zip
01-11
这个转换过程由2FSK调制器完成，如压缩包中的"2FSK调制器.ms14"所示。调制器通常包含一个数字信号处理器（DSP）和一个混频器，通过改变本地振荡器的频率来实现调制。 在Multisim仿真中，"主载波振荡器.ms14"是生成...
MATLAB编程实现2FSK信号的调制与解调（非相干解调）
wwwEERYSTG的博客
09-08 1万+
MATLAB编程实现2FSK信号的调制与解调（非相干解调）
FSK调制解调
热门推荐
gemengxia的博客
03-25 6万+
文章目录1.FSK调制解调原理1.1FSK时域信号1.2相关系数与频谱特性1.3非相干解调1.3.1相乘微分型AFC环解调法1.3.2包络检波解调法1.4相干解调原理1.4.1最佳FSK相干解调器1.4.2易于实现的FSK相干解调器2.FSK的MATLAB仿真2.1FSK信号产生2.2非相干解调FSK2.2相干解调FSK仿真 1.FSK调制解调原理 1.1FSK时域信号 数字频率调制（FM）是利用载波的频率传输信息的一种调制方式，其中最简单的是二进制频移键控（2FSK）。FM分为非连续相位FSK和连续相位F
FSK调制matlab源码程序
FPGA/MATLAB学习教程/源码/项目合作开发
12-31 1万+
clc; clear all; close all; x=[ 1 0 1 0 1 0 1]; % Binary Information bp=.0001; % bit period disp(' Binary information at Trans mitter :'); disp(x); %XX representation o...
用MATLAB编程实现2FSK信号的调制解调（包络检波与相干解调）代码
IoRAGIzAP的博客
11-30 2222
其中，2FSK（二进制频移键控）调制和解调技术是一种常用的调制和解调方式。本文将基于MATLAB编程实现2FSK信号的调制解调，并对不同信噪比下的包络检波与相干解调的误码率性能进行对比分析。综上所述，本文基于MATLAB编程实现了2FSK信号的调制解调，并对不同信噪比下的包络检波与相干解调的误码率性能进行了对比分析。因此，我们需要对不同信噪比下的包络检波与相干解调的误码率性能进行比较。通过对不同信噪比下的误码率性能进行比较，我们可以看到相干解调的性能明显优于包络检波，且在信噪比较高时，误码率近乎为0。
Matlab通信仿真系列——载波频率调制(FSK)仿真
小灰灰的FPGA的博客
11-25 2245
当Δf时1/(2T)的倍数时候，信号是正交的，也就是要保证信号的正交性，必须满足相继频率之间的最小频率间隔是1/(2T)。从图中可以看出载波频率4Hz的两边，分布了对应符号0,1,3,2的频率2.5Hz,3.5Hz,4.5Hz,5.5Hz。最佳检测器，又称为平方律检测器，当传输信号波形等概率时，通过计算平方包络(rk)^2= (rkc)^2+ (rks)θm表示第m个信号由于传输延时而产生的相移，n(t)表示加性带通噪声。Δf是相继两个频率之间的频率间隔，Δf=f(m)-f(m-1)2}的最大值的信号。
MATLAB-Simulink实现的2FSK调制解调仿真研究
- **2FSK调制**：在调制过程中，根据输入的二进制序列，载波在两个预设频率之间切换，生成2FSK信号。 - **2FSK解调**：相干解调是最常见的2FSK解调方法，通过比较接收到的2FSK信号与本地载波的相位，来恢复原始二...
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值