我在插值数据文件时遇到了问题,我已经将其从.csv格式转换为X数组和Y数组,其中X[0]对应于点Y[0]。
我需要在这些值之间进行插值以得到最终的平滑文件。 我正在使用Picoscope输出函数,这意味着每条线在时间上是等间隔的,因此只使用Y值,这就是当您看到我的代码时我尝试以奇怪的方式操作的原因。
需要处理的值的类型为:
X Y
0 0
2.5 0
2.5 12000
7.5 12000
7.5 3000
10 3000
10 6000
11 6625.254
12 7095.154
如果两个相邻的 Y 值相同,它们之间就是一条直线;但当它们变化时,比如从 X = 10 往后,就会成为正弦波。
我一直在查找插值等方程式和其他帖子,但我已经好多年没有做过这种数学了,很遗憾我无法再理解它,因为有两个未知数,我不知道该如何编写 C#程序。
我的数据是:
int a = 0, d = 0, q = 0;
bool up = false;
double times = 0, change = 0, points = 0, pointchange = 0;
double[] newy = new double[8192];
while (a < sizeoffile-1 && d < 8192)
{
Console.Write("...");
if (Y[a] == Y[a+1])//if the 2 values are the same add correct number of lines in to make it last the correct time
{
times = (X[a] - X[a + 1]);//number of repetitions
if ((X[a + 1] - X[a]) > times)//if that was a negative number try the other way around
times = (X[a + 1] - X[a]);//number of repetitions
do
{
newy[d] = Y[a];//add the values to the newy array which replaces y later on in the program
d++;//increment newy position
q++;//reduce number of reps in this loop
}
while (q < times + 1 && d < 8192);
q = 0;//reset reps
}
else if (Y[a] != Y[a + 1])//if the 2 values are not the same interpolate between them
{
change = (Y[a] - Y[a + 1]);//work out difference between the values
up = true;//the waveform is increasing
if ((Y[a + 1] - Y[a]) > change)//if that number was a negative try the other way around
{
change = (Y[a + 1] - Y[a]);//work out difference bwteen values
up = false;//the waveform is decreasing
}
points = (X[a] - X[a + 1]);//work out amount of time between given points
if ((X[a + 1] - X[a]) > points)//if that was a negative try other way around
points = (X[a + 1] - X[a]);///work out amount of time between given points
pointchange = (change / points);//calculate the amount per point in time the value changes by
if (points > 1)//any lower and the values cause errors
{
newy[d] = Y[a];//load first point
d++;
do
{
if (up == true)//
newy[d] = ((newy[d - 1]) + pointchange);
else if (up == false)
newy[d] = ((newy[d - 1]) - pointchange);
d++;
q++;
}
while (q < points + 1 && d < 8192);
q = 0;
}
else if (points != 0 && points > 0)
{
newy[d] = Y[a];//load first point
d++;
}
}
a++;
}
这会生成一个关闭的波形,但仍然非常阶梯状。
所以,有人能看出为什么这不太准确吗? 如何提高其准确性? 或者使用数组的不同方法?
感谢您的关注 :)