随机解调的多频点信号x与伪随机序列经过混频后,被均匀的涂抹到了整个频率轴上,然后经低通滤波,低速均匀采样,最后通过OMP等算法恢复原始信号,整体上是压缩感知求解欠定方程y=Ax的思路。以下给出matlab程序和各个信号的时频域图,更容易理解。

这里的x(n)包含了两个频率:3000Hz和4000Hz,伪随机序列由matlab内部函数randsrc()生成。

首先是绘制频域稀疏信号x(n)与m序列的时域图和频域图

clc;clear all;close all;
%%-------------------------------------------------
Fs = 10000;%采样频率
dt = 1/Fs;%采样时间间隔
L = 1000;
t = (0:L-1)*dt;%生成时间向量
f1=3000;%频率1
f2=4000;%频率2
x1=cos(2*pi*f1*t);
x2=cos(2*pi*f2*t);
x  = x1+x2;

F1 = fft(x)/L;
W = (0:L/2-1)*(Fs/L);%实际频率
subplot(2,1,1);%2行3列,第一个区域
plot(t(1:40),x(1:40),'o');%绘制采样点,只取了一部分区间

figt =(0:400)/100000;
figx = cos(2*pi*f1*figt)+cos(2*pi*f2*figt);
hold on;
plot(figt,figx);
%axis([0 0.1 -2.2 2.2]);
xlabel('时间/s');ylabel('幅值/v');title('频域稀疏信号的时域图');
subplot(2,1,2);%2行1列,第一个区域
stem(W,abs(F1(1:L/2)));%stem用于绘制茎状图
xlabel('频率/Hz');ylabel('幅值');title('频域稀疏信号的频域图');

figure;
pn = randsrc(1,L);
subplot(2,1,1);%2行1列,第一个区域
stairs(t(1:100),pn(1:100));%绘制阶梯图
axis([0 0.01 -1.2 1.2]);%调节坐标轴刻度
xlabel('时间/s');ylabel('幅值/v');title('m序列时域图');

第一幅图片是多频点信号的时域和频域图,由于plot绘图时默认的是直接用直线连接,因此波形不够平滑,为此我绘制了一条非常秘集的线,在线上用圆圈标记出了采样点。

由第二幅图我们可以看到,伪随机序列的频谱非常杂乱。事实上伪随机序列是周期信号,一个周期是TP的话,基波频率为1/TP,周期信号只能做傅里叶级数展开,它的频谱是基波及高次谐波组成的,谐波频率是基波频率的整数倍。

figure;
y = x.*pn;%y为混频信号
subplot(2,1,1);
plot(t,y);
axis([0 0.1 -2.5 2.5]);
xlabel('时间/s');ylabel('幅值/v');title('原信号和m序列相乘后的信号的波形图');
Fy = fft(y)/L;
subplot(2,1,2);
plot(W,abs(Fy(1:L/2)));
xlabel('频率/Hz');ylabel('幅值');title('混频后的信号频谱');

混频后的频谱依然是杂乱无章,但现在频谱的低频区不只是伪随机序列的频谱,它已经包含了原始输入信号的频谱。

