35264

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

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

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

h void min {int ij; double x1x2x3x4; double [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 Введите матрицу n ; fori=0;i 4;i {forj=0;j 4;j {cin c[i][j];}} y0[0]=1; y0[1]=0; y0[2]=0; y0[3]=0; y1[0]=0.0; forj=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];} forj=0;j 4;j...

Украинкский

2013-09-09

90.5 KB

2 чел.

еревко О.В.

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

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

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

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

Хід роботи

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

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

№16

#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<<"Введите матрицу\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<<"Метод Гаусса не применим";

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<<"Метод Гаусса не применим";

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<<"Метод Гаусса не применим";

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;

}

Результат программы

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

Як продовжити рішення, якщо по будь-якому методу отримана система не має


 

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

4040. Социально-этические и гуманистические принципы биологического познания. Генетика и эволюция 108.5 KB
  Введение Современное человечество живет в новом тысячелетии и это заставляет людей уделять большее внимание своему будущему и разумному осмыслению нашего прошлого. Одним из важнейших стала проблема взаимодействия природы и общества. Противоречия меж...
4041. Криоэлектроника. История развития 255.5 KB
  Криоэлектроника. История развития Введение Криогенная (от греческого "криос" — холод, мороз) электроника, или криоэлектроника — направление электроники, охватывающее исследование при криогенных температурах (ниже 120 К) специфических эффек...
4042. Хаос и порядок. Порядок и беспорядок в природе 311 KB
  Этимология понятия «хаос». Соотношение порядка и беспорядка в природе. Хаос, понятие окончательно оформившееся в древнегреческой философии - это трагический образ космического первоединства, начало и конец всего, вечная смерть всего живого и одно...
4043. Методика определения плотности твердого тела 42 KB
  Определение плотности твердого тела. Цель работы Целью данной работы является: Определение плотности твердого тела Приобретение студентами практических навыков в обработке результатов измерений. Рисунок Рабочая формула...
4044. Определение момента инерции диска (стержня) по его крутильным колебаниям 50.5 KB
  Определение момента инерции диска (стержня) по его крутильным колебаниям. Цель работы Целью данной работы является определение момента инерции диска (стержня) относительно оси, проходящей через центр тяжести перпендикулярно плоскости диска...
4045. Определение положения центра тяжести и момента инерции физического маятника 47.5 KB
  Определение положения центра тяжести и момента инерции физического маятника. Цель работы Целью данной работы являются: определение ускорения свободного падения при помощи математического маятника определение положения центра тяжес...
4046. Определение коэффициента внутреннего трения жидкости по методу Стокса 79 KB
  Определение коэффициента внутреннего трения жидкости по методу Стокса. Цель работы. Целью данной работы является определение коэффициента внутреннего трения (вязкости) жидкости. Рисунок. Рабочая формула: Таблица результатов и...
4047. Определение сопротивления с помощью моста Уитстона 85.5 KB
  Определение сопротивления с помощью моста Уитстона I. Цель работы – определение связи между силой тока в цепи гальванического элемента тока и падением напряжением на внешнем участке цепи. Расчет на основании этих данных величины ЭДС, внутреннег...
4048. Исследование электростатических полей методом электролитической ванны 48.5 KB
  Исследование электростатических полей методом электролитической ванны Целью данной работы является изучение электростатических полей, создаваемых заряженными проводниками различной конфигурации, путем моделирования в электролитической ванне Ри...