Comment tracer réponse en fréquence pour une fonction de transfert d'un filtre passe-bande dans Matlab?

je veux écrire un script permettant de tracer un graphique de la fonction de transfert H(f)] pour un filtre passe-bande, |H(f)| à l'encontre de la fréquence et de la phase de H(f) (degrés) à l'encontre de la fréquence, je suis très nouveau à matlab, donc la syntaxe n'est pas à 100%, je suis confondu parce que tout est automatique formaté sous forme de matrice.
ci-dessous mon script:

% RCL circuit: band-pass filter
R=55590; L=0.9571; C=48.811*10.^(-9);      % values of the Resistor and Capacitor
f=(0:60e3); w=2*pi*f;                      % frequency (f) range of measurements
H=(R./(sqrt(R^2+(w*L-(1./(w*C))).^2)));    % Transfer Function

% Magnitude (absolute value) of the transfer function
plot(f,abs(H),'LineWidth',2); grid on; hold on
xlabel('Frequency [Hz]','FontSize',20); ylabel('|H(f)|','FontSize',20)

figure

plot(f,angle(H)*180/pi,'LineWidth',2); grid on; hold on
xlabel('Frequency [Hz]','FontSize',18); 
ylabel('phase of H(f) [degrees]','FontSize',20)

c'est la fonction de transfert de formule im en utilisant
Comment tracer réponse en fréquence pour une fonction de transfert d'un filtre passe-bande dans Matlab?

ci-dessous une autre photo de ce que mes résultats expérimentaux étaient attendus graphique, je ne comprends pas pourquoi MATLAB n'est pas ploting ce que je veux?
Comment tracer réponse en fréquence pour une fonction de transfert d'un filtre passe-bande dans Matlab?

Afin de vous calculer H et de tracer sa valeur absolue par rapport à la fréquence. Quelle est la question?
ma question est: pourquoi n'est-il pas le faire?

InformationsquelleAutor Theo | 2014-01-05

5

Êtes-vous conscient de la bodeplot fonction?

La fonction de transfert pour un simple passe-bande du second ordre est:

Il vous suffit d'insérer vos valeurs dans Matlab tf fonction et tracer avec bodeplot:
```
R = 55590;
L = 0.9571;
C = 48.811*10.^(-9);

% tf creates transfer function object 
sys = tf( [R*C 0] , [L*C R*C 1]);     % [R*C 0] vector of numerator coeffcients
                                      % R*C*s + 0*1
                                      % [L*C R*C 1] vector of denominator coeff.
                                      % L*C*s^2 + R*C*s + 0*1

bodeplot(sys)                         % plot command to plot frequency response
                                      % of magnitude and phase
```
et c'parcelles:

Alternativement, si vous avez besoin de plus de sorties comme l'amplitude et la phase de variables, l'utilisation bode, l'intrigue est égal. Mais bodeplot offre plus de terrain supplémentaire d'options de personnalisation.

Ce qui concerne votre commentaire, que vous avez besoin d'un axe linéaire:

vous suffit d'ajouter les lignes suivantes avant la commande plot:
```
P = bodeoptions;           % handle to plot options
P.MagScale = 'linear';
P.MagUnits = 'abs';
bodeplot(sys,P)            % plot command with custom options
```
de sorte qu'il semble comme suit:

Pour le réglage de l'axe des fréquences limites, utilisation:
```
P.XLim = [1 60e3];
```
ou analogique de l'ampleur:
```
P.YLim = [0 1];
```
Je vous conseille contre un linéaire de la fréquence de l'axe, mais si vous le voulez VRAIMENT, vous pouvez utiliser:
```
P.FreqScale = 'linear';
P.FreqUnits = 'Hz';       % optional
```
Si vous voulez tracer vos données expérimentales avec le graphique ci-dessus, suivez cet exemple.

