17083

Метод Крилова побудови власного багаточлена матриці

Лабораторная работа

Информатика, кибернетика и программирование

Лабораторна робота №18 Тема. Метод Крилова побудови власного багаточлена матриці. Мета. Навчитися знаходити власний багаточлен матриці методом Крилова. Устаткування: лист формату А4 ручка програмне забезпечення Borland C Хід роботи Правила техніки безпеки ...

Украинкский

2013-06-29

66 KB

1 чел.

Лабораторна робота №18

Тема. Метод Крилова побудови власного багаточлена матриці.

Мета. Навчитися знаходити власний багаточлен матриці методом Крилова.

Устаткування: лист формату А4, ручка, програмне забезпечення Borland C++

Хід роботи

  1.  Правила техніки безпеки

Теоретичні дані

Дана  матриця А=  ,

- характеристичний (власний) багаточлен матриці.

- довільний ненульовий вектор.

, , … ,

Система лінійних рівнянь з невідомими  :

або:

яка вирішується одним з методів рішення системи лінійних рівнянь (Метод Гаусса).

Зауваження: якщо система по методу Гаусса не може бути вирішена, то змінити вектор .

№2

 

Листинг программы

#include<iostream.h>

#include<math.h>

void main()

{int i,j;

double x1,x2,x3,x4;

double a[4][5];

double b[4][5];

double c[4][4];

double *y0=new double [4];

double *y1=new double [4];

double *y2=new double [4];

double *y3=new double [4];

double *y4=new double [4];

cout<<"vvedite matrizhu\n";

for(i=0;i<4;i++)

{for(j=0;j<4;j++)

{cin>>c[i][j];}}

y0[0]=1;         y0[1]=0;        y0[2]=0;       y0[3]=0;     y1[0]=0.0;

for(j=0;j<4;j++)

{y1[j]=y0[0]*c[j][0]+y0[1]*c[j][1]+y0[2]*c[j][2]+y0[3]*c[j][3];}

for(j=0;j<4;j++)

{y2[j]=y1[0]*c[j][0]+y1[1]*c[j][1]+y1[2]*c[j][2]+y1[3]*c[j][3];}

for(j=0;j<4;j++)

{y3[j]=y2[0]*c[j][0]+y2[1]*c[j][1]+y2[2]*c[j][2]+y2[3]*c[j][3];}

for(j=0;j<4;j++)

{y4[j]=y3[0]*c[j][0]+y3[1]*c[j][1]+y3[2]*c[j][2]+y3[3]*c[j][3];}

for(i=0;i<4;i++)

{a[i][3]=y0[i];}

for(i=0;i<4;i++)

{a[i][2]=y1[i];}

for(i=0;i<4;i++)

{a[i][1]=y2[i];}

for(i=0;i<4;i++)

{a[i][0]=y3[i];}

for(i=0;i<4;i++)

{a[i][4]=-y4[i];}

if(a[0][0]==0) {cout<<"metod Gaysa ne primenim";

return;}

for(j=0;j<5;j++)

{b[0][j]=a[0][j]/a[0][0];}

for(j=0;j<5;j++)

{b[1][j]=a[1][j];}

for(j=0;j<5;j++)

{a[1][j]=a[1][j]-b[0][j]*b[1][0];}

for(j=0;j<5;j++)

{b[2][j]=a[2][j];}

for(j=0;j<5;j++)

{a[2][j]=a[2][j]-b[0][j]*b[2][0];}

for(j=0;j<5;j++)

{b[3][j]=a[3][j];}

for(j=0;j<5;j++)

{a[3][j]=a[3][j]-b[0][j]*b[3][0];}

if(a[1][1]==0) {cout<<"metod Gaysa ne primenim";

return;}

for(j=0;j<5;j++)

{b[1][j]=a[1][j]/a[1][1];}

for(j=0;j<5;j++)

{b[2][j]=a[2][j];}

for(j=0;j<5;j++)

{a[2][j]=a[2][j]-b[1][j]*b[2][1];}

for(j=0;j<5;j++)

{b[3][j]=a[3][j];}

for(j=0;j<5;j++)

{a[3][j]=a[3][j]-b[1][j]*b[3][1];}

if(a[2][2]==0) {cout<<"metod Gaysa ne primenim";

return;}

for(j=0;j<5;j++)

{b[2][j]=a[2][j]/a[2][2];}

for(j=0;j<5;j++)

{b[3][j]=a[3][j];}

for(j=0;j<5;j++)

{a[3][j]=a[3][j]-b[2][j]*b[3][2];

}

x4=a[3][4]/a[3][3];

x3=b[2][4]-b[2][3]*x4;

x2=b[1][4]-b[1][2]*x3-b[1][3]*x4;

x1=b[0][4]-b[0][1]*x2-b[0][2]*x3-b[0][3]*x4;

cout<<"x1="<<x1<<"\nx2="<<x2<<"\nx3="<<x3<<"\nx4="<<x4;}

