hilbert_function_constraint.m~

来自「一个用EM算法的源程序」· M~ 代码 · 共 23 行

function result = Hilbert_function_constraint(ambientDimension, subspaceDimensions);

load rankTable.mat
subspaceNumber = length(subspaceDimensions);
subspaceCodimentions=ambientDimensions*ones(1,subspaceNumber)-subspaceDimensions;
subspaceCodimentions=sort(subspaceCodimensions);

% Remove Hyperplances
while (subspaceNumber>2) || (subspaceDimensions(1)==1)
    subspaceCodimensions=subspaceCodimensions(2:end);
    subspaceNumber = subspaceNumber - 1;
end

codimensionCombinations=rankTable{ambientDimension,subspaceNumber,1};
charDimensions = rankTable{ambientDimension,subspaceNumber,2};

% Find the right entry
entryFound = false;
entryIndex = 0;
while (entryFound==false) && (entryIndex<size(codimensionCombinations,1))
    entryIndex = entryIndex + 1;
    if codimensionCombinations(entryIndex,:)==
end