۱۳۹۵/۱۱/۰۹, ۰۶:۲۲ ق.ظ
کد:
for i=3:2:m-3
for j=3:2:n-3
((dX(i+1,j+2)+dX(i+1,j))+(dY(i+2,j+1)+dY(i,j+1)))*Pprim(i+1,j+1)=dX(i+1,j+2)*Pprim(i+1,j+3))+(dX(i+1,j)*Pprim(i+1,j-1))+(dY(i+2,j+1)*Pprim(i+3,j+1))+(dY(i,j+1)*Pprim(i-1,j+1))+((us(i+1,j)-us(i+1,j+2)))+(vs(i,j+1)-(vs(i+2,j+1)))% inside control volume
end
end
با سلام
برای حل کردن یک دستگاه معادلات با n مجهول که مقدار n میتواند مقادیر مختلفی باشد به چه صورت است . مثلا برای این مثال تمام متغییر ها موجود میباشد و ماتریس Pprim که یک ماتریس n*n می باشد که فرض میکنیم یه ماتریس صفر باشد وبه جای مقادیری که در بالا ادرس دهی شده مقدار درست قرار بگیرد .
در کل برای دستگاه معادلاتی که تعداد n معادله , n مجهول دارد به چه صورت است . یعنی نوشتن ماتریس ضرایب عملا ممکن نیست و چون بالفرض یک بار این ماتریش 60*60 و یک بار دیگر با تعریف n می شود 50*50 و .....
با تشکر
کد پیاچپی:
for i=3:2:m-3
for j=3:2:n-3
((dX(i+1,j+2)+dX(i+1,j))+(dY(i+2,j+1)+dY(i,j+1)))*Pprim(i+1,j+1)=dX(i+1,j+2)*Pprim(i+1,j+3))+(dX(i+1,j)*Pprim(i+1,j-1))+(dY(i+2,j+1)*Pprim(i+3,j+1))+(dY(i,j+1)*Pprim(i-1,j+1))+((us(i+1,j)-us(i+1,j+2)))+(vs(i,j+1)-(vs(i+2,j+1)))% inside control volume
end
end