通过时域中的逆问题查找参数

1查看(最近30天)
穆罕默德霍奇
穆罕默德霍奇 于2020年5月14日
评论: 沃尔特罗伯森 2020年5月17日
亲爱的沃尔特,
首先,谢谢你花时间并回答我的所有问题。
我在S域,TR中有一个函数。我用durbin将它带入时域中。
该函数有几个参数。我将使用反问题找到一个或几个参数。
我有一个文件的另一个函数,我需要与它匹配我的函数,然后我可以找到我的参数。(最小化)
我担心的是知道我如何在这里使用逆问题吗?
我逃避了durbin,因为它们在S域中。
我决定在我的德林结束时添加几行关于逆问题(最小化),并使用函数名称启动我的Durbin代码[function y = f(x)]。
附加,您可以看到我的tr,f(包括durbin + verewore verse的更多行)和一个垫文件,包括我需要用它最小化的另一个函数的数据。此外,还有一个mfile,名称fint,因为它只是一个函数。
当我使用优化时,我将“位置1中的索引无效。数组指数必须是正整数或逻辑值。”
我经历了旧的姓氏并删除了这个问题。谢谢你。
现在,我有另一个问题,它正在杀了我。
它表示客观函数在初始点处未定义。粉甘群不能继续。“
请注意,我使用fmincon
如果您知道更好的功能,请告诉我。
请找到我的文件,看看内部
再次,谢谢你的时间和答案。
最好的事物
M.H.
1条评论
沃尔特罗伯森
沃尔特罗伯森 20月16日
你没有包括你的f。我可以猜测它调用tr,但我无法猜测您可能进入的s值,或者.mat如何连接到此。

登录评论。

接受答案

穆罕默德霍奇
穆罕默德霍奇 20月16日
亲爱的沃尔特,
请接受我的道歉。我很抱歉。
Attapced,请找到f.m。
如果您需要提出新问题,请询问,请让我这样做。
真挚地
M.H.

更多答案(3)

沃尔特罗伯森
沃尔特罗伯森 20月16日
tthrough = 8.8e-5:1E-9:Tmax;
这成为12001元素的向量
Thrange = tthrough(1):(长度(伸出))/ 100 :(长度(牵伸));
它使用8.8E-5作为第一个元素,增量12001/100 = 120.001至12001,总共100个元素。但 为什么 当为什么使用8.8e-5的时间值作为第一个值,否则将处理 长度 时间向量?
长度(特拉奇)
为了B = 1:长度(特拉奇)
DISP(b)
t =特兰奇(b);
对于第一次迭代,T将是Thrange(1),其是8.8E-5的牵引(1)。
对于第二次迭代,T将为120.001更高,约为120.010088。
Yt(b)=(2 * exp(mu。* t)/ t0)*(s + 0.5 *真实(tr(x,mu)));
MU接近25000.乘以120.010088获得3000102.19489026。exp()的exp()是溢出到inf的6.7 * 10 ^ 1302927。
???为什么你建立牵连只是为了取得其(已知)第一个值及其长度???
我怀疑你的意思是更像
Trange = Linspace(8.8e-5,Tmax,100);

穆罕默德霍奇
穆罕默德霍奇 20月16日
它使用8.8E-5作为第一个元素,增量12001/100 = 120.001至12001,总共100个元素。但 为什么 因为我需要从附加文件找到图15。
为什么使用8.8e-5的时间值作为第一个值,否则将处理 长度 时间向量?
请加载fellah2006.mat,然后你看到我分开了我的特兰奇就像你所看到的那样。实际上,我需要对两个都有相同的阵列。Fellah2006是我从图中获得的实验点,而我的代码有更多的poitn,所以我这样做了。
MU接近25000.乘以120.010088获得3000102.19489026。exp()的exp()是溢出到inf的6.7 * 10 ^ 1302927。
???为什么你建立牵连只是为了取得其(已知)第一个值及其长度???
你是对的。这是我的错误。我想 Trange = Linspace(8.8e-5,Tmax,100)将工作。
谢谢。
此外,我需要在我的论文中承认你。我会在夏天辩护。请让我知道你的名字。
你接受谈论你吗?如果是这样,我可以使用这个名字吗?
什么是最好的说法。
感谢ABND建议
M.H.
1条评论
沃尔特罗伯森
沃尔特罗伯森 2020年5月17日
因为我需要找到附加文件的图15
图15是一段时间递增120的任何东西都没有。图15的X轴高达约10.2 * 10 ^( - 5),不高达12001。
Fellah2006是我从图中获得的实验点
图15在它上面有两条线,并且您将两条线的数据混合在一起在Fellah2006中,导致具有同一时间坐标的多个条目(最后3次是例子。)

登录评论。


穆罕默德霍奇
穆罕默德霍奇 2020年5月17日
谢谢你的关注。
我想我已经解决了你给我的碎片问题的问题。
LINSPACE(8.8E-5,TMAX,100),TMAX = 10E-5。是不是?
M.H.
1条评论
沃尔特罗伯森
沃尔特罗伯森 2020年5月17日
是的,tmax = 10e-5在您的代码中。

登录评论。

标签

s manbetx 845

社区宝藏狩猎

找到Matlab Central中的宝藏,并发现社区如何帮助您!

开始狩猎!