plecs使用C-Script模块实现线性插值算法

提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档


一、定义插值函数

在plecs打开C-Script代码编写模块

double lininterp(double *x, double *y, int vectorLength, double xValue)

*x横坐标向量,c语言中没有现成的向量,故用指针来实现。

*y纵坐标向量。

又因为c语言里传递指针时无法传递指针的长度,所以还需要传输点的个数的数量vectorLength。

最后还需要插值的横坐标xValue。所以就是根据横纵坐标向量和向量长度,得到想要差值横坐标的纵坐标。

double lininterp(double *x, double *y, int vectorLength, double xValue)
{
	int index = getIndex(x, vectorLength, xValue);
	if(index == vectorLength - 1)
		return LINEAR_INTERPOLATION(x[index-1], y[index-1], x[index], y[index], xValue);
	else
		return LINEAR_INTERPOLATION(x[index], y[index], x[index+1], y[index+1], xValue);
}

getIndex这个函数的功能是求出xValue在那两个差值之间,并返回左边那个点的索引值,赋给了index这个变量

除了遇到最右的区块,返回最右边的索引:index == vectorLength - 1。用最右的点与倒数第二个点进行差值。否则会溢出

在这里插入图片描述

宏命令LINEAR_INTERPOLATION实现知道X1,X2,Y1,Y2和X,返回Y

在这里插入图片描述

二、定义函数

1.getIndex

代码如下(示例):

int getIndex(double *x, int vectorLength, double xValue)
{
	int minIndex = 0;
	int maxIndex = vectorLength - 1;
	int middleIndex;
	
	if(xValue < x[minIndex])
		return minIndex;
	else if(xValue > x[maxIndex])
		return maxIndex;
		
	while(1)
	{
		if(maxIndex - minIndex <= 1)
			return minIndex;
		middleIndex = (minIndex + maxIndex)/2;
		if(x[middleIndex] > xValue)
			maxIndex = middleIndex;
		else
			minIndex = middleIndex;
	}
}

2.定义宏命令

代码如下(示例):

#define LINEAR_INTERPOLATION(X1, Y1, X2, Y2, X) ( (Y2-Y1)/(X2-X1)*(X-X1) + Y1 )

#define VECTOR_LENGTH ( ParamDim(0, 0)*ParamDim(0, 1) )

分别定义插值和向量长度。向量长度是行数与列数的乘积

3.指针的获取

setup处
在这里插入图片描述
code declaration部分,获取
在这里插入图片描述
code中start function code部分,分配内存,并赋值
在这里插入图片描述
在仿真结束时释放内存,在terminate function code
在这里插入图片描述

三、返回

调用lininterp函数

Output(0) = lininterp(x, y, VECTOR_LENGTH, fmod(CurrentTime, x[VECTOR_LENGTH - 1]));

总结

提示:这里对文章进行总结:
例如:以上就是今天要讲的内容,本文仅仅简单介绍了pandas的使用,而pandas提供了大量能使我们快速便捷地处理数据的函数和方法。

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值