C++实现整数和复数的一维线性插值，结果与matlab的interp1一样

已于 2022-02-26 17:20:17 修改
阅读量2.3k 收藏 33
点赞数 2
分类专栏： matlab转C 文章标签： c++ matlab
于 2021-08-02 12:51:32 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_47898198/article/details/119322487
版权

C++实现整数和复数的一维线性插值interp1

导言

最近在进行Qt开发，涉及大量的matlab转C的工作，其中包括插值滤波等，这里对matlab的interp1函数进行了研究并用C++进行了重写。
经对比，结果与matlab的interp1()函数生成的插值结果一致。（如果数据不在定义域，则按照线性插值进行插值【与matlab的定义域外插值策略处理结果一致】）

函数需求分析

这里我对工作中调用matlab的interp1()函数的实例为例子进行讲解：

  1. 原matlab程序 ：
	// Ti:已有的时刻的值
	// To:插入的时刻的值
	// I_Id: 已有的时刻对应的y值的实部
	// Q_Id:已有的时刻对应的y值的虚部
   outD=interp1(Ti,I_Id,To,'linear','extrap')+1j*interp1(Ti,Q_Id,To,'linear','extrap');

  1. interp1（）生成的结果
    在这里插入图片描述

  2. C++版interp1（）输入输出参数详情：

参数名介绍
Ti已有的点的x坐标的值
To插入的点的X坐标的值
I_Id已有的点实部的y坐标的值
Q_Id已有的点虚部的y坐标的值
Y插入的点的实部和虚部的Y坐标的值以复数的形式返回
  compoundNumber interp1(vector<double> Ti, vector<double> To, vector<double> I_Id, vector<double> Q_Id) 
  /*
  	参数介绍
  	输入：(数据类型都为double类型的vector数组)
  	请注意：我这里的Ti和To都是排好序的
  		Ti： 已有的点的x坐标的值
  		To： 插入的点的X坐标的值

  		I_Id:已有的点实部的y坐标的值
  		Q_Id:已有的点虚部的y坐标的值

  	公式：Y[i] = y[i-1] + (y[i] - y[i-1]) * (X[i] - x[i-1])

  	输出：插入的点的实部和虚部的Y坐标的值以复数的形式返回

  			数据类型为自定义复数，
  			实部虚部都是一个double类型的vector数组
  */
  1. C++版interp函数生成结果 ：结果表明，生成结果与matlab生成结果一致。
    在这里插入图片描述
    在这里插入图片描述

源码

