es.es.es
یک شنبه 20 بهمن 1392, 20:51 عصر
سلام دوستان من اصلا با متلب برنامه نویسی نکردم یکی از استاد ها برای پروژه یه تیکه کد بهم داده گفته مشکلشو برطرف کنید یه جورایی خودشم بلد نیست چون برنامه نویس نیست این کد با فرمول های خودش میخواد یه نمودار رسم کنه ولی نمودارش درست در نمیاد و یه خط صاف میشه میخوام ببینم مشکلش کجاست
clc
clear
k0=7*10^-7;
h1=10^(-5);
h2=10^(-11);
T(100,1)=300;
t(1)=10^-9;
r(100)=10^(-3);
%-------------------------------------------------------------
for j1=99:1:-1
for n1=1:99
n=n1*h2+h2;
j=j1*h1+h1;
t(n1)=n;
r(j1)=j;
T(j1,1)=300;
T(j1,n1+1)=k0*h2*((2/r(j1)*(T(j1+1,n1)-T(j1,n1))/h1+(T(j1+1,n1)-2*T(j1,n1)+T(j1-1,n1))/h1^2+3*(1/T(j1,n1))*
((T(j1+1,n1)-T(j1,n1))/h1)^2)*T(j1,n1)^3)+T(j1,n1);
end
end
%----------------------------------------------
figure
plot(r,T)
title('soluttion numerical temperature');
legend(numerical,100 mesh points,0);
xlable('raduce r');
ylable('tempreature T');
%-------------------------------------------------------------
figure
plot(t,T)
title('soluttion numerical temperature');
legend(numerical,100 mesh points,0);
xlable('Time t');
ylable('tempreature T');
clc
clear
k0=7*10^-7;
h1=10^(-5);
h2=10^(-11);
T(100,1)=300;
t(1)=10^-9;
r(100)=10^(-3);
%-------------------------------------------------------------
for j1=99:1:-1
for n1=1:99
n=n1*h2+h2;
j=j1*h1+h1;
t(n1)=n;
r(j1)=j;
T(j1,1)=300;
T(j1,n1+1)=k0*h2*((2/r(j1)*(T(j1+1,n1)-T(j1,n1))/h1+(T(j1+1,n1)-2*T(j1,n1)+T(j1-1,n1))/h1^2+3*(1/T(j1,n1))*
((T(j1+1,n1)-T(j1,n1))/h1)^2)*T(j1,n1)^3)+T(j1,n1);
end
end
%----------------------------------------------
figure
plot(r,T)
title('soluttion numerical temperature');
legend(numerical,100 mesh points,0);
xlable('raduce r');
ylable('tempreature T');
%-------------------------------------------------------------
figure
plot(t,T)
title('soluttion numerical temperature');
legend(numerical,100 mesh points,0);
xlable('Time t');
ylable('tempreature T');