VPF::Поправьте/допишите код! - Форум программистов

Tep4k111

Дата 3.12.2013, 11:16 (ссылка)

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

Загрузка ...

Новичок

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

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

Собственно вот условие задачи:
Составить программу и вычислить сумму ряда с точностью e=pow(10,-3), e=pow(10,-6), e=pow(10,-9). Сравнить число итераций для каждого из этих случаев.

user posted image

Вот мой код:

Код

#include <iostream.h>
#include <conio.h>
#include <math.h>
#include <stdio.h>
void main ()
{
float k=0,n,S1,x,S,p,f=1;
printf("Enter x= ");
scanf("%f", &x);
S=0.0;
n=1;   //nachal`noe dlya vichisleniya ryada chlena
while (S>e)
{
for (k=0; k<=n; k++)
S+=pow(1.0,k)*(-1.0)*sin(2.0*k+1.0)*x/2.0*k+1.0; //сама фу-я
n++;
}
cout<<"S="<<S<<endl;
f=1;
while (S>=0)
{
f*=-x*x*(2*k+1)/(2*k*k+5*k+1); // формула Dalamber для вычисления итераций
if (f>e) cout<<"The end"<<endl;
else S+=f;
k++;
}
cout<<"S="<<S<<endl;
cout<<"chislo iteraciy ="<<k;
p=3.14;
if ((x>0.0)&&(x<p/2.0))
S1=0.5*log(tan(x/2.0+p/4.0));
cout<<S1<<endl;
if (S=S1) cout<<"true"<<endl; //сравнение
else cout<<"banan";
getch();
clrscr();
}

Модератор: не забываем пользоваться кнопочкой "Код"

Это сообщение отредактировал(а) bsa - 3.12.2013, 15:18

smalcom

Дата 3.12.2013, 12:53 (ссылка)

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

Загрузка ...

Опытный

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

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

1. используйте тег "код". мусорник же
2. в чём проблема?

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

What DA F...???

feodorv

Дата 3.12.2013, 15:46 (ссылка)

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

Загрузка ...

Эксперт

Профиль
Группа: Комодератор
Сообщений: 2214
Регистрация: 30.7.2011

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

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

if (S=S1) cout<<"true"<<endl; //сравнение
else cout<<"banan";

Числа с плавающей точкой так сравнивать не стоит. У Вас заранее задана точность вычисления, поэтому сравнение должно выглядеть так:

Код


if( fabs( S - S1 ) <= epsilon ) ...

хотя, конечно, не понятно, какая точность задана - относительная или абсолютная...

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

n=1; //nachal`noe dlya vichisleniya ryada chlena

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

for (k=0; k<=n; k++)

Ну и всегда имеете 2 итерации, при которых всё в S складывается и складывается...

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

S+=pow(1.0,k)*(-1.0)*sin(2.0*k+1.0)*x/2.0*k+1.0; //сама фу-я

Смысла возводить 1.0 в степень k никакого нет (а потом умножать на -1), тем более в формуле стоит -1 в степени k (то есть это знакопеременный ряд).
/2.0*k+1.0 компилятором воспринимается несколько иначе, чем Вам хотелось бы, выражение нужно взять в скобки: (2.0*k+1.0)

При чём здесь Даламбер, наверное, не скажет никто, судя по всему, это из какого-то другого задания.

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

while (S>e)

e не только не задана, но даже не прописана в числе переменных.

Никакого сравнения числа итераций нет...

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

float k=0,n

Зачем счётчику цикла быть float???

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

Напильник, велосипед, грабли и костыли - основные инструменты программиста...

borisbn

Дата 3.12.2013, 16:34 (ссылка)

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

Загрузка ...

Эксперт

Профиль
Группа: Завсегдатай
Сообщений: 4875
Регистрация: 6.2.2010
Где: Ростов-на-Дону

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

Цитата(feodorv @ 3.12.2013, 15:46

)

if (S=S1) cout<<"true"<<endl; //сравнение

Это даже не сравнение ))

Tep4k111, пока не отформатируешь по-человечески код, тебе вряд-ли кто-то захочет помогать.

feodorv, пожалуй, исключение, подтверждающее правило smile

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

Женщины отличаются от программистов тем, что у них чары состоят из стрингов

Tep4k111

Дата 3.12.2013, 16:45 (ссылка)

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

Загрузка ...