那么废话不多说，直接上代码

   #include <vector>
   #include <iostream>
   using namespace std;
   const struct comDouble
   {
   	vector<double> real;
   	vector<double> imag;
   
   	double time;
   };
   //复数插值
   comDouble Cominterp1(vector<double> Ti, vector<double> To, vector<INT16> I_Id, vector<INT16> Q_Id) {
   	/*
   	参数介绍
   	输入：(数据类型都为double类型的vector数组)
   	请注意：我这里的Ti和To都是排好序的
   		Ti： 已有的点的x坐标的值
   		To： 插入的点的X坐标的值

   		I_Id:已有的点实部的y坐标的值
   		Q_Id:已有的点虚部的y坐标的值

   	公式：Y[i] = y[i-1] + (y[i] - y[i-1]) * (X[i] - x[i-1])

   	输出：插入的点的实部和虚部的Y坐标的值以复数的形式返回

   			数据类型为自定义复数，
   			实部虚部都是一个double类型的vector数组
   	*/
   	comDouble Y;
   	if (Ti.size() != I_Id.size() || I_Id.size() != Q_Id.size())
   		cout << "采样点个数和信号时刻数量不一致或实虚部数量不一致！" << endl;

   	double maxTi = Ti[Ti.size() - 1];
   	double minTi = Ti[0];
   	double maxTo = To[To.size() - 1];
   	double minTo = To[0];

   	double theK = 0;
   	int before = 1;
   	double real, imag;
   	for (int i = 0; i < To.size(); i++) {
   		待插入点在定义域左边
   		//if (To[i] <= minTi) {
   		//	theK = (I_Id[1] - I_Id[0]) / (Ti[1] - Ti[0]);
   		//	real = theK * (To[i] - Ti[0]) + I_Id[0];
   		//  imag = theK * (To[i] - Ti[0]) + I_Id[0];
   		//	Y.real.push_back(real);
   		//	Y.imag.push_back(imag);
   		//	break;
   		//}
   		//定义域右边
   		if (To[i] >= maxTi) {
   			double theRealK = (I_Id[I_Id.size() - 1] - I_Id[I_Id.size() - 2]) / (Ti[Ti.size() - 1] - Ti[Ti.size() - 2]);
   			double theImagK = (Q_Id[Q_Id.size() - 1] - Q_Id[Q_Id.size() - 2]) / (Ti[Ti.size() - 1] - Ti[Ti.size() - 2]);
   			real = theRealK * (To[i] - Ti[Ti.size() - 1]) + I_Id[I_Id.size() - 1];
   			imag = theImagK * (To[i] - Ti[Ti.size() - 1]) + Q_Id[Q_Id.size() - 1];
   			Y.real.push_back(real);
   			Y.imag.push_back(imag);
   		}
   		else {
   			for (int j = before; j < Ti.size(); j++) {
   				if (To[i] <= Ti[j]) {
   					real = (I_Id[j] - I_Id[j - 1]) / (Ti[j] - Ti[j - 1]) * (To[i] - Ti[j - 1]) + I_Id[j - 1];
   					imag = (Q_Id[j] - Q_Id[j - 1]) / (Ti[j] - Ti[j - 1]) * (To[i] - Ti[j - 1]) + Q_Id[j - 1];
   					Y.real.push_back(real);
   					Y.imag.push_back(imag);

   					before = j;
   					break;
   				}
   			}
   		}
   	}
   	return Y;
   }

   #include <vector>
   #include <iostream>
   using namespace std;
   const struct comDouble
   {
   	vector<double> real;
   	vector<double> imag;
   
   	double time;
   };
   //整数插值
   vector<double> interp1(vector<double> Ti, vector<double> I_Id, vector<double> To)
   {
   	/*
   		参数介绍
   		输入：(数据类型都为double类型的vector数组)
   			Ti： 已有的点的x坐标的值
   			To： 插入的点的X坐标的值

   			I_Id:已有的点实部的y坐标的值

   		公式：Y[i] = y[i-1] + (y[i] - y[i-1]) * (x[i] - x[i-1])
   	*/
   	vector<double> Y;

   	double maxTi = Ti[Ti.size() - 1];
   	//double minTi = Ti[0];
   	//double maxTo = To[To.size() - 1];
   	//double minTo = To[0];

   	double theK = 0;
   	double real = 0.0;
   	int before = 1;
   	for (int i = 0; i < To.size(); i++) {
   		待插入点在定义域左边
   		//if (To[i] <= minTi) {
   		//	theK = (I_Id[1] - I_Id[0]) / (Ti[1] - Ti[0]);
   		//	real = theK * (To[i] - Ti[0]) + I_Id[0];
   		//	Y.push_back(real);
   		//}
   		//定义域右边
   		if (To[i] >= maxTi) {
   			theK = (I_Id[I_Id.size() - 1] - I_Id[I_Id.size() - 2]) / (Ti[Ti.size() - 1] - Ti[Ti.size() - 2]);
   			real = theK * (To[i] - Ti[Ti.size() - 1]) + I_Id[I_Id.size() - 1];
   			Y.push_back(real);
   		}
   		else {
   			for (int j = before; j < Ti.size(); j++) {
   				if (To[i] <= Ti[j]) {
   					real = (I_Id[j] - I_Id[j - 1]) / (Ti[j] - Ti[j - 1]) * (To[i] - Ti[j - 1]) + I_Id[j - 1];
   					Y.push_back(real);

   					before = j;
   					break;
   				}
   			}
   		}
   	}
   	return Y;
   }

