Функция Ханкеля

@AlexChel · Регистрация: 09.11.2016

Студворк — интернет-сервис помощи студентам

Помогите, пожалуйста! Имеется интеграл Макдональда. В него входит функция Ханкеля первого порядка (она же функция Бесселя третьего рода). Как правильно ее записать в Matlab в скрипте? В моем случае интегрирование происходит в диапазоне от минус бесконечности до нуля. Кроме того есть подозрение, что неправильно задаю вектора.

C++
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
clc
n=1
sigma=10.^(-12)
pi=3.14
l=480
len=1500-l
hight1=4
hight2=200
fi0=180-atan(hight1/l)
fi=atan(hight2/len)
ro0=sqrt(hight1.^2+l.^2)
ro=sqrt(hight2.^2+len.^2)
omega=332*10.^6
z=0
z0=0
c=3*10.^8
beta1=pi-fi-fi0
beta2=3*pi-fi-fi0
k=sqrt((omega-4*pi*sigma*i)*omega/c)
alpha1=pi-beta1
alpha2=pi-beta2
R1=sqrt(ro.^2+ro0.^2+(z-z0).^2-2*ro*ro0*cos(alpha1))
R2=sqrt(ro.^2+ro0.^2+(z-z0).^2-2*ro*ro0*cos(alpha2))
x1=2*sqrt(ro*ro0)*sin(beta1/2)/R1
x2=2*sqrt(ro*ro0)*sin(beta2/2)/R2
y1=k*R1
y2=k*R2
M11=int(besselh(n, 1, y1*cosh(E)), E, -inf, x1)
M12=int(besselh(n, 1, y2*cosh(E)), E, -inf, x2)
M1=i*k*M11/2
M2=i*k*M12/2
e1=[0 0 k]
P=[0 0 k]
PE=P*M1-e1*M2

@Nick07 · 09.01.2017, 20:29

Ф-ции Бесселя - встроенные ф-ции матлаба, я с ними работал. Надо поискать встроенные ф-ии Ханкеля или их записать через бесселевские.

@AlexChel · 09.01.2017, 20:56 **[ТС]**

Я находил встроенную и знаю что она есть, но когда пробую прописать, она ошибку выдает. Проблема в том, что в функции у меня есть переменная Е, которая принимает значение только после интегрирования, а Matlab просит, чтобы Е имела значение. Вот что пишет:
Undefined function or variable 'E'.

Error in functionp1 (line 28)
M11=int(besselh(n, 1, y1*cosh(E)), E, -inf, x1)

@Nick07 · 10.01.2017, 08:44

Сообщение от AlexChel

Matlab просит, чтобы Е имела значение

Тогда это вопрос алгоритма, а не Matlab

Я проверял соответствие расчетов ф-ций Бесселя, реализованные в Matlab, фундаментальным таблицам, ошибок не обнаружил.

@AlexChel · 13.01.2017, 20:07 **[ТС]**

Я понимаю, что считает правильно. Меня интересует, как исправить скрипт, чтобы не возникало ошибки. Попробовал уже много вариантов и все равно проблема.
Нашел решение, кому надо - помогу

@Том Ардер · 04.03.2017, 13:36

Сообщение от AlexChel

Нашел решение, кому надо, пишите в личных сообщениях

Комментарий модератора

Правила форума

4.6. Обсуждение вопросов - только в теме на форуме. Приглашения к обсуждению еще где-либо (в том числе и с помощью системы личных сообщений) запрещены, за исключением коммерческих разделов.
4.10. Если вопрос был решен вами самостоятельно, отпишите об этом в своей теме - есть и другие люди, которые столкнутся с той же проблемой, и им поможет ваш ответ.

@Nick07 · 04.03.2017, 20:14

Сообщение от Том Ардер

кому надо, пишите в личных сообщениях

Вы писали свой вопрос не в личные сообщения, а обратились ко всему сообществу. Поэтому логично (да и правилам форума) сообщить всем свое решение.

Новые блоги и статьи Все статьи Все блоги /
Модель микоризы: классовый агентный подход 3 anaschu 06.01.2026 aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114	Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ФедосеевПавел 06.01.2026 Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .	Модель микоризы: классовый агентный подход 2 anaschu 06.01.2026 репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .	Расчёт токов в цепи постоянного тока igorrr37 05.01.2026 / * Дана цепь постоянного тока с сопротивлениями и напряжениями. Надо найти токи в ветвях. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа и решает её. Последовательность действий:. . .
Новый CodeBlocs. Версия 25.03 palva 04.01.2026 Оказывается, недавно вышла новая версия CodeBlocks за номером 25. 03. Когда-то давно я возился с только что вышедшей тогда версией 20. 03. С тех пор я давно снёс всё с компьютера и забыл. Теперь. . .	Модель микоризы: классовый агентный подход anaschu 02.01.2026 Раньше это было два гриба и бактерия. Теперь три гриба, растение. И на уровне агентов добавится между грибами или бактериями взаимодействий. До того я пробовал подход через многомерные массивы,. . .	Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Programma_Boinc 28.12.2025 Советы по крайней бережливости. Внимание, это ОЧЕНЬ длинный пост. Налог на собак: https:/ / ********/ gallery/ V06K53e Финансовый отчет в Excel: https:/ / ********/ gallery/ bKBkQFf Пост отсюда. . .	Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Нашел на реддите интересную статью под названием Anyone know where to get a free Desktop or Laptop? Ниже её машинный перевод. После долгих разбирательств я наконец-то вернула себе. . .

@Nick07 533 / 438 / 47 Регистрация: 17.07.2013 Сообщений: 2,236
	09.01.2017, 20:29
	Ф-ции Бесселя - встроенные ф-ции матлаба, я с ними работал. Надо поискать встроенные ф-ии Ханкеля или их записать через бесселевские. 0

@AlexChel 2 / 2 / 0 Регистрация: 09.11.2016 Сообщений: 11
	09.01.2017, 20:56 [ТС]
	Я находил встроенную и знаю что она есть, но когда пробую прописать, она ошибку выдает. Проблема в том, что в функции у меня есть переменная Е, которая принимает значение только после интегрирования, а Matlab просит, чтобы Е имела значение. Вот что пишет: Undefined function or variable 'E'. Error in functionp1 (line 28) M11=int(besselh(n, 1, y1*cosh(E)), E, -inf, x1) 0

@AlexChel 2 / 2 / 0 Регистрация: 09.11.2016 Сообщений: 11
	13.01.2017, 20:07 [ТС]
	Я понимаю, что считает правильно. Меня интересует, как исправить скрипт, чтобы не возникало ошибки. Попробовал уже много вариантов и все равно проблема. Нашел решение, кому надо - помогу 0