Новичок

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

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

Цитата(feodorv @ 3.12.2013, 15:46)

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

if (S=S1) cout<<"true"<<endl; //сравнение
else cout<<"banan";

Код


if( fabs( S - S1 ) <= epsilon ) ...

хотя, конечно, не понятно, какая точность задана - относительная или абсолютная...

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

n=1; //nachal`noe dlya vichisleniya ryada chlena

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

for (k=0; k<=n; k++)

Ну и всегда имеете 2 итерации, при которых всё в S складывается и складывается...

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

S+=pow(1.0,k)*(-1.0)*sin(2.0*k+1.0)*x/2.0*k+1.0; //сама фу-я

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

while (S>e)

e не только не задана, но даже не прописана в числе переменных.

Никакого сравнения числа итераций нет...

Цитата(Tep4k111 @ 3.12.2013, 12:16

)

float k=0,n

Зачем счётчику цикла быть float???

спасибо за ценный комментарий, но меня всё так интересует, можно ли как-то показать эти итерации? Насчёт float: там целочисленные значения должны быть, мой косяк. Как я понял, Даламбер там вообще не нужен. Переменная "e" была задана через double. Нужно ли так делать, либо лучше вводить? Тогда как ввести степень?

baldina

Дата 3.12.2013, 17:26 (ссылка)

(голосов:1)

Загрузка ...

Эксперт

Профиль
Группа: Завсегдатай
Сообщений: 3433
Регистрация: 5.12.2007
Где: Москва

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

1. надо оформить функцию в виде функции, так удобнее. заодно исправить её, т.к. у вас функция не соответствует картинке

Код


double Series (int k, double x) {
 double arg = 2.0*k+1.0;
 if (k%2) // нечетная степень -1
   x = -x;
 return sin(arg)*x/arg;
}

функцию для расчета точного значения - тоже. и pi при этом взять поточнее, например как результат arccos(-1)

Код


const double pi=acos(-1.0);
double f (double x) {
  return 0.5*log(tan(x/2.0+pi/4.0));
}

2. у вас задание посчитать сумму несколько раз, с заданной точностью, и выяснить число итераций
напрашивается решение - расчет суммы оформить в виде функции. разумно также задать максимальное число итераций, что бы программа всегда завершалась.
Если число итераций совпало с максимальным значением, значит требуемая точность не достигнута.

Код


void sum (
  ,double x 
  ,double y // точное значение 
  ,double epsilon // желаемая точность
) 
{
   double S=0;
   int k;
   int m = ... ; // максимальное число итераций
   for (k=0; k < m; ++k) { 
      S += Series (...);
      if (fabs(S-y) < epsilon) {
          // достигли заданной точности
          break;
      }
   }
   if (k == m)
     cout << "Расчет прерван по числу итераций\n".
   cout << k << ' ' << S << ' ' << y << ' ' << epsilon << endl;
}

3. Вот эту функцию и надо вызывать для различных значений epsilon.

Код


double y = f(x);
sum (x, y, 1e-3); 
sum (x, y, 1e-6); 
sum (x, y, 1e-9);

