Docsity
Docsity

Prepara tus exámenes
Prepara tus exámenes

Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity


Consigue puntos base para descargar
Consigue puntos base para descargar

Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium


Orientación Universidad
Orientación Universidad


Simulación métodos matlab, Guías, Proyectos, Investigaciones de Métodos Numéricos

Solución de problemas, es decir, los métodos que permiten obtener una solución aproximada (en ocasiones exacta) del problema considerado tras realizar algoritmos en matlab

Tipo: Guías, Proyectos, Investigaciones

2020/2021

Subido el 30/03/2022

eugenio-31
eugenio-31 🇨🇴

3 documentos

1 / 2

Toggle sidebar

Esta página no es visible en la vista previa

¡No te pierdas las partes importantes!

bg1
Parte 1
close all ;
clear all ;
clc;
% se definen las variables como simbólicas
syms yt xt D R C1 C2 L iR iC1 iC2 iL vR vC1 vC2 vL
% En MATLAB, para la solución simbólica, todas las ecuaciones se escriben
igualadas a cero.
%------------------------------------------------------------------------
---------
%Ecuaciones de equilibrio
ecn1 = iR-iC1-iL;
ecn2 = iL+iC2;
%Ecuaciones de compatibilidad
ecn3 = xt-vR-vC1;
ecn4 = vC1+vC2-vL;
ecn5 = vC1+vC2-yt;
%Ecuaciones descriptivas
ecn6 = vR-R*iR;
ecn7 = iC1-C1*D*vC1;
ecn8 = iC2-C2*D*vC2;
ecn9 = vL-L*D*iL;
%Se invoca el comando solve para obtener el modelo del circuito
s=solve(ecn1,ecn2,ecn3,ecn4,ecn5,ecn6,ecn7,ecn8,ecn9,vR,vC1,vC2,vL,iR,iC1
,iC2,iL,yt);
%Se extrae la solución para yt
soly=s.yt;
%Se simplifica la solución
y=simplify(soly);
pretty(y)
Parte 2
clear
close all
clc
syms s t
Hs=(7*s^2)/((s^3)+(7*s^2)+((80/29)*s));
pretty(Hs)
H1=(7);
H2=[1 7 80/29];
roots(H1)
roots(H2)
Xs=1/(s+0.2);
Ys=Hs*Xs;
pretty(Ys)
yt=ilaplace(Ys);
pretty(yt)
pf2

Vista previa parcial del texto

¡Descarga Simulación métodos matlab y más Guías, Proyectos, Investigaciones en PDF de Métodos Numéricos solo en Docsity!

Parte 1

close all ; clear all ; clc; % se definen las variables como simbólicas syms yt xt D R C1 C2 L iR iC1 iC2 iL vR vC1 vC2 vL % En MATLAB, para la solución simbólica, todas las ecuaciones se escriben igualadas a cero. %------------------------------------------------------------------------


%Ecuaciones de equilibrio ecn1 = iR-iC1-iL; ecn2 = iL+iC2; %Ecuaciones de compatibilidad ecn3 = xt-vR-vC1; ecn4 = vC1+vC2-vL; ecn5 = vC1+vC2-yt; %Ecuaciones descriptivas ecn6 = vR-RiR; ecn7 = iC1-C1DvC1; ecn8 = iC2-C2DvC2; ecn9 = vL-LD*iL; %Se invoca el comando solve para obtener el modelo del circuito s=solve(ecn1,ecn2,ecn3,ecn4,ecn5,ecn6,ecn7,ecn8,ecn9,vR,vC1,vC2,vL,iR,iC ,iC2,iL,yt); %Se extrae la solución para yt soly=s.yt; %Se simplifica la solución y=simplify(soly); pretty(y)

Parte 2

clear close all clc syms s t Hs=(7s^2)/((s^3)+(7s^2)+((80/29)s)); pretty(Hs) H1=(7); H2=[1 7 80/29]; roots(H1) roots(H2) Xs=1/(s+0.2); Ys=HsXs; pretty(Ys) yt=ilaplace(Ys); pretty(yt)

Parte 3

clear close all clc syms s Hs Xs Ys yu Hs=(7s^2)/((s^3)+(7s^2)+((80/29)s)); Xs=1/(s+0.2); Ys=HsXs; yu=ilaplace(Ys);

Parte 4

clear close all clc; t=0:0.001:10; yt=(1015exp(-(7t)/2).(cosh((29^(1/2)1101^(1/2)t)/ 58)+(19929^(1/2)1101^(1/2).sinh((29^(1/2)1101^(1/2)t)/58))/10643))/ 1014-((1015*exp(-t/5))/1014); plot(t,yt) grid on title('Respuesta y(t) del sistema') xlabel('Tiempo en segundos') ylabel('Amplitud')