IIR Filter mit Mega 16

[Daniel B. (yzdani)]

Hallo,

möchte ein IIR Filter "basteln".
Fpass 100Hz
Fs=1Khz

Der Wert nach dem Filte soll danach auf Port C ausgegeben werden
Umso mehr ich die Frequenz erhöhe sollten dann immer ein kleinerer binär 
wert am Port C anliegen.
Nun ja bei einem Hz blinken die Leds am Port C aber wenn ich die 
Frequenz erhöhe leuchten die Leds dann immer gleich.

Der Sinus hat 1V offset und 1V Amplitude
ADC funktioniert

Was mach ich falsch?
Hat mir jemand ein Lösungsansatz?

Code:

int main(void)
{
  float a[3], b[3], w[3]={0,0,0};
  //koeffizienten
  a[0]=1;
  a[1]=-1.9961;
  a[2]=0.5157;
  b[0]=1;
  b[1]=1.9961;
  b[2]=1;
  float y=0;
  float sample;
  DDRA=0x00;
  DDRC=0xff;
  ADMUX|=(1<<REFS0);
  ADCSRA|=(1<<ADEN)|(1<<ADSC);

    while (1) 
    {
  ADCSRA|=(1<<ADSC);
  while(ADCSRA &(1<<ADSC));
  sample = ADC;  
    w[0]=a[0]*sample-a[1]*w[1]-a[2]*w[2];
    y=b[0]*w[0]+b[1]*w[1]+b[2]*w[2];
    //w[2]=w[1];
    //w[1]=w[0];
    PORTC=y;
    _delay_ms(1);
    }
}

Koeffizienten habe ich mit Matlab (fdatool) ausgerechnet.

[Jim M. (turboj)]

Mit den Auskommentierten 2 Zeilen ist das eine 1-Stufige IIR, da w[1] 
und w[2] immer Null sind.

Höhere Frequenzen würden einen Logic analyser erfordern, menschliche 
Augen sind eher langsam und können zu schnelles Blinken nicht erkennen.

Workaround: Delay verlängern und mit so niedrigen Frequenz testen, wo 
die Augen noch schnell genug fürs LED Blinken sind.

[Matthias S. (Firma: matzetronics) (mschoeldgen)]

Daniel R. schrieb:
> PORTC=y;

Schau dir mal an, wie ein float intern im AVR representiert wird und was 
passiert, wenn man das ohne jedes Casting auf einen 8-Bit Port schreibt. 
Da kann nichts vernünftiges bei passieren. Du kannst mal wenigstens 
einen 16-bit Timer für PWM konfigurieren und den als Ausgang nach einer 
float2short Umsetzung benutzen.
https://gcc.gnu.org/wiki/avr-gcc
https://people.ece.cornell.edu/land/courses/ece4760/Math/Floating_point/

[Peter II (Gast)]

Matthias S. schrieb:
> Schau dir mal an, wie ein float intern im AVR representiert wird

das spielt überhaupt keine rolle.

> wenn man das ohne jedes Casting auf einen 8-Bit Port
> rauswirft.
das ist klar in C definiert.

Ob es das ist was er will, kann ich aber nicht erkennen. Grundsätzlich 
falsch ist es zumindest nicht.

[Matthias S. (Firma: matzetronics) (mschoeldgen)]

Peter II schrieb:
> das ist klar in C definiert.

Definiert ist es, sonst hätte sich der Kompiler geweigert. Aber es ist 
nicht das, was der TE am Ausgang erwartet.

[Peter II (Gast)]

Matthias S. schrieb:
> Definiert ist es, sonst hätte sich der Kompiler geweigert. Aber es ist
> nicht das, was der TE am Ausgang erwartet.

und was soll da der Hinweis auf die interne Darstellung von Float?