INSTITUTO POLITECNICO NACIONAL
ESCUELA SUPERIOR DE INGENIERIA MECANICA Y ELECTRICA
SECCION DE ESTUDIOS DE POSGRADO E INVESTIGACION
MAESTRIA EN INGENIERIA ELECTRONICA
DR. FRANCISCO JAVIER GALLEGOS FUNESCURSO PROPEDEUTICOLENGUAJES DE PROGRAMACIONTel. 5729-6000 ext. 54608, Fax. 5729-6000 ext. 54622, email: [email protected]
OBJETIVOEvaluar y reforzar los conocimientos en lenguajes de programacin de los ingenieros aspirantes al programa de Maestra, para aplicar dichos conocimientos en la solucin de procesos y mtodos matemticos en temas indispensables en la Electrnica y ramas afines.
BIBLIOGRAFIAHebert Schildt, Programacin en Turbo C, Borland-Osbourne, Mc. Graw Hill, 1991.
Brian W. Kernighan, Dennis M. Ritchie, El Lenguaje de Programacin C, Prentice Hall, 1995.
Brian D. Hahn, Essencial MATLAB for Scientists and Engineers, John Wiley and Sons Inc., 1997.
Edward B. Magrab, Computer Integrated Experimentation, Springer-Verlag, 1991.
Steven C. Chapra, Raymond P. Canale, Mtodos Numricos con Aplicaciones en Computadoras Personales, Mc. Graw Hill, 1988.
Rudra Pratap, Getting Started With Matlab: Version 6 : A Quick Introduction for Scientists and Engineers, Oxford University Press, 2002.
Jeffrey Travis, Jim Kring, LabVIEW for Everyone: Graphical Programming Made Easy and Fun, National Instruments Virtual Instrumentation Series, 2006
LENGUAJES DE PROGRAMACION 1. INTRODUCCION.2. LENGUAJE C.3. MATLAB.4. SOLUCION DE PROBLEMAS.5. LABVIEW.
INTRODUCCION
FORTRAN (FORmula TRANslation) 50sCOBOL (COmmon Business-Oriented Language) 50sBASIC (Beginners All-purpose Symbolic Instruction Code) 60sC (Lenguaje de proposito general) Bell Laboratories, AT&T, 60sPascal (multiple Passes) 70s
Mathematica, Wolfram Research, 1980.Maple, University of Waterloo, Canada, 1981, Waterloo Maple Inc. 1988. SPICE, (Simulation Program with Integrated Circuit Emphasis), 1973. MATLAB, (MATrix LABoratory), MathWorks Inc., 1970.LabVIEW, (Laboratory Virtual Instrumentation Engineering Workbench) National Instruments, 1986.
ESTRUCTURA DE UN PROGRAMA EN LENGUAJE C
#include/* BIBLIOTECAS */. . .#include
void main()/* FUNCION PRINCIPAL */{float x,y,a,b;/* VARIABLES LOCALES */clrscr();/* PROGRAMA PRINCIPAL */printf("X = \b" );scanf("%f",&x);printf("Y = \b");scanf("%f",&y);a=x*y;printf(a= %d\n",a);resta( );getch();}
resta( )/* FUNCIONES EXTERNAS */{b=x-y;printf("b= %f\n",b);}
/* Programa #1 */#include #include
void main(){float x;clrscr();printf( "introduce un numero (-10 a 10)");scanf("%f",&x);if (x>10 || x0)printf("el numero es positivo");if(x==0)printf("el numero es cero bisigno");if(x
#include#include#includemain(){int a,b,c,d,e,f,x,y;clrscr();printf("Introduce los datos:\n");printf("\nx= \b");scanf ("%d",&x);printf("\ny= \b ");scanf ("%d",&y);a=x+y;b=x-y;c=x*y;d=x/y;e=pow(x,y);f=sqrt(x);printf("SUMA\n");printf("a= %d\n",a);printf("RESTA\n");printf("b= %d\n",b);printf("MULTIPLICACION\n");printf("c= %d\n",c);printf("DIVISION\n");printf("d= %d\n",d);printf("POTENCIA\n");printf("e= %d\n",e);printf("RAIZ CUADRADA\n");printf("f= %d\n",f);getch();}
MATLAB
FACIL DE APRENDER Y USARPOTENTE, FLEXIBLE Y EXTENSIBLEEXACTO, ROBUSTO Y RAPIDOAMPLIAMENTE UTILIZADO EN INGENIERIA Y CIENCIASUN LENGUAJE RESPALDADO POR UNA COMPAA PROFESIONAL
METODOLOGIA PARA RESOLVER PROBLEMAS DE INGENIERIA
PLANTEAR EL PROBLEMA DESCRIBIR LA INFORMACION DE ENTRADA/SALIDARESOLVER EL PROBLEMA A MANO O CON CALCULADORACREAR UNA SOLUCION MATLABPROBAR LA SOLUCION
MATLAB
INICIALIZACION: Existen 4 mtodos para asignar valores iniciales a MATLAB,
Listas explcitasOperador de dos puntosDel TecladoFunciones de MATLAB
Valores especialespiRepresenta el valor de 3.1415926i, jNmero complejoInfInfinitoNaNNo a numberclockHora (ao, mes, da, hora, minuto, segundo)dateFecha actual (20 jun 96)epsEpsilon, precisin de punto flotanteansValor calculado pero no almacenado en una variable
LISTAS EXPLICITAS
A = [3.5];% El ; al final significa supresin de impresin en pantallaB = [1.5, 3.1];% La coma separa cada valor en la matrizC = [-1,0,0; 1,1,0; 1,-1,0; 0,0,2];% El ; separa los valores de cada fila
% Definiendo una matriz listando cada fila en una lnea aparteC = [-1,0,01,1,01,-1,00,0,2];
F = [1, 52, 64, 197, 42, -42, 55, 82, 22, 109];F = [1, 52, 64, 197, 42, -42, . . .55, 82, 22, 109];% La matriz continua en la siguiente lnea
Ejemplo:
B = [1.5, 3.1];S = [3.0, B];S =
3.0000 1.5000 3.1000
S(2) = -1.0;S =
3.0000 -1.0000 3.1000
S(4) = 5.5;S(8) = 9.5;S =
3.0000 -1.0000 3.1000 5.5000 0 0 0 9.5000
EL OPERADOR DE DOS PUNTOS
% El operador : crea vectores o matrices a partir de una matriz dadaC = [-1,0,0; 1,1,0; 1,-1,0; 0,0,2]; x = C(:,1);y = C(:,2);z = C(:,3);x = y = z =-1 0 0 1 1 0 1 -1 0 0 0 2
w = C(2,:);w = C(4,:);w =w = 1 1 0 0 0 2
% El operador : puede ser usado para generar vectoresH = 1:8;% genera un vector H que tiene los nmeros 1 al 8H =1 2 3 4 5 6 7 8
tiempo = 0.0 : 0.5 : 5.0;% Genera un vector fila llamado tiempo que % contiene los nmeros del 0 al 5 en % incrementos de 0.5tiempo =
0 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000 4.5000 5.0000valores = 10 : -1 : 0;% el incremento tambin puede ser negativovalores = 10 9 8 7 6 5 4 3 2 1 0
Ejemplo:C = [-1,0,0; 1,1,0; 1,-1,0; 0,0,2]; C1 = C(:, 2:3);C2 = C(3:4, 1:2);C1 = C2 =
0 0 1 -1 1 0 0 0 -1 0 0 2x = 0 : 4; y =5 : 5 : 25; [x y] % matriz transpuestaans =
0 5 1 10 2 15 3 20 4 25Ejemplo:clc%limpia pantalla de comandosclear D % limpia variableA = zeros(3); B = ones(3, 2);C = [1 2 3; 4 2 5];D = eye(size(A); A = B =
0 0 0 1 1 0 0 0 1 1 0 0 0 1 1 C =D =
1 2 3 1 0 0 4 2 5 0 1 0 0 0 1
DEL TECLADO
% Entrada de usuarioz = input(Introduzca valores para z en corchetes: );Introduzca valores para z en corchetes: [5.1 6.3 -18.0]z = 5.1000 6.3000 -18.0000
z = input(Indique su nombre: , s);Indique su nombre: VICKY
z =
VICKY
% Salida de datostemp = 78;disp(temp);% imprime variables78 disp(grados F);% imprime textosgrados F
% permite tener mayor control en la salida con formato: %e exponencial, % %f flotante, %g usar el mas corto de e f fprintf(formato, matrices);fprintf(La temperatura es %f grados F\n, temp);La temperatura es 78.000000 grados Ffprintf(La temperatura es %4.1f grados F\n, temp);La temperatura es 78.0 grados F
El punto indica que son operaciones de elemento a elemento
Formatos de exhibicin de nmerosMatlabExhibeEjemploformat shortpor omisin15.2345format long14 decimales15.23453333333333format bank2 decimales15.23format short e4 decimales1.5235e+01format long e15 decimales1.523453333333333e+01
Operaciones en MatlabOperacinEscalaresArreglosSumaa + ba + bRestaa ba bMultiplicacina * ba .* bDivisina / ba ./ bExponenciacina ^ ba .^ b
% Operaciones entre matrices y arreglosC(1) = A(1) * B(1);C(2) = A(2) * B(2);C = A . * B;. . .C(5) = A(5) * B(5);
Ejemplo:A = [2 5 6];B = [2 3 5];C = A .* B;% [4 15 30]C = A ./ B;% [1.0000 1.6667 1.2000]C = A .^ 2;% [4 25 36]D = A .^ B;% [4 125 7776]C = 3.0 .^ A;% [9 243 729]
% Operaciones entre matrices y arreglosd = [1:5; -1:-1:-5];p = d .* 5;q = d .^ 3;
d =
1 2 3 4 5 -1 -2 -3 -4 -5
p =
5 10 15 20 25 -5 -10 -15 -20 -25
q =
1 8 27 64 125 -1 -8 -27 -64 -125% Operaciones entre matrices y arreglosd = [1:3; -1:-1:-3; 1:1:3];e = [1 2 3; 4 5 6; 7 8 9];p = d .* e;q = d * e;
d =e =
1 2 3 1 2 3 -1 -2 -3 4 5 6 1 2 3 7 8 9
p =
1 4 9 -4 -10 -18 7 16 27
q =
30 36 42 -30 -36 -42 30 36 42
FUNCIONES DE MATLABFunciones e instrucciones de seleccin
% La forma general de la instruccin if es la siguienteIf expresin lgica instrucciones . . .end
% if sencillo% if anidado% if - elseif g < 50if g < 50 if intervalo < 1 count = count + 1; count = count + 1; x_int = intervalo / 10; disp(g); disp(g); elseend if b > g x_int = 0.1; b = 0; end end end
% else - elseifIf temp > 100 disp(Demasiado caliente falla equipo)elseif temp > 90 disp(Intervalo operativo normal)elseif temp > 50 disp(Temperatura baja)else disp(Demasiado frio)end
Operadores en MatlabRelacinalesLgicos< menor que~ not mayor que| or>= mayor o igual que== igual~= diferente
CICLOS FOR y WHILEEl ciclo for tiene la siguiente estructura,for ndice = expresin instrucciones . . .endse puede utilizar el operador de dos puntos,for k = inicio : incremento : lmite instrucciones . . .end La expresin para el while es la siguiente,while expresin instrucciones . . .end
Ejemplos:
r =2;if r > 0, vol = (3/4)*pi*r^3;end%---------------------------------If g > 3 | g < 0, a = b; end % si g >3 o g < 0, entonces a = b%---------------------------------If a > 3 & c 3 y c < 0, entonces b =19%---------------------------------r = 2;if r > 3 b = 1;elseif r == 3 b = 2;else b = 0;end%---------------------------------
for r =1:5; vol = (4/3)*pi*r^3; disp([r, vol])end 1.0000 4.1888
2.0000 33.5103
3.0000 113.0973
4.0000 268.0826
5.0000 523.5988%---------------------------------for i = 1:6 for j = 1:20 if j>2*i, break, end endend%--------------------------------r = 0;while r
r =0;while r < 10 r = input(Teclee el radio (-1 para terminar): ); if r < 0, disp(FIN), break, end vol = (3/4)*pi*r^3; fprintf(Volumen = %7.3f \n, vol);end%--------------------------------- %--------------------------------for i = 1:6 v = glacial; x(i) = (i-1) * 0.1; v = [g, l, a, c, i, a, l];End v(4) = c%--------------------------------- %--------------------------------for i = 1:3 % c = a .* b for j = 1:3 c(i,j) = a(i,j) * b(i,j) endend
% Programa completo: ecuacion.m
clc %limpia la pantallaclear all %limpia las variables fprintf('"Menu principal"\n\n');disp('1. Grfica de la Ecuacion de la recta ');disp('2. Grfica de una Ecuacion exponencial ');disp('3. Comportamiento de una seal [D^2 +(R/L)D+(1/LC)]* i(t) ');disp('4.- Salir ');
x=input('Elija la opcion de la operacin que desea realizar: ');
switch x case 1 recta3 case 2 exponencial case 3 Rlc case 4 disp('FIN'); break; end
% recta3.mclc %limpia pantallaclear all %limpia las variables disp('Grafica de la ecuacin de la recta y=mx+b ');
m=input('ingresa el valor de la pendiente ');b=input('ingresa el valor de la b ');xi=input('tecle el valor inicial de x = ');xf=input('tecle el valor final de x = ');
k=1;for x=xi:xf y=(m*x)+b; y1(k)=y; x1(k)=x; k=k+1;end
plot(x1,y1,'-go');xlabel('Eje x');ylabel('Eje y');legend('recta')title('Grfica de la ecuacin de la recta');ecuacion
% exponencial.m
disp('Grafica de una ecuacin exponencial ');a=input('ingresa el valor de la magnitud ');t1=input('tecle el valor inicial de t = ');t2=input('tecle el valor final de t = '); k=1;for t=t1:t2 y=a*exp(-t); y1(k)=y; x1(k)=t; k=k+1;endplot(x1,y1,'-rx');xlabel('t');ylabel('a*e^-t');legend('exponecial')title('Grfica de la ecuacin de la recta');ecuacion
VARIABLES SIMBOLICAS
sym(0.25)ans = 1/4%------------------------syms xf=(x^3 6*x^2 + 11*x 6);pretty(f) 3 2 x - 6 x + 11 x - 6 %------------------------x = sym(x)f = (x-1)*(x-2)*(x-3);collect(f) % agrupa las variablesans = x^3-6*x^2+11*x-6 expand(f) % utiliza identidadesans = x^3-6*x^2+11*x-6
f=x^3 - 6*x^2 + 11*x 6;hornet(f)
ans =
x*(x*(x-6)+11)-6
factor(f)
ans =
(x-1)*(x-2)*(x-3)
%--------------------------------
syms xf = cos(x)^2 + sin(x)^2;simplify(f)ans = 1
SOLUCION DE ECUACIONES
syms xf=solve(x^2+x+1) % de la forma f(x)=0f = -1/2+1/2*i*3^(1/2) -1/2-1/2*i*3^(1/2)pretty(f) [ 1/2] [- 1/2 + 1/2 I 3 ] [ ] [ 1/2] [- 1/2 - 1/2 I 3 ]%----------------------------------------------------------------------------------------------syms a b c xS = a*x^2 + b*x + c;S1=solve(S) % x por defaultS1 = 1/2/a*(-b+(b^2-4*a*c)^(1/2)) 1/2/a*(-b-(b^2-4*a*c)^(1/2))
solve(S,b) % en bans = -(a*x^2+c)/xpretty(S1) [ 2 1/2] [ -b + (b - 4 a c) ] [1/2 ------------------------] [ a ] [ ] [ 2 1/2] [ -b - (b - 4 a c) ] [1/2 ------------------------] [ a ]%----------------------------------------
Ejemplo:
syms x % de la forma f(x)=p(x)s = solve('cos(2*x)+sin(x)=1') s = pi 0 1/6*pi 5/6*pi% representacin flotantedouble(s)
ans = 3.1416 0 0.5236 2.6180
vpa(s, 40) % nmero de digitos ans = 3.141592653589793238462643383279502884197 0. .5235987755982988730771072305465838140329 2.617993877991494365385536152732919070164%Graficacin de variables simbolicas% default -2pi
CALCULO DIFERENCIAL E INTEGRAL
syms x n;f = x^n;diff(f) ans = x^n*n/x
diff(f,n)ans = x^n*log(x)
%------------------------------------syms x tdiff(sin(x^2))ans =
2*cos(x^2)*x
diff(t^6,6)ans =
720diff(S,'v')diff(S,n)diff(S,'v',n)
df = diff(f) % por default usa xdf = diff(f,x)dfa = diff(f,a) % df/dadiff(f,a,2) % d^2 f / da^2 syms a xA = [cos(a*x), sin(a*x); -sin(a*x), cos(a*x)];dy = diff(A)dy = [ -sin(a*x)*a, cos(a*x)*a][ -cos(a*x)*a, -sin(a*x)*a]
ddy = diff(A,a,2) ddy = [ -cos(a*x)*x^2, -sin(a*x)*x^2][ sin(a*x)*x^2, -cos(a*x)*x^2]
LIMITE
syms h n xdc = limit((cos(x+h) cos(x))/h, h,0);
% limit(f)
% limit(f,x,a), limit(f,a)
Ejemplos:syms x a t h;limit(sin(x)/x) % 1limit(1/x,x,0,'right') % inflimit(1/x,x,0,'left') % -inflimit((sin(x+h)-sin(x))/h,h,0) % cos(x)
v = [(1 + a/x)^x, exp(-x)];limit(v,x,inf,'left') % [exp(a), 0]%------------------------------------limit(F,x,a) limit(F,a)limit(F)limit(F,x,a,'right')limit(F,x,a,'left')
INTEGRALES
syms n xint(x^n) % int(x^n, x)ans = x^(n+1)/(n+1)
int(sin(2*x), 0, pi/2);ans =
1
syms t a xf = t*cos(a*x);G = int(f, a, b)G = -t*(-sin(a*b)+sin(a^2))/asyms a t bg = cos(a*t + b);h = int(g, t) h1 = int(g, b) h = h1 = 1/a*sin(a*t+b) sin(a*t+b)
syms x k;f = exp(-(k*x)^2);int(f, -inf, inf);ans = PIECEWISE([csgn(k)/k*pi^(1/2), csgn(k^2) = 1],[Inf, otherwise])1
R = int(S)R = int(S,v)R = int(S,a,b)R = int(S,v,a,b)
SUMATORIAS
syms k XS1 = symsum(1/k^2, 1, Inf)S1 =
1/6*pi^2
S2 = symsum(X^k, k, 0,Inf)S2 = -1/(X-1)
SERIE DE TAYLOR
syms xg = exp(x*sin(x))T = taylor(g, 2, 5)T = exp(5*sin(5))+exp(5*sin(5))*(5*cos(5)+sin(5))*(x-5)
ECUACIONES DIFERENCIALES
syms ydsolve('Dy=1+y^2') % y es la variable dependiente, t es la variable independienteans =tan(t+C1)
y = dsolve('Dy=1+y^2','y(0)=1') % especificando condiciones iniciales y =tan(t+1/4*pi)%-----------------------------------------------------------
x = dsolve('(Dx)^2+x^2=1','x(0)=0')% la solucin de la ecuacin diferencial puede tener ms de una solucin
x = sin(t) -sin(t)
y = dsolve('Dy+4*y = exp(t)','y(0) = 1') y =1/4*exp(t)+exp(-4*t)*(-1/4*exp(t)+1)
%----------------------------------------------------------
y = dsolve('D2y+4*y = exp(2*x)', 'y(0)=0', 'Dy(pi) = 0', 'x') y =-1/2*sin(2*x)*2*(cosh(pi*2)+sinh(pi*2))/(4+2^2)-cos(2*x)/(4+2^2)+1/(4+2^2)*exp(2*x)
TRANSFORMADA DE FOURIER
F(w) = int(f(x)*exp(-i*w*x),x,-inf,inf)
f(x) = 1/(2*pi) * int(F(w)*exp(i*w*x),w,-inf,inf)
%*************************************************syms xf = exp(-x^2); h=fourier(f)h =pi^(1/2)*exp(-1/4*w^2)pretty(h) 1/2 2 pi exp(- 1/4 w )
g=ifourier(h)g =exp(-x^2)pretty(g) 2 exp(-x )
TRANSFORMADA DE LAPLACE
f(s) = int(f(t)*exp(-s*t),0,inf)
f(t) = 1/(2*pi*i)*int(f(s)*exp(s*t),s,c-i*inf,c+i*inf)
%*************************************************
syms tf = t^4;k=laplace(f) k =24/s^5 ilaplace(k)
ans =t^4
TRANSFORMADA Z
F(z) = symsum(f(n)/z^n, n, 0, inf)%*************************************************syms nh=ztrans(2^n)h =1/2*z/(1/2*z-1) simplify(h)ans =z/(z-2)
expand(h)ans =1/2*z/(1/2*z-1)
iztrans(h)ans =2^n
GRAFICACION
plotGrfica linealloglogGrfica logartmicasemilogx Grfica logartmica (x) y lineal (y)semilogy Grfica logartmica (y) y lineal (x)
plot(x,y,'color_linestyle_marker')
EJEMPLO:
x=[0 1 5 10 15 20];lms1=[13.25 13.38 13.89 14.62 15.43 16.73];rom1=[12.78 12.89 13.43 14.12 14.92 16.08];scm1=[16.86 16.71 16.83 16.92 16.81 17.23];a1=[13.16 13.32 13.74 14.12 14.59 15.57]; lms2=[15.47 15.86 16.32 17.06 18.09 19.22];rom2=[15.30 15.52 16.19 16.98 17.95 18.86];scm2=[16.02 16.48 16.43 17.01 17.36 18.21];a2=[15.40 15.86 16.13 16.84 17.39 18.35]; lms3=[20.79 20.54 20.95 22.20 23.42 24.69];rom3=[20.70 21.03 21.40 22.50 23.45 24.48];scm3=[20.88 20.36 20.86 21.24 22.42 23.25];a3=[20.85 20.37 20.89 21.34 22.66 23.61];
plot(x,lms2,'-.rx',x,rom2,'--g+',x,scm2,':ko',x,a2,'-b*',x,lms1,'-.rx',x,rom1,'--g+',x,scm1,':ko',x,a1,'-b*',x,lms3,'-.rx',x,rom3,'--g+',x,scm3,':ko',x,a3,'-b*');xlabel('Impulsive noise percentage');ylabel('MAE values');legend('NLMS-L','ROM','SCM','ANDW');title('Plot MAE for different filters');text(2,22,' \sigma_\epsilon = 0.25')text(2,17,' \sigma_\epsilon = 0.1')text(8,13,' \sigma_\epsilon = 0.05')
EJEMPLO
clc;clear all;
subplot(2,1,1);plot(1,49.93,'or ',1,32.69,'squarek ',2,62.77,'or ',2,35.28,'squarek ',3,52.70,'or ',3,64.35,'squarek ',4,70.53,'or ',4,60.46,'squarek ',5,63.95,'or ',5,50.65,'squarek ');set(gca,'XTickLabel',{'SRBF';' ';'MRBF';' ';'ATMRBF';' ';'MMRBF sc';' ';'MMRBF tb'})ylabel('Efficiency percentage');xlabel('RBF neural networks');grid on;
subplot(2,1,2);plot(1,49.40,'or ',1,67.31,'squarek ',2,37.23,'or ',2,64.72,'squarek ',3,47.30,'or ',3,35.65,'squarek ',4,29.47,'or ',4,39.54,'squarek ',5,36.05,'or ',5,49.35,'squarek ');set(gca,'XTickLabel',{'SRBF';' ';'MRBF';' ';'ATMRBF';' ';'MMRBF sc';' ';'MMRBF tb'})ylabel('Error percentage');xlabel('RBF neural networks');grid on;
EJEMPLO
syms t clcfprintf('AMPLITUD MODULADA\n')A = 10; %amplitudw = 10000 ; %frecuencia portadoram = 0.3 ; %indice de modulacionn = 200 ; %frecuencia naturalv=A*cos(w*t)+((m*A)/2)*(cos(w+n)*t)+((m*A)/2)*(cos(w-n)*t)ezplot(v);grid on
AMPLITUD MODULADA v = 10*cos(10000*t)-100277875666828527/72057594037927936*t
EJEMPLO
clc syms xx = (-pi/1.5)+0.01:0.01:(pi/1.5)-0.01;plot(x,csc(x))title('Grafica de la cosecante');text(0.5,155,'VICKY')text(0.5,130, 'Computacion')xlabel('Funcion de x')ylabel('Amplitud')grid on
clcsyms tr = 100; c = 0.00084; l = 0.23; w = 4000; v = 600*cos(w*t + pi)+ exp(3*t)*cos(w*t); ic = c*diff(v);it = ic;pretty(ic)vr = it*rvc = (1/c)*int(it)vl = l*diff(it)ezplot3(vr,vc,vl,[0,2*pi]) %voltaje totaltitle('Grafica de una ecuacin en 3 dimensiones')xlabel('EJE x'); ylabel('EJE y'); zlabel('EJE z')grid onaxis square 63 84 2016 sin(4000 t) + ---------- exp(3 t) cos(4000 t) - ----- exp(3 t) sin(4000 t) 25000 25 vr =201600*sin(4000*t)+63/250*exp(3*t)*cos(4000*t)-336*exp(3*t)*sin(4000*t) vc =-600*cos(4000*t)+exp(3*t)*cos(4000*t) vl =1854720*cos(4000*t)-7727995653/2500000*exp(3*t)*cos(4000*t)-2898/625*exp(3*t)*sin(4000*t)
syms xclcu=diff (exp(-2*x)-tan(4*x+3));pretty(u)ezplot(u)grid onaxis square 2 -2 exp(-2 x) - 4 - 4 tan(4 x + 3)%Grafica de laplaceclcsyms x s;x=0:0.1:10; %contador s=0.1+i; K=2/s^3*x; % Resultado de laplacestem3(real(K),imag(K),x)hold onplot3(real(K),imag(K),x,'r*') hold offview(-39.5,62)
clcsyms x;a=sin(x);y=laplace(a);ezplot(y)pretty(y) 1 -------- 2 s + 1
clcsyms xu = fourier(sin(x^2))ezplot(u)title('Fourier')text(0,1.9,'VICKY')text(0,1.6,'Computacion')xlabel('Frecuencia')ylabel('Amplitud de la funcion')grid offaxis squareu = 1/2*pi^(1/2)*2^(1/2)*(cos(1/4*w^2)-sin(1/4*w^2))
syms x aclcu = int((sin(x))^2*cos(x))ezplot(u)title('Grafica de la Integal (sin(x))^2*cos(x).')xlabel('Eje X')ylabel('Eje Y')text(0.1,2,'VICKY')text(0.2,2,'Computacin')grid onaxis squareu = 1/3*sin(x)^3
syms tclcfprintf('ONDA ELECTROMAGNETICA COMUN')E = 12; %amplitud electricaB = 20; %amplitud magneticaW = 2000 %frecuencia electricaQ = 1000 %frecuencia magneticat = 0:pi/50:10*pi;plot3(E*sin(W*t),B*cos(Q*t),t)grid onaxis square
t = 0:.000001:.010;y1 = 10*(1-exp(-300*t));y2 = .009*exp(-300*t);[AX,H1,H2] = plotyy(t,y1,t,y2,'plot');title('Respuesta Forzada de un Circuito RC');xlabel('Tiempo que Tarda en Cargarse el Capacitor');text(.0065,8.3, 'Incremento del voltage')set(H1,'LineStyle','--'); %Tipo de linea de la grafiaca 1set(H2,'LineStyle',':'); %Tipo de linea de la grafiaca 1set(get(AX(1),'Ylabel'),'String','Voltage en el Capacitor'); set(get(AX(2),'Ylabel'),'String','Corriente en el capacitor');grid on
CONCLUSIONESESTO ES EL FIN