Вычислить найхитрейшую сумму - Matlab - Страница 2

Зосима · 03.03.2015, 19:38

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

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

чтобы все элементы каждой ячейки каждого массива перемножились, ну например между собой. mass*mass

например MAT{3}.*MAT{1}, но размерности могут не совпасть

@dieORalive · 03.03.2015, 19:42 **[ТС]**

Попробую, спасибо

@dieORalive · 06.03.2015, 14:06 **[ТС]**

Доброго времени суток, очередной вопросик, "модифицированный первый вопрос". На этом мы остановились:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
X = [pi/6 pi/4 pi/3];
a = [ 0.2178+0.0380i 0.7779+0.0516i 0.2796-0.1718i ];
c = [ 0.7844-0.0376i 0.2163-0.0547i 0.7279+0.1828i ];
 
S = zeros(size(X));
for k = 1:numel(X);
    x = X(k);
    b(1) = 1;
    b(2) = 3*cos(x);
    d(1) = cos(x);
    for n = 1:3
        if n > 2
            b(n) = ((2*n-1)/(n-1))*cos(x)*b(n-1)-(n/(n-1))*b(n-2);
        end
        if n > 1
            d(n) = n*cos(x)*b(n)-(n+1)*b(n-1);
        end
        
        S(k) = S(k) + 2*(n+1)/(n*(n+1))*(a(n)*b(n)+c(n)*d(n));
    end
end

Теперь у меня a(n), b(n), c(n) и d(n) вида [1xk cell] (строчка, длинны k), каждая ячейка заполнена в зависимости от Nmax ([Nmin,Nmax] оно растёт). С ростом Nmax кол-во элементов в ячейке растёт, но, например i-ый элемент a(n) имеет такую же размерность (по Nmax), как и b(n), c(n) и d(n). Вобщем вопрос состоит в том, как правильно суммировать такой тип данных. Тоесть, код должен брать a(n), b(n), c(n) и d(n) для k=const(допустим для k=1) Nmax=3, в каждом массиве по 3 элемента, они перемножаются и суммируются, получается одна значение S(k=1), и тд.
Каждый массив считается так, но это не важно, просто для справки:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
a = 0.001:0.001:0.2;
X = 2*pi*a;
 
for k = 1:length(X)
    X1 = X(k); % берётся значение к-го элемента массива Х, с ним идёт цикл по Nmax
    Nmax = round( X1+4*(X1^(1/3))+2 );
    mass = zeros(1,Nmax);
    mass(1) = (sin(X1)/X1)-cos(X1);
    mass(2) = (3/X1)*((sin(X1)/X1)-cos(X1))-sin(X1);
    for n = 3:Nmax
        mass(n) = ((2*(n-1)+1)/X1)*mass(n-1)-mass(n-2);
    end
    MAT{k} = mass; % сохраняем массив в ячейку
end

Найхитрейшая сумма.
У меня пока 2 варианта, либо эти массивы перевести в другой тип данных, либо в матрицу, но тогда при суммировании там надо будет придумывать нули, для успешного перемножения. И второй вариант - наличие в матлабе какой то функции вида: S = modifsum(Nmin, Nmax, datatype) для конкретного типа данных

посоветуте чтото пожалуйста.

Зосима,
Решил перейти к матрицам

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
a = 0.001:0.001:0.2;
X = 2*pi*a;
 
for k = 1:length(X)
    X1 = X(k); % берётся значение к-го элемента массива Х, с ним идёт цикл по Nmax
    Nmax = round( X1+4*(X1^(1/3))+2 );
    %MATT = zeros(3,7);
    mass = zeros(1,Nmax);
    mass(1) = (sin(X1)/X1)-cos(X1);
    mass(2) = (3/X1)*((sin(X1)/X1)-cos(X1))-sin(X1);
    for n = 3:Nmax
        mass(n) = ((2*(n-1)+1)/X1)*mass(n-1)-mass(n-2);
        
    end
    
    MAT{k} = mass; % сохраняем массив в ячейку
    MATT(:,k) = MAT{1,k}.';
end

Код проходит 7 эл-тов цикла для N=3, дальше выдаёт ошибку(Subscripted assignment dimension mismatch.) , т.к. N8=4 и увеличивается вертикальная размерность матрицы на 1, подскажите как исправить эту ошибку, или же как сразу записать результаты в матрицу, не используя тип данных cell, муторный он какой то, и из него нельзя переходить обратно к матрицам, есть там конечно cell2mat, но работает он как то не так как я планировал

Зосима · 14.03.2015, 18:57

К сожалению у менянет сейчас возможности во всем разобраться и дать исчерпывающий ответ

