Почему результат вычислений не зависит от параметра

@valery86 · Регистрация: 08.08.2017

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

Здравствуйте!

Скажите,пожалуйста,почему в этом коде при всех t k получается одинаковое

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
clc;
clear;
t=0;
y=(-10:0.01:10)';
for i=1:length(t)
for j=1:length(y)
    f= @(z,t) (exp(-(z.^2)))./(1+(y(j)+2*sqrt(pi)*z.^2));
    Q(i)=integral(@(z)f(z,t),-inf,inf);
    k(i)=(1/(sqrt(pi)))*Q(i)
end
end
plot(t,k)

Можно ли как-то исправить чтобы с увеличением t (0,1,..,5) уменьшалось к?

@SSC · 14.08.2017, 15:34

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

почему в этом коде при всех t k получается одинаковое

Потому что подинтегральная функция

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

f= @(z,t) (exp(-(z.^2)))./(1+(y(j)+2*sqrt(pi)*z.^2));

от t независит.

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

Можно ли как-то исправить

Измените подинтегральную функцию, для начала проверьте правильно ли Вы ее записали.

@valery86 · 14.08.2017, 16:00 **[ТС]**

SSC,спасибо большое за ответ.

Скажите,пожалуйста,как можно изменить подинт.функцию f= @(z) (exp(-(z.^2)))./(1+(y(j)+2*sqrt(pi)*z).^2) чтобы к убывало с возрастанием т?

@SSC · 14.08.2017, 16:16

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

Matlab M
1
f= @(z,t) (exp(-(z.^2)))./(t(i)+(y(j)+2*sqrt(pi)*z.^2));

PS. Однако при выполнении программы куча предупреждений по недостаточной точности вычисления интеграла

@valery86 · 14.08.2017, 16:48 **[ТС]**

SSC,спасибо,извините еще вопрос,я же задаю много точек y=(-10:0.01:10),т.е 2000 почему программа не выдает эту матрицу,а только 5 значений к?Как можно посмотреть эту матрицу по точкам?

@faraonCar · 14.08.2017, 16:52

Все очень просто
во первых

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

Потому что подынтегральная функция
Сообщение от valery86
f= @(z,t) (exp(-(z.^2)))./(1+(y(j)+2*sqrt(pi)*z.^2));
от t не зависит.

во вторых

Matlab M
1
k(i)=1/(sqrt(pi)*Q(i))

И будет у Вас

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

к убывало с возрастанием т

@valery86 · 14.08.2017, 19:36 **[ТС]**

Действительно,к убывает..Но матрица для каждого т почему-то заполняется одинаковыми значениями т(0)=1,3847е+0,8 т(1)=0.0128..

Есть ли ошибка в коде?Мне просто важно посмотреть на элемент на 1001 шаге,что дает некое понятие о правильности алгоритма.

@faraonCar · 14.08.2017, 19:38

Где вы взяли t тут???

Matlab M
1
(exp(-(z.^2)))./(1+(y(j)+2*sqrt(pi)*z.^2))

@valery86 · 14.08.2017, 19:42 **[ТС]**

Ооо,сорри,вместо корня из Пи т должно быть

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
clc
clear
t=1;
y=(-10:0.01:10)';
for i=1:length(t)
for j=1:length(y)
f=  @(z,t) (exp(-(z.^2)))./(1+(y(j)+2*sqrt(t)*z).^2);
Q(i,j)=integral(@(z)f(z,t),-inf,inf);
k(i,j)=(1/(sqrt(pi)))*Q(i)
end
end
plot(t,k)
data=[t',k]

@faraonCar · 14.08.2017, 19:56

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

только 5 значений к?Как можно посмотреть эту матрицу по точкам?

У Вас размер матрицы k 6 на 2001.

Добавлено через 1 минуту
Проверяйте формулу если что-то не получается.

Новые блоги и статьи Все статьи Все блоги /
Символьное дифференцирование igorrr37 13.02.2026 / * Логарифм записывается как: (x-2)log(x^2+2) - означает логарифм (x^2+2) по основанию (x-2). Унарный минус обозначается как ! в-строка - входное арифметическое выражение в инфиксной(обычной). . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»
SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .

@valery86 0 / 0 / 0 Регистрация: 08.08.2017 Сообщений: 19
	14.08.2017, 16:00 [ТС]
	SSC,спасибо большое за ответ. Скажите,пожалуйста,как можно изменить подинт.функцию f= @(z) (exp(-(z.^2)))./(1+(y(j)+2sqrt(pi)z).^2) чтобы к убывало с возрастанием т? 0

@valery86 0 / 0 / 0 Регистрация: 08.08.2017 Сообщений: 19
	14.08.2017, 16:48 [ТС]
	SSC,спасибо,извините еще вопрос,я же задаю много точек y=(-10:0.01:10),т.е 2000 почему программа не выдает эту матрицу,а только 5 значений к?Как можно посмотреть эту матрицу по точкам? 0

@valery86 0 / 0 / 0 Регистрация: 08.08.2017 Сообщений: 19
	14.08.2017, 19:36 [ТС]
	Действительно,к убывает..Но матрица для каждого т почему-то заполняется одинаковыми значениями т(0)=1,3847е+0,8 т(1)=0.0128.. Есть ли ошибка в коде?Мне просто важно посмотреть на элемент на 1001 шаге,что дает некое понятие о правильности алгоритма. 0