Diagramas de Bode de Un Rotor Simétrico (Alejandra)

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 3

Alejandra Arriola Gómez

30/09/2020

Diagramas de Bode de un rotor simétrico

Respuesta del sistema

Código de MATLAB
clc; clear all;

%Coeficientes de amortiguamiento
c1=0.08333;
c2=0.16;
c3=0.25;
c4=0.375;
c5=0.5;

%Valores de Omega
Omega=0:0.0001:2;

%Respuestas para cada Coeficiente de amortiguamiento

%N=0
R10=(Omega.^0)./(sqrt((1-(Omega.^2)).^2+(2*c1*Omega).^2));
R20=(Omega.^0)./(sqrt((1-(Omega.^2)).^2+(2*c2*Omega).^2));
R30=(Omega.^0)./(sqrt((1-(Omega.^2)).^2+(2*c3*Omega).^2));
R40=(Omega.^0)./(sqrt((1-(Omega.^2)).^2+(2*c4*Omega).^2));
R50=(Omega.^0)./(sqrt((1-(Omega.^2)).^2+(2*c5*Omega).^2));

%N=1
R11=(Omega.^1)./(sqrt((1-(Omega.^2)).^2+(2*c1*Omega).^2));
R21=(Omega.^1)./(sqrt((1-(Omega.^2)).^2+(2*c2*Omega).^2));
R31=(Omega.^1)./(sqrt((1-(Omega.^2)).^2+(2*c3*Omega).^2));
R41=(Omega.^1)./(sqrt((1-(Omega.^2)).^2+(2*c4*Omega).^2));
R51=(Omega.^1)./(sqrt((1-(Omega.^2)).^2+(2*c5*Omega).^2));

%N=2
R12=(Omega.^2)./(sqrt((1-(Omega.^2)).^2+(2*c1*Omega).^2));
R22=(Omega.^2)./(sqrt((1-(Omega.^2)).^2+(2*c2*Omega).^2));
R32=(Omega.^2)./(sqrt((1-(Omega.^2)).^2+(2*c3*Omega).^2));
R42=(Omega.^2)./(sqrt((1-(Omega.^2)).^2+(2*c4*Omega).^2));
R52=(Omega.^2)./(sqrt((1-(Omega.^2)).^2+(2*c5*Omega).^2));

%N=3
R13=(Omega.^3)./(sqrt((1-(Omega.^2)).^2+(2*c1*Omega).^2));
R23=(Omega.^3)./(sqrt((1-(Omega.^2)).^2+(2*c2*Omega).^2));
R33=(Omega.^3)./(sqrt((1-(Omega.^2)).^2+(2*c3*Omega).^2));
R43=(Omega.^3)./(sqrt((1-(Omega.^2)).^2+(2*c4*Omega).^2));
R53=(Omega.^3)./(sqrt((1-(Omega.^2)).^2+(2*c5*Omega).^2));
Alejandra Arriola Gómez
30/09/2020

%N=4
R14=(Omega.^4)./(sqrt((1-(Omega.^2)).^2+(2*c1*Omega).^2));
R24=(Omega.^4)./(sqrt((1-(Omega.^2)).^2+(2*c2*Omega).^2));
R34=(Omega.^4)./(sqrt((1-(Omega.^2)).^2+(2*c3*Omega).^2));
R44=(Omega.^4)./(sqrt((1-(Omega.^2)).^2+(2*c4*Omega).^2));
R54=(Omega.^4)./(sqrt((1-(Omega.^2)).^2+(2*c5*Omega).^2));

figure(1)
subplot(2,2,1)
plot(Omega,R10,Omega,R20,Omega,R30,Omega,R40,Omega,R50);
title('N=0');
xlabel('Frecuencia');ylabel('Amplitud');
getframe;
grid on;
hold on;

subplot(2,2,2)
plot(Omega,R11,Omega,R21,Omega,R31,Omega,R41,Omega,R51);
title('N=1');
xlabel('Frecuencia');ylabel('Amplitud');
getframe;
grid on;
hold on;

subplot(2,2,3)
plot(Omega,R12,Omega,R22,Omega,R32,Omega,R42,Omega,R52);
title('N=2');
xlabel('Frecuencia');ylabel('Amplitud');
getframe;
grid on;
hold on;

subplot(2,2,4)
plot(Omega,R13,Omega,R23,Omega,R33,Omega,R43,Omega,R53);
title('N=3');
xlabel('Frecuencia');ylabel('Amplitud');
getframe;
grid on;
hold on;

figure(2)
plot(Omega,R14,Omega,R24,Omega,R34,Omega,R44,Omega,R54);
title('N=4');
xlabel('Frecuencia');ylabel('Amplitud');
getframe;
grid on;
hold on;
Alejandra Arriola Gómez
30/09/2020

Gráficas de MATLAB

También podría gustarte