|
rfm |

Forum-Newbie
|
 |
Beiträge: 2
|
 |
|
 |
Anmeldedatum: 24.07.08
|
 |
|
 |
Wohnort: ---
|
 |
|
 |
Version: ---
|
 |
|
|
 |
|
Verfasst am: 24.07.2008, 16:11
Titel: IIR Filter
|
 |
Hi Leute,
hab folgendes Problem:
Versuche einen IIR Filter (Ordnung 6) in Matlab zu pogrammieren
Code sieht wie folgt aus:
Für die ersten 7 Werte nehm ich einfach mal an yn=xn weil diese ohnehin falsch wären (Einschwingen) und es vorkommt dass ich dass in der Berechnung x[0] vorkommt, und das nimmt Matlab nicht.
Mein Problem ist dass mein Filter nicht die selben Ergebnisse liefert wie die matlab folgende Matlab Funktion:
Sieht jemand woran das liegen kann?
lg
|
|
|
|
|
Tim |

Forum-Century
|
 |
Beiträge: 140
|
 |
|
 |
Anmeldedatum: 03.11.07
|
 |
|
 |
Wohnort: Stuttgart
|
 |
|
 |
Version: 2011b+aktuellstes Release
|
 |
|
|
 |
|
Verfasst am: 24.07.2008, 17:01
Titel:
|
 |
1. Vor den Zähler-Koeffizienten (b) muss immer ein + sein
2. Für den Nenner-Term darfst du nicht bei a(1) anfangen, sondern bei a(2)*y(k-1)
|
|
|
rfm |
Themenstarter

Forum-Newbie
|
 |
Beiträge: 2
|
 |
|
 |
Anmeldedatum: 24.07.08
|
 |
|
 |
Wohnort: ---
|
 |
|
 |
Version: ---
|
 |
|
|
 |
|
Verfasst am: 25.07.2008, 07:39
Titel:
|
 |
|
 |
|
Super danke,
stimmt das mit a(1) hab ich total übersehen...
Jetzt sieht der code wie folgt aus:
[code]
b=[1 0 -3 0 3 0 -1];
a=[1 1.0794 3.1375 2.0282 2.8856 0.9127 0.7776];
x=[9.8229 -3.6812 -8.4433 6.8455 5.8779 -9.0483 -2.4869 9.9803 -1.2533 -9.5106 4.8175 7.7051 -7.7051 -4.8175 9.5106 1.2533 -9.9803 2.4869 9.0483 -5.8779 -6.8455 8.4433 3.6812 -9.8229 -1.7144e-014 9.8229 -3.6812 -8.4433 6.8455 5.8779 -9.0483 -2.4869 9.9803 -1.2533 -9.5106 4.8175 7.7051 -7.7051 -4.8175 9.5106 1.2533 -9.9803 2.4869 9.0483 -5.8779 -6.8455 8.4433 3.6812 -9.8229 -3.4289e-014 9.8229 -3.6812 -8.4433 6.8455 5.8779 -9.0483 -2.4869 9.9803 -1.2533 -9.5106 4.8175 7.7051 -7.7051 -4.8175 9.5106 1.2533 -9.9803 2.4869 9.0483 -5.8779 -6.8455 8.4433 3.6812 -9.8229 1.9621e-014 9.8229 -3.6812 -8.4433 6.8455 5.8779 -9.0483 -2.4869 9.9803 -1.2533 -9.5106 4.8175 7.7051 -7.7051 -4.8175 9.5106 1.2533 -9.9803 2.4869 9.0483 -5.8779 -6.8455 8.4433 3.6812 -9.8229 -6.8578e-014 9.8229 -3.6812 -8.4433 6.8455 5.8779 -9.0483 -2.4869 9.9803 -1.2533 -9.5106 4.8175 7.7051 -7.7051 -4.8175 9.5106 1.2533 -9.9803 2.4869 9.0483 -5.8779 -6.8455 8.4433 3.6812 -9.8229 -1.5678e-013 9.8229 -3.6812 -8.4433 6.8455 5.8779 -9.0483 -2.4869 9.9803 -1.2533 -9.5106 4.8175 7.7051 -7.7051 -4.8175 9.5106 1.2533 -9.9803 2.4869 9.0483 -5.8779 -6.8455 8.4433 3.6812 -9.8229 3.9242e-014 9.8229 -3.6812 -8.4433 6.8455 5.8779 -9.0483 -2.4869 9.9803 -1.2533 -9.5106 4.8175 7.7051 -7.7051 -4.8175 9.5106 1.2533 -9.9803 2.4869 9.0483 -5.8779 -6.8455 8.4433 3.6812 -9.8229 -3.3317e-013 9.8229 -3.6812 -8.4433 6.8455 5.8779 -9.0483 -2.4869 9.9803 -1.2533 -9.5106 4.8175 7.7051 -7.7051 -4.8175 9.5106 1.2533 -9.9803 2.4869 9.0483 -5.8779 -6.8455 8.4433 3.6812 -9.8229 -1.3716e-013 1 1 1 1 1 1 1 1 1 1 ];
y=zeros(210,1);
for k=1:201
if k<7
y(k)=x(k);
else
y(k)=b(1)*x(k)+b(2)*x(k-1)+b(3)*x(k-2)+b(4)*x(k-3)+b(5)*x(k-4)+b(6)*x(k-5)+b(7)*x(k-6)-a(2)*y(k-1)-a(3)*y(k-2)-a(4)*y(k-3)-a(5)*y(k-4)-a(6)*y(k-5)-a(7)*y(k-6);
k=k+1;
end
end
[/code]
Aber irgendwie kommt immer noch bei weitem nicht das selbe raus wie bei:
[code]
y=filter(b,a,x);
[/code]
Siehst du noch einen Fehler?
lg
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
| RSS
Hosted by:
Copyright © 2007 - 2025
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|