matlab 时域、频域和谱域中信号图

close all;clear all;clc;
data = [];
N=200;%采样点个数
fs = 50; %采样频率
t = 0:1/fs:(N-1)/fs;  %生成信号时间
data = [data csvread('D:\matlab_notebook\PD表\data\data_4s\196.csv')];%导入数据
% plot(data)
x=data(1:200,1);%提取第一列数据 1：N  (1:200,1)
% x=data(1:200,4);%提取第一列数据 1：N  (1:200,1)

%FFT
figure(1);subplot(211);plot(t,x);% 画出原始时域信号
% title('The combined acceleration');%score:2
ylabel('Acc'); %Acceleration
xlabel('Time/s');

signalFFT=abs(fft(x,N));%真实的幅值
Y=2*signalFFT/N;
f=(0:N/2)*(fs/N);
subplot(212);
plot(f,Y(1:N/2+1));
xlim([0 15]);
ylabel('Amp'); %Amplitude
xlabel('Frequency');
title('Spectrum');%Spectrum输入信号的频谱

在这里插入图片描述

%功率谱

%调用periodogram函数
[Pxx,f] = periodogram(x,rectwin(length(x)), N, fs);
subplot(212);
freq = 0:fs/length(x):fs/2;%频率刻度
train_y = 10*log10(Pxx)
for i= 1:size(train_y,2)
    if train_y(1,i)< -40
       train_y(1,i)==-40
    end
end
figure(2);plot(freq,train_y);
grid on;xlim([0 fs/2]);
title('Power spectral density');%用periodogram函数的周期图
xlabel('Frequency');
ylabel('PSD');%功率谱密度/db

在这里插入图片描述

本文链接：https://blog.csdn.net/weixin_50951788/article/details/126366887

智能推荐

Java递归与迭代求斐波那契数列

Fibonacci 数列：指的是这样一个数列：1、1、2、3、5、8、13、21、34……即从第三项开始，每一项等于它的前两项之和。递归程序调用自身的编程技巧称为递归。递归做为一种算法在程序设计语言中广泛应用。一个过程或函数其定义或说明中有直接或间接调用自身的一种方法，它通常把一个大型复杂的问题层层转化为一个与原问题相似的规模较小的问题来求解，递归策略只需少量的...

Android 如何实现时间年月日三级联动呢直接上代码看效果图

2 , 是不是很简单呢？觉的不错的小伙伴可以点个赞。 ...

Redis（三）：Redis的内存淘汰机制与持久化机制

一：Redis的内存淘汰机制 redis 设置过期时间 Redis中有个设置时间过期的功能，即对存储在 redis 数据库中的值可以设置一个过期时间。作为一个缓存数据库，这是非常实用的。如我们一般项目中的 token 或者一些登录信息，尤其是短信验证码都是有时间限制的，按照传统的数据库处理方式，一般都是自己判断过期，这样无疑会严重影响项目性能。我们 set key 的时候，都可以给...

druid连接池配置监控页面

jdbc.properties相关配置 spring-content.xml相关配置 web.xml 相关配置以上三种配置完成后输入localhost：8080/druid/index.jsp,监控页面显示如下图： ...

用基数排序实现单词按字典序排序（包含大写）

输入字符串的处理: 所以字符串中最长单词的长度为mxp 基数排序时将每一个长度小于mxp的字符串后位填补‘A’ 即Cards Cap中mxp=4，Cap–>CapAA 这样所有的字符串就统一长度。代码：...

matlab 时域、频域和谱域中信号图

智能推荐

猜你喜欢