Перейти из форума на сайт.

НовостиФайловые архивы
ПоискАктивные темыТоп лист
ПравилаКто в on-line?
Вход Забыли пароль? Первый раз на этом сайте? Регистрация
Компьютерный форум Ru.Board » Компьютеры » Прикладное программирование » Вопросы программирования на FORTRAN (ФОРТРАН)

Модерирует : ShIvADeSt

 Версия для печати • ПодписатьсяДобавить в закладки
Страницы: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329

Открыть новую тему     Написать ответ в эту тему

akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Обсуждаются все вопросы, связанные с программированием на ФОРТРАН, как общего так и конкретного характера.
Постарайтесь дать как можно больше информации о возникшей проблеме -- это в конце концов в ваших же интересах чтобы вам помогли...

прежде чем просить помощи в задании
платное решение задач

ресурсы этого топика
ссылка на подборку ресурсов, собранных посетителями этого форума
 
то, чем мы решили поделиться
ссылка на страничку программ etc собственного изготовления, которыми любезно делятся наши форумчане


если вам вдруг не отвечают или ответ вас не устраивает
и вообще полезно прочитать всем спрашивающим
 
просьба к пишущим и отвечающим все большие листинги оформлять тегом more
и отключать графические смайлики при размещении фортран-кода

Всего записей: 24056 | Зарегистр. 06-12-2002 | Отправлено: 18:11 14-01-2007 | Исправлено: akaGM, 09:47 01-03-2020
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Igorr
посмотри её определение через бесселЯ...

Всего записей: 24056 | Зарегистр. 06-12-2002 | Отправлено: 21:42 11-09-2009
Igorr

Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
akaGM
Ее определение я знаю - сходящийся ряд произведений степенной и бесселя. Самому проссумировать пока не дошли руки. Мне нужен либо готовый алгоритм через аппроксимацию с независящим числом членов ряда от значений аргументов, либо библиотека.

Всего записей: 2003 | Зарегистр. 01-05-2002 | Отправлено: 22:00 11-09-2009
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Igorr

Цитата:
Мне нужен либо готовый алгоритм через аппроксимацию

а через аппроксимацию самих бесселей?

Цитата:
произведений степенной и бесселя

скорее степенно-гиперболической, для двух-то переменных...
 

Цитата:
либо библиотека.

это вряд ли, самому придётся...

Всего записей: 24056 | Зарегистр. 06-12-2002 | Отправлено: 14:00 14-09-2009 | Исправлено: akaGM, 19:00 14-09-2009
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Igorr
посм сюда, может пригодится
 stacks.iop.org/Met/40/S5

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 16:50 14-09-2009
Igorr

Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
akaGM

Цитата:
а через аппроксимацию самих бесселей?

Я рассматривал разные варианты аппроксимаций. Аппроксимация бесселей усложняет проблему, т.к. они выражаются через сложный полином при |аргументе| > 3. Самый первый кандидат - это аппроксимация "exp(z2)", но тогда мне понадобятся функции Ломмеля smu,nu(z); это уже не критично, т.к. они выражаются через 1F2. Этот вариант я думаю использовать, если не найду Umu(z,p).
 
terminat0r

Цитата:
stacks.iop.org/Met/40/S5

Спасибо; эта статья часто попадается при поиске, но меня пересылают на www.iop.org, а там денег хотят, аж $30 - не позволительно для меня.

Всего записей: 2003 | Зарегистр. 01-05-2002 | Отправлено: 17:31 14-09-2009
Andrew10

Advanced Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Igorr
 
Статью можно попросить здесь:
http://forum.ru-board.com/topic.cgi?forum=93&topic=3234&start=80
 
 
Собственно выяснилось, что у меня самого есть доступ к этой статье, так что держи:
 
_http://rapidshare.de/files/48341234/me3102.pdf.html
 
 

Всего записей: 780 | Зарегистр. 26-02-2005 | Отправлено: 17:44 14-09-2009 | Исправлено: Andrew10, 17:52 14-09-2009
Igorr

Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Andrew10
Большое спасибо; уже скачал; сейчас смотрю.
 
Доп.
Предложенный метод (замены одной проблемы другой, причем, решаемой не во всей области значений аргументов) мне показался менее приемлемым, чем через аппроксимацию экспоненты. Еще раз спасибо за статью.
 
