我这个程序红色部分显示的是输入参数数目不足,求大神给个建议,解决一下。
% params = enkf(datafilename,numEns,numsteps,paramsFile,data,outflowVesselIDs,N)
%
% this code uses an Ensemble Kalman Filter to optimize the windkessel
% parameters associated with the CoW code
%
% datafilename: the output filename
% numEns: number of ensemble members to be used
% numsteps: number of timesteps to be taken (total time = numsteps*.05)
% paramsFile: designates the parameter file to be used to setup the CoW
% data: NEED TO FILL THIS IN********************
% N: number of nodes per vessel
% if omitted, it defaults to four
%
% params: optimized parameter values (laca,raca,lmca,rmca,lpca,rpca)
%
% NOTE: AS WRITTEN THIS CODE ONLY WORKS WITH 6 OUTFLOW DATA SETS!!
%
% junctions file and inflow files are hard-coded as:
% junctionsFile.txt
% postprocesseddata/BasilarInflow.txt
% postprocesseddata/LCarotidInflow.txt
% postprocesseddata/RCarotidInflow.txt
%
% data files are hard-coded as:
% postprocesseddata/LACAOutflow.txt
% postprocesseddata/RACAOutflow.txt
% postprocesseddata/LMCAOutflow.txt
% postprocesseddata/RMCAOutflow.txt
% postprocesseddata/LPCAOutflow.txt
% postprocesseddata/RPCAOutflow.txt
%
% File Created: 04/08/08
% Last Edited: 04/09/08
function params = enkf(datafilename,numEns,numsteps,paramsFile,data,N)
if ~exist('data','var')
data = [{'postprocesseddata/LACAOutflow.txt',16}
{'postprocesseddata/RACAOutflow.txt',15}
{'postprocesseddata/LMCAOutflow.txt',11}
{'postprocesseddata/RMCAOutflow.txt',10}
{'postprocesseddata/LPCAOutflow.txt',5}
{'postprocesseddata/RPCAOutflow.txt',4}];
end
if ~exist('N','var')
N = 4;
end
numDatasets = length(data);
if ~exist('outflowVesselIDs','var')
outflowVesselIDs = [16 15 11 10 5 4];
end
outflowLocations = outflowVesselIDs*N-1;
dataLocations = [data{:,2}]*N-1;
disp(sprintf('Running Ensemble Kalman Filter (EnKF) with %d ensemble members over %d timesteps',numEns,numsteps))
%filenames etc.
junctionsFile = 'junctionsFile.txt';
basFile = 'postprocesseddata/BasilarInflow.txt';