Расчет элементов матрицы и циклическая нумерация - Matlab - Страница 2

Зосима · 11.01.2015, 21:13

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

Маленький Лунат, солнце, а может в формуле ошибка? или данных(M,K)? :sratch:

@Маленький Лунат · 19.01.2015, 11:54 **[ТС]**

а можно ли еще в этот цикл добавить проверку найденных переменных? Если U = 0, то присвоить ему значение предыдущего элемента, а если и предыдущий ноль, то пре-предыдущего

Зосима · 19.01.2015, 15:49

Сообщение от Маленький Лунат

Если U = 0, то присвоить ему значение предыдущего элемента

не совсем понял, что ты хочешь

У нас есть 3 штуковины:
U(:,i+1), U(:,i), U(:,i-1)
что из них нужно сравнивать?

@Маленький Лунат · 19.01.2015, 16:00 **[ТС]**

U(:,i+1) вот эту штуковину)

Зосима · 19.01.2015, 16:25

Маленький Лунат, имеешь ввиду такой вариант?

Matlab M

clear, clc
 
M = 5+rand(54);
K = diag(2*ones(1,54))+diag(ones(1,51),3)+diag(ones(1,51),-3);
E = eye(54);
 
n = 1e4;
U = zeros(54,n); % заготовка из нулей
t = 0:0.001:(n-1)*0.001; % массив времени
for i = 2:n-1
    q = repmat([sin(7*i*t(i));0;0],18,1); % повторяем вектор 18 раз
    W = M\(-M*q - K*U(:,i))*t(i)^2 + E*(2*U(:,i) - U(:,i-1));
    if W==0;
        U(:,i+1) = U(:,i);
    else
        U(:,i+1) = W;
    end
end
plot(t,U)

но загвоздка в том, что W нулем не становится, поэтому все без изменений

@Маленький Лунат · 19.01.2015, 21:10 **[ТС]**

как же не становится нулем?!..а на графике

Matlab M

1	plot (U, 'DisplayName','U', 'YDataSource', 'U'); figure(gcf)

видно,что идет синусоида и определенное расстояние, которое идет по нулевой оси..

Зосима · 20.01.2015, 09:36

Маленький Лунат, ну дык у тебя там 54 кривые

не исключено, что где-то нули (ибо вектор q их содержит)
но мне думается не в нулях загвоздка

система почему-то неустойчива

@Маленький Лунат 0 / 0 / 0 Регистрация: 08.12.2014 Сообщений: 16
	19.01.2015, 11:54 [ТС]	22
	а можно ли еще в этот цикл добавить проверку найденных переменных? Если U = 0, то присвоить ему значение предыдущего элемента, а если и предыдущий ноль, то пре-предыдущего 0

Зосима 5243 / 3571 / 379 Регистрация: 02.04.2012 Сообщений: 6,473 Записей в блоге: 17
	19.01.2015, 15:49	23
	Сообщение от Маленький Лунат Если U = 0, то присвоить ему значение предыдущего элемента не совсем понял, что ты хочешь У нас есть 3 штуковины: U(:,i+1), U(:,i), U(:,i-1) что из них нужно сравнивать? 0

Зосима 5243 / 3571 / 379 Регистрация: 02.04.2012 Сообщений: 6,473 Записей в блоге: 17
	20.01.2015, 09:36	27
	Маленький Лунат, ну дык у тебя там 54 кривые не исключено, что где-то нули (ибо вектор q их содержит) но мне думается не в нулях загвоздка система почему-то неустойчива 0

Зосима 5243 / 3571 / 379 Регистрация: 02.04.2012 Сообщений: 6,473 Записей в блоге: 17
	11.01.2015, 21:13	21
	Маленький Лунат, солнце, а может в формуле ошибка? или данных(M,K)? :sratch: 0

@Маленький Лунат 0 / 0 / 0 Регистрация: 08.12.2014 Сообщений: 16
	19.01.2015, 16:00 [ТС]	24
	U(:,i+1) вот эту штуковину) 0