Відповідь:

Приклад:

Відповідь:

Висновок: я навчилася знаходити власний багаточлен матриці методом Крилова


 

А также другие работы, которые могут Вас заинтересовать

76701. Нормы права и его структура 180 KB
  На известном этапе развития общества возникает потребность упорядочивания существующих в нем отношений. Таким образом выделяются: нормыобычаи нормыморали нормыправа и корпоративные нормы. Здесь права и обязанности совпадают это является главным отличием его от других социальных норм.
76702. Роль интегрирования PR при создании и поддержании имиджа бренда 86.5 KB
  Актуальность реферата заключается в том, что подобно человеку товар имеет свое лицо, которое может обеспечить ему рыночный успех или стать причиной полного провала. Лицо товара - это некий уникальный образ, существующий в сознании потребителя.
76703. Организация раздаточных 653.92 KB
  Нечеткая работа раздаточной может привести к снижению качества готовых блюд и ухудшает обслуживание посетителей. Это повышает их ответственность за качество оформление и правильный выход блюд. На предприятиях с самообслуживанием когда торговый зал работает продолжительное время блюда отпускают раздатчицы.
76704. Оздоровительная физическая культура 30.99 KB
  Физическая культура часть образа жизни человека система специальных упражнений и спортивной деятельности направленная на развитие его физических и духовных сил. Физическая культура как часть общей культуры направлена на гармоническое развитие всех природных сущностных сил...
76705. Печёночная долька. Железистый и сосудистый варианты строения 971.38 KB
  Структура и функции печени взаимосвязанные вероятностные сущности. Представления о структуре печени это предмет анатомии печени и гистологии печени. Представления о функциях печени это предмет физиологии печени и биохимии печени.
76706. ПЕРВЫЙ ЭТАП СЕСТРИНСКОГО ПРОЦЕССА 423.5 KB
  Она может быть получена от самого пациента членов его семьи свидетелей случившегося самой медицинской сестры её коллег. Для этого в первую очередь необходимо определить проблемы пациента связанные с актуальным или потенциальным состоянием его здоровья.
76707. Los Angeles 33.76 KB
  Los Angeles officially the City of Los Angeles, often known by its initials L.A., is the most populous city in the U.S. state of California and the second-most populous in the United States, after New York City, with a population at the 2010 United States Census of 3,792,621.
76708. Biography of Winnie the Pooh 29.95 KB
  Like many other characters from Milne, Winnie was given the name of one of the real toys Christopher Robin (1920-1996), son of the writer. In turn, the teddy bear Winnie the Pooh was named for the mother bear named Winnipeg (Winnie) contained in the 1920 at the London Zoo.
76709. Развитие и современное состояние автомобилизации 658.64 KB
  Значение автомобильной промышленности и перспективы её развития определяются тем, какое место занимает автотранспорт в транспортно-энергетической инфраструктуре, и его общей ролью в национальной экономике той или иной страны.