Thermometer mit NTC: Berechnung der Temperatur funktioniert nicht

Sebastian D. (--sebi--) 24.10.2010 15:47 #1909893

Lesenswert?

•

Ich möchte ein Thermometer mit einem NTC bauen.
Aber die Berechnung der Temperatur funktioniert nicht. Ich möchte den 
ADC Wert ermitteln(klappt), Die Temperatur daraus ausrechnen (klappt 
nicht), und die Temperatur per RS232 an meinen Laptop schicken(klappt).
Bei meinem Laptop kommt immer nur "25.01" oder "-273,15" an.
Am ADC ist der NTC und ein 10k Pullup angeschlossen.
Das ist der Code:

#define F_CPU 1000000
#include <stdio.h>
//#include <avr/interrupt.h>
//#include "General.h"
#include "RS232.c"
#include <util/delay.h>
#include <math.h>
//#include "RS232.h"
//#include "Delay.h"
uint16_t ADC_Read(uint8_t channel)
    ADMUX = (ADMUX & ~(0x1F)) | (channel & 0x1F);
  ADCSRA |= (1<<ADSC);            // eine Wandlung "single conversion"
  while (ADCSRA & (1<<ADSC) );    // auf Abschluss der Konvertierung warten
    return ADCW;                    // ADC auslesen und zurückgeben
int main (void)
    RS232_Init ();
    ADMUX = (1<<REFS0);      // VCC als Referenz
      ADCSRA = (1<<ADPS1) | (1<<ADPS0);     // Frequenzvorteiler
      ADCSRA |= (1<<ADEN);                  // ADC aktivieren
    printf ("Hello world...\n");
    char itoabuf[20];
    char dtostrfbuf[20];
    int raw;
    float temp;
    while (1)
        _delay_ms(100);    
        _delay_ms(100);
        _delay_ms(100);
        _delay_ms(100);
        _delay_ms(100);
        _delay_ms(100);
        _delay_ms(100);
        _delay_ms(100);
        _delay_ms(100);
        _delay_ms(100);
        raw=ADC_Read(0);
        temp=(1282088/(4300+ln(raw/(1023-raw))*298.16))-273.15;
 //Hier wird die Temperatur berechnet
        dtostrf(temp,3,1,dtostrfbuf);
 //float zu char[]
        printf(dtostrfbuf);
        printf("\n");
        itoa(raw,itoabuf,10);
        printf(itoabuf);
        printf("\n");

Antwort Zitieren Bearbeiten Thread verschieben Thread sperren Anmeldepflicht aktivieren Anpinnen Thread löschen Thread mit anderem zusammenführen Beitrag melden

eProfi Gast 24.10.2010 15:55 #1909905

Lesenswert?

•

▲
▼

Muss man, damit ln funktioniert, nicht noch etwas anderes als math.h 
einbinden / einstellen ?

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Jemand Gast 24.10.2010 16:43 #1909980

Lesenswert?

•

▲
▼

Ist "ln" nicht der Logarithmus und über welchen Temperaturbereich soll 
die Messung sein ?
Welche Werte nimmt "raw" an, also bevor es in die Temperatur umgerechnet 
wird.
Einfacher als diese Floatberechnung wäre erstmal eine Tabelle, die für 
einen bestimmten ADC-Wert den entsprechenden Temperaturwert hat.

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Sebastian D. (--sebi--) 24.10.2010 16:58 #1910004

Lesenswert?

•

▲
▼

raw nimmt den richtigen Wert an. z.B. 528. Wenn ich darauf die 
Temperatur mit dem Taschenrechner berechne, komme ich auch auf die 
richtige Temperatur.

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Stefan W. (wswbln) 24.10.2010 17:24 #1910037

Lesenswert?

•

▲
▼

Bist Du sicher, dass die Lib eine funktionierende ln()-Funktion hat?

Ich habe das damals so gemacht und es funktionietrte:

#define V_ref  5.0      // 5V
// NTC-thermistor defines
#define Beta  3977.000    // the Beta-value for the NTC sensor
#define R_th  10000.00    // the thermistor resistance @ 25 degree Celsius
#define T_amb  298.0000    // the temperature in Kelvin for 25 degree Celsius
#define T_zero  273.0000    // the temperature in Kelvin for 0 degree Celsius
    ADCSRA |= (1<<ADSC);        //Start conversion
  while (ADCSRA & (1<<ADSC)) {}    //wait until converstion completed
    V_ADC = (ADC * V_ref)/1024;        // calculate the voltage over the NTC
  Temperature = (Beta / ((log(V_ADC/(V_ref-V_ADC))/log(exp(1))) + (Beta/T_amb))) - T_zero;

Antwort Zitieren Bearbeiten Löschen Beitrag melden

MaWin Gast 24.10.2010 17:30 #1910041

Lesenswert?

•

▲
▼

> temp=(1282088/(4300+ln(raw/(1023-raw))*298.16))-273.15;


Schreib 1023.0 und überlege mal wirklich ob floating point so schlau auf 
einem uC ist.

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Pete K. (pete77) 24.10.2010 19:42 #1910198

Lesenswert?

•

▲
▼

Hast Du auch dem Linker gesagt, dass er die Floating-Point Bibliothek 
mit einbinden soll?



Sebastian Deppisch schrieb:
> #include "RS232.c"

