提取采样时间点答案

【问题标题】：Extracting sampled Time Points提取采样时间点
【发布时间】：2014-05-01 00:43:53
【问题描述】：

我有一个 matlab 曲线，我想从中绘制并找到 17 个不同时间样本的浓度值

以下是我想从中提取 17 个不同时间点的浓度值的曲线

以下是以分钟为单位的时间点

t = 0,0.25,0.50,1,1.5,2,3,4,9,14,19,24,29,34,39,44,49. minutes samples

以下是我为绘制上图而编写的函数

function c_t = output_function_constrainedK2(t, a1, a2, a3,b1,b2,b3,td, tmax,k1,k2,k3)



K_1   = (k1*k2)/(k2+k3);
K_2   = (k1*k3)/(k2+k3);
DV_free= k1/(k2+k3);


c_t = zeros(size(t));

ind = (t > td) & (t < tmax);


c_t(ind)= conv(((t(ind) - td) ./ (tmax - td) * (a1 + a2 + a3)),(K_1*exp(-(k2+k3)*t(ind)+K_2)),'same');

ind = (t >= tmax);


c_t(ind)= conv((a1 * exp(-b1 * (t(ind) - tmax))+ a2 * exp(-b2 * (t(ind) - tmax))) + a3 * exp(-b3 * (t(ind) - tmax)),(K_1*exp(-(k2+k3)*t(ind)+K_2)),'same');



plot(t,c_t);
axis([0 50  0 1400]);
xlabel('Time[mins]');
ylabel('concentration [Mbq]');
title('Model :Constrained K2');
end

如果可能的话，请给我一些想法，我可以如何改变上述函数，以便我可以得出上述 17 个不同时间点的浓度值

以下是我用来绘制曲线的输入值

 output_function_constrainedK2(0:0.1:50,2501,18500,65000,0.5,0.7,0.3,3,8,0.014,0.051,0.07)

【问题讨论】：

为什么不在调用函数时将您想要的时间点作为t 输入参数输入？
@David：对不起，我是 matlab 的初学者，我所拥有的是时间 t=50 分钟，这是扫描时间，我不可能更改该值，是的，我可以简单地输入不同的 t 值和调用该函数，但根据我目前的工作指南，它是不可接受的

标签： signal-processing matlab sampling

【解决方案1】：

这将为您提供所需时间点的浓度值。您必须将其放在 output_function_constrainedK2 函数中，以便您可以访问变量 t 和 c_t。

T=[0 0.25 0.50 1 1.5 2 3 4 9 14 19 24 29 34 39 44 49];
concentration=interp1(t,c_t,T)

【讨论】：

尊敬的先生，让我更准确地重新表述我的问题，在您在回答中所说的完全正确之前，我已经完成了相同的过程来获取不同时间间隔的值......我需要编写一个功能或代码来完成完全相同的工作，我坚信你的方式是绝对正确的，但作为我项目的一部分，我必须编写一个函数或代码来从曲线中提取浓度值在不同的时间点。这是我的问题
我编辑了我的答案，但我仍然不确定您真正要问的是什么，或者您是否尝试过解决这个问题。
问题在于对数据进行采样，如果我使用给定的输入值运行原始函数，我会获得 501 个采样点，现在我想创建 17 个不同的时间帧（基于帧的信号处理方法） ...而且您的工作方法是绝对正确的，但是根据项目，我必须编写一个“for循环”来提取不同的时间样本..您解决它的方式没有错，我也这样做了，但是当我处理 for 循环时，我遇到了错误，所以我发布了这个问题。如果可能的话，请建议我如何通过 for 循环解决同样的问题
为此提出一个新问题，举一个非常小的问题示例，例如您的真实问题，即小而独立。您现在提出的问题与您在此处发布的问题完全不同。
尊敬的先生，我正在做作业，所以根据规则我必须使用 for 循环提取浓度值，for 循环必须在索引 ind 之前，但这只是示例分配的一小部分，我已经像你一样接近了这个样本。如果可能的话，请告诉我如何从上面的代码中优化 K3 的值，初始化值 k3=0.07