سلامبا سلام خدمت دوستان عزیز
میخواستم ببینم کسی میتونه راجب به انتقال حرارت یک بعدی در استوانه و کره با کد نویسی متلب کمکم کنه.؟
1- شرایط مرزی دو طرف جداره ها ره ثابت فرض کنید.
2- ضخامت استوانه 2 میلی متر شعاع داخلیش 300 میلیمتره و ضخامت کره هم 6 میلیمتره و شعاع داخلیش هم 220 میلیمتره
3- انتقال حرارت یک بعدی و فقط در جهت شعاعی است.
من خودم که برنامه نویسیش کردم دماهام منفی در میان
اگه کسی کد نوشتشو داره یا میتونه کمک کنه ممنون می شم.
function y = smpl4plt(x)
if x<1
y = x^3;
elseif x<3
y = x;
else
y = -x^2 + 8*x - 12;
end% main_smpl4plt
x = -2:0.1:5;
n = length(x);
y = zeros(1,n);
for k = 1:n
y(k) = smpl4plt(x(k));
end
plot(x,y)
xlabel('x'), ylabel('y')
set(gca,'XTick',[-2,1,3,5],'YTick',[1,3])
axis([-2,5,-8,4.4])
grid
n = length(x);
y = zeros(1,n);
>> x=0:0.00001:1;
y=x.^3;
plot(x,y,'b')
hold on
x=1:0.00001:2;
w=sin(x);
plot(x,w,'r')
hold on
x=2:0.00001:3;
z=cos(x);
plot(x,z,'g')
سلام...کد سختی نبود،یه خورده وقت گذاشتم و درست شد،اصلا هم نیاز به روش خاصی و انجام گزاره شرطی هم نبود...
تایع X^3 و sinx و cosx رو در بازه های مختلف رسم کردم و نمودارش هم به شرح زیر هست..نکتش فقط نوشتن Hold on هست.بفرمایید:
کد:>> x=0:0.00001:1; y=x.^3; plot(x,y,'b') hold on x=1:0.00001:2; w=sin(x); plot(x,w,'r') hold on x=2:0.00001:3; z=cos(x); plot(x,z,'g')
دقیقا همینو داخل پنجره Comman window کپی و پیست نمایید.
در ابتدای هر دستور که میخواین غیرفعال شه شیفت رو نگه دارید و کلید 5 رو بزنین.
یاعلی
ممنون
i/20ii;
بچه ها من این کد رو برای نشون دادن خطوط اطراف استوانه چرخان نوشتم.ولی اجرا نمیشه.میشه کمکم کنین
x=0:0.01
y=pi:0.01:2*pi;
s=0:1:3;r=roots(s);s=-r.*sin(x)+sim(x)./r-log(r)./2.*pi;
polar(x,r1);
hold on;
polar(y,r1);
hold on;
polar(x,r2);
hold on;
polar(y,r2);
hold on;
polar(x,r3);
hold on;
polar(x,r4);
i;
