Я имею дело с проблемой DAE и для точности я намерен использовать схему дискретизации с 3 точками коллокации (скажем, я использую 10 конечных элементов).Затем я хочу использовать только решения на границах конечных элементов для вычисления данной метрики, потому что я считаю, что только там у меня есть измерения для сравнения.
Измерение с учетом всех точек: def _FIMmatrix(m, i, j):
return m.gFIM[i,j] == sum(sum(m.sen[k,1,i,t]*m.sen[k,1,j,t] for t in m.t) \
for k in m.e)/len(m.t)/len(m.e)
m.FIMmatrix = Constraint(m.p, m.p, rule=_FIMmatrix)
Любые советы приветствуются.Я подумал о стратегии, в которой локальный счетчик также должен быть равен mod (k, m) eq 1, где mod - это модуль, а m - количество точек коллокации, но не удалось найти оператор модуля в Pyomo.
Васко Луна