matlab + digitaler regler

[tobias hofer (Gast)]

Hallo

Gerne würde ich einmal wissen wie man am besten einen digitalen PI 
Regelalgorithmus in Maltlab testen kann.

Und zware möchte ich den programmierten PI-Regler mit der entsprechenden 
Strecke in Matlab simulieren um die Regelparameter einstellen zu können.

Ideal wäre es wenn man für die Strecke nicht extra eine DGL austellen 
müsste und die diskret durchrechnen sondern die Matlab Befehle wie TF() 
usw brauchen könnte.

Hat jemand von euch schoneinmal so etwas realisiert?

[gabi (Gast)]

ja, ich!

[tobias hofer (Gast)]

Na ja freut mich zu höhren.

Kannst du mir hier ein paar tips geben wie du das realisiert hast?
oder ein entsprechendes matlab file?

gruss

[gabi (Gast)]

mit matlab-workspace ->simulink

[tobias hofer (Gast)]

Hallo Gabi

Ehrlich gesagt eine nicht sehr hilfreiche Antwort!
Ich verstehe nicht so recht was du mit diesen Antworten bezweckst. Wenn 
du keine Auskunft geben willst, kannst dann musst du ja nicht zwingend 
Antworten.

Ich kenne Matlab und ich kenne Simulink. Mir ist aber nicht klar wie ich 
in Simulink einen Regel-Algorithmus einbinden kann, das war ja meine 
Frage. Dazu wäre ich froh ein paar tips oder bsp. zu bekommen.

Gruss Tobias

[gabi (Gast)]

ok... das ist ein konkretes problem.
wie lautet dein regel-algorithums?
um was für eine regelstrecke handelt es sich?

[tobias hofer (Gast)]

>wie lautet dein regel-algorithums?
normaler PI regler

>um was für eine regelstrecke handelt es sich?
um anzufangen und einarbeitung in matlab einfache PT1 strecke

[gabi (Gast)]

nice... schau mal hier
Beitrag "Re: MATLAB (digitaler Regler)"
hilft dir das schon weiter?

[tobias hofer (Gast)]

Danke für den Link. Was ich jetzt noch vermisse wie binde ich meinen 
algorithmus den ich ja programmiere in das simulink model(anstelle der 
simulink blöcke) ein?

Es gibt sicherlich eine möglichkeit einen Simulink Block selbst zu 
erstellen der dann meine programmierte funktion aufruft.

[gabi (Gast)]

im einfachsten fall -> transfer funktion... wo du dann zähler, nenner 
nach doppelklick eingeben kannst (num, den)!

[tobias hofer (Gast)]

meine idee:

in matlab programmiere ich den pi-regler gleich wie ich den in meinem 
zielsystem prozessor porgrammieren würde.
das ist ja eigentlich der sinn der ganzen geschichte.

in simulink arbeite ich mit den blöcken transferfkt. für die strecke 
oszilloskop für die ausgabe usw. damit muss ich in matlab weniger 
programmieren.

kann man diese module zusammenschalten wie ich mir das vorstelle? wenn 
ja wie?

[gabi (Gast)]

du kannst mit gedrückter linker maustaste von den ausgängen der blöcke 
verbindungen zu den eingängen anderer blöche ziehen (abzweigungen mit 
zusätzlich gedrückter strg-taste)... kannst dir dann beliebige 
regelkreise aufbauen.
im einfachten fall:
sollwert -> summierstelle -> regler -> strecke -> rückkopplung auf 
summierstelle (vorzeichen beachten)...

[tobias hofer (Gast)]

wie kann ich selbst einen block programmieren?
bzw. wie kann ich dann die verbindung von simulink zu meinem 
programmcode herstellen?

[Thomas (Gast)]

schau dir mal den "Embedded Matlab Function" Block an.

[tobias hofer (Gast)]

@Thomas

Das ist genau was ich gesucht habe, funktioniert bestens.

Besten Dank

[tobias hofer (Gast)]

Hallo

Nun bin ich ein bisschen am Üben. Das erste Beispiel Programm sollte ein
einfaches RC-Glied sein das ich wie folgt in die Funktion PT1 
programmiert habe.

function y = pt1(u)

dt = 0.001;
tau = 0.02;

uc_1 = 0;

for k=1:length(u)
    y(k) = 1/tau*(u(k)-uc_1)*dt + uc_1;
    uc_1 = y(k);    
end

uc_1 sollte in dieser funktion eingentlich gespeichert werden bis zum 
nächsten aufruf von Simulink. Das funktioniert jedoch noch nicht so 
richtig.

Hat jemand eine Idee wie ich das richtig realisieren kann?

Im Anhang das Simulink Model.

[Thomas (Gast)]

persistent uc_1;
if isempty(uc_1)
     uc_1 = 0;
end

deinen Code solltest du aber auch nochmal überarbeiten, in deiner 
Embedded Matlab Function sind u und y Skalare.

[tobias hofer (Gast)]

Hallo Thomas

Besten dank funtioniert, der Code für alles diese Interessiert
sieht folgendermassen aus:

function y = pt1(u)

    dt = 0.001;
    tau = 0.02;

    persistent uc_1;
    if isempty(uc_1)
         uc_1 = 0;
    end
    
    y = 2*1/tau*(u-uc_1)*dt + uc_1;
    uc_1 = y;