Двумерное нестационарное уравнение теплопроводности методом конечных разностей

@gar · Регистрация: 09.01.2015

Author24 — интернет-сервис помощи студентам

Здравствуйте, необходимо решить по неявной схеме. В решении использую метод прогонки.

Алгоритм переведен из кода в матлабе. Никак не могу заставить его работать.

Помогите люди добрые и умные.

@VSI · 09.01.2015, 17:26

gar, в Mathcad'е при создании расчетного документа Вы ОБЯЗАНЫ следовать ОСНОВНОМУ правилу: размещение вычислительных блоков в документе (согласно алгоритму Вашего расчета) должно осуществляться СЛЕВА НАПРАВО и СВЕРХУ ВНИЗ...

@gar · 09.01.2015, 17:34 **[ТС]**

Спасибо за ответ, исправил.

Не подскажите в чем ошибка?

@VSI · 09.01.2015, 17:43

gar, то, что на розовом фоне, вызывает "непонятки"... Еще: вычисление а0 и в0 происходит ОДИНАКОВО?... и т.д.
P.S. Иногда проще самому с 0 написать, чем разгребать чужой код...

@Krasme · 09.01.2015, 17:48

gar, или матрица, или вектор, что-то одно должно быть..

Кликните здесь для просмотра всего текста

Двумерное нестационарное уравнение теплопроводности методом конечных разностей

@gar · 09.01.2015, 18:06 **[ТС]**

Привожу исходный код в матлабе

Matlab M

for i = 1:1:N
        for   j = 1:1:M
                y(i,j) = u(i,j);
        end;
......
 
for k = 1:K1
    for j = 1:N    
        x(j,: ) = f_matmod2(N, a, b, c, y(j,: ));

Как в маткаде отправить вектор из матрицы? Или нужно выводить поэлементно?

@Krasme · 09.01.2015, 18:13

gar, в маткаде можно выделить столбец матрицы.... y^<j>
или строку, предварительно транспонировав.... (y^T)^<j>

@gar · 09.01.2015, 18:28 **[ТС]**

Krasme, Спасибо, но все равно ругается:

Эта ошибка возникает из за разного количества элементов в векторах?

@Krasme · 09.01.2015, 18:33

Сообщение от gar

Эта ошибка возникает из за разного количества элементов в векторах?

да, из-за этого

@gar · 09.01.2015, 19:40 **[ТС]**

Krasme,
Поставил везде одинаковое ранжирование, ошибка не уходит

@Krasme · 09.01.2015, 19:48

индексы массивов не могут быть отрицательными

Кликните здесь для просмотра всего текста

@gar · 09.01.2015, 20:05 **[ТС]**

Знаю, перепробывал все варианты и в конце с психу поставил везде нули.

Если ставить отсчет с 2, то получаю ошибку деления на ноль.

Не пойму где конкретно не совпадает количество элементов. Можно ли как то вывести промежуточный результат вычисления?

@Krasme · 09.01.2015, 20:27

желтое - временно вырезано и отключено вычисление...
видны размерности все ваших параметров...

Кликните здесь для просмотра всего текста

ничего, что матрица с - другого размера?

@gar · 10.01.2015, 16:49 **[ТС]**

Вроде разобрался с количеством элементов, результата не выдает, ждал больше 1,5 часа.

Нид хэлп

@Krasme · 10.01.2015, 17:05

gar, какая у вас версия?
у меня mathcad 15 M030.. все работает, дает хоть какой-то результат, правильно ли, я не знаю..

Кликните здесь для просмотра всего текста

@Symon · 10.01.2015, 17:59

Сообщение от gar

for i = 1:1:N
for j = 1:1:M
y(i,j) = u(i,j)

1. Ваша функция tdma заканчивается строкой X_i. Это значит, результатом работы программы является одно число, равное последней координате вектора X. Думаю, что вам надо получить вектор Х, а не одно число. Тогда уберите индекс в этой строке.
2. В другой функции U я заметил неэффективный двойной цикл (см. приведенную цитату). Не проще ли записать просто y<-u.
3. В других деталях не разбирался.