随机解调-多频点信号与伪随机序列混频相关推荐

  1. 随机解调系统的伪随机序列m序列的产生及频谱图

    这是我压缩感知分栏的第一篇文章,我将压缩感知.模拟信息转化领域中的一些概念.公式推导.程序写在博客里,一方面能够记录我的学习历程,另一方面可以分享给别人. 随机解调是压缩感知理论的一种实际应用,它针对 ...

  2. 基于EWT的单通道时频混合信号的分离研究

    基于EWT的单通道时频混合信号的分离研究 主要探讨基于经验小波分解(EWT)的数据自驱动分解方法是否适用于时频混叠信号,这里通过实验数据与分类效果对该方法进行研究和探讨. 首先定义混合信号,为了方便对 ...

  3. 频谱细化matlab程序,分享FFT频谱细化程序(处理单频点信号)

    function resultCorrect=spectrumcorrectfftftmethod(inputDate,fStart,fEnd,fDelt,windowType,fs) %功能:采用F ...

  4. 【随机共振】基于随机共振的高灵敏度GPS信号捕获算法

    算法结构图 首先input输入到SR 这个语句意思就是每次输入一段数据进行处理,因为这个系统是实时反馈的系统,所以在处理的时候,必须是一段段处理,而不能和原来一样一次性进行处理,所以我们每次处理的数据 ...

  5. matlab进行fm调制与解调,基于matlab的fm信号调制与解调.doc

    基于matlab的fm信号调制与解调.doc 调制就是将基带信号的频谱搬移到信道通带中或者其中的某个频段上的过程,而解调是将信道中来的频带信号恢复为基带信号的反过程.调制的目的是把要传输的模拟信号或数 ...

  6. 信号的调制与解调matlab仿真,基于MATLAB对信号调制与解调的仿真

    基于MATLAB对信号调制与解调的仿真 2.3 2PSK的基本原理和调制解调实现 数字相移键控,记作PSK(Phase shift-keying ),二进制相位键控记作2PSK. 2PSK就是根据基带 ...

  7. matlab实现幅度调制,利用matlab实现信号幅度的调制与解调钟媛

    利用matlab实现信号幅度的调制与解调钟媛 1<MATLAB 语言>课程论文利用 MATLAB 实现信号幅度的调制与解调 姓名: 钟媛学号:12010245219专业:电子信息工程班级: ...

  8. 基于FPGA的AM信号调制与解调详细步骤

    详细程序及原理参考原文: FPGA学习(二)--实现AM信号调制与解调 - 子木的文章 - 知乎 https://zhuanlan.zhihu.com/p/37203478 使用的软件是Vivado ...

  9. MATLAB用相干解调DSB信号,AM DSB信号解调的MATLAB实现

    AM &DSB信号解调的MATLAB实现 一.实验目的 掌握AM&DSB解调器电路原理: 掌握AM&DSB解调仿真原理. 了解MATLAB软件,学习应用MATLAB软件的仿真技 ...

最新文章

  1. 028_CSS外边距
  2. php进程守护进程,php 多进程实现守护进程的实例代码
  3. java web后台_java web 后台那些事
  4. C语言指针怎么存二维数组,C语言怎么用指针代替二维数组
  5. js轮询导致服务器瘫痪_面试官:讲一下什么是负载均衡,什么是轮询策略随机策略哈希策略...
  6. drools 7.x 决策表使用
  7. django migrate 抛出异常:ValueError: Found wrong number (0) of constraints for …
  8. Linux下patch打补丁命令
  9. Matlab模拟四旋翼飞行器PID控制仿真
  10. IDEA设置鼠标滚轮控制缩放大小
  11. php获得视频文件扩展名,php 获取文件扩展名的 n 种方法
  12. 计算机笔记本硬盘,笔记本取证之--笔记本硬盘拆卸
  13. 精选教程:来亲手开发个问答社区,干掉知乎!
  14. 【opencv学习笔记】2用摄像头识别指定颜色
  15. 学习 stm32(TTL)串口通信控制16路舵机控制板(维特智能)
  16. C++——百分制成绩转五分制成绩
  17. 速写在网站建设中的重要地位
  18. GitHub下载 无法分配请求的地址_Hexo+Github--搭建个人博客(一)准备工作amp;amp;环境搭建
  19. iOS拍照和录制视频 摄像头使用
  20. 个人电子邮箱如何在网页进行登录?

热门文章

  1. Java基础案例3-4:学生投票系统
  2. UnityShdaer笔记第二课笔记-贴花
  3. ROS:启动管理工具launch
  4. AI赋能医疗大健康产业 千亿市场触手可及
  5. C/C++/Linux「学习必备」珍藏多年的学习网站分享 | 工欲善其事必先利其器
  6. “端到端”是什么意思
  7. 计算机初始化,电脑初始化是什么意思?
  8. 动态规划01背包问题之跳跃点解法
  9. 可可金融COO李劳:启蒙时期3类优质数字资产的逻辑与研判
  10. matlab课程数字图像处理设计题目,数字图像处理课程设计题目和要求_2013