Btw., C-Programme werden nicht included

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Peter Gast 25.10.2010 12:05 #1910859

Lesenswert?

•

▲
▼

Moin,

ich kenn mich nich so genau mit dem AVR aus, aber die math.h lib hat 
meines Wissens keine "ln" Funktion, sondern nur eine "log" bzw. "log10" 
Funktion. Das entspricht auch so in etwa dem Original der Formel aus der 
Lego-NXT Sammlung. Also mal ausprobieren:
=(1282088/(4300+LOG(RAW/(1023-RAW))*298,16))-273,15

ABER, hast Du auch den gleichen NTC-Widerstand ??

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Karl H. (kbuchegg) 25.10.2010 12:11 #1910871

Lesenswert?

•

▲
▼

Peter schrieb:
> Moin,
>
> ich kenn mich nich so genau mit dem AVR aus, aber die math.h lib hat
> meines Wissens keine "ln" Funktion, sondern nur eine "log" bzw. "log10"

log IST der ln

  raw / (1023-raw)


einfach mal mit Ganzzahlarithmetik ausprobieren, was da rauskommt.
Oder um es so zu sagen: Wenn er auf dem Taschenrechner diesen 
Teilausdruck hier ausrechnet, dann nach der Division alle Kommastellen 
wegwerfen. Dann sieht sein Ergebnis auf dem Taschenrechner nämlich auch 
plötzlich ganz anders aus. Zb wird alles für raw kleiner 512 zu einer 
schnöden 0 (was für den ln dann wieder gar nicht gut ist)

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Peter Gast 25.10.2010 12:29 #1910894

Lesenswert?

•

▲
▼

Karl heinz Buchegger schrieb:
> log IST der ln

Danke für den Hinweis. Der OP hat aber "ln" in seiner Berechnung. Das 
sollte doch eher schief gehen, oder ?
Der Hinweis mit der Ganzzahlarithmetik ist deshalb interessant, weil 
jetzt auch klar wird, warum "MaWin" vorgeschlagen hat, mit 1023.0 zu 
rechnen. Dadurch wird der gesamte Ausdruck eine Floatberechnung. Oder 
sehe ich das falsch ?

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Karl H. (kbuchegg) 25.10.2010 12:30 #1910897

Lesenswert?

•

▲
▼

Peter schrieb:
> Karl heinz Buchegger schrieb:
>> log IST der ln
>
> Danke für den Hinweis. Der OP hat aber "ln" in seiner Berechnung. Das
> sollte doch eher schief gehen, oder ?


Nochmal:
In Programmiersprachen ist die Funktion log der Logarithmus zur Basis e 
(vulgo der 'natürliche Logarithmus', vulgo ln). Der Logarithmus zur 
Basis 10 schreibt sich meistens log10.


> jetzt auch klar wird, warum "MaWin" vorgeschlagen hat, mit 1023.0 zu
> rechnen. Dadurch wird der gesamte Ausdruck eine Floatberechnung. Oder
> sehe ich das falsch ?

Das siehst du richtig.

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Sebastian D. (--sebi--) 25.10.2010 13:22 #1910971

Lesenswert?

•

▲
▼

Danke. Ich werde es mal mit X.0 probieren

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Karl H. (kbuchegg) 25.10.2010 13:23 #1910973

Lesenswert?

•

▲
▼

Sebastian Deppisch schrieb:
> Danke. Ich werde es mal mit X.0 probieren

Und das hier lesen

http://www.mikrocontroller.net/articles/FAQ#Datentypen_in_Operationen

und darüber nachdenken, was da alles daraus folgt.

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Sebastian D. (--sebi--) 30.11.2010 17:46 #1957241

Lesenswert?

•

▲
▼

Also es funktioniert jetzt.
Vielen Dank.
Und µC Programmieren klappt bei mir jetzt auch unter Linux.
@Peter
Woher hast du die Formel?
Die ist doch von meiner Homepage, und die Formel hab ich so erstellt.
http://mylegopage.my.funpic.de/index.php?page=selbstbau&sub=ntc-nxt
Und da es meine ist, hab ich sogar nicht nur den gleiche, sondern auch 
den selben NTC ;)

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Reinhard Kern Gast 01.12.2010 00:56 #1957764

Lesenswert?

•

▲
▼

Sebastian Deppisch schrieb:
> Woher hast du die Formel?
> Die ist doch von meiner Homepage, und die Formel hab ich so erstellt.

Was immer du damit sagen willst - es gibt eine Formel für das Verhalten 
eines NTC, die steht in Wikipedia und in jedem einschlägigen Datenblatt. 
Die muss man für die Anzeige nach T auflösen, und wie das halt in der 
Mathematik so ist, kommt dabei IMMER das Gleiche heraus. U.a. deswegen 
gibt es auch keine Patente auf mathematische Formeln, jedenfalls bisher 
nicht.

Gruss Reinhard

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Sebastian D. (--sebi--) 01.12.2010 16:59 #1958485

Lesenswert?

•

▲
▼

Das ist mir schon klar, aber 
"=(1282088/(4300+LOG(RAW/(1023-RAW))*298,16))-273,15"(so wie Peter es 
geschrieben hat) Habe ich so auf meiner Homepage stehen.

Antwort Zitieren Bearbeiten Löschen Beitrag melden

Thermometer mit NTC: Berechnung der Temperatur funktioniert nicht

Antwort schreiben