世纪电源网社区logo
社区
Datasheet
标题
返回顶部
未解决

非线性曲线反推求助

[复制链接]
查看: 905 |回复: 26
1
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-4-29 16:20:57
10问答币
最近遇到一个数学问题需要反向推导曲线(或函数),具体问题是这样:
有一未知曲线f1(t),已知曲线f2(t)=|f1(t)-f1(t-τ)|,其中延迟时间τ是已知的,可否由已知的曲线f2(t)反推出曲线f1?

比如上图绿色的f2(t)曲线是已知的,f1(t)曲线的延迟时间τ也是已知的,如何反向绘制出曲线f1?(只需推出两虚线内的相对解)
附上曲线文件 反推曲线.rar (21.47 KB, 下载次数: 2)
收藏收藏3
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-4-29 17:02:26
 
有无可能呢?
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-4-29 18:50:29
 
思考了一下感觉延迟时间τ跟△t性质有些像所以对f1进行了求导运算发现求导结果可以同f2完全重合起来

反过来对f2积分应该就能反推出f1。
gnaijnaoul
  • 积分:1292
  • |
  • 主题:3
  • |
  • 帖子:141
积分:1292
LV6
高级工程师
  • 2019-4-30 08:48:15
 
信号延迟相减其实就是微分信号了,只是实际使用中如果有噪声,会有噪声放大效应。
可以去参考一下韩京清先生跟踪微分器(TD)试试。
另外你的f2是取得绝对值,那f1 岂不是有两种可能。

boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-4-30 09:04:47
 
非常感谢!实际应用中是存在噪声的,我去学习一下。
x轴的周期是已知的可以在此处加个判断来选择需要的那个结果。
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-4-30 13:51:29
 
看了一些资料,这种跟踪微分器(TD)应该是针对数字微分器的一种实用算法,同模拟微(积)分器一样理想的一阶微(积)分器是没办法使用的。
用Saber仿真软件做了下验证:

                            仿真电路

                           仿真结果
如果用的是一阶理想微分器对应上图bode图中的红色线,噪声频率越高增益越大,在实际电路中输出几乎都被噪声淹没。
采用“二阶”微分器后高频段的增益是可控的,选择恰当的参数可以得到较理想的微分结果。
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-4-30 10:52:25
 
重合只因f1的gradient很小。
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-4-30 11:29:27
 
一楼的Mathcad文件不是很准确把一个变量写成了常量,改正后的文件 反推曲线.rar (24.77 KB, 下载次数: 1)
重新验证后感觉是跟延迟时间τ有关,如果τ太大重合性会有点偏差

greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-4-30 13:51:43
 
是不是已知f1是一个三角加一个正弦,只是不知它们的权?
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-4-30 14:02:50
 
只是接近这样的趋势,实际是类似于二极管的那种非线性曲线而且受工作环境影响会发生变化。
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-4-30 14:29:58
 
即不是两个已知函数的叠加,是的话就好办。
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-4-30 16:19:45
 
我这里写了个推导过程可以用来验证延迟相减近似等于微分结果

求图中黑点处的斜率,这里用△y/△t近似的方式计算,如果原曲线线性度比较好△t(或者τ)可以取的大一些,如果线性度不好△t就不能太大否则结果偏差大。
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-4-30 16:57:41
 
是的, gradient 变化小,△t大小影响不大。如果f1不需要很exact,把它拟成一个Polynomial ,则容易多了。
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-4-30 21:03:38
 
f1可能需要很exact,其实f1是个中间计算过程,在已知τ1和两个f2的条件下最终要求解的是τ2值。

如上图τ1和τ1下的f2是已知的,τ2下的f2是已知的,我目前的想法是先求出f1再算τ2,不知您是否有好的求解方法?
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-5-1 10:23:21
 
到底f2是个什么函数? Polynomial, Transcendental, Trigonometric ? 可微分,可积分?
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-5-1 10:30:54
 
f2是个实测曲线理论上来自于f1(t)-f1(t-τ),具体属于什么函数说不太清楚,f1又是个不确定的非线性方程。
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-5-1 10:56:59
 
那起码是一组数据,既然这样,在目标时间范围内,把f2拟成一个polynomial 或什么的函数,可行么?

曾想过一个很烦的方法,在频域里弄,令f2和f1的Complex Fourier Coeff 分别是C2k和C1k,τ time delay,那么有C2k=(1-e^-j2πkτ/T)*C1k,得C1k,得f1 。
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-5-1 12:21:59
  • 倒数10
 
有一种非均匀采样的算法跟你的这种有些相似也是通过时域-频域-时域变换来分析不过需要高速采集卡和电脑,我是想找出一种简单点的低成本的方案。
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-5-1 16:14:45
  • 倒数9
 
能看看实际的f2数据或曲线吗,有点好奇。
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-5-1 20:22:49
  • 倒数8
 

f1和f2都是频率的函数,上图是f1分别为非线性和线性扫频时得到的f2结果,下面的绿色曲线是做的频-压转换(鉴频器或锁相环)。
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-5-2 11:37:55
  • 倒数7
 
和之前的波形大异其趣。怎么我觉得绿线是f2,f1是红线的modulating signal .
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-5-2 13:52:41
  • 倒数6
 
是的,可能我没表述清,红线绿线都是f2一个是真实的一个是其频率特性。包括之前的分析都是对频率特性的分析。
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-5-2 15:43:29
  • 倒数5
 
这样的话,f1 会是个 Sine function ?
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-5-2 17:16:47
  • 倒数4
 
是的,f1是个正弦波 实际为激光频率
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
  • 2019-5-2 18:22:37
  • 倒数3
 
测发射回波,测距离? 积分f2会是个什么样的 wave ?
boy59
  • boy59
  • 离线
  • LV10
  • 总工程师
  • 积分:16430
  • |
  • 主题:118
  • |
  • 帖子:2779
积分:16430
LV10
总工程师
  • 2019-5-3 09:37:45
  • 倒数2
 
是用来测距的,一般都是先将频率信号转化为电压信号再做处理,如果直接处理真实的频率波形还没操作过。

图中绿色曲线为积分f2(经F-V转换)后的结果。
greendot
  • 积分:31625
  • |
  • 主题:0
  • |
  • 帖子:8701
积分:31625
LV12
专家
最新回复
  • 2019-5-3 12:16:29
  • 倒数1
 
(之前这里说的,不要理会了)
热门技术、经典电源设计资源推荐

世纪电源网总部

地 址:天津市南开区黄河道大通大厦8层

电 话:400-022-5587

传 真:(022)27690960

邮 编:300110

E-mail:21dy#21dianyuan.com(#换成@)

世纪电源网分部

广 东:(0755)82437996 /(138 2356 2357)

北 京:(010)69525295 /(15901552591)

上 海:(021)24200688 /(13585599008)

香 港:HK(852)92121212

China(86)15220029145

网站简介 | 网站帮助 | 意见反馈 | 联系我们 | 广告服务 | 法律声明 | 友情链接 | 清除Cookie | 小黑屋 | 不良信息举报 | 网站举报

Copyright 2008-2024 21dianyuan.com All Rights Reserved    备案许可证号为:津ICP备10002348号-2   津公网安备 12010402000296号