Если основная задача - рассчет суммы, то с матрицами и ячейками можно не заморачиваться

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
clear, clc
global a c
a = 0.001:0.001:0.2;
X = 2*pi*a;
 
S = zeros(size(X)); 
for k = 1:length(X)
    X1 = X(k); % берётся значение к-го элемента массива Х, с ним идёт цикл по Nmax
    Nmax = round( X1+4*(X1^(1/3))+2 );
  
    for n = 1:Nmax
        if n==1
           b(n) = cos(pi*X1);
           d(n) = (sin(X1)/X1)-cos(X1);
        elseif n==2
           b(n) = 1-sqrt(sin(X1);
           d(n) = (3/X1)*((sin(X1)/X1)-cos(X1))-sin(X1);
        else
           b(n) = cos(X1) + b(n-1) - sin(b(n-2));
           d(n) = ((2*(n-1)+1)/X1)*d(n-1)-d(n-2);
        end
        S(k) = S(k) + (2*n+1)/(n*(n+1))*( a(n)*b(n) + c(n)*d(n) );
    end
end
plot(X, real(S), X, imag(S) )

@dieORalive · 17.03.2015, 12:20 **[ТС]**

Зосима, Добрый день

Поразмыслил я и перевёл в матрицы все данные, не знаю вот только как считаются такие вещи внутри суммы, как a(n)*b(n), где a и b это матрицы разных размерностей - a(n,x), b(n,y)

Зосима · 17.03.2015, 12:43

dieORalive, тут уже и y затесался?

было бы проще, если бы иметь данные (формулы) для всех переменных

@dieORalive · 17.03.2015, 16:36 **[ТС]**

да там похожие формулы просто, я уж там их сам модифицировал, я конечно могу, если интересно скинуть полный код в личку

Сообщение от Зосима

Зосима · 18.03.2015, 10:28

кидай

Новые блоги и статьи Все статьи Все блоги /
Когда логика программы не спасает от человеческих ошибок Maks 18.08.2026 В последнее время всё чаще и чаще сталкиваюсь с таким явлением, как абсолютная невнимательность (или глупость) пользователей. Проявляется это чаще всего на работе в коллективе. Допустим, человек с. . .	Лето уходит kumehtar 17.08.2026	Мысли в слух kumehtar 17.08.2026 Забавно, насколько сейчас стала доступна информация. Например о магии, духовном развитии, медитациях, и других подобных направлениях, ранее зачастую тайных, передаваемых от учителя к ученику. Хотя. . .	Перемещение строк из ТЧ в другой документ с учетом текущего пробега Maks 17.08.2026 Реализация из решения ниже выполнена на примере нетипового документа "Автозапчасти", с ТЧ "Шины". За основу взят алгоритм отсюда: https:/ / www. cyberforum. ru/ blogs/ 359708/ 10838. html Задача: . . .
Саморегулирующийся социальный контракт для сервера cross-section. Hrethgir 14.08.2026 С кодом конечно таких глубоких размышлений пока не было, впрочем я уже привык к алгоритмизации. Суть предмета записи: снова в диалоге с нейросетью (я взял пока себе ник для учётки админа - Rector). . . .	Часы электронные Uhbif79 12.08.2026 Выкладываю программу часов. Программа позволяет: 1. Использовать системное время и дату, 2. Есть возможность вводить время и дату вручную. 3. Реализованы 2 будильника: начало и конец рабочего дня. . . .	Часы с будильником на основе класса QLCDNumber Uhbif79 12.08.2026 Всем добрый день, выкладываю программу часов с будильником на основе класса QLCDNumber. Здесь я пробовал самостоятельно создавал классы, впервые столкнулся с видимостью переменной одного класса из. . .	Установка MinGW GCC 16.2 и CMake 8Observer8 10.08.2026 VK Видео: https:/ / vkvideo. ru/ video-240781534_456239017 YouTube: eY5-5PyI9NM Текстовая версия

@dieORalive 0 / 0 / 0 Регистрация: 13.09.2014 Сообщений: 21
	03.03.2015, 19:42 [ТС]
	Попробую, спасибо 0

@dieORalive 0 / 0 / 0 Регистрация: 13.09.2014 Сообщений: 21
	17.03.2015, 12:20 [ТС]
	Зосима, Добрый день Поразмыслил я и перевёл в матрицы все данные, не знаю вот только как считаются такие вещи внутри суммы, как a(n)*b(n), где a и b это матрицы разных размерностей - a(n,x), b(n,y) 0

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	17.03.2015, 12:43
	dieORalive, тут уже и y затесался? было бы проще, если бы иметь данные (формулы) для всех переменных 0

Зосима 5246 / 3574 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	18.03.2015, 10:28
	кидай 0