VPF::[Численные методы] Монте-Карло

MagicPRO

Дата 24.5.2007, 09:58 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

Профиль
Группа: Участник
Сообщений: 200
Регистрация: 4.10.2006
Где: Kremencity

Репутация: нет
Всего: нет

Народк кто знает ссылку или может у кого-то есть пример решение системы линейных уравнений методом Монте-Карло???Заранее благодарен

--------------------

smile=iwashere

Alexeis

Дата 24.5.2007, 10:17 (ссылка)

(нет голосов)

Загрузка ...

Амеба

Профиль
Группа: Админ
Сообщений: 11743
Регистрация: 12.10.2005
Где: Зеленоград

Репутация: 23
Всего: 459

Модератор: не забывайте в заголовке указывать язык программирования, если это не язык программирования то укажите область из которой эта задача.

Заголовок исправил.

--------------------

Vit вечная память.

Обсуждение действий администрации форума производятся только в этом форуме

гениальность идеи состоит в том, что ее невозможно придумать

Misha_White

Дата 24.5.2007, 14:13 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

Профиль
Группа: Участник
Сообщений: 207
Регистрация: 25.4.2007
Где: Кодинск

Репутация: 1
Всего: 1

MagicPRO,
Блин, знакомое че0то... А как этот метод подругому называется?

--------------------

Я не волшебник... я только учусь...

Maka6er

Дата 27.5.2007, 11:06 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

Профиль
Группа: Участник
Сообщений: 150
Регистрация: 27.2.2007
Где: Киев

Репутация: нет
Всего: 3

Вот сайт мого препода по численным методам
yzorin.ho.com.ua
Только там лекции на Украинском языке, но все примеры есть.

MagicPRO

Дата 28.5.2007, 23:08 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

Профиль
Группа: Участник
Сообщений: 200
Регистрация: 4.10.2006
Где: Kremencity

Репутация: нет
Всего: нет

Misha White, этот метод по другому называется метод статистически-вероятностное моделирование! У меня РГР на эту тему!

--------------------

smile=iwashere

MagicPRO

Дата 1.6.2007, 08:38 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

Профиль
Группа: Участник
Сообщений: 200
Регистрация: 4.10.2006
Где: Kremencity

Репутация: нет
Всего: нет

Вот переделал код с Бэйсика с книги по Численным методам, но он чето матерится

Код


#include<iostream.h>
#include<conio.h>
#include<stdlib.h>
#include<math.h>
int sgn(float);
int main()
{ clrscr();
  cout<<"\t\t\t Monte-Carlo line equatations =)\n\n";
  cout<<" Enter N: "; int N; cin>>N;
  float **A=new float*[N]; float **B=new float*[N];
  float *W=new float[N];
//vvod
  cout<<" Enter coefficients of system:\n";
  for(int i=0;i<N;i++)
    { A[i]=new float[N+1]; B[i]=new float[N];
      cin>>A[i][i]; B[i][i]=fabs(A[i][i]);
      for(int j=1;j<N;j++)

      { cin>>A[i][j]; B[i][j]=B[i][j-i]+fabs(A[i][j]); }
      cin>>A[i][N]; }
  for(i=0;i<N;i++) W[i]=A[i][N]/(1-B[i][N-1]);
////
  cout<<" Enter B: "; int b; cin>>b;
  cout<<" Enter M: "; int m; cin>>m;
  int t=1, y=0, s;
  float c; int V;
M: s=b; y=1;
if (t>m)  goto P;
X: c=random(101)*.01;
   for(int j=N-1;j>=0;j--)
     if(c>B[s][j])
       { if(j==N-1) { t++; y+=V*W[s]; goto M; }
       else { V=V*sgn(A[s][j+1]); s=j+1; goto X;}
       }
  V=V*sgn(A[s][0]); s=0; goto X;
  float x=y/m;
P:  cout<<"Koren = "<<x;
  getch();
  for(i=0;i<=N;i++) delete[] A[i];
  for(i=0;i<N;i++) delete[] B[i];
  delete[] A; delete[] W;
  return 1;
}

int sgn(float t)
{ if(t==0) return 0;
  if(t>0) return 1;
  return -1; }

может кто-то знает в чем проблема????

--------------------

smile=iwashere

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »