当前位置: 首页 > news >正文

Matlab信号处理3:fft(快速傅里叶变换)标准使用方式

Fs = 1000;      % 采样频率
T = 1/Fs;       % 采样周期:0.001s
L = 1500;       % 信号长度
t = (0:L-1)*T;  % 时间向量.  时间向量从0开始递增,0s~1.499sS = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);  % 模拟原信号
X = S + 2*randn(size(t));                  % 加入干扰信号subplot(221);
plot(t,S);
xlabel("时间(s)");
ylabel("幅值");
title("原信号波形");subplot(222);
plot(t,X);
xlabel("时间(s)");
ylabel("幅值");
title("加入干扰信号波形");Y=fft(X);
len = length(Y);     % Len = L
P2=abs(Y/L);         % 计算幅值并进行能量还原% 751个,P1提取P2中左半部分的第1~751个
P1 = P2(1:L/2+1);  
% 绘制单边谱,所以除了首尾元素,其余信号幅度翻倍
P1(2:end-1) = 2*P1(2:end-1);   
%根据奈奎斯特定理,此处采样频率Fs为1000,因此频域分析的信号频率最大为500
f = Fs*(0:(L/2))/L;  % 0 ~ 500Hzsubplot(223);
plot(f,P1);
xlabel("频率(Hz)");
ylabel("幅值");
title("含有干扰信号的傅里叶变换");

运行效果:

http://www.lryc.cn/news/164218.html

相关文章:

  • Python|合并两个字典的几种方法
  • ElementUI浅尝辄止24:Message 消息提示
  • 让照片动起来的软件,轻松制作照片动效
  • 【图解RabbitMQ-7】图解RabbitMQ五种队列模型(简单模型、工作模型、发布订阅模型、路由模型、主题模型)及代码实现
  • Linux命令200例:write用于向特定用户或特定终端发送信息
  • javaee spring整合mybatis spring帮我们创建dao层
  • 修改Tomcat的默认端口号
  • Open3D Ransac拟合空间直线(python详细过程版)
  • 题目:2729.判断一个数是否迷人
  • 微服务模式:服务发现模式
  • 9.4 数据库 TCP
  • 普通用户使用spark的client无法更新Ranger策略
  • Git超详细教程
  • C++ 回调函数
  • xilinx FPGA IOB约束使用以及注意事项
  • 如何统计iOS产品不同渠道的下载量?
  • 大模型学习
  • Redis原理:IntSet
  • 【已解决】Splunk 8.2.X 升级ES 后红色报警
  • 香橙派使用外设驱动库wiringOP 配合定时器来驱动舵机
  • C++学习笔记--函数重载(2)
  • 代码随想录算法训练营Day56 || ● 583. 两个字符串的删除操作 ● 72. 编辑距离
  • chrome_elf.dll丢失怎么办?修复chrome_elf.dll文件的方法
  • 代码随想录32|738.单调递增的数字,968.监控二叉树,56. 合并区间
  • BIO NIO AIO演变
  • JVM GC垃圾回收
  • 【数据结构】队列知识点总结--定义;基本操作;队列的顺序实现;链式存储;双端队列;循环队列
  • 嵌入式学习之链表
  • 静态代理和动态代理笔记
  • [SM6225][Android13]user版本默认允许root和remount