Телескопы покупают здесь


A A A A Автор Тема: Непрерывная функция относительной яркости  (Прочитано 536 раз)

0 Пользователей и 1 Гость просматривают эту тему.

Оффлайн gree3xАвтор темы

  • Новичок
  • *
  • Сообщений: 42
  • Благодарностей: 0
  • 10^10^122
    • Сообщения от gree3x
Здравствуйте, дорогие форумчане! С наступившим НГ и наступающим Рождеством вас! Снова обращаюсь к вам за помощью. Возможно кто-то из вас хорошо знает матанализ и может подсказать на счёт функции. Проблема в следующем: относительная яркость наблюдаемых объектов равна (D/d)^2. Где D- диаметр объектива, а d-диаметр выходного значка, что эквивалентно квадрату от увеличения. Для наблюдения дипская (особенно тусклых объектов неба) стараются использовать т.н. равнозрачковые увеличения, при которых d=6мм. Т.е. световой поток будет полностью поглощен глазом по всей площади зрачка глаза. При больших увеличениях d<6мм яркость пятна будет уменьшаться вследствие сокращения поля зрения, а при меньших d>6мм аппертура будет диафрагмироваться и световой поток также будет падать. Вопрос в следующем: существует ли такая непрерывная функция яркости, которая достигает максимального значения в точке, где d=6мм и равна D/36, а в остальных значениях убывает экспоненциально (?) до нуля в обе стороны? Ломаю голову давно уже, может есть идеи как найти такую функцию? Заранее спасибо))
СТФ_Мираж7,Celestron_SkyMaster_25x100

Оффлайн Toth

  • *****
  • Сообщений: 2 603
  • Благодарностей: 175
    • Сообщения от Toth
существует ли такая непрерывная функция яркости, которая достигает максимального значения в точке, где d=6мм и равна D/36, а в остальных значениях убывает экспоненциально (?) до нуля в обе стороны?
Наверное, не только непрерывная, но и гладкая. Такая наверное не годится
y=exp(x-6) при x<=6
y=exp(6-x) при x>6   ?
А иначе - чтоб прямо экспоненциально, не получится.
А если без точной exp-зависимости, то таких функций много, хотя бы https://ru.wikipedia.org/wiki/Гауссова_функция при μ = 6


Оффлайн gree3xАвтор темы

  • Новичок
  • *
  • Сообщений: 42
  • Благодарностей: 0
  • 10^10^122
    • Сообщения от gree3x
А если без точной exp-зависимости, то таких функций много, хотя бы https://ru.wikipedia.org/wiki/Гауссова_функция при μ = 6
Использовал функцию распределения Гаусса. Импирически привёл значение основания для уровня 0.5 для случая d>6 имеет значение (72)^1/2. Функция имеет вид 1*1.12^(-(x-6)^2).
Это справедливо если функция симметрична, но в левой части на уровне 0.5 она равна (18)^1/2.
Нужна несимметричная функция.
СТФ_Мираж7,Celestron_SkyMaster_25x100

Оффлайн Toth

  • *****
  • Сообщений: 2 603
  • Благодарностей: 175
    • Сообщения от Toth
Ну, тогда типа y = x^a*exp(-bx),
коэфф. a,b - подобрать.
Или просто нарисовать нужный график, а потом по точкам подобрать нужный полином с помощью Пафнутия Львовича.

Оффлайн gree3xАвтор темы

  • Новичок
  • *
  • Сообщений: 42
  • Благодарностей: 0
  • 10^10^122
    • Сообщения от gree3x
Пафнутия Львовича
Можно по подробнее?
СТФ_Мираж7,Celestron_SkyMaster_25x100

Оффлайн Toth

  • *****
  • Сообщений: 2 603
  • Благодарностей: 175
    • Сообщения от Toth
Я про это - https://ru.wikipedia.org/wiki/Многочлены_Чебышёва

Но там довольно нудный алгоритм.
Я как-то писал прогу для расчета. Вроде работает.
Можете просто выдать значения X, Y , я подберу формулу многочлена.

Оффлайн gree3xАвтор темы

  • Новичок
  • *
  • Сообщений: 42
  • Благодарностей: 0
  • 10^10^122
    • Сообщения от gree3x
Можете просто выдать значения X, Y , я подберу формулу многочлена

Левая часть до максимума описывается уравнением (d/dr)^2 где d-диаметр выходного зрачка, а dr=6 (равнозрачковый диаметр в мм) и правая часть функции после максимума убывает и имеет вид (dr/d)^2  значения X,Y я рассчитывал по ним.
СТФ_Мираж7,Celestron_SkyMaster_25x100

Оффлайн Toth

  • *****
  • Сообщений: 2 603
  • Благодарностей: 175
    • Сообщения от Toth
Так это получается не гладкая функция. У ее производной - перепад в 6.
Тут полиномы не помогут. То есть можно, но получится 100500 порядка полином, чтоб хорошо подходил.
А почему нельзя кусочно-гладкую, то есть если x< 6 .. если x>6 ..  ?

Оффлайн gree3xАвтор темы

  • Новичок
  • *
  • Сообщений: 42
  • Благодарностей: 0
  • 10^10^122
    • Сообщения от gree3x
А почему нельзя кусочно-гладкую, то есть если x< 6 .. если x>6 ..  ?
Задать функцию интервально не проблема, но сама задача состоит в том как получить полное уравнение функции для всех значений и не получить бесконечность..
СТФ_Мираж7,Celestron_SkyMaster_25x100

Оффлайн Toth

  • *****
  • Сообщений: 2 603
  • Благодарностей: 175
    • Сообщения от Toth
Может как то так подобрать -

Оффлайн gree3xАвтор темы

  • Новичок
  • *
  • Сообщений: 42
  • Благодарностей: 0
  • 10^10^122
    • Сообщения от gree3x
Может как то так подобрать -
У меня по этой формуле строит нечто иное
СТФ_Мираж7,Celestron_SkyMaster_25x100

Оффлайн Toth

  • *****
  • Сообщений: 2 603
  • Благодарностей: 175
    • Сообщения от Toth
Странно.
Вот эта 0.222*x^2*EXP(-0.35*x)

Оффлайн gree3xАвтор темы

  • Новичок
  • *
  • Сообщений: 42
  • Благодарностей: 0
  • 10^10^122
    • Сообщения от gree3x
Вот эта 0.222*x^2*EXP(-0.35*x)
В каком приложении строили? У меня Grapher Free
СТФ_Мираж7,Celestron_SkyMaster_25x100

Оффлайн Toth

  • *****
  • Сообщений: 2 603
  • Благодарностей: 175
    • Сообщения от Toth
меня Grapher Free
У меня Excel.
И в Паскаль то же самое получается.
y:=0.222*x*x*Exp(-0.35*x);

Оффлайн gree3xАвтор темы

  • Новичок
  • *
  • Сообщений: 42
  • Благодарностей: 0
  • 10^10^122
    • Сообщения от gree3x
У меня Excel
Хорошо. Спасибо. Видимо проблема в андройдовских программах. Доберусь до компа проверю.
СТФ_Мираж7,Celestron_SkyMaster_25x100