Ingetdera.. De koefficienterna definierar en oscillator med frekvens motsvarande hög respektive låg frekvens. Närmast kan det ses som en DFT med enbart två buckets (en centrerad på F+delta och en på F-delta).AndersG skrev:Jo, jag vet. Har tom tentat det en gång. Vad jag var ute efter var närmast vilken typ av filter jag får med koefficienterna uträknade som ovan?
Således ett rektangulärt?Du får ett FIR-filter ja, men inget av de fyra grundtyperna.
Nej, det stämmer. Läste slarvigt. Så hur blir de då ett filter, vilket det bevisligen blir när jag multiplicerar det samplade indatat med cosinus för Fc?Det stämmer inte med ett sinc-filter heller.. du har ingen division med x. Koefficienterna har alla amplitud 1.
Du kan se det lite som att du först mixar ner din önskade signal till DC, och därefter tar ett medelvärde över alla tapparna. Den komplexa varianten är egentligen renare och bättre och enklare att förstå, tycker jagVad jag egentligen undrar över är vilken typ av filter jag får med:Kod: Markera allt
cos(2*pi*(F-Fdelta/2)/Fs*(0:N-1));
Kod: Markera allt
exp(j*2*pi*(F-Fdelta/2)/Fs*(0:N-1));
Kod: Markera allt
//
// Calculate Basic exponential smoothing for hi
//
resultACCA = __builtin_clr(); // clear out accumulator
// square signal
resultACCA = __builtin_mpy(outputSignal[0], outputSignal[0], NULL, NULL, 0, NULL, NULL, 0);
ftemp1 = __builtin_sac(resultACCA, 0); //
resultACCA = __builtin_clr(); // clear out accumulator
// multiply by weighthing factor
resultACCA = __builtin_mpy(ftemp1, f2, NULL, NULL, 0, NULL, NULL, 0);
ftemp1 = __builtin_sac(resultACCA, 0); // Weighted new value
//spi_send((ftemp1^0x8000)>>4);
resultACCA = __builtin_clr(); // clear out accumulator
// multiply by weighthing factor
resultACCA = __builtin_mpy(hi_ave, f1, NULL, NULL, 0, NULL, NULL, 0);
ftemp2 = __builtin_sac(resultACCA, 0); // weighted old value
//
//hi_ave = ftemp1 + ftemp22;
resultACCA = __builtin_clr(); // clear out accumulator
resultACCA = __builtin_add(resultACCA, ftemp1, 0);
resultACCA = __builtin_add(resultACCA, ftemp2, 0);
hi_ave = __builtin_sac(resultACCA, 0);
//
//spi_send((hi_ave^0x8000)>>4);
//
// Calculate Basic exponential smoothing for lo
//
resultACCB = __builtin_clr(); // clear out accumulator
// square signal
resultACCB = __builtin_mpy(outputSignal2[0], outputSignal2[0], NULL, NULL, 0, NULL, NULL, 0);
ftemp1 = __builtin_sac(resultACCB, 0); //
resultACCB = __builtin_clr(); // clear out accumulator
// multiply by weighthing factor
resultACCB = __builtin_mpy(ftemp1, f2, NULL, NULL, 0, NULL, NULL, 0);
ftemp1 = __builtin_sac(resultACCB, 0); // Weighted new value
//spi_send((ftemp1^0x8000)>>4);
resultACCB = __builtin_clr(); // clear out accumulator
// multiply by weighthing factor
resultACCB = __builtin_mpy(lo_ave, f1, NULL, NULL, 0, NULL, NULL, 0);
ftemp2 = __builtin_sac(resultACCB, 0); // weighted old value
//
//hi_ave = ftemp1 + ftemp22;
resultACCB = __builtin_clr(); // clear out accumulator
resultACCB = __builtin_add(resultACCB, ftemp1, 0);
resultACCB = __builtin_add(resultACCB, ftemp2, 0);
lo_ave = __builtin_sac(resultACCB, 0);Variablerna cnt och cnt2 räknas upp vid varje inkommande sample från ADC. Blir litet förvirrande i texten med "sample" som dels kan avse sample an inkommande vågform och dels att jag samplar differensen mellan hög och låg. Uppifrån och ned:case SITOR_NEXTBIT: // Rest of bits..
//
// Measure time efter last sample and check on change
//
if(value != last_value)
{
adj = (cnt2 - 40)/4;
if(adj>4) adj = 4;
if(adj<-4) adj = -4;
//spi_send3(adj);
}
last_value = value;
if(cnt>=(80+adj))
{
cnt = 0;
cnt2 = 0;
sitor_value = value;
sitor_state = SITOR_NEXTBIT;
if(value)
{
sitorbyte = sitorbyte | mask;
LATBbits.LATB1 = 1; // Wiggle för att visa en byte
LATBbits.LATB1 = 0;
}
mask = mask <<1;
break;
}
break;
Men som sagt var: Hur äter man en elefant? Jo en bit i taget, så här är bitarna:
