GPOPS-II Interpolation Function Handle returned by OptimTraj · MatthewPeterKelly/OptimTraj#20

(4 评论) (0 反应) (0 负责人)MATLAB (704 star) (219 fork)batch import

enhancementhelp wantedquestion

描述

In gpopsWrapper.m, OptimTraj provides an interpolation function handle for the GPOPS-II solution using pchip, as per the code below. Is pchip the most accurate way to interpolate the GPOPS-II solution?

soln.interp.state = @(t)( interp1(tSoln',xSoln',t','pchip',nan)' ); soln.interp.control = @(t)( interp1(tSoln',uSoln',t','pchip',nan)' );

The formula for the GPOPS-II state is given in equation (25) of the GPOPS-II paper:

Patterson, Michael A., and Anil V. Rao. "GPOPS-II: A MATLAB software for solving multiple-phase optimal control problems using hp-adaptive Gaussian quadrature collocation methods and sparse nonlinear programming." ACM Transactions on Mathematical Software (TOMS) 41.1 (2014): 1.

http://dl.acm.org/citation.cfm?id=2558904

贡献者指南

技术栈: 无
领域: data
议题类型: research
难度: 3
预计时间: half day
活动状态: stale
清晰度: clear
前置要求: MATLABnumerical interpolationGPOPS II paper
新手友好度: 30
研究方向: The issue questions the accuracy of pchip interpolation used in gpopsWrapper.m. A contributor should review equation (25) in the GPOPS II paper (Patterson & Rao, 2014) to determine if a more accurate interpolation method (e.g., spline or direct evaluation of polynomials) should be used. The file gpopsWrapper.m contains the relevant interpolation code. The contributor could implement and test alternative interpolation methods, comparing results against known solutions.