С Новым годом! Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
[O]Clic[K]
1 / 1 / 0
Регистрация: 28.03.2012
Сообщений: 55
#1

Вычисление косинуса разложением в ряд (C++ и Asm) - C++

06.06.2014, 02:20. Просмотров 614. Ответов 10
Метки нет (Все метки)

В общем такая ситуация:
Нужно посчитать cos используя разложение в Ряд
Вот мои наброски
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
#include <stdio.h>
#include <math.h>
 
extern "C"
{ 
    float _asmtest(float) ;
}
 
 
int main() 
{
    float x=0;  
    printf("x: ");
    scanf("%f",&x);
    
    printf("C cos = %f\n", cos(x));
    printf("My cos = %f\n", _asmtest(x));
    
    return 0;
}
Сам же косинус
Assembler
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
36
37
38
39
40
41
42
.386
.model flat, c
.code
PUBLIC _asmtest
_asmtest PROC near source:dword
    
    fld1
    fld source 
    fld1
    
ls: fmul st[0],st[1]            
    fmul st[0],st[1]
    inc cntr
    fidiv cntr
    inc cntr
    fidiv cntr
    fchs
    fadd st[2],st[0]
    fchs
 
    fmul st,st[1]
    fmul st,st[1]
    inc cntr
    fidiv cntr
    inc cntr
    fidiv cntr
    fadd st[2],st[0]
    
    cmp cntr,10
    jb ls
    
    fstp temp
    fstp temp
                
   
                
ret
_asmtest endp 
.data
temp    dq  0
cntr    dw  0
end
Упорно не могу понять, где у меня косяк
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
06.06.2014, 02:20
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Вычисление косинуса разложением в ряд (C++ и Asm) (C++):

Вычисление косинуса через ряд - C++
Здравствуйте =) Пытаюсь посчитать косинус через ряд Тейлора, написал программку, но не могу понять, как она считает (результаты не сходятся...

Вычисление ln(1 + x) разложением в ряд - C++
Всем доброго время суток бьюсь никак не могу въехать как написать программу вычисления ln(1 + x) разложением в ряд в книге дана формула: ...

Вычисление ch(x) разложением в ряд - C++
Ребят, не получается решить задание, помогите, пожалуйста. Задачу и формулы надо переписывать! Редактор формул внизу страницы. В...

Вычисление функций разложением в ряд - C++
Доброго времени суток ))) Прошу помочь начинающему пользователю С++ :-[ Я бы хотела разобраться как выполнить данное задание в коде С++...

Вычисление функции разложением в ряд - C++
Доброго времени суток,народ...Помогите найти ошибку:help: #include &lt;iostream&gt; #include &lt;clocale&gt; #include &lt;cmath&gt; using...

Вычисление функции,разложением ее в ряд - C++
Помогите, пожалуйста, в решении задачи. НИЧЕГО НЕ ПОЛУЧАЕТСЯ Написать программу нахождения суммы ряда с заданной точностью ε. ...

10
castaway
Эксперт С++
4916 / 3024 / 370
Регистрация: 10.11.2010
Сообщений: 11,081
Записей в блоге: 10
Завершенные тесты: 1
06.06.2014, 08:46 #2
А формула разложения в ряд как выглядит?
0
[O]Clic[K]
1 / 1 / 0
Регистрация: 28.03.2012
Сообщений: 55
06.06.2014, 10:16  [ТС] #3
Вычисление косинуса разложением в ряд (C++ и Asm)
0
castaway
Эксперт С++
4916 / 3024 / 370
Регистрация: 10.11.2010
Сообщений: 11,081
Записей в блоге: 10
Завершенные тесты: 1
06.06.2014, 11:14 #4
cntr надо обнулять каждый раз при вызове функции, но это ладно...
После 4 идёт 6 или 8?
Почему степень всегда одна считается?
0
[O]Clic[K]
1 / 1 / 0
Регистрация: 28.03.2012
Сообщений: 55
06.06.2014, 15:45  [ТС] #5
После 4 идёт 6

Добавлено через 4 часа 23 минуты
В общем, в промежутке от -3 до 6 косинус считает нормально, после же числа резко увеличиваются в размерах, с чем это связанно и как это исправить?
0
azbest
41 / 41 / 8
Регистрация: 12.03.2013
Сообщений: 148
06.06.2014, 15:57 #6
это известная проблема компьютерных вычислений в "плохо обусловленых задачах", просто для определения значения функции убирай период 2*pi, значение для которого считаеш должно быть между [ 0..2*pi )
0
[O]Clic[K]
1 / 1 / 0
Регистрация: 28.03.2012
Сообщений: 55
06.06.2014, 16:13  [ТС] #7
Как я понял период надо убирать перед тем как передать значение ассемблеру, но что-то не догоню как его убрать
0
azbest
41 / 41 / 8
Регистрация: 12.03.2013
Сообщений: 148
06.06.2014, 16:26 #8
да до передачи "упростить"
C++
1
2
3
4
5
#define DPI 4*acos(0.0)
...............
x-=((int)(x/DPI))*DPI;
printf("C cos = %f\n", cos(x));
......................
Добавлено через 5 минут
правда такой вариан даст диапазон (-2*pi; 2*pi]
можно еще поставить проверку после "упрощения"
C++
1
    if (x<0) x+=DPI;
тогда уже наверняка.
0
[O]Clic[K]
1 / 1 / 0
Регистрация: 28.03.2012
Сообщений: 55
06.06.2014, 16:34  [ТС] #9
Что-то не помогает...
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
#include <stdio.h>
#include <math.h>
#define DPI 4*acos(0.0)
 
extern "C"
{ 
    float _asmtest(float) ;
}
 
 
int main() 
{
    float x=0;  
    printf("x: ");
    scanf("%f",&x);
    x-=((int)(x/DPI))*DPI;
    if (x<0) x+=DPI;
    printf("C cos = %f\n", cos(x));
    printf("My cos = %f\n", _asmtest(x));
    
    return 0;
}
0
azbest
41 / 41 / 8
Регистрация: 12.03.2013
Сообщений: 148
06.06.2014, 16:42 #10
Если предыдущий вариант считал нормально для (-3 6) тогда проверь какое значение х когда передаеш его в функцию. код я писал не проверяя, можеш что-то забыл))

