/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/ /* [ Created with wxMaxima version 0.8.7 ] */ /* [wxMaxima: title start ] Circuit RLC Paral.lel [wxMaxima: title end ] */ /* [wxMaxima: comment start ] Josep Pujadas i Jubany (jpujades a gmail punt com) 04-abril-2011 [wxMaxima: comment end ] */ /* [wxMaxima: comment start ] Els textos no porten accents perque es vegin be tant en sistemes operatius nous (UTF-8) com en antics (ISO-8859-15). NO FEU SERVIR DONCS ACCENTS A LES EXPLICACIONS. [wxMaxima: comment end ] */ /* [wxMaxima: comment start ] Fes: - Cell - Evaluate All Cells - Celda - Evaluar todas las celdas per a veure els resultats. [wxMaxima: comment end ] */ /* [wxMaxima: comment start ] Dades, amb unitats del SI (Sistema Internacional). E [V], [Volt] f [Hz], [Hertz] R [Ohm] L [H], [Henry] C [F], [Faraday] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ E:48$ f:50$ R:3$ L:10*10^(-3)$ C:150*10^(-6)$ /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Calcul d'omega, [rad] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ %omega:2*%pi*f$ float(%omega); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Calcul de XL (reactancia inductiva) i XC (reactancia capacitiva), [Ohm] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ XL:%omega*L$ float(XL); if C # 0 then XC:1/(%omega*C) else XC:0$ float(XC); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Calcul de G (conductancia), [S] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ if R = 0 then G:0 else G:1/R$ float(G); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Calcul de BL (susceptancia inductiva) i BC (susceptancia capacitiva), [S] Calcul de B (susceptancia), [S] RECORDA: - B negativa -> Circuit inductiu (predomina la bobina). - B positiva -> Circuit capacitiu (predomina el condensador). - B zero -----> Circuit ressonant, si hi ha bobina i condensador. [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ if XL = 0 then BL:0 else BL:1/XL$ float(BL); if XC = 0 then BC:0 else BC:1/XC$ float(BC); B:-BL+BC$ float(B); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Calcul de Y (admitancia total), [S] L'angle s'obte en radians (per a les sinusoidals). Per veure'l millor es dona tambe en graus sexagesimals. RECORDA: - Angle negatiu -> Circuit inductiu (predomina la bobina). - Angle positiu -> Circuit capacitiu (predomina el condensador). - Angle cero ----> Circuit ressonant, si hi ha bobina i condensador. [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Y:sqrt(G^2+B^2)$ float(Y); if G = 0 then (if B > 0 then PHI:%pi/2 else (if B < 0 then PHI:-%pi/2 else PHI:0)) else PHI:atan(B/G)$ float(PHI); float(PHI*180/%pi); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Calcul de Z (impedancia total), [Ohm] L'angle s'obte en radians (per a les sinusoidals). Per veure'l millor es dona tambe en graus sexagesimals. RECORDA: - Angle positiu -> Circuit inductiu (predomina la bobina). - Angle negatiu -> Circuit capacitiu (predomina el condensador). - Angle cero ----> Circuit ressonant, si hi ha bobina i condensador. [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ if Y = 0 then Z:0 else Z:1/Y$ float(Z); float(-PHI); float(-PHI*180/%pi); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Calcul de Emax (f.e.m. maxima d'alimentacio), [V] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Emax:E*sqrt(2)$ float(Emax); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Sinusoidal de la f.e.m. d'alimentacio, [V] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ e(t):=Emax*sin(%omega*t); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Sinusoidal de la intensitat del circuit, [A] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Imax:Emax*Y$ float(Imax); i(t):=Imax*sin((%omega*t)+PHI); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Sinusoidals de les intensitats a R, L i C, [A] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Irmax:Emax*G$ float(Irmax); ir(t):=Irmax*sin((%omega*t)); Ilmax:Emax*BL$ float(Ilmax); il(t):=Ilmax*sin((%omega*t)-(%pi/2)); Icmax:Emax*BC$ float(Icmax); ic(t):=Icmax*sin((%omega*t)+(%pi/2)); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Representacio grafica de les sinusoidals. La tensio es dibuixa sempre amb una amplitud la meitat de la de la intensitat. [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ Ymax:max(Imax,Irmax,Ilmax,Icmax)$ plot2d([(Y/2)*e(t),i(t),ir(t),il(t),ic(t)], [t,-0.25/f,1.25/f], [y,-Ymax-5,Ymax+5], [style, [lines,3],[lines,3],[lines,1],[lines,1],[lines,1]], [legend, "(Y/2)*e(t)","i(t)", "ir(t)", "il(t)", "ic(t)"], [ylabel,"A,V"], [plot_format, gnuplot], [gnuplot_preamble, "set grid;"])$ /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Frequencia de ressonancia, [Hz] [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ if L = 0 or C = 0 then fr:0 else fr:1/(2*%pi*(sqrt(L*C)))$ float(fr); /* [wxMaxima: input end ] */ /* [wxMaxima: comment start ] Comprovacio grafica de les intensitats. Les sinusoidals i(t) i iii(t) son la mateixa funcio. Nomes es veu una funcio perque l'altra queda amagada a sota. [wxMaxima: comment end ] */ /* [wxMaxima: input start ] */ iii(t):=ir(t)+il(t)+ic(t)$ plot2d([i(t),iii(t)], [t,-0.25/f,1.25/f], [y,-Imax-5,Imax+5],[legend, "i(t)","iii(t)"], [ylabel,"A"], [plot_format, gnuplot], [gnuplot_preamble, "set grid;"])$ /* [wxMaxima: input end ] */ /* Maxima can't load/batch files which end with a comment! */ "Created with wxMaxima"$