domingo, 29 de abril de 2012

Visualizacion continuada de la posicion del mouse

Hola,
me gustaría saber como puedo visualizar mediante un "edit text" la posición del puntero, cada vez que haga clic con el mouse sobre la grafica en forma continuada.
José Jeremías Caballero
Servicios de Programación en MatLab
https://programacioncaballero.blogspot.com/
Lima - Perú

Para mayor informacion, contactarse a jjcc94@hotmail.com

miércoles, 18 de abril de 2012

Ecuaciones de Fitzhugh - Nagumo

Considerando las ecuaciones de Fitzhugh-Nagumo:
dv/dt= -v*(v-a)*(v-1)-w
dw/dt=b*(v-c*w)
siendo a, b y c parámetros que cumplen:
0<a<1, b>1 y c>0.
La variable v representa el potencial de acción y la variable w modela los canales de sodio y potasio.
Ilustrar (utilizando Matlab) la existencia de condiciones iniciales (v,0), cuyas correspondientes órbitas se aproximan a equilibrios distintos. Se supone que a=1/4, b=3/2 y c=10.
Lo que me sucede, es que no entiendo como meter las condiciones iniciales que me pide en matlab. Sé representar las órbitas para los puntos de equilibrio pero sin condiciones iniciales y por mas que he probado, no he conseguido nada.

function ode45caballero8
global a b c
fprintf('CONDICIONES: 0<a<1, b>1 , c>0\n')
fprintf('\n EJEMPLO: a=0.25  b=1.5  c=10  I=[0 1]  y0=[-4 0]\n\n')
a=input('Ingrese el valor de a: ');
b=input('Ingrese el valor de b: ');
c=input('Ingrese el valor de c: ');
I=input('Ingrese el intervalo I=');
y0=input('Ingrese ondición inicial y0=');
[t,y]=ode45(@potencial_modelado1,[I(1) I(2)],[y0(1) y0(2)]);
figure(gcf)
plot(t,y(:,1), t,y(:,2));
end
function dydt=potencial_modelado1(t,y)
global a b c
dydt=zeros(2,1);
dydt(1)=-y(1)*(y(1)-a)*(y(1)-1)-y(2);  %dv/dt= -v*(v-a)*(v-1)-w
dydt(2)=b*(y(1)-c*y(2));               %dw/dt=b*(v-c*w)
end
EJECUCION EN MATLAB
>> ode45caballero8
CONDICIONES: 0<a<1, b>1 , c>0
EJEMPLO: a=0.25    b=1.5    c=10  I=[0 1]    y0=[-4 0]
Ingrese el valor de a: 0.25
Ingrese el valor de b: 1.5
Ingrese el valor de c: 3
Ingrese el intervalo I=[0 1]
Ingrese condicion inicial y0=[1 2]
José Jeremías Caballero
Servicios de Programación en MatLab
https://programacioncaballero.blogspot.com/
Lima - Perú

Para mayor informacion, contactarse a jjcc94@hotmail.com

jueves, 12 de abril de 2012

Punto de interseccion de exponencial y una recta

Hallar el punto de intersecion de una exponecial  y una recta.
clear all
figure(1)
figure(gcf)
syms  x
y1=2;
y2=x+1;
figure(gcf)
plot([-5 5],[ 2 2],[-5 5],[-5 5]+1)
hold on
X=solve(y1==y2);
plot(X,2,'*r')
hold off

clear all
figure(2)
figure(gcf)
syms  x
y1=exp(x)-20;
fplot(char(y1),[3.5 4.5],'r')
hold on
y2=5*x+20;
fplot(char(y2),[3.5 4.5],'b')
x=fzero(char(y1-y2),[3.5 4.5]);
F=inline(char(y2),'x');
plot(x,F(x),'*m')
hold off


 José Jeremías Caballero
Servicios de Programación en MatLab
https://programacioncaballero.blogspot.com/
Lima - Perú

Para mayor informacion, contactarse a jjcc94@hotmail.com

jueves, 5 de abril de 2012

Ecuaciones Diferenciales ode45

Estoy trabajando con el comando ode45, el cual requiere la asignación de una función en un archivo .m. Me encuentro definiendo dicha función f.m donde aparece la ecuación diferencial a resolver por ode45, pero ahora necesito que en esta ecuación aparezca el valor de una variable de entrada definida en el Command Window principal, como por ejemplo 'a'.
function ode45caballero5
global a
a=input('ingrese el valor de a:');
% Ejemplo
% y'=a+2t+y 1<=t<=5
% y(1)=0.5
[T,Y]=ode45(@f,[1,5],0.5);
plot(T,Y)
function dydt=f(t,y)
global a
dydt=a+2*t+y;
EJECUCION EN MATLAB
>> ode45caballero5
ingrese el valor de a:10



function ode45caballero6
global a
a=input('ingrese el valor de a:');
[t,x]=ode45(@fint2,[0,5],[-4,0]);
figure(gcf)
plot(t,x(:,1));
end

function dydt=fint2(t,x)
global a
dydt=[x(2);a-x(2)];
end
EJECUCION EN MATLAB
>> ode45caballero6

ingrese el valor de a:6

José Jeremías Caballero
Servicios de Programación en MatLab
Lima - Perú
Para mayor informacion, contactarse a jjcc94@hotmail.com

