2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему
 
 matlab, интеграл с параметром, график
Сообщение13.02.2006, 18:18 


13/02/06
4
вопрос по matlab
подскажите, пожалуйста, с помощью каких функций можно нарисовать график такой вот штуки
f(t)=\int_0^t \eta(t-p) exp (-x^2/4p) dp
интеграл не берется явно
спасибо

 Профиль  
                  
 
 
Сообщение13.02.2006, 18:35 
Экс-модератор
Аватара пользователя


23/12/05
12046
x - в данном случае константа?

 Профиль  
                  
 
 
Сообщение13.02.2006, 18:40 


13/02/06
4
здесь можно считать константой
итогом будет график функции двух переменных (t,x).

 Профиль  
                  
 
 
Сообщение13.02.2006, 18:41 
Экс-модератор
Аватара пользователя


23/12/05
12046
Дайте требуемый диапазон t и х, и я построю Вам этот график

 Профиль  
                  
 
 
Сообщение13.02.2006, 18:51 


13/02/06
4
x = 0:20
t = 0:10
интересен не сколько график, а как ;)
спасибо

 Профиль  
                  
 
 Re:
Сообщение14.02.2006, 10:31 
Экс-модератор
Аватара пользователя


23/12/05
12046
Вечерком нарисовал я эту функцию... но - хотел спать, поэтому:
а) Не перепроверял - доверяй, но проверяй :wink:
б) Не писал комментарии: если что-то будет непонятно - добро пожаловать.
Выбросил я эту $\eta$ - потом домножить на константу можно всегда, и получил такую функцию для вычисления интеграла:
Код:
function integ_par=integ_par(x,t,number_of_points)
fxt=zeros(1,number_of_points);
integ_par=zeros(length(t),length(x));
for ii=1:length(t)
    p=0:t(ii)/number_of_points:t(ii);
    for jj=1:length(x)
        if ~isempty(p)   
            fxt=(t(ii)-p(2:end)).*exp(-0.25*x(jj)^2./p(2:end));
            integ_par(ii,jj)=0.5*sum((fxt(1:end-1)+fxt(2:end))).*t(ii)/number_of_points;
        end
    end
end

В данном случае вычисление ведется методом трапеций.

Вызов функции и построение графика осуществляется например так:
Код:
x=0:.1:20;
t=0:.1:10;
integ=integ_par(x,t,50);
[tt,xx]=meshgrid(t,x);
mesh(t,x,integ)

В результате имеем:
Изображение

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 6 ] 

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group