y(i) ist für Matlab nur eine mitteloptimale Bezeichnung, weil i*i = -1 (imaginäre Einheit). Das gleiche gilt für j. Wenn man mal nicht aufpasst, passieren da interessante Dinge.
Zitat:
was sind hier genau meine filterkonstanten?
Die musst Du dir schon selbst ausdenken - in der Funktion sind die nicht enthalten.
Hast Du denn nähere Vorgaben, was das für ein Filter sein soll?
der Befehl filter funktioniert auch mit Übertragungsfunktionen siehe Code
Code:
data = [random('Poisson',1:6,1,6)]';
n = 1:1:length(data);
% Anzahl der Datenpunkte für gleitenden Mittelwert
windowSize = 3;
% PT1 Zeitkonstante und Verstärkung
T = 0.01;
P = 1;
figure() hold on
plot(n,data,'r') plot(n,filter(ones(1,windowSize)/windowSize,1,data),'b') plot(n,filter([T 1], [1], data),'c')
leged('Daten','Filter gleitender Mittelwert','Filter PT1 mit T = 0.01',1) hold off
hey, danke. aus der diskussion konnte ich so einiges herausziehen. suche einen lowpass filter - also tiefpass filter. mach mich jetz noch mal an die arbeit.
gruß, lure
nur warum wird mein beispiel instabil wenn ich num und denum vertausche??
Weil Du einmal G(s) = 1/(0.01s + 1) (PT1) und das andere mal G(s) = (0.01s + 1)/1 = 0.01s + 1 (D-Irgendwas) als Übertragungsfunktion definierst. Die verdrehte Fassung hat eine Verstärkung, die ab 1/0.01s = 100rad/s mit steigender Frequenz ansteigt. Also ein Hoch-Verstärker - niedrige Frequenzen passieren unverstärkt.
Zitat:
funktioniert filter doch nicht so wie ich mir das denke?
Doch, ich glaube schon. Nur darf man bei einer Übertragungsfunktion nicht nach belieben Zähler und Nenner vertauschen - aber das hattest Du ja auch schon angemerkt (und im Beispiel vergessen, dich dran zu halten).
clear clc closeall % data = [random('Poisson',1:100,1,100)]';
data = zeros(5,1);
data = [data; ones(40,1)];
dt = 10E-6;
t = (0:1:length(data)-1).*dt;
% Anzahl der Datenpunkte für gleitenden Mittelwert
windowSize = 3;
% PT1 Zeitkonstante und Verstärkung
T = 0.05;
D = 1;
P = 1;
num = [P];
den = [T 1];
FILTER = tf(num,den)
figure() hold on
plot(t,data,'r') plot(t,filter(ones(1,windowSize)/windowSize,1,data'),'b') plot(t,filter(den,num, data'),'c') ylim([01.2]) legend('Daten','Filter gleitender Mittelwert','Filter mit Übertragungsfunktion',4) grid on
hold off
clear clc closeall % data = [random('Poisson',1:100,1,100)]';
data = zeros(5,1);
data = [data; ones(40,1)];
dt = 10E-6;
t = (0:1:length(data)-1).*dt;
% Anzahl der Datenpunkte für gleitenden Mittelwert
windowSize = 3;
% PT1 Zeitkonstante und Verstärkung
T = 0.05;
D = 1;
P = 1;
num = [P];
den = [T 1];
FILTER = tf(num,den)
figure() hold on
plot(t,data,'r') plot(t,filter(ones(1,windowSize)/windowSize,1,data'),'b') plot(t,filter(num,den, data'),'c') ylim([01.2]) legend('Daten','Filter gleitender Mittelwert','Filter mit Übertragungsfunktion',4) grid on
hold off
und nun sage mir
sieht auch nur eins der Bilder nach einer Sprungantwort eines PT1 Gliedes aus?
war ich zu doof das PT1 richtig zu definieren?
funktioniert filter so einfach gar nicht?
sind meine Testdaten irgendwie falsch?
ok
man muss natürlich das Filter noch mit der Abtastzeit diskretisieren!!
Code:
clear clc closeall % data = [random('Poisson',1:100,1,100)]';
data = zeros(5,1);
data = [data; ones(40,1)];
dt = 10E-3;
t = (0:1:length(data)-1).*dt;
% Anzahl der Datenpunkte für gleitenden Mittelwert
windowSize = 3;
% PT1 Zeitkonstante und Verstärkung
T = 0.05;
D = 1;
P = 1;
num = [P];
den = [T 1];
FILTER = tf(num,den)
c2d(FILTER,dt)
figure() hold on
plot(t,data,'r') plot(t,filter(ones(1,windowSize)/windowSize,1,data'),'b') plot(t,filter(0.1813,[1-0.8187], data'),'c') ylim([01.2]) legend('Daten','Filter gleitender Mittelwert','Filter mit Übertragungsfunktion',4) grid on
hold off
jetzt wäre es nur noch schön zu wissen wie man per Befehl auf numerator und denumerator einer Übertragungsfunktion kommt...
also
[num, den] = "Befehl"(c2d(FILTER,dt))
das muss man dann eben aus der Hilfe zurecht suchen...
puhh
nun ist meine kleine Matlab-Welt wieder heil *G*
wieder eine Menge gelernt beim Versuch Fragen im Forum zu beantworten ;-)
Verschoben: 19.03.2009, 14:04 Uhr von Bijick Von Programmierung nach Signalverarbeitung
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
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.
RSS