L'estimation de la fonction de transfert à partir d'une entrée de données de sortie

J'ai deux colonnes de la matrice de données obtenues à partir de NI-DAQ. La première colonne représente les données de sortie d'un moteur de groupe électrogène (avec chauffeur) et la deuxième colonne est la donnée d'entrée (onde carrée). Je veux trouver la fonction de transfert à l'aide de tfest sans Simulink. Est-il possible? J'ai Système d'Identification de la boîte à outils.

Comment puis-je joindre un .tapis de fichier à ce poste? Mes données sur gist https://gist.github.com/anonymous/6484844

OriginalL'auteur Matt1993 | 2013-09-08

  1. 7

    Vous avez déjà eu la bonne idée, je ne sais pas où vous avez obtenu coincé.
    Voici le code qui permet de résoudre votre problème, je l'ai testé, il fonctionne très bien. Sachez qu'une simple entrée probablement vous obtient de meilleurs résultats, ce qui signifie par exemple une seule étape au lieu d'une onde carrée.

    % load your data (text file with your two columns)
load data.txt;

% sample index, reducing of input to get single step instead of square wave
x = 1:1:length(data);
data(x > 2900,:)=[];
x(x > 2900)=[];

% plot data
figure(1)
plot(x,data(:,1)); hold on
plot(x,data(:,2)); hold off

% prepare data for tftest, 100 is a random chosen sampling time
tfdata = iddata(data(:,1),data(:,2),100);

% estimate system, factor 5 -> number of poles (variable as desired)
sys = tfest(tfdata,5);

% plot step response (factor 5 comes from input)
figure(2)
step(5*sys)

    Edit: la sortie pour le nombre de pôles np=5:

    sys =


  From input "u1" to output "y1":
    -3.337e-05 s^4 + 1.326e-07 s^3 + 1.639e-11 s^2 + 1.221e-14 s + 1.064e-19
  ----------------------------------------------------------------------------
  s^5 + 0.005287 s^4 + 1.516e-06 s^3 + 4.517e-10 s^2 + 2.896e-14 s + 2.037e-19

    Edit#2: dans votre question précédente vous demandez s'il serait préférable d'utiliser idfrd ou iddata - mais pensez-vous réellement avoir la réponse en fréquence des données? En fait, je ne pense pas qu'il devrait faire une grande différence, si vous choisissez le nombre de pôles assez élevé. Juste essayer ce que ça fonctionne le mieux pour vous. La solution est la même.

    +1 pour la douleur dans le ... de télécharger les données et de les aider quelqu'un d'autre
    merci 😉 j'ai attendu d'avoir le même problème assez rapidement, afin que je puisse l'essayer aussi maintenant.
    Je m'attendais à (np=1): sys = 19.71 / ( s + 7.8125 ). Pourquoi les coefficients de si petits lors de l'utilisation de MATLAB tfest?
    C'est certainement pas la fonction de transfert basée sur vos données en sortie, il pourrait être la fonction de transfert de votre entrée de données détaillées. Vous pouvez mélangé à quelque chose?
    Merci @thewaywewalk je vais maintenant la conception de PI (ou PID) contrôleur pour obtenir zéro-ss erreur pour l'étape d'entrée pour cette plante. Vous avez dit que vous vous attendez à avoir le même problème assez rapidement - avez-vous trouvé un bon blog ou page du forum pour montrer comment la conception de PI à partir de sys?

    OriginalL'auteur thewaywewalk