Фортран для профессионалов
Излагаются средства математической библиотеки IMSL, входящей в состав профессиональных версий Фортрана фирм Microsoft и Compaq, позволяющие анализировать линейные модели и выполнять матричные вычисления. Приводятся алгоритмы, положенные в основу процедур библиотеки, а также варианты их программных реализаций. Поясняются все используемые при описании процедур библиотеки понятия. Излагаемый материал иллюстрируется большим числом примеров.Предназначено для научных работников, инженеров, преподавателей, студентов и аспирантов вузов.
часть 2
Содержание
Содержание книги "Фортран для профессионалов : математическая библиотека IMSL"
Отрывок из книги
О. В. Бартеньев. Фортран. Математическая библиотека IMSL. (Часть 2) 40 nz - размер векторов z и zhat. На входе, когда ipad = 0, nz ≥ nx + ny - 1; когда ipad = 1, nz ≥ 2α3β5γ, где 2α3β5γ - наименьшее число, большее или рав-ное nx + ny - 1; α, β и γ - неотрицательные целые числа. На выходе nz со-держит значение, использованное RCONV. z - вектор, содержащий свертку векторов x и y. zhat - вектор, содержащий дискретное преобразование Фурье вектора z. Комментарий. Если размер вектора z меньше необходимой величины, то возникнет завершающая ошибка, сообщающая приемлемый размер nz. Например: *** FATAL ERROR 1 from RCONV. The length of the vector Z must be at *** least 10 while NZ = 8 is given, where NZ is based on NX = 8 and *** NY = 3. Stop - Program terminated. Описание: Пусть nx - размер x, а ny - размер y. Если вычисляется циклическая свертка (ipad = 0), то nz ≥ nx + ny - 1 и более короткий вектор расширяется до размера более длинного за счет добавления завершающих нулей. Затем вычисляются элементы вектора z: ∑=+−=znjjjiiyxz11, i = 1, 2, …, nz, где индекс вектора x - это MODULO(j, nz). Метод вычисления zi основан на том, что дискретное преобразование Фурье циклической свертки двух наборов чисел (векторов) равно произве-дению дискретных преобразований Фурье этих векторов: nnnyxz=, где ∑=−−π−=zznmnnmimnezz1/)1)(1(2. После вычисления произведения дискретных преобразований Фурье векторов x и y выполняется обратное преобразование Фурье этого произве-дения и таким образом находится искомая свертка. Важно, чтобы величина nz равнялась произведению небольших простых чисел. Тогда вычислитель-ная сложность будет пропорциональна nzlognz. Заметим, что в случае вещественных векторов выполняются не ком-плексные, а вещественные преобразования Фурье, что приводит к сущест-венному сокращению вычислительных затрат и используемой памяти.
Внимание!
При обнаружении неточностей или ошибок в описании книги "Фортран для профессионалов : математическая библиотека IMSL (автор О. Бартеньев)", просим Вас отправить сообщение на почту help@directmedia.ru. Благодарим!
и мы свяжемся с вами в течение 15 минут
за оставленную заявку