Произведение матриц

Обсуждение вопросов, связанных с Partial Differential Equations Toolbox и Comsol MultiPhysics (Femlab)Анализ данных и статистика в MATLAB. Проектирование систем управления в MATLAB\Simulink

Модератор: Admin

Solu
Пользователь
Сообщения: 1
Зарегистрирован: Чт фев 01, 2018 2:41 am

Произведение матриц

Сообщение Solu » Чт фев 01, 2018 3:00 am

Добрый день. Никак не могу понять, где кроется ошибка

Код: Выделить всё

clc 
clear

syms phi B A p q s r t v i j k l m n a b

epsl(1,1)=0;
epsl(2,2)=0;
epsl(1,2)=1;
epsl(2,1)=-1;

for i=1:2
    for j=1:2
        A(i,j)=symsum(symsum(((kroneckerDelta(i,p)*kroneckerDelta(j,q)+kroneckerDelta(i,p)*epsl(j,q))*cos(phi)+...
            (epsl(i,p)*kroneckerDelta(j,q)+epsl(i,p)*epsl(j,q))*sin(phi))*B(p,q), p, 1, 2), q, 1, 2)
    end
end

Код: Выделить всё

Error using sym/subsindex (line 836)
Invalid indexing or function definition. Indexing must follow MATLAB indexing. Function arguments must be symbolic variables, and
function body must be sym expression.

Error in Untitled2 (line 13)
        A(i,j)=sum(sum(((kroneckerDelta(i,p)*kroneckerDelta(j,q)+kroneckerDelta(i,p)*epsl(j,q))*cos(phi)+...

пробовала и sum и symsum, но ничего не получается. нужно найти произведение матриц, где суммирование идет, как указано, и получить матрицу А. Пробовала в MuPAD то же самое. Сперва все выписал, но во второй раз стал игнорировать индексы p, q. непонятно. И как можно в Matlab и MuPAD задать символы Леви-Чивиты?