35273

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

Практическая работа

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

h void min { int klj; double [3][3]b[3][3]y0[3]y1[3]y2[3]y3[3]y4[3]yn1yn2yn3yn4Sum1Sum2Sum3Sum4; double x1x2x3x4d0d1d2d3102030213132; cout Vvedite mtritsy endl; fork=0;k =3;k { forl=0;l =3;l cin b[k][l]; } cout Vvedite nylevou vektor endl; fork=0;k =3;k cin y0[k]; fork=0;k =3;k { forl=0;l =3;l { yn1=b[k][l]y0[l]; Sum1=Sum1yn1; } y1[k]=Sum1; } fork=0;k =3;k { forl=0;l =3;l { yn2=b[k][l]y1[l]; Sum2=Sum2yn2; } y2[k]=Sum2; } fork=0;k =3;k {...

Украинкский

2013-09-09

36 KB

0 чел.

еревко О.В.

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

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

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

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

Хід роботи

  1.  Правила техніки безпеки
  2.  Теоретичні дані
  3.  Індивідуальне завдання.

Знайти власний багаточлен матриці А по методу Крилова, використовуючи метод Гаусса.

№16

#include <iostream.h>

void main()

{

int k,l,j;

double a[3][3],b[3][3],y0[3],y1[3],y2[3],y3[3],y4[3],yn1,yn2,yn3,yn4,Sum1,Sum2,Sum3,Sum4;

double x1,x2,x3,x4,d0,d1,d2,d3,a10,a20,a30,a21,a31,a32;

cout<<"Vvedite matritsy A"<<endl;

for(k=0;k<=3;k++)

{

 for(l=0;l<=3;l++)

 cin>>b[k][l];

}

cout<<"Vvedite nylevou vektor"<<endl;

for(k=0;k<=3;k++)

cin>>y0[k];

for(k=0;k<=3;k++)

{

 for(l=0;l<=3;l++)

 {

  yn1=b[k][l]*y0[l];

  Sum1=Sum1+yn1;

 }

 y1[k]=Sum1;

}

for(k=0;k<=3;k++)

{

 for(l=0;l<=3;l++)

 {

  yn2=b[k][l]*y1[l];

  Sum2=Sum2+yn2;

 }

 y2[k]=Sum2;

}

for(k=0;k<=3;k++)

{

 for(l=0;l<=3;l++)

 {

  yn3=b[k][l]*y2[l];

  Sum3=Sum3+yn3;

 }

 y3[k]=Sum3;

}

for(k=0;k<=3;k++)

{

 for(l=0;l<=3;l++)

 {

  yn4=b[k][l]*y3[l];

  Sum4=Sum4+yn4;

 }

 y4[k]=Sum4;

}

cout<<"Koordinati vektora y1:"<<endl;

for(k=0;k<=3;k++)

cout<<" "<<y1[k];

cout<<endl;

cout<<"Koordinati vektora y2:"<<endl;

for(k=0;k<=3;k++)

cout<<" "<<y2[k];

cout<<endl;

cout<<"Koordinati vektora y3:"<<endl;

for(k=0;k<=3;k++)

cout<<" "<<y3[k];

cout<<endl;

cout<<"Koordinati vektora y4:"<<endl;

for(k=0;k<=3;k++)

cout<<" "<<y4[k];

a[0][0]=y3[0];a[0][1]=y2[0];a[0][2]=y1[0];a[0][3]=y0[0];a[0][4]=-y4[0];

a[1][0]=y3[1];a[1][1]=y2[1];a[1][2]=y1[1];a[1][3]=y0[1];a[1][4]=-y4[1];

a[2][0]=y3[2];a[2][1]=y2[2];a[2][2]=y1[2];a[2][3]=y0[2];a[2][4]=-y4[2];

a[3][0]=y3[3];a[3][1]=y2[3];a[3][2]=y1[3];a[3][3]=y0[3];a[3][4]=-y4[3];

d0=a[0][0];

if(d0==0)cout<<"Metod Gayssa ne premenim"<<endl;return;

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

{a[0][j]=a[0][j]/d0;

}

a10=a[1][0];

a20=a[2][0];

a30=a[3][0];

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

{a[1][j]=a[1][j]-a[0][j]*a10;

a[2][j]=a[2][j]-a[0][j]*a20;

a[3][j]=a[3][j]-a[0][j]*a30;

}

d1=a[1][1];

if(d1==0)cout<<"Metod Gayssa ne premenim"<<endl;return;

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

{a[1][j]=a[1][j]/d1;

}

a21=a[2][1];

a31=a[3][1];

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

{a[2][j]=a[2][j]-a[1][j]*a21;

a[3][j]=a[3][j]-a[1][j]*a31;

}

d2=a[2][2];

if(d2==0)cout<<"Metod Gayssa ne premenim"<<endl;return;

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

{a[2][j]=a[2][j]/d2;

}

a32=a[3][2];

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

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

}

d3=a[3][3];

if(d3==0)cout<<"Metod Gayssa ne premenim"<<endl;return;

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

{a[3][j]=a[3][j]/d3;

}

x4=a[3][4];

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

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

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

cout<<"p1="<<x1;

cout<<"p2="<<x2;

cout<<"p3="<<x3;

cout<<"p4="<<x4;

}

  1.  Контрольні питання
  •  Дайте визначення власного багаточлена матриці.
  •  Сформулюйте задачу знаходження власного багаточлена матриці по методу Крилова.
  •  Яка ідея методу Крилова знаходження власного багаточлена?
  •  Які допоміжні методи використовуються в методі Крилова?
  •  Як продовжити рішення, якщо по будь-якому методу отримана система не має єдиного рішення?

3


 

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

1253. Дільниця виготовлення гачкових колекторів електродвигунів приводу підлоготерів 403.5 KB
  Призначення та галузь застосування електродвигунів постійного струму малої потужності. Опис конструкції колектора електродвигуна малої потужності. Технічні умови на виготовлення колектора. Заходи по вдосконаленню діючого технологічного процесу. Нормування технологічного процесу виготовлення гачкового колектора.
1254. Расчет затрат участка по ремонту приборов системы питания АТП на 500 автомобилей ПАЗ-2305 123.5 KB
  Классификация предприятий автомобильного транспорта. Малые предприятия и индивидуальное предпринимательство на АТП. Расчет фонда заработной платы рабочих. Расчет накладных расходов, связанных с общей организацией цехового производства.
1255. Проектирование двухэтажного жилого здания 164 KB
  Объемные строительные системы, имеющие надземную и (или) подземную части, включающие в себя помещения, сети инженерно-технического обеспечения и системы инженерно-технического обеспечения. Здание в плане имеет неправильную форму . Здание двухэтажное с высотой этажа3,3м. Размеры здания в осях 14,3х13,5м. Отметка подошвы фундамента -1,85 м. Грунтовые воды отсутствуют.
1256. Создание предприятия ИП Британ 73.5 KB
  Наименование предприятия, его специализация форма собственности и адрес. Схема действующего гаражно-технологического процесса, его описание. Описание методов организации технологических процессов ТО и ремонта ПС. Ведомость оборудования и оснастки на предприятии по форме.
1257. Соціально-психологічне дослідження ревнощів чоловiкiв 115.5 KB
  Науково-теоретичні підходи до розуміння та дослідження почуття ревнощів чоловіків. Літературно-філософські підході до розуміння природи ревнощів. Методика емпірико-практичного дослідження почуття ревнощів у чоловіків. Емпірико-практичне дослідження особливостей почуття ревнощів у чоловіків. Експериментальне дослідження співвідношення ревнощів з іншими почуттями та емоціями.
1258. Технология изготовления детали подъемника (вала редуктора) 91.5 KB
  Тихоходный вал редуктора, массой 11 кг., является телом вращения. Вал считается жёстким, поскольку отношение длины к максимальному диаметру меньше 15. Выбор заготовки и метод её изготовления. Определение групп допусков на исходные размеры заготовки.
1259. Разработка технических требований к сборочной единице 349.5 KB
  Техническое описание сборочной единицы. Выбор степени точности зубчатой передачи и класса точности подшипников. Выбор посадок шлицевых соединений и резьбовых соединений. Выбор исполнительных размеров калибров и их расчет. Выбор параметров точности зубчатого колеса.
1260. Анализ рынка. Абсолютная доля рынка по сегментам 375 KB
  Абсолютная доля рынка по сегментам. Уровень овладения рынком или горизонтальное проникновение. Уровень приверженности марке. Модель формирования позиции отдельного потребителя по отношению к продукту. Средняя важность атрибута.
1261. Економічне обґрунтування діяльності підприємства 497 KB
  Дані методичні вказівки покликані допомогти студентам напряму Архітектура, Реставрація творів мистецтва і Дизайн (рівень бакалавра) в написанні, а головне у виконанні тих певних обчислень, які необхідні для визначення вартості підприємницького проекту та обґрунтуванні виробничої програми.