MATLAB的答案

如何绘制相图FFT和代码时间历史down-sampled信号?

47个视图(30天)
奥马尔·阿里
奥马尔·阿里 2023年5月6日23:00
评论道: 明星黾 2023年5月8日在47
跑:
你好,
到目前为止,我已经完成了这个代码,给了我一个频率和振幅图。我ussure如何从这得到相图或时间的历史由给定速率采样信号等我nwould感激援助。
data = 0.829668000000000 [0
0.00833300000000000 - 0.829801000000000
0.0166670000000000 - 0.828810000000000
0.0250000000000000 - 0.829008000000000
0.0333330000000000 - 0.828352000000000
0.0416670000000000 - 0.827073000000000
0.0500000000000000 - 0.826807000000000
0.0583330000000000 - 0.826828000000000
0.0666670000000000 - 0.826857000000000
0.0750000000000000 - 0.827056000000000
0.0833330000000000 - 0.827139000000000
0.0916670000000000 - 0.828689000000000)
data = 12×2
0 0.8297 0.0083 0.8298 0.0167 0.8288 0.0250 0.8290 0.0333 0.8284 0.0417 0.8271 0.0500 0.8268 0.0583 0.8268 0.0667 0.8269 0.0750 0.8271
t =数据(:1);
s =数据(:,2);
t = 0.0083;
Fs = 1 / Ts;
Fn = f / 2;
L = 4000;
smean =意味着(s);
Fts = fft (s-smean) / L;
阵线= linspace(0, 1,修复(L / 2) + 1) * Fn;
4 = 1:元素个数(艘);
图
情节(阵线,abs (Fts (Iv)) * 2)
指数超过数组元素的数量。指标不得超过12。
网格
包含(“频率”)
ylabel (“振幅”)

在回答这个问题。

接受的答案

明星黾
明星黾 2023年5月7日2时38分
跑:
也许这- - - - - -
data = 0.829668000000000 [0
0.00833300000000000 - 0.829801000000000
0.0166670000000000 - 0.828810000000000
0.0250000000000000 - 0.829008000000000
0.0333330000000000 - 0.828352000000000
0.0416670000000000 - 0.827073000000000
0.0500000000000000 - 0.826807000000000
0.0583330000000000 - 0.826828000000000
0.0666670000000000 - 0.826857000000000
0.0750000000000000 - 0.827056000000000
0.0833330000000000 - 0.827139000000000
0.0916670000000000 - 0.828689000000000)
data = 12×2
0 0.8297 0.0083 0.8298 0.0167 0.8288 0.0250 0.8290 0.0333 0.8284 0.0417 0.8271 0.0500 0.8268 0.0583 0.8268 0.0667 0.8269 0.0750 0.8271
t =数据(:1);
s =数据(:,2);
t = 0.0083;
Fs = 1 / Ts;
Fn = f / 2;
L =大小(数据,1);
smean =意味着(s);
NFFT = 2 ^ nextpow2 (L);
FTs = fft ((s-smean)。*损害(L)) / L;
阵线= linspace (0, 1, NFFT / 2 + 1) * Fn;
4 = 1:元素个数(艘);
图
次要情节(2,1,1)
情节(阵线,abs (FTs (Iv)) * 2)
网格
ylabel (“振幅”)
次要情节(2,1,2)
情节(阵线,rad2deg(角(FTs (Iv))))
网格
包含(“频率”)
ylabel (“阶段(°)”)
也可以考虑使用 打开 弧度的阶段,为: “rad2deg(打开(角(FTs (Iv))))” 。(我在原始代码做了一些改进。)
。
4评论
显示3年长的评论
明星黾
明星黾 2023年5月8日在47
像往常一样,这是我的荣幸!
相反,我也会这么做
F = 8
Fs = 30
(y、x) =黄油(5 F / (Fs / 2)“低”)
inputSignal = s;
outSignal = filtfilt (y, x, inputSignal);
图
情节(t, outSignal)
没有理由包括向量的时间信号过滤,所以我取消了它。
但是如果你想连接两个 列 向量,使用方括号 [a, b] (尽管的逗号是可选的)如果你想连接两个 行 向量 [一个;b] 使用分号来创建一个第二行(垂直连接它们)。在所有情况下,vectos必须有相同的行和列维度。
实际上,我将设计不同的滤波器(我更喜欢椭圆滤波器,因为他们是计算效率更高),然而设计巴特沃斯滤波器我将开始 buttord 函数,那么 黄油 然而,生产zero-pole-gain输出,然后使用 zp2sos 函数产生second-order-section表示,和电话 filtfilt 为:
outSignal = filtfilt (sos, g, inputSignal);
过滤将更有效和过滤器将总是与这些变化是稳定的。
。

登录置评。

更多的答案(0)

类别

找到更多的在数字滤波在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!