Проблема остается.

Всего записей: 2003 | Зарегистр. 01-05-2002 | Отправлено: 18:25 14-09-2009 | Исправлено: Igorr, 18:59 14-09-2009
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Igorr
http://ifile.it/xosm1qh
вот еще одна статтья, того же автора по сходной тематике

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 18:48 14-09-2009
Igorr

Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
terminat0r
Цитата:
еще одна статтья
Подход в этой статье аналогичен подходу в предыдущей; и, к сожалению, не решает моей проблемы. Спасибо за статью.

Всего записей: 2003 | Зарегистр. 01-05-2002 | Отправлено: 19:45 14-09-2009
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Igorr
Эти функции  -частный случай гипергеометрической? тогда
http://www.cpc.cs.qub.ac.uk/summaries/ACPA_v1_0.html
(проверено, работает, но для некоторых аргументов надо использовать формулы аналитического продолжения для pFq)
и скомбинировать в Lommel
--
Я бы на вашем месте искал в пакетах, где она может использоваться, например в Geant4 или что-то похожее

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 15:55 15-09-2009 | Исправлено: terminat0r, 15:58 15-09-2009
Igorr

Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
terminat0r

Цитата:
Эти функции  -частный случай гипергеометрической?
Нет, конечно - тогда бы и вопроса не было. По крайней мере я такого представления - Umu(z,y) = f(z,y,pFq(a...;b...;f0(z,y))) - не встречал.
 

Цитата:
Я бы на вашем месте искал в пакетах, где она может использоваться, например в Geant4 или что-то похожее
Я и ищу, только не знаю таких пакетов. В моей памяти есть смутное ощущение, что какие-то функции Ломмеля были в фортрановской библиотеке БЭСМ-6. В теме "БЭСМ-6, Эльбрус, ЕС" http://forum.ru-board.com/topic.cgi?forum=35&bm=1&topic=19594#1 я спросил про исходники библиотек БЭСМ-6 и CERN, но пока тихо. А в Geant4 этой функции нет.

Всего записей: 2003 | Зарегистр. 01-05-2002 | Отправлено: 17:41 15-09-2009
terminat0r



Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
насколько я понимаю например с этой статьи, проблем никаких нет, можно просто считать в лоб,
там заданы даже условия когда можно обрывать ряд
http://ifile.it/ksw10dr

Всего записей: 2084 | Зарегистр. 31-03-2002 | Отправлено: 00:28 16-09-2009
Igorr

Silver Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
terminat0r
Цитата:
насколько я понимаю например с этой статьи, проблем никаких нет, можно просто считать в лоб,  
там заданы даже условия когда можно обрывать ряд  
Спасибо за статью. Повторюсь: вопрос не в количестве членов ряда, определяющего функцию (их число оценивается элементарно, зная поведение бесселя в зависимости от индекса), а в том, что нужен
Цитата:
Igorr
Цитата:
готовый алгоритм через аппроксимацию с независящим числом членов ряда от значений аргументов
Что-то типа как sin(x)/x=1+a1x2+...+a5x10, или через многочлены Чебышева (как у Ю.Люка "Спец. мат. функции и их аппроксимации"), или что-то аналогичное.

Всего записей: 2003 | Зарегистр. 01-05-2002 | Отправлено: 01:22 16-09-2009
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
кто-нибудь занимался оптимизацией с использованием генетических алгоритмов?

Всего записей: 24056 | Зарегистр. 06-12-2002 | Отправлено: 11:58 18-09-2009
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
тогда такие два вопроса:
 
1. кто-нибудь глобальной оптимизацией "увлекался"?
 
2. кто-нибудь сталкивался с целевой функцией, кот. является узкощелевой (типа
многомерных дельта-функций), на которых отказываются работать стандартные
МНК'ашные методы?

Всего записей: 24056 | Зарегистр. 06-12-2002 | Отправлено: 17:46 27-09-2009
sergeymd21

Newbie
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
дана к примеру матрица:
 
a  b  c  d  e  f
1, 0, 0, 1, 0, 1, a
0, 1, 1, 0, 0, 1, b
0, 1, 1, 0, 0, 0, c
1, 0, 0, 1, 0, 1, d
0, 0, 0, 0, 1, 0, e
1, 1, 0, 1, 0, 1  f
 
