请教一个问题(内详)
12-12
一个线性系统,系统函数分别是 H(s)/ h(t)
当频率采样间隔(df)太大的时候,T = 1/df 过于小, h(t)在t = T的时候还没有到零(该系统的冲击响应时间太长了。。。。,我要用h(t)做卷积,T较小的时候不单单是h(t)没到零,而是基本不对。)
但是如果 df太小,我又需要计算过多的频域点(fc/dt, fc是H(s)的带宽),每计算一个频域点都非常的费时
现在的问题非常的尴尬,我觉得较大的df就已经能够很好的反映H(s)的频域特征,但是其反变换所得的h(t)却与实际相差甚远(还是由于1/dt太小)
我现在想这样解决,非常山寨,不知道有没有相应的理论支持,我用大的df计算H(s)的值,然后再两个相邻的频域采样点之间线性插值,用插值之后的H(s)反变换求h(t).这样既计算较少的频点,又可以得到较长的冲击响应时间(1/df).
我想知道有没有相应的方法解决我遇到的这个问题,我这种山寨的方法有没有相应的理论?
=========================
可能的解决方案,采用较大的df,不能太大,至少能把峰抓住
用pade逼近或是vector-fitting 去拟合H(s)得到可物理实现的形式H1(s)
对H1(s)求解析的反变换,得到h1(t),这样就差不多可以了。
※ FROM: 59.66.211]
※ 来源:·水木社区 newsmth.net·[FROM: 211.99.222]
当频率采样间隔(df)太大的时候,T = 1/df 过于小, h(t)在t = T的时候还没有到零(该系统的冲击响应时间太长了。。。。,我要用h(t)做卷积,T较小的时候不单单是h(t)没到零,而是基本不对。)
但是如果 df太小,我又需要计算过多的频域点(fc/dt, fc是H(s)的带宽),每计算一个频域点都非常的费时
现在的问题非常的尴尬,我觉得较大的df就已经能够很好的反映H(s)的频域特征,但是其反变换所得的h(t)却与实际相差甚远(还是由于1/dt太小)
我现在想这样解决,非常山寨,不知道有没有相应的理论支持,我用大的df计算H(s)的值,然后再两个相邻的频域采样点之间线性插值,用插值之后的H(s)反变换求h(t).这样既计算较少的频点,又可以得到较长的冲击响应时间(1/df).
我想知道有没有相应的方法解决我遇到的这个问题,我这种山寨的方法有没有相应的理论?
=========================
可能的解决方案,采用较大的df,不能太大,至少能把峰抓住
用pade逼近或是vector-fitting 去拟合H(s)得到可物理实现的形式H1(s)
对H1(s)求解析的反变换,得到h1(t),这样就差不多可以了。
※ FROM: 59.66.211]
※ 来源:·水木社区 newsmth.net·[FROM: 211.99.222]
你的名字很酷,炸眼看着熟悉,仔细一个还别有洞天
理论依据不足。
因为你的h(t)本身的响应持续时间就很长,这本身就说明了H(w)的带宽很窄,或者说在频域上H(w)存在比较陡峭的突变点。例如理想矩形滤波器的时域响应就很长,是因为H(w)在其在截至频率处是跳变的。那么,频域的采样间隔df就必须很小,才能采集到这些跳变的地方。因此必然就有采样时间T很长。
的确是这样,H(w)频带很窄,就是在低 中 高频段有几个谐振点。
问题是我的输入含有较多的高频分量,必须要有H(w)的高频相应,晕啊
如果H(w)主要是离散的分布在频谱上的若干的频点上的,那么可以考虑减少时间轴上的采样率,虽然这样会使H(w)在整个频带上重叠,但是只要那几个主要的谐振峰没有重叠在一起就可以。从理论上讲,带通信号的采样率只要是其带宽的两倍就够了。
射频专业培训教程推荐