, , .
. , . , , .
. .
, . , . . :
1) ( );
2) ( );
3) ( ).
, , - .
MatLab, . (. 10).
g='squareg'; %
b='squareb1'; %
c = input('Koiffisient teploprovodnosti ');
a=0;
f=1;
d=1;
:
[p,e,t]=initmesh(g);
:
u0=zeros(size(p,2),1);
ix=find(sqrt(p(1,:).^2+p(2,:).^2)<0.4);
u0(ix)=ones(size(ix));
:
nframes=20;
tlist=linspace(0,0.1,nframes);
:
u1=parabolic(u0,tlist,b,p,e,t,c,a,f,d);
:
x=linspace(-1,1,31);y=x;
[unused,tn,a2,a3]=tri2grid(p,t,u0,x,y);
:
newplot;
Mv = moviein(nframes);
umax=max(max(u1));
umin=min(min(u1));
for j=1:nframes,...
u=tri2grid(p,t,u1(:,j),tn,a2,a3);i=find(isnan(u));u(i)=zeros(size(i));...
surf(x,y,u);caxis([umin umax]);colormap(cool),...
axis([-1 1 -1 1 0 1]);...
Mv(:,j) = getframe;...
u=tri2grid(p,t,u1(:,j),tn,a2,a3);i=find(isnan(u));u(i)=zeros(size(i));...
surf(x,y,u);caxis([umin umax]);colormap(cool),...
axis([-1 1 -1 1 0 1]);...
Mv(:,j) = getframe;...
end
:
movie(Mv,10)
, , .
|
|
, , , .
(, , ). , F, W.
:
( );
(-);
, (-, - . .);
(, .).
: , , , (, ).
815%, 8085% ( ).
. 4.
4.
, /3 | ||
2500 | (, , .) | |
22002500 | , | |
18002200 | , | |
( ) | ||
5001800 | , , | - |
() | ||
500 |