Решение задач математической физики в системе MatLab. Рындин Е.А - 28 стр.

UptoLike

ВУЗ: 

ТТИ ЮФУ | Таганрог

Составители: 

Рубрика: 

% Преобразование вектора-строки значений искомой функции
% в узлах координатной сетки в матрицу размерности n x m x s,
% удобную для представления результатов в графическом виде
for l=1:s
for i=1:n
for j=1:m
U(i,j,l)=u(n*m*(l-1)+m*(i-1)+j);
end
end
end
% Построение графика искомой функции U(x,y,t)
for l=1:s
figure
surf(y,x,U(:,:,l))
xlabel('y','FontSize',13)
ylabel('x','FontSize',13)
zlabel('U','FontSize',13)
grid on
colormap('cool')
axis([min(y) max(y) min(x) max(x) min(min(min(U))) max(max(max(U)))])
pause(0.1)
M(l)=getframe;
end
for l=s+1:2*s-2
figure
surf(y,x,U(:,:,2*s-l))
xlabel('y','FontSize',13)
ylabel('x','FontSize',13)
zlabel('U','FontSize',13)
grid on
colormap('cool')
axis([min(y) max(y) min(x) max(x) min(min(min(U))) max(max(max(U)))])
pause(0.1)
M(l)=getframe;
end
% Отображение волнового процесса в динамическом режиме
figure
ans=1;
while ans==1
movie(M,10,10)
ans=menu('Повторить просмотр результатов?','ДА','НЕТ');
end
Приведенное выше текстовое описание функции сохраняется, как и в
предыдущем примере, в виде m-файла с именем f_wave2d.m. В отличие от
m-файла, приведенного выше, все переменные и массивы в данном случае по
умолчанию локальны, т.е. после выполнения функции они автоматически
удаляются из оперативной памяти, за исключением возвращаемых функцией
28

Страницы