性能提升（vector转用动态数组，这里以复数来举例）

   #include <vector>
   #include <iostream>
   using namespace std;
   struct comDoubleC
   {
   	double *real;
   	double *imag;
   	double time;
   	double Rsize;
   	double Isize;
   
   	comDoubleC(int reallen, int imaglen) {
   		this->real = new double[reallen];
   		this->imag = new double[imaglen];
   		this->time = 0;
   		this->Rsize = 0;
   		this->Isize = 0;
   	}
   };
   //复数插值
   comDoubleC Cominterp1(
           double Ti[],
           int Tilen,
           double To[],
           int Tolen,
           INT16 I_Id[],
           INT16 Q_Id[],
   		comDoubleC bh_IQ) {
       /*
       参数介绍
       输入：(数据类型都为double动态数组)
       请注意：我这里的Ti和To都是排好序的
           Ti： 已有的点的x坐标的值
           To： 插入的点的X坐标的值

           I_Id:已有的点实部的y坐标的值
           Q_Id:已有的点虚部的y坐标的值

   		bh_IQ：结果---在函数外管理内存


       公式：Y[i] = y[i-1] + (y[i] - y[i-1]) * (X[i] - x[i-1])

       输出：插入的点的实部和虚部的Y坐标的值以复数的形式返回

               数据类型为自定义复数，
               实部虚部都是一个double类型的动态数组
       */
       double real;
       double imag;
       double theRealK;
       double theImagK;
       double maxTi = Ti[Tilen - 1];
       int before = 1;

       for (int i = 0; i < Tolen; i++) {
           //待插入点在定义域右边
           if (To[i] >= maxTi) {
               theRealK = (I_Id[Tilen - 1] - I_Id[Tilen - 2]) / (Ti[Tilen - 1] - Ti[Tilen - 2]);
               theImagK = (Q_Id[Tilen - 1] - Q_Id[Tilen - 2]) / (Ti[Tilen - 1] - Ti[Tilen - 2]);
               real = theRealK * (To[i] - Ti[Tilen - 1]) + I_Id[Tilen - 1];
               imag = theImagK * (To[i] - Ti[Tilen - 1]) + Q_Id[Tilen - 1];
               bh_IQ.real[i] = real;
               bh_IQ.imag[i] = imag;
           }
           else {
               for (int j = before; j < Tilen; j++) {
                   if (To[i] <= Ti[j]) {
                       real = (I_Id[j] - I_Id[j - 1]) / (Ti[j] - Ti[j - 1]) * (To[i] - Ti[j - 1]) + I_Id[j - 1];
                       imag = (Q_Id[j] - Q_Id[j - 1]) / (Ti[j] - Ti[j - 1]) * (To[i] - Ti[j - 1]) + Q_Id[j - 1];
   					bh_IQ.real[i] = real;
   					bh_IQ.imag[i] = imag;

                       before = j;
                       break;
                   }
               }
           }
       }
       return bh_IQ;
   }

注意

该文章仅个人学习使用，欢迎大家一起交流学习

