Diskrete Daten aus Übertragungsfunktion erstellen

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Peter K.

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 30.04.2012, 10:20 Titel: Diskrete Daten aus Übertragungsfunktion erstellen

Schönen guten Tag zusammen,

Ich möchte aus einer Übertragungsfunktion mit gegebenen Inputdaten die korrespondierenden Outputdaten ermitteln.

Im Moment gehe ich so vor:

Code:

for i=a_den+1:samples
y(i)=0;
for f=0:a_num-1
y(i)=y(i)+num(a_num-f)*u(i-f);
end
for f=0:a_den-1
y(i)=y(i)-den(a_den-f)*y(i-f);
end
end

Funktion ohne Link?

wobei, a_num = Anzahl der Zählerkoeffizienten und a_den = Anzahl der Nennerkoeffizienten.
Dieser Codeteil soll möglichst flexible auf die gewünschte Ordnung der Übertragungsfunktion eingestellt sein. Nun möchte ich lediglich wissen ob dies bisher so korrekt ist, da ich das Ergebnis (sprich der Ausgang) schlecht überprüfen kann.

MfG

Flashmaster

Forum-Anfänger


	Beiträge: 27

	Anmeldedatum: 02.05.12

	Wohnort: ---

	Version: ---

Verfasst am: 02.05.2012, 10:41 Titel:

Schau dir mal die lsim-Funktion an Wink

DSP

Forum-Meister


	Beiträge: 2.117

	Anmeldedatum: 28.02.11

	Wohnort: ---

	Version: R2014b

Verfasst am: 02.05.2012, 11:25 Titel:

Wie du in diesem Vergleich sehen kannst, macht dein Code nicht das was du erwartest.

Code:

Ts = 1/1000;
samples = 1000;
u = ones(samples,1);
y = zeros(samples,1);

% Vergleich mit einfachem Mittelwertfilter
windowSize = 9;
num = ones(1,windowSize)/windowSize;
den = [1, zeros(1,windowSize-1)];
y_filt = filter(num,den,u);

a_num =length(num);
a_den =length(den);

% Zeitachse
t = 0:Ts:(samples-1)*Ts;

for i=a_den+1:samples
%y(i)=0;
for f=0:a_num-1
y(i)=y(i)+num(a_num-f)*u(i-f);
end
for f=0:a_den-1
y(i)=y(i)-den(a_den-f)*y(i-f);
end
end

% Darstellung der Sprungantwort der Filter
plot(t,y,t,y_filt,'r.--');
grid on;

Funktion ohne Link?

Die Filter Funktion ist im Gegensatz zu lsim eine Standardfkt. und man benötigt keine Toolbox. Sie realisiert genau das, was dein Code machen soll. Im Übrigen könntest du die beiden inneren for-Schleifen zu einer zusammen fassen.


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	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 - 2024 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.