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;}

Відповідь:

Приклад:

Відповідь:

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


 

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

1596. Клинические исследования молочной железы 18.95 KB
  Молочная железа, вымя коровы - железистый орган, состоящий из 4 четвертей; каждая из них внизу заканчивается соском. Вымя осматривают сзади и сбоку; при этом обращают внимание на его форму, сохранность волосяного покрова, цвет кожи; выявляют повреждения, кожные заболевания или их следы
1597. Конструкция искусственных вагин используемых в ветеринарии 19.23 KB
  Искусственная вагина - прибор, состоящий из цилиндра, изготовленного из металла, резины или эбонита, в просвет которого вставлена эластичная резиновая трубка. Искусственная вагина для быка имеет цилиндр из толстой резины и эластичную камеру, концы которой завернуты на концы цилиндра.
1598. Контроль за животными в послеродовый период. Ранняя акушерско-гинекологическая диспансеризация на фермах 21.01 KB
  Необходимость этого обусловлена тем, что организм животного, ослабленный беременностью и родами, в этот период необходимо в рацион включать легкопереваримые корма, богатые витаминами и минеральными солями, и скармливать эти корма нужно небольшими порциями.
1599. Кормление, содержание и эксплуатация производителей 21.45 KB
  Кормление производителей по энергетическому уровню должно быть таким, чтобы они постоянно находились в заводской кондиции. Во избежание ожирения или, напротив, снижения упитанности их ежемесячно взвешивают.
1600. Краткая история развития ветеринарного акушерства, гинекологии и биотехники размножения животных. Роль и достижения отечественных ученых 21.79 KB
  Элементы ветеринарного акушерства зарождались в то время, когда человек начал приручать животных, и возникла необходимость в простейших приемах родовспоможения. Постепенно накапливался акушерский опыт у определенных людей, появились так называемые знахари, которые передавали свои знания другим поколениям.
1601. Макроскопическая оценка качества спермы: объем, цвет, консистенция, запах 18.5 KB
  Доброкачественная сперма должна содержать достаточное количество живых, устойчивых во внешней среде и способных принять участие в оплодотворении спермиев.
1602. Маститы. Распространение и экономический ущерб 20.23 KB
  Мастит – воспаление молочной железы, развивающиеся как следствие воздействия механических, термических, химических и биологических факторов.
1603. Методика и степень разбавления спермы 19.18 KB
  Для разбавления спермы используют специальные среды, создающие следующие преимущества: получается большой объем спермы, что имеет важное практическое значение для интенсивного использования ценных производителей, в сперме, смешанной с некоторыми разбавителями, спермии переходят из активного в анабиотическое состояние.
1604. Методы определения течки, охоты, полового и овуляции у разных видов животных 20.38 KB
  Течка — выделение прозрачной слизи из половых органов самки. Характеризуется покраснением и набуханием слизистой оболочки преддверия влагалища и влагалищной части шейки матки.