4. Если введенное значение x не лежит в диапазоне ]0..pi/2[, то и ряд считать смысла нет, надо попросить пользователя вернуть правильное значение

Это сообщение отредактировал(а) baldina - 3.12.2013, 17:53

Tep4k111

Дата 4.12.2013, 01:12 (ссылка)

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

Загрузка ...

Новичок

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

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

Цитата(baldina @ 3.12.2013, 17:26)

Код


double Series (int k, double x) {
 double arg = 2.0*k+1.0;
 if (k%2) // нечетная степень -1
   x = -x;
 return sin(arg)*x/arg;
}

функцию для расчета точного значения - тоже. и pi при этом взять поточнее, например как результат arccos(-1)

Код


const double pi=acos(-1.0);
double f (double x) {
  return 0.5*log(tan(x/2.0+pi/4.0));
}

Код


void sum (
  ,double x 
  ,double y // точное значение 
  ,double epsilon // желаемая точность
) 
{
   double S=0;
   int k;
   int m = ... ; // максимальное число итераций
   for (k=0; k < m; ++k) { 
      S += Series (...);
      if (fabs(S-y) < epsilon) {
          // достигли заданной точности
          break;
      }
   }
   if (k == m)
     cout << "Расчет прерван по числу итераций\n".
   cout << k << ' ' << S << ' ' << y << ' ' << epsilon << endl;
}

3. Вот эту функцию и надо вызывать для различных значений epsilon.

Код


double y = f(x);
sum (x, y, 1e-3); 
sum (x, y, 1e-6); 
sum (x, y, 1e-9);

4. Если введенное значение x не лежит в диапазоне ]0..pi/2[, то и ряд считать смысла нет, надо попросить пользователя вернуть правильное значение

спасибо!

Tep4k111

Дата 6.12.2013, 22:49 (ссылка)

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

Загрузка ...

Новичок

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

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

Цитата(Tep4k111 @ 4.12.2013, 01:12)

Цитата(baldina @ 3.12.2013, 17:26)

Код


double Series (int k, double x) {
 double arg = 2.0*k+1.0;
 if (k%2) // нечетная степень -1
   x = -x;
 return sin(arg)*x/arg;
}

функцию для расчета точного значения - тоже. и pi при этом взять поточнее, например как результат arccos(-1)

Код


const double pi=acos(-1.0);
double f (double x) {
  return 0.5*log(tan(x/2.0+pi/4.0));
}

Код


void sum (
  ,double x 
  ,double y // точное значение 
  ,double epsilon // желаемая точность
) 
{
   double S=0;
   int k;
   int m = ... ; // максимальное число итераций
   for (k=0; k < m; ++k) { 
      S += Series (...);
      if (fabs(S-y) < epsilon) {
          // достигли заданной точности
          break;
      }
   }
   if (k == m)
     cout << "Расчет прерван по числу итераций\n".
   cout << k << ' ' << S << ' ' << y << ' ' << epsilon << endl;
}

3. Вот эту функцию и надо вызывать для различных значений epsilon.

Код


double y = f(x);
sum (x, y, 1e-3); 
sum (x, y, 1e-6); 
sum (x, y, 1e-9);

спасибо!

Вот что получилось... укажите на недочёты.

Код


#include <iostream.h>
#include <conio.h>
#include <math.h>
#include <stdio.h>
void main ()
{
int x;
cout<<"Enter the number x =";
cin>>x;
int iter=0;
 int iter2=0;
int iter3=0;  
 double S=0.0, arg, F;
double S2=0.0;
 double S3=0.0;
double pi=3.14;
 double e=pow(10,-3);
double e2=pow(10,-6);
 double e3=pow(10,-9);
double F4=0.5*log(tan(x/2.0+pi/4.0));

      for (int k=0; fabs(S-F4)>e; k++)
     {
      arg=2.0*k+1.0;
      double F=pow(-1.0,k)*sin(arg)*x/arg;
      S+=F;
      iter+=1;
     }

cout<<"\niter ="<<iter<<"\nS ="<<S<<"\nF2="<<F4;

     for (int k2=0; fabs(S2-F4)>e2; k2++)
    {
     arg=2.0*k2+1.0;
     double F2=pow(-1.0,k2)*sin(arg)*x/arg;
     S2+=F2;
     iter2+=1;
    }

cout<<"\niter2 ="<<iter2<<"\nS2 ="<<S2<<"\nF2="<<F4;

    for (int k3=0; fabs(S3-F4)>e3; k3++)
       {
    arg=2.0*k3+1.0;
    double F3=pow(-1.0,k3)*sin(arg)*x/arg;
    S3+=F3;
    iter3+=1;
       }

cout<<"iter3 ="<<iter3<<"\nS3 ="<<S3<<"\nF2="<<F4;
getch();
clrscr();
}

baldina

Дата 6.12.2013, 23:31 (ссылка)

(голосов:2)

Загрузка ...

Эксперт

Профиль
Группа: Завсегдатай
Сообщений: 3433
Регистрация: 5.12.2007
Где: Москва

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

код не отформатирован, вырви глаз
контроля числа итераций нет, возможно зацикливание (ряд кстати сходится медленно)
бессмысленно считаете с точностью 1e-9, задавая pi с точностью 1e-2
зачем числа 1e-3, 1e-6, 1e-9 вычислять при помощи pow()?
переменные iter1..3 лишние, у вас есть k
контроль значения х отсутствует

про copy-paste я уж не говорю, видимо функции С++ вы ниасилили

Добавлено через 59 секунд
ну вобщем на троечку, если у преподавателя настроение хорошее smile

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