Anuncio

Colapsar
No hay ningún anuncio todavía.

cálculo simbólico, sustituir datos

Colapsar
X
 
  • Filtro
  • Hora
  • Mostrar
Borrar todo
nuevos mensajes

  • Matlab cálculo simbólico, sustituir datos

    Buenos días, estoy intentando sacar la incertidumbre del campo magnético mediante matlab. He conseguido hacerlo bien con un copia pega de la ventana de comandos al script de las derivadas, pero intento hacerlo todo en un script.
    Lo que intento, y pregunto, es como hacer para que después de hacer el cálculo de las derivadas parciales las variables simbólicas se sustituyan por su valor correspondiente.
    Lo que he intentado pero no funciona es:

    syms u i r d

    A =diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),i);
    B=diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),r);
    C=diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),d);

    i=double(50)
    Di=2;
    r=double(0.03);
    Dr=0.001;
    d=double(0.01);
    Dd=0.001;
    u=double(4*pi*10^(-7));

    incert=sqrt((A*Di)^2+(B*Dr)^2+(C*Dd)^2)

    Muchas gracias

  • #2
    Re: cálculo simbólico, sustituir datos

    Ya encontre la solución, por si a alguien le interesa:

    %%%% CALCULO DE INCERTIDUMBRES DE FORMA SENCILLA EN MATLAB
    %%%% ejemplo hecho con el campo magnético en el eje de la espira


    syms u i r d

    A=diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),i);
    B=diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),r);
    C=diff(u*i*r^2/(2*sqrt((d^2+r^2)^3)),d);


    Di=2; %incertidumbres sistematicas
    Dr=0.001;
    Dd=0.001;

    %%%%%%%% I r d u
    X=subs(A,{i,r,d,u},{50,0.03,0.02,4*pi*10^(-7)});%cuidado con los paréntesis y los corchetes
    Y=subs(B,{i,r,d,u},{50,0.03,0.02,4*pi*10^(-7)});
    Z=subs(C,{i,r,d,u},{50,0.03,0.02,4*pi*10^(-7)});

    incert=sqrt((X*Di)^2+(Y*Dr)^2+(Z*Dd)^2)
    Última edición por draykaler; 02/06/2012, 01:06:19.

    Comentario

    Contenido relacionado

    Colapsar

    Trabajando...
    X