Решение уравнения теплопроводности в MatLab

[Kafari]

Спасибо огромное, теперь у меня, кажется, получилось то что нужно! По крайней мере выглядит очень похоже на правду
Делаю так вот:

function [  ] = Temp( T1,T2,T0,A,L )

s=figure;
step=0.1;
t=0;
h=1;
Ti=T0*ones(L,1);
Tk=T0*ones(L,1);
colormap(bone);

while ishandle(s)
    Ti(1)=T1;
    Ti(L)=T2;
    Tk(1)=T1;
    Tk(L)=T2;
    for x=2:h:(L-1)
        Tk(x)=(A*step/h^2)*(Ti(x-1)+Ti(x+1)-2*Ti(x)) + Ti(x);
    end
    image(Ti);
    Ti=Tk;
    t=t+step;
    pause(step);
end

end

Теперь задали реализовать это же на С++... Но, думаю, теперь уже справлюсь))

[panika_]

Kafari
Ты написала программку?

[Kafari]

Да. Она двумя сообщениями выше

[panika_]

какая я внимательная....
а пояснить можно?)... (и что какими буквами обозначено)

[Kafari]

Где-то я уже это писала, но в общем так:
step шаг по времени, t - время, h - шаг по координате, Ti и Tk - два массива температур от координаты в моменты времени i и i+1.

Но я еще раз говорю, моя программа не рисует график. Она просто визуализирует процесс температуропроводности. А как рисовать график, я, если честно, не очень себе представляю - температура ведь зависит от координаты и от времени.. Это у тебя должен получиться анимированный график? Как-то это не наглядно...

[ewert]

Kafari в сообщении #369041 писал(а):

А как рисовать график, я, если честно, не очень себе представляю - температура ведь зависит от координаты и от времени.. Это у тебя должен получиться анимированный график?

Накопите в двумерном массиве значения на нескольких временных слоях -- и выведите тот массив на картинку простейшей функцией mesh. Достаточно красивенько выйдет.

[ohotnik6666]

полностью ввожу код приведенный выше, и выдает ошибку: ??? function [ ] = Temp( T1,T2,T0,A,L )
|
Error: Function definitions are not permitted at the prompt or in scripts.
в чем дело?

[Maik2013]

Kafari Здравствуй, не подскажешь пожалуйста, Программа расчетов как выглядит,
думаю это код сохраняем в M-файл. Очень нуждаюсь.

Используется синтаксис Matlab M

function [  ] = Temp( T1,T2,T0,A,L )
s=figure;
step=0.1;
t=0;
h=1;
Ti=T0*ones(L,1);
Tk=T0*ones(L,1);
colormap(bone);
while ishandle(s)
    Ti(1)=T1;
    Ti(L)=T2;
    Tk(1)=T1;
    Tk(L)=T2;
    for x=2:h:(L-1)
        Tk(x)=(A*step/h^2)*(Ti(x-1)+Ti(x+1)-2*Ti(x)) + Ti(x);
    end
    image(Ti);
    Ti=Tk;
    t=t+step;
    pause(step);
end
end