Тут нужно найти начальное покрытие ну то есть оно ищется так: построчно просматривается матрица и где есть единица составляется пара из строки и столбца
 
a,adf - 1-ая строка и 1-ый и 4-ый и 6-ый столбцы
 
Вот оно найдено: C0={(a,adf),(b,bcf),(c,bc),(d,adf),(e,e),(f,abdf)}
 
ВНИМАНИЕ ВОПРОС:как бы мне это оформит в c++ , какой тип данных ,ну то есть я могу сделать массив arr[5][5] и туда это все помещать,но там будут оставаться не заполненные области т.к. эл-ов не равное кол-во.  
 
Мне просто потом нужно работать с этими группами искать их пресечение, объединение и т.д  
             
 
Добавлено:
изменяюсь я ошибся форумом это в c++ ,удалите!!

Всего записей: 3 | Зарегистр. 10-06-2009 | Отправлено: 19:29 28-09-2009
karakurt2



Advanced Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
akaGM, я постоянно "увлекаюсь" глобальной оптимизацией.  
Но, к сожалению, реализую алгоритмы самостоятельно, без использования функций математических библиотек. Например, реализовывал simulated annealing для решения задачи об оптимальном расположении меток на карте (np-полная задача). Генетические алгоритмы изучал в студенчестве, но применять в реальной практике не приходилось.
 
PS: на алгоритмические вопросы Вам смогут ответить быстрее и квалифицированнее на dxdy.ru

Всего записей: 733 | Зарегистр. 06-12-2003 | Отправлено: 20:56 05-10-2009 | Исправлено: karakurt2, 20:57 05-10-2009
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
karakurt2

Цитата:
Но, к сожалению, реализую алгоритмы самостоятельно

и о чём же здесь сожалеть?
 
можно примеры, хотя бы на пальцАх? и каковы критерии, что достигаемые экстремумы -- глобальные?
 
меня интересует многомерная параметрическая оптимизация, при этом функционал -- сильно осциллирующая функция с большим количествои локальных экстремумов
 

Цитата:
Например, реализовывал simulated annealing

какое совпадение...
только что освоил и прикрутил себе SA
сначала смотрел ASA, но счёл её слишком переусложнённой, потом набрёл на что-то более простое...

Всего записей: 24056 | Зарегистр. 06-12-2002 | Отправлено: 12:58 06-10-2009 | Исправлено: akaGM, 13:09 06-10-2009
karakurt2



Advanced Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
akaGM
 
Вот несколько статей, которые у меня остались с тех пор. Кроме того, этот метод описывается в Numerical Recipes. Ещё я встречал его применение в геостатистике для построения корреляции каротажных кривых скважин.
 
Насколько я понимаю, достижение глобального минимума функционала не гарантируется. Если это можно было бы легко проверить, терялся бы смысл самого существования оптимизационных алгоритмов.
 
PS: спасибо за ссылку на интересный сайт.

Всего записей: 733 | Зарегистр. 06-12-2003 | Отправлено: 18:23 06-10-2009
akaGM

Platinum Member
Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
karakurt2

Цитата:
Насколько я понимаю, достижение глобального минимума функционала не гарантируется. Если это можно было бы легко проверить...

ну так как _ты_ это проверяешь?
в свете
Цитата:
реализую алгоритмы самостоятельно

 
я вот, например, исследую сечения функционала, плюс одно время ковырялся с развёртками пространства решений на одномерное...

Цитата:
Кроме того, этот метод описывается в Numerical Recipes

да, я это знаю, да и не вызвал SA особых проблем...
 
не за что...
тебе спасибо за
Цитата:
несколько статей

Всего записей: 24056 | Зарегистр. 06-12-2002 | Отправлено: 18:37 06-10-2009 | Исправлено: akaGM, 19:34 06-10-2009
Открыть новую тему     Написать ответ в эту тему

Страницы: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329

Компьютерный форум Ru.Board » Компьютеры » Прикладное программирование » Вопросы программирования на FORTRAN (ФОРТРАН)


Реклама на форуме Ru.Board.

Powered by Ikonboard "v2.1.7b" © 2000 Ikonboard.com
Modified by Ru.B0ard
© Ru.B0ard 2000-2024

BitCoin: 1NGG1chHtUvrtEqjeerQCKDMUi6S6CG4iC

Рейтинг.ru