从两个信号域计算频率响应时间在matlab

67(30天)

显示旧的评论

亚历克斯Dashevsky 2018年1月2日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/375211-calculate-frequency-response-from-two-signal-on-domain-time-in-matlab

编辑: Nghi阮2020年9月4日

你好,我在时间域数组样本的输入和输出信号。向量是2002个样本的大小。我想找到频率响应。我写[h, w] = freqz (y、x, 2002);图()图(abs (h))但我非常奇怪的图。我想看到过滤器。我要计算频率响应的智能手机(音频过滤器)我做错了什么?

2的评论
显示1年长的评论藏1年长的评论

亚历克斯Dashevsky 2018年1月2日

编辑:亚历克斯Dashevsky 2018年1月2日

x.mat
y.mat

x是y的输入和输出

[h, w] = freqz (x, y, 2002);情节(abs (h))

登录置评。

在回答这个问题。

答案(1)

明星黾 2018年1月2日

2
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/375211-calculate-frequency-response-from-two-signal-on-domain-time-in-matlab answer_298388

你不能得到任何有用的信息从你的数据,因为噪音。然而,使用正确的采样频率(你没有状态,或给一次向量),这将给你的所有信息可以得到从你的数据:

                             D1 =负载(x.mat);
                            
                             D2 =负载(y.mat);
                            
                             x = D1.x;
                            
                             y = D2.y;
                            
                             L =元素个数(x);%长度的信号
                            
                             Fs = 1;%采样频率
                            
                             Fn = f / 2;%奈奎斯特频率
                            
                             FTx = fft (x) / L;
                            
                             FTy = fft (y) / L;
                            
                             阵线= linspace(0, 1,修复(L / 2) + 1) * Fn;%频率向量
                            
                             4 = 1:长度(艘);%指数向量
                            
                             TFyx = FTy. / FTx;
                            
                             图(1)
                            
                             次要情节(2,1,1)
                            
                             semilogx(阵线,20 * log10 (abs (TFyx (Iv))))
                            
                             ylabel (“| H (f) | (dB)”)
                            
                             次要情节(2,1,2)
                            
                             semilogx(阵线,角(TFyx (Iv)) * 180 /π)
                            
                             ylabel (“阶段(°)”)
                            
                             包含(的频率(赫兹))

注意,您只能使用一个频率,所以如果你想在频率、频率响应的振幅 “y” 的振幅 “x” 。如果你想要的频率响应,输入必须是一个扫描信号(如啁啾)或宽带噪声(实质上是一个脉冲响应)。当你这样做,我的代码将为宽带工作输入信号产生宽带输出信号。

16条评论
显示15个年长的评论隐藏15年长的评论

明星黾 2018年1月2日

是的。我计算复杂的传递函数为:

TFyx = FTy. / FTx;

当我绘制数据在时域中,有两个信号的不同的振幅 180° 的阶段。

如果你想做回波分析, cceps (链接), rceps (链接)或相关功能可能是你想要的。然而,当我做了分析,我认为没有回音:

                                  L =元素个数(x);%长度的信号
                                 
                                  Fs = 4.8 e + 4;%采样频率
                                 
                                  Fn = f / 2;%奈奎斯特频率
                                 
                                  t = linspace (0, L, L) / Fs;%的时间向量
                                 
                                  Ts =意味着(diff (t));%采样间隔
                                 
                                  FTx = fft (x) / L;
                                 
                                  FTy = fft (y) / L;
                                 
                                  阵线= linspace(0, 1,修复(L / 2) + 1) * Fn;%频率向量
                                 
                                  4 = 1:长度(艘);%指数向量
                                 
                                  TFyx = FTy. / FTx;
                                 
                                  图(1)
                                 
                                  次要情节(2,1,1)
                                 
                                  semilogx(阵线,20 * log10 (abs (TFyx (Iv))))
                                 
                                  ylabel (“| H (f) | (dB)”)
                                 
                                  次要情节(2,1,2)
                                 
                                  semilogx(阵线,角(TFyx (Iv)) * 180 /π)
                                 
                                  ylabel (“阶段(°)”)
                                 
                                  包含(的频率(赫兹))
                                 
                                  图(2)
                                 
                                  次要情节(2,1,1)
                                 
                                  semilogx(阵线,20 * log10 (abs (FTx (Iv))))
                                 
                                  ylabel (“rm \ itF \ \ {x (t) \}’)
                                 
                                  网格
                                 
                                  次要情节(2,1,2)
                                 
                                  semilogx(阵线,20 * log10 (abs (FTy (Iv))))
                                 
                                  ylabel (“rm \ itF \ \ {y (t) \}’)
                                 
                                  包含(的频率(赫兹))
                                 
                                  网格
                                 
                                  图(3)
                                 
                                  情节(t t, x,, y)
                                 
                                  网格
                                 
                                  rc = rceps (y);
                                 
                                  cc = cceps (y);
                                 
                                  图(4)
                                 
                                  次要情节(2,1,1)
                                 
                                  情节(t, rc)
                                 
                                  轴([min (t) max (t) -0.1 + 0.2])
                                 
                                  标题(“真正的倒频谱”)
                                 
                                  次要情节(2,1,2)
                                 
                                  情节(t, cc)
                                 
                                  轴([min (t) max (t) 2 + 2])
                                 
                                  标题(复杂的倒频谱的)

明星黾 2018年1月2日

不。你会得到一个 H (f) 。

这条线计算复杂的传递函数:

TFyx = FTy. / FTx;

这些线条图传递函数级 “H (f)” 在上次要情节 ,较低的阶段次要情节 :

                                  图(1)
                                 
                                  次要情节(2,1,1)
                                 
                                  semilogx(阵线,20 * log10 (abs (TFyx (Iv))))
                                 
                                  ylabel (“| H (f) | (dB)”)
                                 
                                  次要情节(2,1,2)
                                 
                                  semilogx(阵线,角(TFyx (Iv)) * 180 /π)
                                 
                                  ylabel (“阶段(°)”)
                                 
                                  包含(的频率(赫兹))

复杂的传递函数定义的傅里叶变换输出(这里我假设 “y” )除以输入的傅里叶变换(我认为 “x” )。你可以计算出大小(使用腹肌功能(使用)和阶段角从复杂的传递函数函数),正如我在我的代码完成。

如果你想计算系统的频率响应,您必须提供一个广泛的频率输入,并记录系统的输入和输出在同一时间。然后你可以计算传递函数的时域数据通过傅里叶变换。

Nghi阮 2020年9月4日

编辑:Nghi阮 2020年9月4日

你好明星黾,

可能我有一个问题关于x轴的单位吗?为什么在赫兹,但不是在rad / s ?

非常感谢。

Nghi

登录置评。

在回答这个问题。

类别

信号处理信号处理工具箱转换、相关性和建模转换 Cepstral分析

信号处理 DSP系统工具箱滤波器设计与分析滤波器的分析

找到更多的在Cepstral分析在帮助中心和文件交换

社区寻宝

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

开始狩猎!

从两个信号域计算频率响应时间在matlab

2的评论
显示1年长的评论藏1年长的评论

答案(1)

16条评论
显示15个年长的评论隐藏15年长的评论

另请参阅

类别

标签

社区寻宝

从两个信号域计算频率响应时间在matlab

2的评论 显示1年长的评论藏1年长的评论

答案(1)

16条评论 显示15个年长的评论隐藏15年长的评论

另请参阅

类别

标签

社区寻宝

2的评论
显示1年长的评论藏1年长的评论

16条评论
显示15个年长的评论隐藏15年长的评论