реклама на сайте
подробности

 
 
> Логарифм, Вычесление логарифма
vov136
сообщение Aug 31 2004, 03:24
Сообщение #1





Группа: Новичок
Сообщений: 1
Регистрация: 31-08-04
Пользователь №: 569



Помогите пожалуста можно ли зашить в МК, и в какой, вычисление десятичного логарифма? Заранее благодарен.
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
Starick
сообщение Mar 30 2006, 14:31
Сообщение #2


Участник
*

Группа: Участник
Сообщений: 34
Регистрация: 4-09-05
Пользователь №: 8 211



Цитата(-Tумблер- @ Nov 5 2004, 17:19) *
Цитата(vov136 @ Aug 31 2004, 07:24 AM)
Помогите пожалуста можно ли зашить в МК, и в какой, вычисление десятичного логарифма? Заранее благодарен.


В библиотеках С-компилеров просто должна быть
готовая функция.
:P

Ага...только при использовании этой функции флэш забивается чуть ли не на половину (У меги 8535)!!! Я уже молчу о времени выполнения данной функции.....быстрее было бы на калькуляторе просчитать blink.gif
Так что советую два выхода:
- либо численные методы, как было сказано выше.
- либо ставь сопроцессор, который бы тебе считал логарифм, сооствественно с вытекающими последствиями
Go to the top of the page
 
+Quote Post
defunct
сообщение Mar 31 2006, 01:10
Сообщение #3


кекс
******

Группа: Свой
Сообщений: 3 825
Регистрация: 17-12-05
Из: Киев
Пользователь №: 12 326



Цитата(Starick @ Mar 30 2006, 16:31) *
Ага...только при использовании этой функции флэш забивается чуть ли не на половину (У меги 8535)!!! Я уже молчу о времени выполнения данной функции.....быстрее было бы на калькуляторе просчитать blink.gif
Так что советую два выхода:
- либо численные методы, как было сказано выше.
- либо ставь сопроцессор, который бы тебе считал логарифм, сооствественно с вытекающими последствиями


Что-то не сходится.. Ну не может функция расчета логарифма занимать пол-флеша у mega8535. Одно из двух либо она будет работать очень быстро (если столько памяти занимает) либо вы что-то преврали.
В общем стандартное решение:
#include <math.h>

y = log10(x);
Go to the top of the page
 
+Quote Post
radamateur
сообщение Apr 2 2006, 18:15
Сообщение #4





Группа: Новичок
Сообщений: 3
Регистрация: 4-01-06
Пользователь №: 12 830



Цитата
Что-то не сходится.. Ну не может функция расчета логарифма занимать пол-флеша у mega8535. Одно из двух либо она будет работать очень быстро (если столько памяти занимает) либо вы что-то преврали.
В общем стандартное решение:
#include <math.h>

y = log10(x);


Вообще-то используя строку :

#include <math.h>

кроме логарифмической функции попадает куча ненужной информации,
что приводит к быстрому расходу памяти в микроконтроллере, оссобенно в tiny и classic.
Предлагаю использовать примерно такую функцию для всех x>0 ( вместо log (x) log_x (x),
а вместо log10 (x) log10_x (x) ):

double log_x (double );

double log_x (double x){
int m= ; /* (m-длина ряда, целое число) присвойте m значение , чем больше длина ряда,
тем точнее значение логарифма, но тем дольше идет его расчет */

double b,c,d,e,f;
b=(x-1)/(x+1);
с=b*b;
d=1/b;
f=0;
e=-1;
while (m>0){
d*=c;
e+=2;
f+=(d/e);
m--;
}
f*=2;
return f;
}

А для вычисления десятичного логарифма важно помнить след. соотношение:
lg (x) = ln (x) / ln (10 ) , а точнее это будет выглядеть так :

log10_x (x) = log_x (x)/log_x (10);
Go to the top of the page
 
+Quote Post
radamateur
сообщение Apr 4 2006, 17:13
Сообщение #5





Группа: Новичок
Сообщений: 3
Регистрация: 4-01-06
Пользователь №: 12 830



А вот для десятичного логарифма и к тому-же работаяющая гораздо быстрее.
Для достаточно больших x, требуеся достаточно большая длина ряда. В нижеприведенной функции
x преобразуется к виду 1<x<10, а далее вычисляется десятичный логарифм. Точность L= 10^(-6),
если нужна большая - увеличьте длину ряда m исходя из соотношения :

(9/11)^(2*m-1) < L/(100*(2*m-1))

В данном случае L=10^(-6) т.е. m=36


double log10_x (double );

double log10_x (double x){
int m= 36; // m-длина ряда, целое число

double a,b,c,d,e,f,g,i,j;

if (x==1) return 0;

a=x;
i=1;
if (x<1){
a=1/x;
i=-1;
}

j=0;
while (a>10){
a/=10;
j++;
}

if (a==10){
j++;
return i*j;
}

b=(a-1)/(a+1);
с=b*b;
d=1/b;
f=0;
e=-1;
while (m>0){
d*=c;
e+=2;
g=d/e;
f+=g;
m--;
}
f*=2;
f*=0.43429448; // для точности 10^(-6)
// округлите 1/ln(10) до 7-го знака, для 10^(-7) до 8-го и т.д..
f+=j;
f*=i;

return f;
}
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- vov136   Логарифм   Aug 31 2004, 03:24
- - KRS   Да в любой можно запихать вопрос в точности и ско...   Aug 31 2004, 08:34
- - Himer   Насколько помню из математики логарифм через ряд в...   Aug 31 2004, 09:00
- - Kot   Если много лишней памяти и мало времени - табличн...   Sep 24 2004, 01:18
- - -Tумблер-   Цитата(vov136 @ Aug 31 2004, 07:24 AM)Помогит...   Nov 5 2004, 14:19
- - impuls-v   Флеш забивается потому что ты используеш всю библи...   Apr 2 2006, 11:53
- - AndyBig   ЦитатаФлеш забивается потому что ты используеш всю...   Apr 2 2006, 18:11
- - vitus_strom   Я делал на основе знакопостоянного CORDIC'a та...   Apr 5 2006, 08:16
- - Rst7   А при чем тут CORDIC (быстрый поворот вектора) к л...   Apr 5 2006, 09:27
- - vitus_strom   при том что это семейство алгоритмов, в которое вх...   Apr 5 2006, 11:09


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 2nd August 2025 - 05:52
Рейтинг@Mail.ru


Страница сгенерированна за 0.01403 секунд с 7
ELECTRONIX ©2004-2016