Надо выбрать компоненты из списка, указать их в файле и запустить. Там в README.txt написано. Главное формат фортрановский соблюдать - то есть если I5 - то это число должно занять вместе с пробелами пять символов.
Зарегистрирован: 14.04.2005 Сообщения: 1278 Откуда: Россия
Добавлено: 04 Мар 2006 [08:40] Заголовок сообщения:
Вот мой алгоритм для расчета орбиты многоступеньчатой ракеты. Параметры установлены для нашей сверхдешевой ракеты - конкурента Falcon-1
Я даже попытался учесть сопротивление воздуха. Правда это справедливо до 300 м/сек. При сверхзвуковой скорости сопротивление воздуха рассчитывается по формуле k*x - где x-число Маха, k- коэффициент, зависящий от высоты. Желающие могут добавить эту формулу и посмотреть,что выйдет.
Цитата:
// Traektr0View.cpp : implementation of the CTraektr0View class
// для MFC
#include "cmath"
/////////////////////////////////////////////////////////////////////////////
double gm0;
double a;
double R;//Радиус орбиты
double R0;//Радиус Земли
double vy;
double vx;
double dvy;
double dvx;
double v;//скорость ракеты
double alfa;//угол полета
double dl;
double dt;
double dR;
double dbeta;
double beta;//угловой путь ракеты
double x;
double y;
double c;
double pi;
int i0a;
int j0a;
int fl;
char znm2[20]={0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
double mn;//масса ракеты
double m1;//масса газов, выбрасываемых за секунду
double v1;//удельный импульс
double kx1;
double kx2;
double gama;//угол положения ракеты
double vn;
double v2;
double nu;//вязкость воздуха
double r;//лобовой радиус ракеты
CTraektr0View::CTraektr0View()
{
// TODO: add construction code here
R=6371000;R0=6371000;
c=6371000/150;
vx=0;dvx=0;dvy=0;vy=0;dt=0.01;
v1=2620;
fl=100;
r=3.2;
alfa=0;pi=3.141592;
gm0=6.672*5.978*10000000000000;
v=0;//sqrt(pow(vy,2)+pow(vx,2));
}
void CTraektr0View::OnDraw(CDC* pDC)
{
CTraektr0Doc* pDoc = GetDocument();
ASSERT_VALID(pDoc);
// TODO: add draw code for native data here
CClientDC ClientDC(this);
//if(fl==10)goto a1;
a=gm0/pow(R,2);
if(alfa<pi)alfa=acos(vx/v);
else asin(vy/v);
dl=vx*dt;
dR=vy*dt;
R=sqrt(pow((R+dR),2)+pow(dl,2));
dbeta=asin(dl/R);
Вы не можете начинать темы Вы не можете отвечать на сообщения Вы не можете редактировать свои сообщения Вы не можете удалять свои сообщения Вы не можете голосовать в опросах