Об РСН.

Не так давно наткнулся на радиостанцию "Русская Служба Новостей" - 107.0 FM. (www.rusnovosti.ru) Радио в моём вкусе - музыки нет, только новости и всякие ток-шоу.
Есть, например, Михаил Леонтьев (программа "Диагноз"). Его всегда интересно послушать, но по отношению к звонящим в студию он ведёт себя совершенно по-хамски, если говорят что-то против Путина. Но это уж часть имиджа.
В четверг послушал программу "Вечерний приём" с Гончаровым. Довольно развязный дядя ("Мне нафиг не нужен такой чиновник, которого нужно совестить, чтоб он наконец начал работать"), но зато он так душевно общается со звонящими и темы у него интересные (например - жизнь инвалидов в России и "скука как двигатель всякого безкультурия и бездуховности"). Ещё один интересный момент - задавать вопросы и реплики можно через его ЖЖ (ahav10.livejournal.com).
Сегодня была программа Гастроном. Два довольно вальяжных человека (ведущая и гость) очень вкусно рассуждали о постной еде. Советовали есть чечевицу (оказывается, она содержит все минералы, бывает разных сортов, а её английское название происходит от слова пост) и хомус(арабское блюдо, паста из гороха и кунжутных семян с лимоном). Призывали поститься и думать о душе. Не ожидал такое услышать на радиостанции с явной заточкой на богатеньких.

TeX notes #1

В pdf можно включать файлы в формате jpg и png. А я-то мучался с векторизацией...

В TeX - e двойные русские кавычки ёлочкой (« ») делаются с помощью лигатур (если подключен babel): << и >>.

Ещё нашёл запись в блоге у jetxee про создание pdf на русском: http://sovety.blogspot.com/2007/07/pdf-latex.html

Этот удивительный Autotrace

Понадобилась мне картинка для презентации - иллюстрация теории центральных мест Кристаллера. Нашёл вордовский документ с чем-то похожим (правда зачем-то использовали какие-то стёбанные шрифты, типа шрифтов punk в TeX-е).

Задача первая - как вытянуть эту картинку из doc-a? C этим справился быстро, помогла заметка в Компьютерре (рубрика RTFM).
Надо просто сохранить в формате HTML. Все картинки сохраняются в gif-файлы (а не в png, как в заметке говорилось).
Теперь надо вставить картинку в tex-документ. Т.е. надо перевести в ps/eps, а лучше сразу в pdf.
Ну, думаю -чего проще. Открываем в XnView, сохраняем в ps. Ага, щаз. Получили чёрный прямоугольник.

Обращаемся к Интернету. И тут оказалось, что ps - векторный формат :0 ! И для перевода из битмапа в векторную картинку необходима автотрассировка. Это умеют делать Adobe Illustrator, Bitmap Tracer и многие другие проги - вот здесь
http://www.masternewmedia.org/ru/ digital_imaging/ autotracing_how_to_
convert_bitmaps_20050819.htm

длинный список таких прог. Я выбрал единственную бесплатную прогу из него - AutoTrace ( http://autotrace.sourceforge.net/ ). На вход она gif не берёт, конвертнул в bmp. Скачал в бинарниках (неохота компилять исходники). Оказалась это такая консольная утилита - один экзешник и с ним штук тридцать dll. Нда, думаю, оригинальная поставка.

Но ещё более оригинальным оказался формат командной строки.
Казалось бы чего естественнее
autotrace file.bmp file.ps
Нет - too many files. Если оставить autotrace file.gif - получаем кучу цифер, выданных в stdout (такого цирка ещё не видел я). Вызов с ключом -help прояснил ситуацию, однако при этом система выкинула ошибку о попытке записи куда-то не туда. Такую же ошибку выдаёт при попытке трассировки в ps (даже с правильными ключами). Сыровата прога, даа..

В итоге удалось конвертировать мой c_theory.bmp в pdf и eps (в цветном виде). Делалось это так:
autotrace -output-format eps -output-file c.eps c_theory.bmp
autotrace -output-format pdf -output-file c.pdf c_theory.bmp

Есть ещё надстройка над Autotrace с GUI на Java. Называется Delineate. Надо бы и её попробовать.

Расчёты

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

Делаем расчёт с растяжением (hx *= 2) и сжатием (hx *= 0.5).

В первом случае в центре образуется область с малой температурой - значение в центре растёт ОЧЧЕНЬ медленно по сравнению с максимумом решения:

time = 0 center = 1.225234 sol_max = 1.859148

time = 0.5417 center = 2.689374 sol_max = 14.994066

time = 0.5810 center = 3.058009 sol_max = 120.199303

time = 0.5836 center = 3.088310 sol_max = 485.715073

При сжатии функция очень быстро вырождается в первую радиально-симметричную сф.

Неискажённая функция долго держит свою форму (досчитал до sol_max = 6000), хотя значение в центре всё равно сильно отстаёт.

Полезные команды Matlab

Список файлов в папке: files = dir('dirname'). В files сохраняется список файлов в данной директории виде массива структур с полями: name, date, bytes (размер файла в байтах), isdirectiry(1, если name - папка).

Сохранение рисунка из figure в графический файл - print. Можно выбирать формат файла (bmp, eps, jpeg, много других) и регулировать разрешение.
C помощью разных опций эта команда также может печатать графики, сохранять их в clipboard и сохранять всё содержимое figure.

Число обусловленности матрицы (отношение максимального с.з. к минимальному) - cond. Некая оценка этого числа - rcond. Если матрица хорошо обусловлена, то rcond близко к 1, если плохо - к 0. Как сказано в хелпе "rcond более эффективна, но менее надёжна, чем cond".

Перед отъездом

Завтра еду в Ялту на NPDE2007. Поехали!! По приезде необходимо сделать программу расчёта тепловых структур нормальным консольным приложением - задавать файл с начальной функцией, время, до которого считать и файл для выдачи результата. Или без времени - CountTillBlowUp, каждую десятую по времени выдаём дамп.
А потом с помощью скрипта скармливать целые директории с нач функциями. Хотя в таком подходе есть здесь есть свои трудности...