MATLAB的答案

高度扭曲的系统的传递函数;影响锤测试和系统的频率响应

15的观点(30天)
丹尼尔
丹尼尔 2019年9月4日
编辑: 丹尼尔2019年9月5日
你好,每个人!
两周了,我处理传递函数估计。我目前研究的系统是高度扭曲的高频率由于其结构动力学。我分析它的频率响应和真的越来越混乱的2、5 kHz。
图- H.jpg
我想弥补这个系统传递函数的使用。我有一些输入和输出数据(附件你会发现图形说明输入和输出信号)。它通过影响锤进行了测试。削减长话短说:输入-锤冲击,系统的输出-响应(振荡)。
Data.jpg
我试图处理它与系统辨识工具箱。到目前为止我是:
2 d.png
你能提供给我任何暗示这个函数在时域怎么走吗?我的工作的目的是复制输入信号不失真。
非常感谢您的任何反馈。

在回答这个问题。

答案(1)

拉吉
拉吉 2019年9月4日
不知道为什么你的麻烦发这么多细节。据我所知,一行的问题是,你有一个传递函数S域和您要将其转换为时间域。您可以使用 ilaplace 命令。
注意:你应该为了使用这个符号数学工具箱。
7评论
显示 6年长的评论
丹尼尔
丹尼尔 2019年9月5日
编辑:丹尼尔 2019年9月5日
关闭所有
清晰的所有
clc
t = linspace (0, 200000);%采样率:100千赫;t -时间[s]
t = t ';
%导入数据,减少样本的数量
输入= importdata (“C: \用户Lechowicz \ \电脑\ input.txt”,”);
输出= importdata (“C: \用户Lechowicz \ \电脑\ output.txt”,”);
%画出数据
图(“名字”,“原始数据”,“NumberTitle”,“上”);
次要情节(2,1,1)
情节(t,输入)
标题(输入信号的);
包含(“时间[s]”);
ylabel (“迫使[N]”);
次要情节(2,1,2)
情节(t,输出)
标题(的输出信号);
包含(“时间[s]”);
ylabel (“迫使[N]”);
%传递函数估计和转换到时间域
data = iddata(输出,输入,0.00001);
sys =特遣部队(数据、3、2);
(num,窝)= tfdata (sys,“v”);
信谊年代
tf_s = poly2sym ((num), s) / poly2sym((穴),s);
%得到传递函数为时间域
time_domain_tf = ilaplace (tf_s);
time_domain_tf = vpa (time_domain_tf);
%建立时域矩阵的传递函数
为我= (linspace (0, 200000))”
tf_time_domain_matrix = 10.753107173906815697416094377728 . * exp(-3.327484297190415992653272529398。*我)- 840.46923468767926446235542117068 . * exp(-246.19890346395301270080891184263 *我)。* cos(19126.216965822094847808546864243。*我)- 14175.620532280572917636521940787。* exp(-246.19890346395301270080891184263 *我)。*罪(19126.216965822094847808546864243 . *);
结束
%我几乎确定我在这里搞砸了
output_corrected = conv (tf_time_domain_matrix、输出);
%画出数据
图(“名字”,“复制输入信号”,“NumberTitle”,“上”);
情节(linspace(0、2、长度(output_corrected)), output_corrected);
标题(“复制输出信号”);
包含(“时间[s]”);
ylabel (“迫使[N]”);
我放弃了。你能指出我的错误吗?
据我理解,纠正我的输出数据我应该得到类似的输入数据,即单一的影响。我不能实现它的坚果。
我也附上我的数据。
总之,越来越好!非常感谢你的帮助。

登录置评。


翻译的