@Krasme 6833 / 4893 / 2066 Регистрация: 02.02.2014 Сообщений: 13,057
	09.01.2015, 19:48	11
	индексы массивов не могут быть отрицательными Кликните здесь для просмотра всего текста 0

@gar 0 / 0 / 0 Регистрация: 09.01.2015 Сообщений: 8
	09.01.2015, 20:05 [ТС]	12
	Знаю, перепробывал все варианты и в конце с психу поставил везде нули. Если ставить отсчет с 2, то получаю ошибку деления на ноль. Не пойму где конкретно не совпадает количество элементов. Можно ли как то вывести промежуточный результат вычисления? Изображения 0

@VSI Модератор $Эксперт по математике/физике$ 5246 / 4032 / 1386 Регистрация: 30.07.2012 Сообщений: 12,331
	09.01.2015, 17:26	2
	gar, в Mathcad'е при создании расчетного документа Вы ОБЯЗАНЫ следовать ОСНОВНОМУ правилу: размещение вычислительных блоков в документе (согласно алгоритму Вашего расчета) должно осуществляться СЛЕВА НАПРАВО и СВЕРХУ ВНИЗ... 0

@VSI Модератор $Эксперт по математике/физике$ 5246 / 4032 / 1386 Регистрация: 30.07.2012 Сообщений: 12,331
	09.01.2015, 17:43	4
	gar, то, что на розовом фоне, вызывает "непонятки"... Еще: вычисление а0 и в0 происходит ОДИНАКОВО?... и т.д. P.S. Иногда проще самому с 0 написать, чем разгребать чужой код... 0

@Krasme 6833 / 4893 / 2066 Регистрация: 02.02.2014 Сообщений: 13,057
	09.01.2015, 17:48	5
	gar, или матрица, или вектор, что-то одно должно быть.. Кликните здесь для просмотра всего текста 0

@Krasme 6833 / 4893 / 2066 Регистрация: 02.02.2014 Сообщений: 13,057
	09.01.2015, 18:13	7
	gar, в маткаде можно выделить столбец матрицы.... y^<j> или строку, предварительно транспонировав.... (y^T)^<j> 1

@gar 0 / 0 / 0 Регистрация: 09.01.2015 Сообщений: 8
	09.01.2015, 18:28 [ТС]	8
	Krasme, Спасибо, но все равно ругается: Эта ошибка возникает из за разного количества элементов в векторах? Миниатюры 0

@Krasme 6833 / 4893 / 2066 Регистрация: 02.02.2014 Сообщений: 13,057
	09.01.2015, 18:33	9
	Сообщение от gar Эта ошибка возникает из за разного количества элементов в векторах? да, из-за этого 0

@Krasme 6833 / 4893 / 2066 Регистрация: 02.02.2014 Сообщений: 13,057
	09.01.2015, 20:27	13
	желтое - временно вырезано и отключено вычисление... видны размерности все ваших параметров... Кликните здесь для просмотра всего текста ничего, что матрица с - другого размера? 1

@Krasme 6833 / 4893 / 2066 Регистрация: 02.02.2014 Сообщений: 13,057
	10.01.2015, 17:05	15
	gar, какая у вас версия? у меня mathcad 15 M030.. все работает, дает хоть какой-то результат, правильно ли, я не знаю.. Кликните здесь для просмотра всего текста 1

@Symon $Эксперт по математике/физике$ 2615 / 2229 / 684 Регистрация: 29.09.2012 Сообщений: 4,578 Записей в блоге: 13
	10.01.2015, 17:59	16
	Сообщение от gar for i = 1:1:N for j = 1:1:M y(i,j) = u(i,j) 1. Ваша функция tdma заканчивается строкой X_i. Это значит, результатом работы программы является одно число, равное последней координате вектора X. Думаю, что вам надо получить вектор Х, а не одно число. Тогда уберите индекс в этой строке. 2. В другой функции U я заметил неэффективный двойной цикл (см. приведенную цитату). Не проще ли записать просто y<-u. 3. В других деталях не разбирался. 1