и еще для точности поменяй float на double
0
[O]Clic[K]
1 / 1 / 0
Регистрация: 28.03.2012
Сообщений: 55
06.06.2014, 18:16  [ТС] #11
Короче как-то не так от периода надо освободиться, но учитывая мои сегодняшние жесткие тупняки и очень посредственное знания программирования, я что-то не догоняю

Добавлено через 18 минут
все, вроде бы сообразил и вроде бы даже работает, работает правда немного не так
Использовал
C++
1
x = modf(DPI, &x)
Добавлено через 55 минут
Наконец то прога работает более-менее сносно, но почему то для больших цифр в районе +-1млн результат отличается от калькулятора Винды, причем довольно ощутимо, даже знак не совпадает) Может быть кто-нибудь подскажет почему
C++
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
#include <stdio.h>
#include <math.h>
#define DPI 4*acos(0.0)
 
extern "C"
{ 
    float _asmtest(float) ;
}
 
 
int main() 
{
    float x=0;
    float i=0;
    printf("x: ");
    scanf("%f",&x);
    if (x>0){
    while (fabs(x)>=DPI){
            x=x-DPI;}
    }
    if (x<0) {
        while (fabs(x)>=DPI){
        x=x+DPI;}
    }
    //x = x-modf(DPI, &x);
    //printf("x_new: %f\n", x);
    printf("C cos = %f\n", cos(x));
    printf("My cos = %f\n", _asmtest(x));
    
    return 0;
}
0
06.06.2014, 18:16
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
06.06.2014, 18:16
Привет! Вот еще темы с ответами:

Вычисление arcsin разложением в ряд - C++
здравствуйте,помогите пожалуйста с темой Вычислить с заданной точностью значение функции , используя ее разложение в ряд:

Вычисление функции разложением в ряд - C++
Здравствуйте, помогите написать программу пожалуйста. лекцию прослушал и не понимаю как. мой вариант на 1й фотографии нужно сделать по...

Вычисление функции разложением в ряд - C++
1. Составить алгоритм и написать программу вычисления следующих сумм. 2. Составить программу вычисления суммы первых n=10; 20; 30 членов...

Вычисление Sinx/x разложением в ряд - C++
Всем доброго время суток. Вторые сутки бьюсь, задали написать на C++ программу вычисление Sinx/x разложением в ряд: Sinx/x = 1 – (x^2)/3!...


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

Или воспользуйтесь поиском по форуму:
11
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.