2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Плохо обусловленная система
Сообщение13.10.2016, 21:25 
Какие могут быть причины, что в плохо обусловленной СЛАУ погрешность в правой части не даёт больших погрешностей решения?
Например, есть у меня матрица с большим числом обусловленности (порядка $10^4$). Но, как бы я не меняла правую часть, относительная погрешность решения возрастает не более, чем в $10$ раз по сравнению с относительной погрешностью правой части. Почему?

 
 
 
 Re: Плохо обусловленная система
Сообщение13.10.2016, 21:36 
Аватара пользователя
А как Вы меняете правую часть?

 
 
 
 Re: Плохо обусловленная система
Сообщение13.10.2016, 21:49 
Евгений Машеров в сообщении #1159577 писал(а):
А как Вы меняете правую часть?

Сначала изменяла каждую координату на несколько %. Потом уже как-то произвольно меняла, но рост не больше, чем в 10 раз.

 
 
 
 Re: Плохо обусловленная система
Сообщение14.10.2016, 07:26 
Аватара пользователя
Ну, вот у нас задача $Ax=b$
Представляя матрицу A через сингулярное разложение $A=S\Lambda C$, где S и С ортогональные матрицы, а $\Lambda$ диагональная. Число обусловленности равно отношению максимального и минимального сингулярных чисел.
Решение есть $x=C^T\Lambda^{-1}S^Tb$
Разлагая возмущение вектора b по базису, образованному столбцами матрицы S, видим, что основное влияние на возмущение решения x оказывают компоненты, направленные вдоль столбцов, соответствующих малым сингулярным числам и изменение обратно пропорционально величине сингулярных чисел.
Далее могу лишь предположить, что среди сингулярных чисел матрицы одно очень мало, а остальные существенно больше его, что даёт высокое число обусловленности, однако если возмущение почти ортогонально столбцу, соответствующему этому сингулярному числу, его действие на решение определяется прочими сингулярными числами, не малыми. В Ваших экспериментах чисто случайно могла быть "почти ортогональность" этого рода. То есть подобрать малое возмущение b, дающее сильное возмущение x, можно, но подбирая специально, соответственно структуре матрицы А.

 
 
 [ Сообщений: 4 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group