plotting when variables are conditional

2 次查看(过去 30 天)
t = logspace(-3,13,100)
for j = 1:100
q = inv((At * Cdi * A) + (t(j) * Ht * H)) * At * Cdi * do
c1(j) = transpose((A * q) - do) * Cdi * ((A * q) - do)
end
m = inv((At * Cdi * A) + (t("C1 = N") * Ht * H)) * At * Cdi * do
depth = 1:50
plot(m,depth)
I want to plot the solution m againt depth
but for the the value of t when C1 is closest to N (i.e 20), and this value changes for each RUN
i believe an if-statement might help me achieve this, pls help
thanks
  2 个评论
Walter Roberson
Walter Roberson 2019-11-7
It looks unlikely to me that c1 will happen to be bit-for-bit identical to any integers.
Will you be wanting to find the index at which C1 is closest to N? Would you be wanting to find the last index at which C1 does not exceed N? The first index at which C1 is at least as high as N ? Do you want to interpolate the index at which c1 exactly equals N, and use that to interpolate a the corresponding t value?
Prince Igweze
Prince Igweze 2019-11-7
yea, ill like to use the value at which C1 is closest to N, and use the corresponding t value

请先登录,再进行评论。

采纳的回答

Walter Roberson
Walter Roberson 2019-11-7
[~, idx] = min(abs(c1 - N));
t(idx)
  2 个评论
Prince Igweze
Prince Igweze 2019-11-7
thanks but how do i implement this in the plot
Walter Roberson
Walter Roberson 2019-11-7
[~, idx] = min(abs(c1 - N));
m = inv((At * Cdi * A) + (t(idx) * Ht * H)) * At * Cdi * do;
It seems unlikely to me that you want depth to be your y axes. It seems more likely that you would want
plot(depth, m)

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Logical 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by