martes, 3 de abril de 2012

Impresion de una matriz cuadrada desde exterior al interior

Realizar un programa en MATLAB en el que se muestre la matriz nxn (matriz cuadrada), compuesta por los digitos desde 1 hasta nxn, comenzado desde el exterior hasta el interior de la matriz, en sentido de rotación horaria. El tamaño n es el resultado de un valor de entrada de un número entero
Programa en matlab
 
function  ordenada_desde_exterior_una_matriz1(n)
j1=1;     f1=n; 
j2=2;     f2=n;
j3=n-1;   f3=1;
j4=n-1;   f4=2;
s=0;
opcion=1;
while s<n^2
switch (opcion)
    case 1
          for i=j1:f1, s=s+1;  b(j1,i)=s, pause(0.4), end
          j1=j1+1; f1=f1-1;   opcion=2;
       
    case 2 
         for i=j2:f2, s=s+1;   b(i,f2)=s, pause(0.4), end
         j2=j2+1; f2=f2-1;    opcion=3;
         
   case 3 
        for i=j3:-1:f3, s=s+1;  b(j3+1,i)=s,  pause(0.4),end
        j3=j3-1; f3=f3+1;    opcion=4;
         
   case 4
        for i=j4:-1:f4,  s=s+1;  b(i,f4-1)=s, pause(0.4),end
        j4=j4-1;   f4=f4+1;   opcion=1;      
    end
end
 b
EJECUCION DEL PROGRAMA

>> ordenada_desde_exterior_una_matriz(4)
b =

     1     2     3     4
    12    13    14     5
    11    16    15     6
    10     9     8     7


>> ordenada_desde_exterior_una_matriz(6)
b =
     1     2     3     4     5     6
    20    21    22    23    24     7
    19    32    33    34    25     8
    18    31    36    35    26     9
    17    30    29    28    27    10
    16    15    14    13    12    11
José Jeremías Caballero
Servicios de Programación en MatLab
Lima - Perú

Para mayor informacion, contactarse a jjcc94@hotmail.com

lunes, 2 de abril de 2012

Importación y exportacion de datos a Excel

Hola, me gustaría crear una matriz que guarde estos datos
 Inicial Primario X Y Ancho Altura Visible
  A     1     10    20    30    40    on
  A     2    10    20    30    40    on
  A     3    10     20    30    40   on
  B     1     10      30    30   40   off
 ¿Se podría hacer?
 Luego quiero utilizar estos valores para un GUI
 Un saludo y gracias a quien me pueda echar una mano

 function varargout = crear_matrices_strings2_jose_caballero(varargin)
  clc
  gui_Singleton = 1;
  gui_State = struct('gui_Name',       mfilename, ...
                     'gui_Singleton',  gui_Singleton, ...
                     'gui_OpeningFcn', @crear_matrices_strings2_jose_caballero_OpeningFcn, ...
                     'gui_OutputFcn',  @crear_matrices_strings2_jose_caballero_OutputFcn, ...
                     'gui_LayoutFcn',  [] , ...
                     'gui_Callback',   []);
  if nargin && ischar(varargin{1})
      gui_State.gui_Callback = str2func(varargin{1});
  end
 
if nargout
      [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
  else
      gui_mainfcn(gui_State, varargin{:});
  end
 
function crear_matrices_strings2_jose_caballero_OpeningFcn(hObject, eventdata, handles, varargin)
handles.output = hObject;
guidata(hObject, handles);
 
function varargout = crear_matrices_strings2_jose_caballero_OutputFcn(hObject, eventdata, handles)
  a(1)=handles.text1; a(2)=handles.text2; a(3)=handles.text3;
  a(4)=handles.text4; a(5)=handles.text5; a(6)=handles.text6;
  a(7)=handles.text7; a(8)=handles.text8; a(9)=handles.text9;
  a(10)=handles.text10; a(11)=handles.text11; a(12)=handles.text12;
  a(13)=handles.text13; a(14)=handles.text14;
  for i=1:5
       set(a(i),'string',' ')
  end
 
%Impresion desde un archivo de texto
  fid =fopen('crear_matrices_strings1_jose_caballero.txt');
  datos=textscan(fid,'%s %s %s %s %s %s %s');
  fclose(fid);
  for i=1:7
       set(a(i),'String',datos{i}(1:end))
  end
  columnas=char(65:90);
%exportacion a excel el archivo de texto con matlab
for i=1:7
   xlswrite('crear_matrices_strings1_jose_caballero.xlsx',datos{i}(1:end),1,columnas(i));
  end
  %importacion desde excel hacia matlab
  [num,txt,raw]=xlsread('crear_matrices_strings1_jose_caballero.xlsx');
  for i=8:14
       set(a(i),'String',raw(:,i-7))
  end
  handles.output = hObject;
  guidata(hObject, handles);
  varargout{1} = handles.output;
 
function pushbutton1_Callback(hObject, eventdata, handles)
  web http://www.lawebdelprogramador.com/foros/Matlab/  -browser
 
function pushbutton2_Callback(hObject, eventdata, handles)
  web http://www.matlabcaballero.blogspot.com/  -browser
José Jeremías Caballero
Servicios de Programación en MatLab
Lima - Perú

Para mayor informacion, contactarse a jjcc94@hotmail.com