Форум программистов [Powered by Invision Power Board]

Код

#include<iostream>
#include<math.h>
#include<fstream>
double a=16,b=0.54,c=5,u1=5,u2=65; 
float uin[500],uout[500],t,dt,uid=0,uod=0;
int i=0;
using namespace std;

void peremennie()
    {
    cout<<"\n\n\nОбратите внимания на формулы!\n";
    cout<<"Пока 0<=tu<1 расчёт происходит по формуле u=a*(1-exp(-b*ttemp)), когда 1<=t<2 по u=a*(1-exp(-b))*exp(-c*(ttemp-1))\n";
    cout<<"* U - напряжение; t - время; a,b,c=константы\n\n\n";
        cout<<"Введите a\n";
        cin>>a;
        cout<<"Введите b\n";
        cin>>b;
        cout<<"Введите c\n";
        cin>>c;
        cout<<"Введите u1\n";
        cin>>u1;
        cout<<"Введите u2\n";
        cin>>u2;
}
void standart()
    {    cout<<"\n\n\nОбратите внимание на формулы!\n";
        cout<<"Пока 0<=tu<1 расчёт происходит по формуле u=a*(1-exp(-b*ttemp)), когда 1<=t<2 по u=a*(1-exp(-b))*exp(-c*(ttemp-1))\n";
        cout<<"* U - напряжение; t - время; a,b,c=константы\n\n\n";
        cout<<"a = "<<a<<"\nb = "<<b<<"\nc = "<<c<<"\nu1 = "<<u1<<"\nu2= "<<u2<<"\n\n\n";
    }

float nachf(float ttemp)
  {
  float u;
  if(ttemp<=1) u=a*(1-exp(-b*ttemp));
  else u=a*(1-exp(-b))*exp(-c*(ttemp-1));
  return u;
 } 

void rash(int n)
  {float ud;
   i=0;uid=0;uod=0;ud=0;
   dt=2./(n-1);
   for(t=0;t<=2;t+=dt)
    {
         uin[i]=nachf(t);
          uid+=uin[i]*uin[i];    
        if (uin[i]<=u1)
              uout[i]=(u2/u1*uin[i]);
        else 
            uout[i]=u2;
          uod+=uout[i]*uout[i];
          i++;
    }
   uid=sqrt(uid/(n-1));
   uod=sqrt(uod/(n-1));
  }

void stemp(int n)
  {
   t=0;dt=2./(n-1);
   for(i=0;i<n;i++)
     {
      cout<<"Время:"<<t<<"c  "<<"Напряжение на входе= "<<uin[i]<<"В  "<<"Напряжение на выходе= "<<uout[i]<<"В  \n";t+=dt;
     }
    cout<<"Действующее входящее напряжение:"<<uid<<"B\n"<<"Действующее выходящее напряжение:"<<uod<<"B\n";
  }

void save(int n)
  {
   ofstream fd("Report");
   t=0;dt=2./(n-1);
   fd<<"Действующее входящее напряжение:"<<uid<<"B\n"<<"Действующее выходящее напряжение:"<<uod<<"B\n";
   for(i=0;i<n;i++)
     {
      fd<<"Время:"<<t<<"c  "<<"Uвх= "<<uin[i]<<"В  "<<"Uвых= "<<uout[i]<<"В  \n";
      t+=dt;
     }
  }
void start()
  {
     cout<<"\n\n\nДобро пожаловать!\n1. Просмотреть текущие значения для расчётов\n2. Изменить значения\n3. Произвести расчёты\n4. Просмотреть текущие результаты расчётов\n5. Записать расчёты в файл\n6. Выход из программы\n";
  }

int main()
{
 int n,temp;
 char quit;
 start(i);
 while((quit=getchar())!='6')
  {
   switch(quit)
 {    
    case '1':
        {
        standart();
        break;}    
    case '2':
        {
                peremennie();
                break;
                }
       case '3':{
        cout<<"Сейчас расчёт производится со следущими значениями:";
        rash();
                cout<<"Введите количество точек (от 3 до 500):";cin>>n;
                if(n<=500&&n>=3)
                 {
                  temp=n;
                  Calc(n);
                 }
                else{cout<<"Вы ввели неверное значение\n";n=temp;}
                break;
               }
      case '4':{
                stemp(n);
                break;
               }
      case '5':{
                save(i);
        cout<<"Сохранение прошло успешно! \n\n\n";
                break;
               }
      default:start();
    }
  }
}