اوه اوه چرا اینکار عجیب رو کردید؟!!! این کار که حرفش رو زدید (منظورم توضیحتونه، نه برنامه اش) معنیش اینه که شما اصلاً نمی دونید روش عددی حل دستگاههای IVP چیکار میکنن.
برای اینکه از متلب خواستید یک بردار سطری 165درایه ای رو به عنوان سطر اول، و یک بردار سطری 235درایه ای رو به عنوان سطر دوم در ماتریس B بریزه؛ روشنه که همچین کاری غیر ممکنه.
احتمالاً ناوارد بوده.
دومی ضرب و تقسیم اسکالر با برداره؛ واضحه که نیازی به عملیات درایه به درایه نیست. کسی که مسلط به متلب باشه، دومی رو به کار نمی بره.
دو راه بیشتر ندارید؛ یا باید تابع مورد نظر در پوشه جاری باشه، یا باید آدرس اون پوشه رو (با استفاده از set path یا add path یا pathdef) به مسیرهای متلب اضافه کنید.
x = [1; 2; 3; 8; 9; 10];
y = [9; 8; 2; 3; 1; 4];
z = [2; 8; 6; 4; 9; 7];
[xi, yi] = meshgrid(x, y);
f = scatteredInterpolant(x, y, z);
zi = f(xi, yi);
surf(xi, yi, zi)
البته کلی تنظیمات داره؛ که من همه رو پیش فرض در نظر گرفتم.
داری اذیت می کنی دیگه!!! یعنی یک بردار که 6 تا درایه داشته باشه. نکنه "بردار" و "درایه" نمی دونی به چی میگن؟!!!
در ضمن اون مشکل همون موقع حل شد؛ یعنی برای بردار 6درایه ای جواب میده.
با رنگ آبی در بالا، نظر دادم.
سؤال من که باز هم فراموش کردید جواب بدید، این بود که صورت مسأله به صورت ریاضی چیه؟ یعنی همونی که از شما خواستن که حل کنید، چیه؟
خب پس تمرکزتون روی اصطلاح غلطیه؛ توی حل عددی معادلات دیفرانسیل کوپل و غیرکوپل اصلاً مطرح نمی شه؛ چون رویه های حل عددی اصلاً براشون مهم نیست که یک دستگاه معادلات دیفرانسیل کوپل باشند یا نه. توی حل تحلیلی (حل معادلات با دست) کوپل بودن و نبودن مطرح میشه و حل معادلات کوپل (با دست) سخت تره.
در مورد...
بله متلب برای حلشون تابع داره. از این سؤالتون پیداست که هنوز نمی دونید دنبال چی هستید. در مبحث حل عددی معادلات دیفرانسیل معمولی، سه نوع مسأله مطرح میشه: IVP, BVP, DAE
با توجه به اون سؤال قبلیتون که در مورد جرم و فنر پرسیده بودید، حدس من اینه که شما می خواید یک مسأله IVP رو حل کنید. در این صورت...
برای اینکه جواب سؤال 1 رو بگیری؛ دستور زیر و نتیجه اش رو ببین.
>> logspace(1,5,5)
ans =
10 100 1000 10000 100000
>>
در مورد سؤال 2، فلسفه خاصی نداره؛ اگه ربطی به یک مسأله فیزیکی نداشته باشه، فقط یک انتخاب بد می تونه باشه.