傅里叶变换—信号观测时长和频率分辨率
傅里叶变换—信号观测时长和频率分辨率
前言
一、开始验证
附录
MATLAB代码:
%% observe time T and frequency resolution verified
f1 = 1; %the frequency of child signal 1
f2 = 1.5; %the frequency of child signal 2 Unit Hz
Fs = 50; %sampling rate , Unit Hz
delta_t = 1/Fs;
% t = 0:delta_t:1-delta_t; %the total time length is 1s, contains 1 peroid of singal 1 and 1.5 peiod of singal 2
t = 0:delta_t:2-delta_t; %the total time length is 2s, contains 2 peroid of singal 1 and 3 peiod of singal 2
s1 = sin(2*pi*f1.*t); %child signal 1
s2 = sin(2*pi*f2.*t); %child signal 2
s = s1 + s2; %signal;
N = length(t); %signal length and frequency axis length;
fft_s1 = fft(s1,N); %fft of child signal 1
fft_s2 = fft(s2,N); %fft of child signal 2
fft_s = fft(s,N); %fft of signal
delta_f = Fs/N; %resolution of frequncy
axis_f = (0:N-1)*delta_f; %frequency axis
figure(1);
subplot(2,2,1)
plot(t,s1); %plot signal 1
hold on
plot(t,s2,'-'); %plot signal 2
xlabel('time s');
ylabel('signal');
title('time domain wave');
legend('child signal 1','child signal 2');
subplot(2,2,3)
plot(axis_f, abs(fft_s1)); %spectrum of signal 1;
hold on
plot(axis_f, abs(fft_s2),'-'); %spectrum of signal 2;
xlabel('frequency Hz');
ylabel('fft result');
title('observation time 1s');
legend('child signal 1','child signal 2');subplot(2,2,2)
plot(t,s); %plot signal
xlabel('time s');
ylabel('signal');
title('time domain wave');
subplot(2,2,4)
plot(axis_f, abs(fft_s)); %spectrum of signal 1;
xlabel('frequency Hz');
ylabel('fft result');
title('observation time 1s');% %separately observe the spectrum of signal that not full period
% f1 = 1; %the frequency of child signal 1
% f2 = 1.5; %the frequency of child signal 2 Unit Hz
% Fs = 50; %sampling rate , Unit Hz
% delta_t = 1/Fs;
% % t = 0:delta_t:1-delta_t; %the total time length is 1s, contains 1 peroid of singal 1 and 1.5 peiod of singal 2
% t = 0:delta_t:2-delta_t; %the total time length is 2s, contains 2 peroid of singal 1 and 3 peiod of singal 2
% s1 = sin(2*pi*f1.*t); %child signal 1
% s2 = sin(2*pi*f2.*t); %child signal 2
% s = s1 + s2; %signal;
% N = 500000; %fft number in order to observe more detailed
% fft_s1 = fft(s1,N); %fft of child signal 1
% fft_s2 = fft(s2,N); %fft of child signal 2
% fft_s = fft(s,N); %fft of signal
% delta_f = Fs/N; %resolution of frequncy
% axis_f = (0:N-1)*delta_f; %frequency axis
% figure(1);
% subplot(2,2,1)
% plot(t,s1); %plot signal 1
% hold on
% plot(t,s2,'-'); %plot signal 2
% xlabel('time s');
% ylabel('signal');
% title('time domain wave');
% legend('child signal 1','child signal 2');
% subplot(2,2,3)
% plot(axis_f, abs(fft_s1)); %spectrum of signal 1;
% hold on
% plot(axis_f, abs(fft_s2),'-'); %spectrum of signal 2;
% xlabel('frequency Hz');
% ylabel('fft result');
% title('observation time 1s');
% legend('child signal 1','child signal 2');
%
%
% subplot(2,2,2)
% plot(t,s); %plot signal
% xlabel('time s');
% ylabel('signal');
% title('time domain wave');
% subplot(2,2,4)
% plot(axis_f, abs(fft_s)); %spectrum of signal 1;
% xlabel('frequency Hz');
% ylabel('fft result');
% title('observation time 1s');
傅里叶变换—信号观测时长和频率分辨率相关推荐
- python判定固定时长固定频率的音频是否连续
我有一段时长大约为3.34 s 频率为1khz 的连续音频.该音频输入一个设备,经转换输出后,我想判断该音频是否还连续. 好了,这就是我的需求,看起来比较棘手,但我们稍稍了解音频采样率,码流,音频格式 ...
- Matlab FFT变换细节(信号采样频率,FFT变换点数,频率分辨率)
问题: 在做深度学习的故障诊断中,发现代码直接将原始信号fft之后直接将实频域信号输入网络中进行诊断,虽说效果比较不错95% 但因为输入的是双边谱且频率范围远超故障特征频率同时由于单个样本的点数只有1 ...
- matlab 采样点数,信号频率、采样频率、频率分辨率以及FFT信号补零
采样点数,信号频率.采样频率.采样点数的区分 包含matlab代码讲解示例 清晰明了 采样点数,信号频率.采样频率.采样点数 首先,频率指的是物质在单位时间内完成周期性变化的次数叫做频率,常用f表示. ...
- 关于频谱分析中两个重要指标:频率分辨率和时间分辨率的理解及计算
当今最常见时频分析方法主要有四种,分别是基于短时傅立叶变换法,基于小波变换法,Choi-Williams分布法和Hilbert-Hang变换法,经实验测得Hilbert-Huang具有最高的频率分辨率 ...
- 时间分辨率、频率分辨率
时间分辨率 时域图横坐标上的最小时间间隔.其大小代表着对模拟信号采样后得到的离散信号中点的密集程度. 数学公式: T s = 1 F s {T_s} = \frac{1}{{{F_s}}} Ts=F ...
- 【转载】采样频率、采样点数、频率分辨率
原作者: 辰忆书阁 原链接:http://www.360doc.com/myfiles.aspx?reg=1&app=1&type=3 1.频率分辨率的2种解释 解释一:频率分辨率可以 ...
- ffmpeg获取视频时长和分辨率
文章目录 ffmpeg获取视频文件时长和分辨率 获取视频时长 获取视频分辨率 ffmpeg获取视频文件时长和分辨率 获取视频时长 import os import subprocess# video_ ...
- 【频域分析及处理】1. 数字角频率w、模拟角频率Ω、频率分辨率
数字角频率.模拟角频率的 由来 已知有某实际存在的原始 模拟信号: y ( t ) = s i n ( Ω t ) = s i n ( 2 π f t ) y(t)=sin(Ωt)=sin(2\pi ...
- matlab混叠现象与频率分辨率,连续时间信号频谱分析研究及MATLAB实现
0.引言在信号处理过程中,频域分析方法往往比时域分析方法更方便和有效.对于确知连续时间信号,其频域分析可以通过连续时间傅里叶变换来进行,但是,这样计算出来的结果仍然是连续函数,计算机不能直接加以处理. ...
最新文章
- 【组队学习】【31期】青少年编程(Scratch 四级)
- php arrayaccess 二维,php的ArrayAccess(数组式访问接口)
- python - 栈与队列(只有代码)
- sqlserver2000 php5.4,PHP添加Mssql/sqlserver2000扩展,适用于Php5.2/Php5.3/Php5.4
- php小偷程序原理分析2
- 将数据归一化到任意区间范围的方法
- Mongodb存储特性与内部原理
- AndroidStudio_安卓原生开发_Android开发中界面调试很别扭? 设置应用屏宽屏高_应用大小_design_width_in_dp---Android原生开发工作笔记140
- 人生 第一页HTML
- python怎么解微分方程组_python能解微分方程吗
- CSS实现动态文字效果
- 使用poi解析word转html,并处理word中图片
- Java动态代理机制原理详解(JDK 和CGLIB,Javassist,ASM)
- 弹性伸缩系统的技术演进与落地实践
- 直接拿来用!最火的前端开源项目
- 分享如何将多个pdf文件合并成一个pdf?
- c语言常用函数库 c++常用函数库
- 抖音表白程序Python版,明人不说暗话,我喜欢你
- 万字详文告诉你如何做 Code Review!
- cisp题库700道(带答案)