“PEDAGOGS”
international research journal ISSN:
2181-4027
_SJIF:
4.995
www.pedagoglar.uz
Volume-8, Issue-4, April - 2022
253
2. M.file da f1 ni saqlab olinadi.
function
y=f1(x,t)
y=x-t;
end
3. M.file da u0 ni saqlab olinadi.
function
y=u0(x)
y=x-x^2;
end
4. M.file da g1 ni saqlab olinadi.
function
z=g1(t)
z=t;
end
5. M.file da g2 ni saqlab olinadi.
function
y=g2(t)
y=t;
end
2. M.file da pdt_Oshkormas_teskari_matrissa
ni saqlab olinadi. Shu faylni saqlagandan keyin
dasturni ishga tushiramiz. Iloji bo`lsa, bitta papkaga hammasini saqlashimiz lozim. Dasturni ishga
tushirishda aniq harakat qilish kerak.
function
U=pdt_Oshkormas_teskari_matrissa(u0,f,c1,c2,a,b,n,m)
h=a/(n-1); tau=b/(m-1); r=tau/h^2; U=zeros(n,m); A=zeros(n,n); U(1,1:m)=c1; U(n,1:m)=c2;
for
j=1:m
for
i=1:n
f1(i,j)=feval(f,h*(i-1),tau*(j-1));
end
end
for
i=1:n
U(i,1)=feval(
'u0'
,h*(i-1));
end
A(1,1)=1;A(n,n)=1;
for
i=2:n-1
A(i,i-1)=-r;A(i,i)=1+2*r;A(i,i+1)=-r
end
“PEDAGOGS”
international research journal ISSN:
2181-4027
_SJIF:
4.995
www.pedagoglar.uz
Volume-8, Issue-4, April - 2022
254
%AU=D ni yaratish va yechish
D(1)=feval(
'g1'
,0);D(n)=feval(
'g2'
,b);
for
j=2:m
for
i=2:n-1
D(i)=U(i-1,j)+tau*f1(i,j);
end
V=inv(A);X=V*D';U(1:n,j)=X';
end
U=U';
% Command Window oynasiga quyidagi kodni yozish lozim
% [U]=pdt_Oshkormas_teskari_matrissa('u0','f',0,0,1,0.05,10,10)
% Grafigi mesh(U)
3. Natija chiqarish bo`limi.
Command Window oynasiga quyidagilarni yozish lozim
>> [U]=pdt_Oshkormas_teskari_matrissa('u0','f',0,0,1,0.05,10,10)
A =
1.0000 0 0 0 0 0 0 0 0 0
-0.4500 1.9000 -0.4500 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 1.0000
A =
1.0000 0 0 0 0 0 0 0 0 0
-0.4500 1.9000 -0.4500 0 0 0 0 0 0 0
0 -0.4500 1.9000 -0.4500 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0
|