Форум программистов, компьютерный форум, киберфорум
Pascal (Паскаль)
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.80/5: Рейтинг темы: голосов - 5, средняя оценка - 4.80
12 / 8 / 6
Регистрация: 10.07.2011
Сообщений: 374
Записей в блоге: 1

Непонятная ошибка в точности вычисления arccos

31.10.2013, 14:53. Показов 971. Ответов 5
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Написал код, но он вычисляет arccos с погрешностью:
Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
var
    n:integer;
    x, s, t, e: real;
begin
    e := 0.000001;
    Readln(x);
    s := x;
    t := x;
    n := 0;
    if abs (x)  < 1 then begin // -1 <= x <= 1
        while (abs(t) > e ) do begin
      inc(n);
      t := t * sqr(x) * (2*n-1)/((2*n)*(2*n+1));
            s := s + t;
        end;
        s := Pi/2 - s;
        writeln('arccos (', x,') = ', s);
        writeln(Arctan(Sqrt(1-sqr(x))/x)  );
    end else
        writeln('Error'); 
end.
Нашел на форуме практически такой же, но в нем все ок:
Pascal
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
program Project1;
{$APPTYPE CONSOLE}
 
uses
   SysUtils, Math;
 
function MyArcSin(X : double) : double;
const
   eps = 0.00001;
var
   i : integer;
   v, curr : double;
begin
   i := 0;
   Result := x;
   curr := x;
   repeat
      inc(i);
      curr := curr * sqr(x) * (2*i - 1) / 2 / i;
      v := curr / (2*i + 1);
      Result := Result + v;
   until v < eps;
end;
 
const
   eps = 0.000001;
var
   x : double;
begin
   write('X (-1 .. 1) = '); readln(x);
   // Тут можешь проверять, входит ли введенное значение в интервал -1 .. 1
   writeln(MyArcSin(x):10:6);
   writeln(arcsin(x):10:6); // Проверка
   readln;
end.
Где теряется точность в моем коде?
0
cpp_developer
Эксперт
20123 / 5690 / 1417
Регистрация: 09.04.2010
Сообщений: 22,546
Блог
31.10.2013, 14:53
Ответы с готовыми решениями:

Arccos. Ошибка.
Чувствую себя полным идиотом. Вот задача: y= ArcCos (ArcCos x ) - ArcCos x^3. Вот решение: program zadan2; uses crt; var ...

Ошибка при использовании arccos
uses crt; var x : integer; y : real; BEGIN clrscr; writeln('Введите X '); readln(x);

Не точности вычисления C++
double N=1.0100000000000000; N=(N-floor(N))*10; cout&lt;&lt;N; Выводит: 0.10000000000000009 - что ни есть правильно (добавляет 9 в...

5
 Аватар для Новичок
1682 / 1098 / 489
Регистрация: 17.07.2012
Сообщений: 5,360
31.10.2013, 15:20
Цитата Сообщение от Bohes2013 Посмотреть сообщение
погрешностью
А погрешность большая?
0
12 / 8 / 6
Регистрация: 10.07.2011
Сообщений: 374
Записей в блоге: 1
31.10.2013, 15:32  [ТС]
Если запустите мой код, то первой строкой выведется подсчитанное моей программой число, а второй--с помощью функции tan. Во втором же коде, эти числа совпадают

Добавлено через 1 минуту
Вот, конкретный пример: Ввожу 0.8. Вывод: arccos (0.8) = 0.685462993461563
0.643501108793284

Добавлено через 3 минуты
Самое интересное, что формула в обеих программах одна и та же, да и переменные аналогичны
0
696 / 570 / 414
Регистрация: 31.03.2013
Сообщений: 1,029
31.10.2013, 16:05
Bohes2013,
Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
const eps = 0.00001;
function MyArcSin(x:double):double;
var i:integer;
   curr:double;
begin
  i:= 0;
  Result:=x;
  curr:=x;
  repeat
    inc(i);
    curr := curr*x*x*(2*i-1)*(2*i)/(4*i*i);
    Result:= Result + curr/(2*i+1);
  until curr/(2*i + 1)<eps;
end;
 
var x:double;
begin
  repeat
    write('|x|<1, x = '); readln(x);
  until abs(x)<1;
  writeln(MyArcSin(x):10:6);
  writeln(arcsin(x):10:6);
  readln;
end.
0
12 / 8 / 6
Регистрация: 10.07.2011
Сообщений: 374
Записей в блоге: 1
01.11.2013, 02:07  [ТС]
Цитата Сообщение от Samhain Посмотреть сообщение
Bohes2013,
Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
const eps = 0.00001;
function MyArcSin(x:double):double;
var i:integer;
   curr:double;
begin
  i:= 0;
  Result:=x;
  curr:=x;
  repeat
    inc(i);
    curr := curr*x*x*(2*i-1)*(2*i)/(4*i*i);
    Result:= Result + curr/(2*i+1);
  until curr/(2*i + 1)<eps;
end;
 
var x:double;
begin
  repeat
    write('|x|<1, x = '); readln(x);
  until abs(x)<1;
  writeln(MyArcSin(x):10:6);
  writeln(arcsin(x):10:6);
  readln;
end.
Код, который вы модифицировали, уже имеел высокую точность (после модификации она понизилась). Мне нужны были рекомендации или объяснения по исправлению моего кода. Собственно, я нашел в чем тут обстоит дело, поэтому тема закрыта.
0
 Аватар для Новичок
1682 / 1098 / 489
Регистрация: 17.07.2012
Сообщений: 5,360
01.11.2013, 12:49
Bohes2013,а вы не могли бы написать,а то я с рядами особо не дружу.Буду благодареню
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
raxper
Эксперт
30234 / 6612 / 1498
Регистрация: 28.12.2010
Сообщений: 21,154
Блог
01.11.2013, 12:49
Помогаю со студенческими работами здесь

Вычисления до заданной точности
Исправьте ошибку пожалуйста! #include &lt;stdio.h&gt; #include &lt;math.h&gt; #include &lt;conio.h&gt; int fact(int x) { if(x&lt;=1) return 1;...

Определение точности вычисления интеграла
Задание: Определить точность вычисления значения интеграла (ε - ?) взяв сумму первых трех членов разложения: ...

Проверить вычисления на переполнение и потерю точности
Надо написать функции арифметических действий с проверкой переполнения,потери и генерировать исключение если оно случается. Вообщем то...

Различия в точности вычисления в одно/многопоточном (OMP) режиме
Различия в точности, скорости вычисления в одно/многопоточном (OMP) режиме Привет. Стала задача оптимизировать софт и ускорить за...

Вычислить сумму ряда. Вычисления заканчивать при достижении заданной степени точности
Помогите пожалуйста доработать программу, в программировании новичок. #include &quot;stdafx.h&quot; #include &lt;iostream&gt; ...


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

Или воспользуйтесь поиском по форуму:
6
Ответ Создать тему
Новые блоги и статьи
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++
8Observer8 10.03.2026
Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:
Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE
8Observer8 10.03.2026
Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .
Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки
Neotwalker 08.03.2026
Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .
Модульная разработка через nuget packages
DevAlt 07.03.2026
Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .
Модульный подход на примере F#
DevAlt 06.03.2026
В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .
Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование
8Observer8 05.03.2026
Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .
SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D
8Observer8 04.03.2026
Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .
SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий
8Observer8 02.03.2026
Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru