Форум программистов, компьютерный форум CyberForum.ru

Интерполяция в канонической форме - C++

Восстановить пароль Регистрация
 
StalkerSmerch
1 / 1 / 2
Регистрация: 04.03.2012
Сообщений: 159
12.06.2012, 00:11     Интерполяция в канонической форме #1
Где то ошибка, вот только где. Просто ответ не выводит

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
#include <iostream>
#include <stdlib.h>
#include <iostream>
#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <time.h>
#define N 3
#define razm 3
 
using namespace std;
 
void m_e(int k,double mas[][N+1],int n,int otv[]);
 
int main(void)
{
  setlocale(LC_ALL, "Russian");
  double mas[N][N+1];
  double x[N];
  int otv[N];
  int i=0,j=0,k=0,n=0;
  system("cls");
   float massx[razm];
 float massy[razm];
// float mas[razm][razm+1];
//  int i, j, k, n=0;
 printf("Введите x и y %i точек:\n",razm);
 for(i=0;i<razm;i++)
 {
     scanf("%f%f",&massx[i],&massy[i]);
 }
 
for(i=0;i<razm;i++)
{
    for(j=0;j<razm;j++)
    mas[i][j]= (pow(massx[i],j));
    mas[i][razm]=massy[i];
}
  for(i=0;i<razm;i++)
 {
     for(j=0;j<razm+1;j++)
     printf("%f ",mas[i][j]);
     printf("\n");
 }
 
 
  for (i=0;i<n+1;i++)
    otv[i]=i;
  for (k=0;k<n;k++)
  {
    m_e(k,mas,n,otv );
    if (fabs(mas[k][k])<0.0001)
    {
      printf( "нет единственного решения" );
      return (0);
    }
    for (j=n;j>=k;j--)
      mas[k][j]/=mas[k][k];
    for (i=k+1;i<n;i++)
      for (j=n;j>=k;j--)
        mas[i][j]-=mas[k][j]*mas[i][k];
  }
  printf ("Система имеет единственное решение\n");
  for (i=0;i<n;i++)
    x[i]=mas[i][n];
  for (i=n-2;i>=0;i--)
    for (j=i+1;j<n;j++)
      x[i]-=x[j]*mas[i][j];
 
  printf( "Ответ:\n" );
  for (i=0;i<n;i++)
    for (j=0;j<n;j++)
      if (i==otv[j])
    {
        printf("%f\n", x[j] );
        break;
      }
  return (0);
}
 
void m_e(int k,double mas[][N+1],int n,int otv[])
{
  int i,j,i_max=k,j_max=k;
  double temp;
  for (i=k;i<n;i++)
    for (j=k;j<n;j++)
      if ( fabs(mas[i_max][j_max])<fabs(mas[i][j]))
      {
        i_max=i;
        j_max=j;
      }
  for (j=k;j<n+1;j++)
  {
    temp=mas[k][j];
    mas[k][j]=mas[i_max] [j];
    mas[i_max][j]=temp;
  }
  for (i=0; i<n;i++)
  {
    temp=mas[i][k];
    mas[i][k]=mas[i][j_max];
    mas[i][j_max]=temp;
  }
  i=otv[k];
  otv[k]=otv[j_max];
  otv[j_max]=i;
  }
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
12.06.2012, 00:11     Интерполяция в канонической форме
Посмотрите здесь:

C++ Интерполяция Лагранжа
C++ интерполяция по Лагранжу
C++ Многомерная интерполяция
C++ Лагранжевая интерполяция
C++ линейная интерполяция
C++ Интерполяция изображения c++
C++ Интерполяция Лагранжа
C++ Линейная интерполяция
C++ Класс должен соответствовать Ортодоксально канонической форме класса
C++ Квадратичная интерполяция
C++ Квадратичная интерполяция
Тригонометрическая интерполяция C++

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
grizlik78
Эксперт С++
 Аватар для grizlik78
1887 / 1419 / 103
Регистрация: 29.05.2011
Сообщений: 2,967
12.06.2012, 02:11     Интерполяция в канонической форме #2
Там где вывод, циклы до n, а n равно нулю. Вот и не выводится ничего.
StalkerSmerch
1 / 1 / 2
Регистрация: 04.03.2012
Сообщений: 159
14.06.2012, 20:36  [ТС]     Интерполяция в канонической форме #3
grizlik78, Спасибо.
Yandex
Объявления
14.06.2012, 20:36     Интерполяция в канонической форме
Ответ Создать тему
Опции темы

Текущее время: 23:46. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru