Подякувати Студентському архіву довільною сумою
Admin
26.02.2023 12:38
Дякуємо, що користуєтесь нашим архівом!
ВИКОРИСТАННЯ ФУНКЦІОНАЛЬНОЇ ДЕКОМПОЗИЦІЇ ДЛЯ РОЗВ’ЯЗКУ ОБЧИСЛЮВАЛЬНИХ ЗАДАЧ
Інформація про навчальний заклад
ВУЗ:
Національний університет Львівська політехніка
Інститут:
Не вказано
Факультет:
КН
Кафедра:
Кафедра ЕОМ
Інформація про роботу
Рік:
2011
Тип роботи:
Лабораторна робота
Предмет:
Паралельні та розподілені обчислення
Частина тексту файла (без зображень, графіків і формул):
Мiнiстерство освiти і науки, молоді та спорту України
Національний університет “Львівська політехніка”
Кафедра ЕОМ
Лабораторна робота №1
з диципліни: «Паралельні і розподілені обчислення»
на тему: «ВИКОРИСТАННЯ ФУНКЦІОНАЛЬНОЇ ДЕКОМПОЗИЦІЇ ДЛЯ РОЗВ’ЯЗКУ ОБЧИСЛЮВАЛЬНИХ ЗАДАЧ»
МЕТА РОБОТИ. Вивчити методи декомпозицій задач. Набути навиків розв’язування задач з використанням функціональної декомпозиції.
Завдання:
Вираз, який слід обрахувати, заданий наступним чином:
x = (y2y1Y32)(y1’y1) + y1y2’ + Y3 + y2’Y3y1Y3;
При чому елементи визначаються згідно правил:
y1=A×b, де bi=18/(i+18)2 , і=1,2,...n
y2=A1(b1-c1);
Y3 = A2B2-A2C2, де Cij=18/(i+2j)
Послідовність виконання.
1. Аналіз завдання.
Для заданого виразу вхідними даними є:
розмірність матриць – n;
матриці ;
вектори-стовпці .
Ці параметри повинні вводитися з клавіатури, або генеруватися випадковим чином (крім розмірності). При чому, елементи всіх матриць та векторів є цілими додатними числами, більшими за нуль.
Вектор-стовпець та матриця обраховуються, виходячи з уведеної розмірності, зауважимо, що значення їх елементів завжди менші одиниці і різко спадають зі збільшенням розмірності.
Наприклад для n=3, значення вектора-стовпця будуть становити: b1=18/(1+18)2=0.0499; b2=18/(2+18)2=0.0450; b3=18/(3+18)2=0.0408,а значення матриці , відповідно:
C11=18/(1+2*1)=6.0 C12=18/(1+2*2)=3.6 C13=18/(1+2*3)=2.5714
C21=18/(2+2*1)=4.5 C22=18/(2+2*2)=3.0 C23=18/(2+2*3)=2.25
C31=18/(3+2*1)=3.6 C32=18/(3+2*2)=2.5714 C33=18/(3+2*3)=2.
При утворенні враховуємо, що результатом множення матриці А на вектор-стовпець b є вектор-стовпець, елементи якого будуть раціональними числами(тобто матимуть значущу дробову частину).
При утворенні враховуємо, що результатом віднімання двох векторів-стовпців є вектор-стовпець, елементи якого можуть бути меншими за нуль цілими числами. Далі, при множенні цілочисельної додатної матриці А1 на результат віднімання, отримаємо вектор-стовпець з цілочисельними елементами довільного знаку.
При утворенні враховуємо, що присутні лише операції віднімання та множення, а тому вихідний результат може бути додатнім або відємним і завжди матиме значущу дробову частину.
Таким чином, згідно поставленої задачі, в обчисленні загального виразу приймають участь три різні елементи – два вектори стовпці та матриця .
Перший доданок загального виразу містить шість множники – два транспоновані вектор-стовпці y1,(тобто вектор-рядки), вектор-стовпець , вектор-стовпець y2 та матрицю піднесену до квадрату. Оскільки, згідно правил матричних обчислень, добуток не є комутативною операцією, всі множення слід виконувати в тій послідовності, яка задана, вищий пріоритет мають тільки множники які знаходяться в дужках, вони виконуються першими. Результатом множення рядка на стовпець є число, числа на стовпець є стовпець, а матриці на матрицю – матриця, а число на матрицю дасть матрицю. Тому, в загальному, перший доданок буде матрицею.
Другий доданок складається з двох множників вектора-стовпця y1 і транспонованого вектора стовпця y2 (вектора рядка). В результаті множення стовпця на рядок отримаємо матрицю.
Третій доданок складається з однієї матриці Y3.
Четвертий доданок складається з 4 множників, де 2 множники матриця Y3 вектор-стовпець y1 і транспонований вектор стовпець y2 (Вектор-рядок). При множенні вектор-рядок на матрицю отримаємо вектор-рядок, при множенні вектора-рядок на вектор-стовпець отримаємо скаляр і при множенні скаляра на матрицю отримаємо матрицю.
Таким чином, з попереднього випливає, що остаточний результат є матриця, елементи якої можуть бути як додатними так і від’ємними і завжди мають дробову частину.
2. Декомпозиція задачі.
Однозначно, всі обчислення безпосередньо залежать від розмірності даних, тому найперше, слід забезпечити ввід змінної n, що визначає цю розмірність. Далі, можна паралельно виконувати обчислення значень вектора b та матриці С2, оскільки вони незалежні від інших параметрів. Крім того, на тому ж рівні декомпозиції слід визначати вхідні дані, тобто вводити з клавіатури, або генерувати випадковим чином матриці та вектори-стовпці . Наступний рівень декомпозиції – це знаходження елементів виразу. Значення залежить від введеної матриці А та обрахованого вектора b. Значення залежить від введеної А1 та різниці векторів b1 і c1, тому знайти його можна лише після обчислення (b1 - c1). Зауважимо, що множення на константу не є окремою операцією, як і транспонування векторів. Аналогічно, знаходимо . Подальша декомпозиція відбувається згідно заданої послідовності операцій та врахування залежностей отриманих на кожному рівні даних. Повна схема декомпозиції обчислення заданого виразу приведена нижче.
3. Об’єднання частин виразу проведено безпосередньо у схемі декомпозиції, оскільки воно однозначно визначається порядком обчислень.
4. Для написання програми, що ілюструє процес обчислення виразу згідно розробленої схеми декомпозиції, слід врахувати наступні особливості:
Дані, що вводяться з клавіатури, або генеруються випадковим чином повинні бути цілими числами, більшими за нуль. Результати проміжних обрахунків будуть містити дробову частину, виняток складає лише елемент y2. Тому слід обирати відповідний тип даних(багатобайтний).
Оскільки на окремо взятому рівні декомпозиції обраховуються незалежні частини підвиразів, функції що їх програмно реалізують повинні викликатися псевдо-одночасно.(Тобто перш ніж виконувати будь яку функцію третього рівня, слід виконати всі функції другого рівня і т.д.)
Схема декомпозиції обчислення виразу
5. Текст програми
n = input('n = ');
r = input('Щоб задати масиви вручну введіть 1, рандомом 2: ');
if (r == 2)
A = randint(n,n,10);
mess = 'A = ';
dlmwrite('1.txt', mess, 'delimiter', '');
dlmwrite('1.txt', A, '-append','delimiter', '\t', 'precision',15);
A1 = randint(n,n,10);
mess = 'A1 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', A1, '-append','delimiter', '\t', 'precision',15);
b1 = randint(n,1,10);
mess = 'b1 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', b1, '-append','delimiter', '\t', 'precision',15);
c1 = randint(n,1,10);
mess = 'c1 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', c1, '-append','delimiter', '\t', 'precision',15);
A2 = randint(n,n,10);
mess = 'A2 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', A2, '-append','delimiter', '\t', 'precision',15);
B2 = randint(n,n,10);
mess = 'B2 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', B2, '-append','delimiter', '\t', 'precision',15);
else
A = randint(n,n,10);
A1 = randint(n,n,10);
b = randint(n,1,10);
b1 = randint(n,1,10);
c1 = randint(n,1,10);
A2 = randint(n,n,10);
B2 = randint(n,n,10);
C2 = randint(n,n,10);
for i = 1:n
for j = 1:n
A(i,j) = input('A[i,j] = ');
end
end
for i = 1:n
for j = 1:n
A1(i,j) = input('A1[i,j] = ');
end
end
for i = 1:n
b1(i) = input('b1[i] = ');
end
b1 = b1';
for i = 1:n
c1(i) = input('c1[i] = ');
end
c1 = c1';
for i = 1:n
for j = 1:n
A2(i,j) = input('A2[i,j] = ');
end
end
for i = 1:n
for j = 1:n
B2(i,j) = input('B2[i,j] = ');
end
end
end;
for i = 1:n
b(i) = 18/(i+18)^2;
end
mess = 'b = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', b', '-append','delimiter', '\t', 'precision',15);
y1 = A*b';
mess = 'y1 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', y1, '-append');
temp1 = b1-c1;
y2 = A1*temp1;
mess = 'y2 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', y2, '-append');
for i = 1:n
for j = 1:n
C2(i,j) = 18/(i+2*j);
end
end
mess = 'C2 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', C2, '-append','delimiter', '\t', 'precision',15);
Temp1 = A2*B2;
mess = 'Y3 = ';
Temp2 = A2*C2;
Y3 = Temp1 - Temp2;
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', Y3, '-append','delimiter','\t');
t1 = y1'*y1;
mess = '(y1"*y1) = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t1, '-append','delimiter', '\t', 'precision',15);
t2 = y2*y1';
mess = 'y2*y1 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t2, '-append','delimiter', '\t', 'precision',15);
t3 = Y3^2;
mess = 'Y3^2 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t3, '-append','delimiter', '\t', 'precision',15);
t4 = t2*t3;
mess = '(y2*y1*Y3^2) = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t4, '-append','delimiter', '\t', 'precision',15);
t5 = t4*t1;
mess = '(y2*y1*Y3^2)(y1"*y1) = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t5, '-append','delimiter', '\t', 'precision',15);
t6 = y1*y2';
mess = 'y1*y2" = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t6, '-append','delimiter', '\t', 'precision',15);
t7 = y2'*Y3;
mess = 'y2"*Y3 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t7, '-append','delimiter', '\t', 'precision',15);
t8 = t7*y1;
mess = 'y2"*Y3*y1 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t8, '-append','delimiter', '\t', 'precision',15);
t9 = t8*Y3;
mess = 'y2"*Y3*y1*Y3';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t9, '-append','delimiter', '\t', 'precision',15);
t10 = t5 + t6;
mess = '(y2*y1*Y3^2)(y1"*y1) + y1*y2" = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t10, '-append','delimiter', '\t', 'precision',15);
t11 = t10+Y3;
mess = '(y2*y1*Y3^2)(y1"*y1) + y1*y2" + Y3 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', t11, '-append','delimiter', '\t', 'precision',15);
x = t11+t9;
mess = 'x = (y2*y1*Y3^2)(y1"*y1) + y1*y2" + Y3 + y2"*Y3*y1*Y3 = ';
dlmwrite('1.txt', mess,'-append', 'delimiter', '');
dlmwrite('1.txt', x, '-append','delimiter', '\t', 'precision',15);
6. Результати виконання програми
x = (y2*y1*Y3^2)(y1"*y1) + y1*y2" + Y3 + y2"*Y3*y1*Y3 =
-434197.79050664 -478021.665828007 -15923.6024351779
-137150.519147917 43514.6935823115 -15361.1347316156
-330288.387510253 -321856.839011331 -35483.0658862551
Висновки
В роботі використано паралелізм на рівні підзадач, оскільки передбачається, що кожен блок зі схеми декомпозиції є реалізований у виді функції. Це є середньоблоковий паралелізм. Обмін даними відбувається через використання спільних змінних. Присутня залежність даних між різними рівнями декомпозиції, але в межах одного рівня її немає. Є залежність за керуванням, оскільки послідовність обчислювального процесу наперед однозначно відома. Залежність за ресурсами та вводом/виводом може бути визначена лише у відношенні до певної обчислювальної системи.
28.03.2013 13:03-
Ділись своїми роботами та отримуй миттєві бонуси!
0%
Оголошення від адміністратора