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


 

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

48366. Термічний аналіз. Температура та методи її вимірювання 961.5 KB
  Температура та методи її вимірювання Питання що розглядаються: Визначення температури. Шкали температур. Методи вимірювання температури.
48367. Технологические основы машиностроения 2.48 MB
  Формы организации технологических процессов механической механичной обработки Основные факторы которые влияют на точность обработки Способы обработки металлов резанием и элементы режима резания О прогрессивных методах обработки
48368. СУЧАСНА УКРАЇНСЬКА ЛІТЕРАТУРНА МОВА. СТИЛІ МОВИ 742 KB
  Виникнення української мови. Мова як соціальне явище: а соціальна сутність мови; б варіювання мови в суспільстві; в взаємодія мов в багатоетнічному суспільстві; г мова і культура. Загальна структура мови.
48369. Бухгалтерський облік у торговельних підприємствах 747 KB
  Торговельна діяльність це ініціативна самостійна діяльність юридичних осіб і громадян по здійсненню купівлі та продажу товарів з метою одержання прибутку. Роздрібна торгівля це продаж товарів кінцевому споживачеві в розмірах що відповідають його щоденним некомерційним потребам а сам товар як правило придбавається на оптових базах або в оптових продавців. Оптова торгівля це продаж торговельними підприємствами товарів іншим підприємствам і організаціям роздрібної торгівлі ресторанного господарства суб'єктам підприємництва. Магазини...
48370. Цели автоматизации конструкторско-технологической подготовки производства 13.61 KB
  В своем историческом развитии САПР ТП постепенно расширяли арсенал своих средств. Однако подавляющее большинство САПР ТП в том числе и ныне существующих не способны поддерживать автоматизацию принятия решений в процессе проектирования на основе технологических знаний. История развития САПР ТП показала бесперспективность алгоритмического подхода.
48371. Основные требования к применяемым комплексам конструкторско-технологической подготовки производства 12.93 KB
  Основные требования к применяемым комплексам конструкторско-технологической подготовки производства Системы комплекса должны быть взаимосвязаны что позволяет постепенно наращивать количество рабочих мест и объединять их в единый комплекс. Минимальную стоимость комплекса и обеспечение всех необходимых предприятию функции. Внедряемый комплекс САПР должен удовлетворять следующим критериям: Соответствие функциональных возможностей этапам основного и вспомогательного производства на конкретном предприятии Возможность безбумажного обмена...
48373. Подготовка конструкторской документации с использованием трехмерных моделей деталей 12.58 KB
  После построения 3Dмодели детали или сборки либо непосредственно в ходе построения конструктор может получить ее чертеж избежав таким образом рутинного создания видов средствами плоского черчения. Плоский чертеж будет создан автоматически и с абсолютной точностью независимо от сложности модели. В Компас3D объемные модели и плоские чертежи ассоциативны между собой. Компас3D располагает мощными средствами редактирования модели которые позволяют задавать параметрические связи и ассоциации как между отдельными элементами деталей так и...
48374. Подготовка конструкторской документации с использованием трехмерных моделей сборок 12.6 KB
  После построения 3Dмодели детали или сборки либо непосредственно в ходе построения конструктор может получить ее чертеж избежав таким образом рутинного создания видов средствами плоского черчения. Плоский чертеж будет создан автоматически и с абсолютной точностью независимо от сложности модели. В Компас3DV объемные модели и плоские чертежи ассоциативны между собой. Компас3D располагает мощными средствами редактирования модели которые позволяют задавать параметрические связи и ассоциации как между отдельными элементами деталей так...