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
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?
采纳的回答
Walter Roberson
2019-11-7
[~, idx] = min(abs(c1 - N));
t(idx)
2 个评论
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 Center 和 File Exchange 中查找有关 Logical 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!