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.  Контрольні питання
  •  Дайте визначення власного багаточлена матриці.
  •  Сформулюйте задачу знаходження власного багаточлена матриці по методу Крилова.
  •  Яка ідея методу Крилова знаходження власного багаточлена?
  •  Які допоміжні методи використовуються в методі Крилова?

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


 

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

73218. Дифракция света 219.5 KB
  Дифракционные решетки и их применение. Основное свойство дифракционной решетки - способность раскладывать падающий на неё пучок света по длинам волн поэтому она используется в качестве диспергирующего элемента в спектральных приборах.
73219. Магнетики 279.5 KB
  В настоящее время установлено, что все материальные тела в той или иной мере обладают магнитными свойствами. Ими обладают не только макроскопические тела, но и отдельные атомы, молекулы, а также электроны и атомные ядра. Все вещества (кроме вакуума) можно назвать магнетиками.
73220. Управление персоналом в системе современного менеджмента 117.5 KB
  Та его часть которая состоит в штате организации называется кадрами. С социально производственной точки зрения иными словами по той роли которую работники играют в организации и в процессе достижения ее целей те же термины называют уже совершенно иные понятия.
73221. Екологічна безпека як складова національної безпеки Основні поняття та завдання, що вирішує дисципліна 105 KB
  Екологічна безпека як складова національної безпеки Основні поняття та завдання що вирішує дисципліна. Сутність обєкти субєкти екологічної безпеки Мета задачі та структура курсу його місце та значення у підготовці фахівців з екології. Становлення та розвиток екологічної безпеки. Критерії та ознаки екологічної безпеки.
73222. Оплата труда на предприятиях 1.43 MB
  Организация оплаты труда определяется 3мя взаимосвязанными факторами: 1 нормированием труда которое определяет основу для расчета ЗП 2 тарифной системой которая включает: тарифную ставку тарифную сетку и тарифный квалификационный справочник и определяет правильность оценки количества и качества труда 3 формами и системами оплаты труда устанавливающими порядок расчета ЗП по категориям сотрудников. До начала экономических реформ регулирование заработной платы в нашей стране осуществлялось в централизованном порядке: все тарифные...
73223. Себестоимость продукции 3.47 MB
  Себестоимость продукции является нижним пределом цены а также выполняет много других функций. Себестоимость транспортной продукции это текущие затраты предприятия на производство и реализацию продукции выраженные в денежной форме. Производственная себестоимость затраты связанные с производством продукции включая общепроизводственные и общехозяйственные расходы.
73224. Управление персоналом на предприятиях 804.02 KB
  Рабочие Служащие это категория работников в деятельности которых преобладают те или иные виды умственного труда. разрабатывают и внедряют в производство новые виды продукции технологические процессы технические и экономические нормативы формы и методы организации производства и труда; обеспечивают производство необходимой документацией материалами оснасткой ремонтным и другим обслуживанием; проводят анализ показателей деятельности предприятия разрабатывают и реализуют мероприятия по их улучшению. В формуле отражается одна из норм...
73225. Введение в экономику предприятий транспорта 1.3 MB
  Материальное производство включает: Промышленность; Сельское и лесное хозяйство; Грузовой транспорт; Связь обслуживающая материальное производство; Строительство; Торговля; Общественное питание; Информационно-вычислительное обслуживание и др. В процессе работы транспорта не создается новая продукция а осуществляется перемещение грузов пассажиров с одного пункта в другой. В процессе перевозок грузов с ними не происходят не физические не химические изменения а происходит лишь изменение их местонахождения по отношению к...
73226. Основные фонды предприятий транспорта 899.56 KB
  Классификация основных фондов предприятия По роли в производственном процессе: Активные: машины и оборудование силовые машины измерительные прибор вычислительная техника инструмент Пассивные: здания строения сооружения передаточные устройства транспорт пх инвентарь. Для производственных предприятий в том числе и предприятий транспорта типичный состав основных фондов следующий: Здания и строения административные здания производственные строения склады и т. Племенной скот Соотношение отдельных групп основных фондов в их общем...