complexify
Замените ureal атомы по суммам ureal и ucomplex (или ultidyn) атомы
Синтаксис
MC = complexify(M,alpha) MC = complexify(M,alpha,'ultidyn')
Описание
Область команды complexify заменяет ureal атомы с суммами ureal и ucomplex атомы с использованием usubs. Опционально сумма может состоять из ureal и ultidyn атом.
complexify используется для улучшения обусловленности вычислений устойчивой устойчивости (robstab) для ситуаций, когда преимущественно бывают ureal неопределенные элементы.
MC = complexify(M,alpha) результаты в каждом ureal атом в MC имеющий ту же Name и NominalValue как соответствующий ureal атом в М. Если Range - область значений значений одной ureal атом от M, затем область значений соответствующего атома мочевины в MC является
[Range(1)+alpha*diff(Range)/2 Range(2)-alpha*diff(Range)/2]
Чистый эффект заключается в том, что та же реальная область значений покрыта реальной и сложной неопределенностью. Действительная область значений параметра уменьшается на равные величины на каждом конце, и alpha представляет (в относительном смысле) уменьшение общей области значений. The ucomplex атом добавит это сокращение области значений назад в MC, но как мяч с реальными и мнимыми частями.
The ucomplex атом имеет NominalValue 0 и Radius равно alpha*diff(Range). Его имя - имя оригинала ureal атом, добавленный символами '_cmpxfy'.
MC = complexify(M,alpha,'ultidyn') то же самое, кроме ограниченного усилением ultidyn вместо ucomplex используются атомы атомы. The ultidyn атом имеет свою Bound равно alpha*diff(Range).
Примеры
Комплексный неопределенный параметр
Чтобы проиллюстрировать комплексификацию, создайте неопределенный действительный параметр, приведите его к неопределенной матрице и примените 10% комплексификацию.
a = umat(ureal('a',2.25,'Range',[1.5 3])); b = complexify(a,.1); as = usample(a,200); bs = usample(b,4000);
Сделайте график поля точек значений, которые может взять комплексная матрица (скаляр), а также значений исходного неопределенного действительного параметра.
plot(real(bs(:)),imag(bs(:)),'.',real(as(:)),imag(as(:)),'r.') axis([1 3.5 -0.2 0.2])
