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

Відповідь:

Приклад:

Відповідь:

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


 

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

52450. Использование компьютерных технологий в исследовательской деятельности 29 KB
  Исследовательские работы учащихся в наши дни часто опираются на возможности которые даёт компьютер. Подчас использование компьютерных технологий сводится к представлению результатов исследования ребёнок создаёт презентацию своей работы. На самом деле презентация это самый последний этап работы и самый сложный из всех. Сложный этап определение целей и постановка задач исследовательской работы.
52451. Environmental Situation And Climate In The World 146.5 KB
  Procedure Stges Chirmn Guests Students Notes Introduction Wrmingup Speking Summing up Home ssignment Ldies nd gentlemen Let me open our conference. There re representties of the Ecologicl orgniztions from meric Gret Britin Ukrine the Climte Experts the Professor of Oxford University nd Students of this University t our conference. Mr Simon nd Mr Winter their interesting story of mking the cquintnce of the pretty girl nd invittion her to our conference.
52455. Національна кухня 11.66 MB
  CDplyer crds. Pudding is typicl English dish. Stek nd fries risotto but we lso do severl English dishes s well. Were rgentinin people surprised when they herd tht n English chef ws going to open resturnt here Yes they were very I think people do not usully expect the English to be good cooks.
52456. National cuisine 532 KB
  Tody we continue to discuss different questions nd fcts bout mels nd food. ll living things need food becuse it is the bsic fuel of our life. Wht food is very good for us Which of the foods re not so good for us nd wht food cn even cuse illnesses Write list of 15 best nd 15 worst foods. Nowdys in spoken English you my often come cross such phrse s junk food.
52457. Cultural Development of Young People 45 KB
  Der friends Tody t the lesson well try to revel the notion âCulturl development of young peopleâ.: I think Ill be doctor becuse its n importnt nd necessry profession which helps people to keep in good helth. Pupils brinstorm their ides nd come to the conclusion tht â culturl developmentâ includes the following things: eduction upbringing socil surroundings trditions the wy we tret other people the wy we work nd study the wy we dress...
52458. We and the English-speaking world. Comparing cultures 139.5 KB
  Guests: investors from Germny representtives of school dminist rtion correspondents from locl mss medi techers of the lyceum. Im gld to introduce you our guests from Germny Peter Freundlieb nd Julius Zitsek. nd I hope our guests from Germny will help us wont they The guests introduce themselves nd tell some words bout their country nd their mission in Ukrine. I think they will explin us the reson: boy nd girl worn in ntionl Ukrinin costumes hnd in the Germn guests bred nd slt on the ntionl Ukrinin towel nd spek bout Ukrinin...