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

最新推荐文章于 2024-02-27 10:33:30 发布
最新推荐文章于 2024-02-27 10:33:30 发布
阅读量8.7k 收藏 204
点赞数 36
文章标签： 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
关注
  • 36
    点赞
  • 踩
  • 204
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
matlab_2ask、2fsk、2psk仿真_
10-02
波形仿真，已通过测试，其中test文件是的，其他文件是其它类型可以参考
【数字信号调制】基于matlab二进制数字频率调制（2FSK）【含Matlab源码 998期】.zip
11-29
CSDN海神之光上传的代码均可运行，亲测可用，直接替换数据即可，适合小白； 1、代码压缩包内容 主函数：fsk2.m； 调用函数：其他m文件；无需运行 运行结果效果图； 2、代码运行版本 Matlab 2019b；若运行有误，根据提示修改；若不会，私信博主； 3、运行操作步骤 步骤一：将所有文件放到Matlab的当前文件夹中； 步骤二：双击打开fsk2.m文件； 步骤三：点击运行，等程序运行完得到结果； 4、仿真咨询 如需其他服务，可私信博主或扫描博客文章底部QQ名片； 4.1 博客或资源的完整代码提供 4.2 期刊或参考文献复现 4.3 Matlab程序定制 4.4 科研合作 功率谱估计： 故障诊断分析： 雷达通信：雷达LFM、MIMO、成像、定位、干扰、检测、信号分析、脉冲压缩 滤波估计：SOC估计 目标定位：WSN定位、滤波跟踪、目标定位 生物电信号：肌电信号EMG、脑电信号EEG、心电信号ECG 通信系统：DOA估计、编码译码、变分模态分解、管道泄漏、滤波器、数字信号处理+传输+分析+去噪、数字信号调制、误码率、信号估计、DTMF、信号检测识别融合
【MATLAB源码-第149期】基于MATLAB的2ASK,2FSK,2PSK,2DPSK等相干解调仿真，输出各节点波形。
Matlab程序猿i的博客
02-27 1846
操作环境：1、算法描述2ASK（二进制幅移键控）、2FSK（二进制频移键控）、2PSK（二进制相移键控）和2DPSK（二进制差分相移键控）是数字调制技术中的基本调制方式，它们在无线通信、数据传输等领域有着广泛的应用。相干解调是这些调制方式中一个重要的解调技术，它要求接收端的本地振荡器与发送端的载波同步，以实现信号的高效、准确解调。接下来，我将详细介绍这四种调制技术及其相干解调原理。
matlab程序 地震 相干噪声_无网格局部Petrov-Galerkin法及程序设计（15）——代码：计算RPIM形函数...
weixin_39609457的博客
01-04 350
参考资料G.R.Liu. 2010. Mesh Free Methods: Moving Beyond the Finite Element Method. Second Edition.G.R.Liu. and Y.T.Gu. 2005. An Introduction to Meshfree Methods and Their Programming张雄.2003.无网格方法数值实现Matla...
FPGA 数字信号处理之 FSK 调制、解调的实现与仿真基于 verilog + ise + modelsim + matlab （保姆级）
大杂烩
05-30 6224
试验说明： 笔者复现FSK调制、解调的过程中，在网上查了很多很多相关文章，至少在本文之前，没发现能在参考一篇博文的情况下能完整复现实验结果的。总之，全乎的基本没有。 本笔记（文）的目标是让入门者或者说小白在参考本文的基础上，能一步一步完整复现实验结果，少走弯路，尽力节省时间，因为时间很宝贵。 另：本实验只是复现参考博文的结果，有关概念、原理相关知识请参考原文。 主要参考博文链接如下： https://bestfpga.blog.csdn.net/article/details/81166187 实验目标：
matlab对信号加噪代码,Matlab给信号加噪声
weixin_42561249的博客
03-17 8362
在信号处理中经常需要把噪声叠加到信号上去，在叠加噪声时往往需要满足一定的信噪比，这样产生二个问题，其一噪声是否按指定的信噪比叠加，其二怎么样检验带噪信号中信噪比满足指定的信噪比。在MATLAB中可以用randn产生均值为0方差为1的正态分布白噪声，但在任意长度下x=randn(1,N)，x不一定是均值为0方差为1(有些小小的偏差)，这样对后续的计算会产生影响。在这里提供3个函数用于按一定的信噪比把...
2fsk调制解调MATLAB仿真代码
09-19
2fsk信号的相干解调MATLAB代码。此为一般课程所能用到，通过对2fsk信号的产生进行模拟，以及信号在信道中传播，进过滤波器等产生等。画出图像
通信原理与MATLAB(七):2FSK的调制解调
热门推荐
菜yuan~的博客
12-27 2万+
通信原理和MATLAB
【Matlab】通信系统实验：二进制数字调制与仿真
codeforce111的博客
06-19 771
利用Matlab实现对二进制数字调制与解调的仿真，包括以下三部分：二进制数字相移键控（2PSK），二进制数字振幅键控（2ASK），二进制数字频移键控（2FSK）。title('2ASK相干解调波形');title('2FSK相干解调波形');title('2ASK已调信号波形')
2ASK和2FSK调制解调仿真
Quiet_Superman的博客
01-25 1万+
实验一2ASK和2FSK调制解调仿真 一、实验目的 1．熟悉2ASK和2FSK调制解调原理。 2．掌握编写2ASK和2FSK调制解调程序的要点。 3．掌握使用Matlab调制解调仿真的要点。 二、实验内容 1．根据2ASK和2FSK调制解调原理，设计源程序代码。 2．通过Matlab软件仿真给定信号的调制波形。 3. 对比给定信号的理论调制波形和仿真解制波形。 三、实验原理 1．2ASK 二进制振幅键控（2ASK）信号码元为： ...
MATLAB编程实现2FSK信号的调制与解调（非相干解调）
wwwEERYSTG的博客
09-08 1万+
MATLAB编程实现2FSK信号的调制与解调（非相干解调）
FSK调制解调的matlab、Verilog、VHDL代码
03-19
这是FSK二进制频移键控的仿真和FPGA实现的代码，用matlab做了连续、非连续fsk的仿真，用quartus进行verilog和VHDL实现，最后用modelsim进行仿真。
用Matlab编程仿真2FSK调制和相干解调系统.doc
06-03
用Matlab编程仿真2FSK调制和相干解调系统 1.2、设计要求 （1）信源为八位二进制随机信号，比特率为50bit/s，载波频率f1=200Hz,f2=100Hz； （2）画出信源，2FSK已调信号，解调器乘以相干载波后的信号，经过低通...
2FSK.zip_2fsk matlab
07-15
2FSK仿真演示系统，可以任意输入码字序列，和载波频率，动态产生调制波形。可以清楚的看到整个调制解调过程的时域及频域的波形。
FSK 扩频仿真：FSK 调制信号的扩频-matlab开发
05-31
这是带有 PN 序列发生器的 FSK 调制扩展。 它是一个模拟文件，说明了扩频通信的影响。 我对合作很感兴趣。 [email protected]
A毕业设计：MATLAB的FSK和PSK调制与解调系统
最新发布
03-17
1. **FSK调制**： - 将数字信号转换为序列信号，每个符号映射到一个频率。 - 在特定频率上生成正弦波或载波信号。 - 对载波信号进行调制，将其频率调整为代表0和1的两种不同频率。 2. **FSK解调**： - 接收调制...
一种基于软件无线电的FSK调制设计与实现
01-13
用MATLAB仿真了数据产生、映射、高斯滤波、相位积分、相位-幅度转换和正交调制等FSK调制的链路级设计，用XC5VSX35T、AD9779等芯片完成了PN序列生成、滤波因子卷积、正余弦寻址查找、混频求和等FSK调制的电路实现，...
瑞利信道下卷积码
lin_yulin的博客
03-21 1162
%2FSK调制在瑞利信道下，卷积码 clear all SNR=-10:10; N=30000; %消息比特个数 berjj=zeros(1,length(SNR)); L=7; %约束长度 trel=poly2trellis(L,[171 133]); tblen=6*L; SymbolRate=2000; %符号速率...
matlab 2ask相干解调误码率太低,2ASK 2PSK 2FSK在不同信噪比下三种调制方式误码性能MATLAB仿真...
weixin_39604280的博客
03-16 2516
clear,clc;L=10000;t=0:0.01*pi:2*pi;SNR=-4:0.5:10;snr=SNR/10;for ii=1:length(snr)Sigma(ii)=1/sqrt(2*(10^snr(ii)));%标准差psk_error1=0; fsk_error1=0; ask_error1=0; %误码数初始值为0psk_error2=0; fsk_error2=0; as...
2FSK调制的MATLAB实现
05-12
2FSK调制是一种基于频移键控技术的数字调制方式。在MATLAB中，可以通过以下步骤实现2FSK调制： 1. 确定调制信号的基本参数，包括信号带宽、调制频率、采样频率等。 2. 生成调制信号的基带信号，可以使用sawtooth、sine、square等函数实现，其中sine函数可以用于产生正弦波，square函数可以用于产生方波。 3. 生成调制信号，使用fskmod函数进行2FSK调制，该函数需要输入基带信号、调制频率、采样频率等参数。 4. 对调制信号进行波形显示和频谱分析，可以使用plot和fft等函数实现。 以下是一个简单的2FSK调制代码示例： ``` % 设置调制信号参数 fc1 = 100; % 第一信号频率 fc2 = 200; % 第二信号频率 fs = 1000; % 采样频率 t = 0:1/fs:1; % 时间序列 m1 = sin(2*pi*fc1*t); % 第一基带信号 m2 = sin(2*pi*fc2*t); % 第二基带信号 % 2FSK调制 f_dev = 50; % 频偏 mod_signal = fskmod([m1', m2'], [fc1, fc2], fs, f_dev); % 波形显示 figure; plot(t, mod_signal); % 频谱分析 figure; freq_axis = linspace(-fs/2, fs/2, length(mod_signal)); mod_spectrum = fftshift(abs(fft(mod_signal))); plot(freq_axis, mod_spectrum); ``` 在上述代码中，首先通过设置调制信号的参数，生成两个基带信号m1和m2，然后使用fskmod函数对其进行2FSK调制，产生调制信号mod_signal。最后使用plot和fft函数分别对调制信号进行波形显示和频谱分析。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值