diff options
Diffstat (limited to 'misdata.sci')
| -rw-r--r-- | misdata.sci | 41 |
1 files changed, 41 insertions, 0 deletions
diff --git a/misdata.sci b/misdata.sci new file mode 100644 index 0000000..c1ced50 --- /dev/null +++ b/misdata.sci @@ -0,0 +1,41 @@ +function varargout = misdata(varargin) + [lhs,rhs] = argn(0) +//------------------------------------------------------------------------------ +// checking the number of inputs + if rhs <> 1 then + error(msprintf(gettext("%s:Wrong number of input arguments.\n"),"misdata")) + end +//------------------------------------------------------------------------------ + ioData = varargin(1) + if typeof(ioData) <> "iddata" then + error(msprintf(gettext("%s: Plant input data must be ""iddata"" type. "),"misdata")) + end + inputMat = ioData.InputData;inputMat = linearINTRP(inputMat,abs(ioData.Ts));ioData.InputData = inputMat; + outputMat = ioData.OutputData;outputMat = linearINTRP(outputMat,abs(ioData.Ts));ioData.OutputData = outputMat; + varargout(1) = ioData +endfunction + +function varargout = linearINTRP(matData,Ts) + // looking for overall nan values + nanData = isnan(matData);nanIndex = find(nanData == %T) + if ~size(nanIndex,'*') then + varargout(1) = matData + else + tempMat = [] + matSize = size(matData,'r') + // looking for nan in each column + for ii = 1:size(matData,'c') + nanData = isnan(matData(:,ii));nanIndex = find(nanData == %T); + if ~size(nanData,'*') then + tempMat = [tempMat matData(,ii)] + else + timeData = (linspace(1*Ts,matSize*Ts , matSize))'; + nanMat = isnan(matData(:,ii)); + data = matData(:,ii) + data(nanMat) = interp1(timeData(~nanMat), data(~nanMat), timeData(nanMat)); + tempMat = [tempMat data] + end + end + varargout(1) = tempMat + end +endfunction |