Utilisation bode d'obtenir également des données formatées à partir de votre fonction de transfert, à l'instar de vos données expérimentales et l'utilisation semilogx à tracer.
```
freqVec = logspace(-1, 3, 5000);
[mag, phs] = bode(sys, freqVec * (2*pi));
mag = db(mag(:));
phs = phs(:);
figure;
subplot(211)
semilogx(freqVec, mag); hold on
semilogx(freqVec, experimentalDataMagnitude); hold off
grid on
title('System Bode Plot')
ylabel('Magnitude (dB)')
subplot(212)
semilogx(freqVec, phs); hold on
semilogx(freqVec, experimentalDataPhase); hold off
grid on
ylabel('Phase (deg)')
xlabel('Frequency (Hz)')
```
- non, je n'étais pas, mais ce n'est pas vraiment la chose im après, j'ai besoin d'axe linéaire parce que mes graphiques doivent être semblables à mes données expérimentales.
- Un linéaire de l'axe des fréquences est très rare, pour un linéaire de grandeur ont un coup d'oeil à mon montage.
- merci beaucoup, mais pouvez-vous commenter votre script comme il y a beaucoup de nouvelles fonctions, et où avez-vous d'entrée de la fonction de transfert de l'équation? ou fait-il partie d'une fonction prédéfinie? merci
- J'ai ajouté quelques commentaires, je ne sais pas ce que je doit faire des commentaires, en fait, tout est très clair. Ont également un regard sur le documentaire pour les fonctions particulières (je l'ai déjà lié à certains). L'exemple en bas, juste contient quelques fonctions de base pour le tracé logarithmique des données. Doit être facile à comprendre. Veuillez également envisager d'accepter la réponse (cliquez sur la case sur la gauche, près de ma question) - que je pense que c'entièrement répondu à votre question. Merci 🙂
- je pense que j'ai besoin pour être plus précis, je me sens bombardés d'informations et de fonctions, désolé si je n'était pas clair, je voulais dire en commentaire ce que la fonction logspace n'? est freqvec une variable? ce n'phs (:), est-ce une fonction de la bibliothèque standard? ce genre de chose? merci:)
- freqVec est une loi logarithmique espacés vectoriel engendré par logspace. Il est utilisé pour la fréquence comme un exemple. mag et phs sont les vecteurs de données pour l'amplitude et la phase, généré par le tracé de bode de la fonction basée sur freqVec. Ces 3 variables décrivent la réponse en fréquence de votre fonction de transfert sys. Le côlon opérateur (:) fait juste une colonne de vecteur. Vos données expérimentales devraient également avoir des valeurs pour la fréquence, l'amplitude et la phase. Vous la tracez ensemble avec semilogx en utilisant hold on pour garder la parcelle avant.
- pouvez-vous s'il vous plaît expliquez-moi pourquoi l'axe des x a toujours l'air logarithmique? :S
- parce que vous utilisez semilogx à l'intrigue, c'est l'idée. 😉 Sinon utiliser plot. Mais sérieusement, une fréquence linéaire de l'axe est de la connerie...
- je ne sais pas pourquoi à l'aide d'un axe linéaire, j'ai juste dit de le faire de cette façon. désolé, j'ai l'impression que je n'étais pas au clair sur l'importance de l'axe linéaire chose. jetez un oeil à mon montage, pouvez-vous encore de l'aide?
- Je veux dire que je vous ai dit avant de savoir comment vous obtenez un linéaire de l'axe des fréquences, l'utilisation de la bodeoptions pour définir P.FreqScale = 'linear' et utiliser la normale plot de fonctionner dans tous les autres cas.
- j'ai essayé le 'linéaire', chose qu'il n'a pas l'air de fonctionner, les graphiques s'est avéré rien de semblable à ce que j'attendais (à partir de l'exemple de graphique dans mon édition), elle ressemblait à ce que mon script d'origine drew :s
- avez-vous réglé les limites de l'axe avec P.XLim = [lowerBound upperBound]?
- yep viens de vérifier
- Je pense que je ne peux pas vous aider davantage. Je vous ai donné toutes les informations dont vous avez besoin. Vous avez besoin de vérifier les résultats intermédiaires, étape par étape, qui sera également vous aider à comprendre ce qui se passe réellement.
- -_- mes valeurs des composants étaient des ordres de grandeur de ce qui l'aurait été
InformationsquelleAutor thewaywewalk

Vous devez vous connecter pour publier un commentaire.