宇智波盆
关注
  • 2
    点赞
  • 踩
  • 33
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
MATLAB插值代码(一维插值interp1和二维插值interp2)
06-21
MATLAB插值代码(一维插值interp1和二维插值interp2) 本资源为一维插值interp1和二维插值interp2的使用方法代码, 为完整代码, 可直接运行 % yi = interp1(x,y,xi,method); % 已知样本点坐标x,y，求xi处的函数值yi，插值方法是method % method有以下几种： % 'nearest'邻近点插值 % 'linear'线性插值（默认） % 'spline'三次样条函数插值 % 'cubic'三次函数插值 % 常用的是'spline'和'cubic'
【无标题】c++ 实现 interpolate.interp1d插值
weixin_42731150的博客
10-13 984
c++ 实现interpolate.interp1d
matlab interpolation,matlab interpolation有什么作用
weixin_30009115的博客
03-20 1733
满意答案玫瑰园喜乐蒂推荐于 2013.07.26采纳率：56%等级：10已帮助：1021人展开全部 关于MATLAB 插值(Interpolation)2008-11-02 10:43插值：在已知数据之间计算估计值的过程。1． 一维插值(1D Interpolation)由interp1实现，用多项式技术计算插值点。Yi=interp1(x,y,xi,method)y—函数值矢量， x...
C++实现matlab中的interp1和interp2插值
nneerr123的专栏
02-09 7614
头文件interpfun.h#ifndef INTERPFUN_H#define INTERPFUN_H#include"math.h"#include&lt;stdio.h&gt;#include &lt;string.h&gt; #ifndef SafeDeleteVec#define SafeDeleteVec(X) { if((X)) delete (X); (X)=NULL;}#endi...
MATLAB2016笔记（七）：数据分析
qq_52441682的博客
08-04 1950
一、多项式及其函数 （一）多项式的表达式和创建 （二）多项式求根、由根创建多项式 （三）多项式的四则运算 （四）多项式的导数、积分、计算函数值 （五）多项式运算函数及操作指令 （六）有理多项式（部分分式展开，求导） 二、数据插值（近似表达式预测插值点） （一）概述：插值问题 （二）一维插值 （三）二维插值 三、函数极限 （一）极限概念 （二）求极限的函数 四、函数数值积分 （一）数值积分问题的数学 （二）一元函数的数值积分 （三）多重数值积分...
C++实现MATLAB函数interp
09-06
interp函数的实现需要调用filter函数（MATLAB内置），为此，我特地实现了自己的filter函数，可以从我上传的资源里面找到
Interp1 c++实现
weixin_42827083的博客
05-27 2647
在网上找了一下，有是有但是我下载下来用的时候结果不对。想修改一下但是搞得迷迷糊糊的，就干脆写了一个，不过只有最简单的线性插值的实现，新手可以直接拿过去用。也希望老鸟也可以不吝赐教，提高一下效率或者优化下结构。 cpp文件 //--------------------------------------------------------------------------- #pragma hdrstop #include "Interpfun.h" //------------------------
一维插值C++算法（三种方法）
weixin_37928884的博客
01-28 1331
一维插值C++算法（三种方法）
线性插值用c语言程序源代码,实现快速线性插值的实现
weixin_33277564的博客
05-24 1567
Fast Linear InterpolationF = lininterp1f(X,Y,XI,Ydefault) returns the value of the 1-D function Y at the points XI using linear interpolation. Length(F)=length(XI). The vector X specifies the coordina...
7种插值算法的c++代码实现
09-13
7种插值算法的c++代码实现，1 拉格朗日插值(POLINT) 2 有理函数插值(RATINT) 3 三次样条插值(SPLINE(二阶导数值)->SPLINT(函数值)) 4 有序表的检索法(LOCATE(二分法), HUNT(关联法)) 5 插值多项式(POLCOE(n2), POLCOF(n3)) 6 二元拉格朗日插值(POLIN2) 7 双三次样条插值(SPLIE2)
Interp1_C工程源码_interp1c源码_interp1_
10-02
与matlab结果一致的以为插值函数，结果完全一致。
基于C++的双线性插值算法设计与实现
05-12
基于C++的双线性插值算法设计与实现
三次样条差值c语言实现interp.zip
11-01
三次样条差值c语言实现，解方程部分对系数矩阵分屏内存进行了优化，适合大量数据使用
更快的一维线性插值：interp1qr：与内置 MATLAB 函数“interp1q”功能相同，但速度快 3 倍。-matlab开发
05-30
更快的一维线性插值：'interp1qr' 按照公式yi = y1 +（y2-y1）/（x2-x1）*（xi-x1），使用'x'和'y'对'xi'点进行一维线性插值，得到'yi'。 变量： - 'x' 是列向量 [mx 1]，单调递增。 -“ y”是矩阵[mxn]，对应于...
matlab最临近插值法代码-interp1:受MATLAB启发的一维数据插值
05-28
一维数据插值。 导入此功能 节点模块 运行npm install interp1 将一个导入添加到npm包中， import interp1 from 'interp1'; 然后，您可以在代码中使用该函数。 CDN 将以下脚本标记[removed][removed]放在index.html...
matlab中二维插值函数interp2的使用详解
09-17
主要介绍了matlab中二维插值函数interp2的使用详解，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
插值算法与matlab代码 - 副本.zip_interp1的代码_skyvi2_spline_三次样条插值_插值 matlab
07-15
MATLAB中的插值函数为interp1，其调用格式为： yi= interp1(x,y,xi,'method') 其中x，y为插值点，yi为在被插值点xi处的插值结果；x,y为向量， 'method'表示采用的插值方法，MATLAB提供的插值方法有几种： 'method'...
C++实现复数的滤波filter
qq_47898198的博客
08-04 4035
C++实现复数的滤波filter导言函数需求分析源码注意 导言 最近在进行Qt开发，涉及大量的matlab转C的工作，其中包括插值滤波等，这里对matlab的filter函数进行了研究并用C++进行了重写。 经对比，结果与matlab的filter()函数生成的滤波结果一致。 函数需求分析 这里我对工作中调用matlab的filter()函数的实例为例子进行讲解： 原matlab程序 ： temp1_i=filter(lbf,1,real(x_d(1:lbn+LEN))); temp1_q=
matlab一维线性插值
最新发布
06-01
在Matlab中进行一维线性插值可以使用interp1函数。该函数的语法如下： ``` Vq = interp1(X,V,Xq) ``` 其中，X是原始数据的横坐标向量，V是原始数据的纵坐标向量，Xq是要进行插值的横坐标向量，Vq是插值后的纵坐标...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

宇智波盆

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值