Matlab - programa edo representar grafica

hola miren tengo que representar esta funcion con matlab y stoy bloqueado. he escrito el siugiente codigo pero me he kedado ataskado, si alguien me dijese como continuar o si lo que he escrito esta bien. estaria muy agradecido

el probelema me da la edo siguiente d2y/ dx2 = k/H con H=10, k=1, y(0)=0 i y’(0)=a


pause(3)
disp(`utilitzem la formula donada que integrada ens queda: ')

pause(3)
y = dsolve (‘k/H=D2y’, ‘x’)

pause(3)
disp(`si ara apliquem les condicions: ')

pause(3)
y = dsolve (‘(k/H=D2y’,’Dy(0)=a’,’y(0)=0’, ‘x’)

pause(3)
disp(`posant un exemple concret per a K=1, H=10')
y = dsolve (‘(1/10=D2y’,’Dy(0)=a’,’y(0)=0’, ‘x’)

disp('A PARTIR DE AQUI ES CUANDO DUDO Y NO SE COMO CONTINUAR YA QUE ESTO QUE HE ESCRITO CREO K NO ESTA BIEN')

z=0;
ezplot(z,0,a)
hold on

pause(2)
ezplot(y,0,a)
hold off

Hola Anton.
Lo que quieres hacer me parece que esto: Como vas a resolver un PVI de segundo grado con condiciones iniciales, entonces para eso utilizas la funcion dsolve() de matlab.
Ahora para graficar la solucion del PVI necesitas tener un dominio,
por ejemplo:
D2y/d2x=f(y,x) xo<=x<=xf
y(xo)=yo
y'(xo)=a

Entonces esta solucion de este PVI debe graficarse en el dominio [xo xf]
Tambien he visto que quieres hacer variar tu 'a', en el programa que he hecho lo puedes variar a tu gusto. Tambien los valores de k y H.

%==================================================


clear
%Hola.
%Los valores de H1,k1, a1, xf puedes variar a tu gusto.
% valores variables: H1,k1, a1,xf
% x0<=X<=xf
% D^2y/d^2x=k1/H1 con y(x0)=0, Dy(xf)=a1
H1=10; k1=1; a1=5; x0=0; xf=1;
k=num2str(k1);
H=num2str(H1);
a=num2str(a1);
y=dsolve(['D2y=',k,'/',H],['y(0)=0,Dy(0)=',a],'x');
f=inline(y);
X=linspace(x0,xf,20);
figure(gcf)
plot(X,f(X))
grid
%================================

EJECUCION
>> sistemasdsolve1
y te saldrá una grafica en el primer plano. El programa lo he guardado con el nombre de sistemasdsolve1.m, por eso cuando lo ejecuto en la ventana de comandos lo llamo con ese nombre, pero sin la extension.



Saludos.
JOSE JEREMIAS CABALLERO
ASESOR DE PROYECTOS CON MATLAB
PROFESOR DE METODOS NUMERICOS CON MATLAB
PROGRAMADOR EN MATLAB
[email protected]

muxisimas gracias!!! la verdad es que pinta muy bien, las funciones son un poco dificieles para mi pero kreo k buscando entendere cuales as utiliizado. muxas gracias

hola Anton.
Me alegra saber:
Te dejo lo que hace cada una de las funciones en forma concisa.

clear : Sirve para borrar variables de la memoria de matlab
% :Sirve para poner comentarios al programa.
num2str : convierte numeros a cadenas
dsolve : sirve para resolver un PVI con condiciones iniciales
inline(y) : sirve para crear un archivo function para y para poder evaluarlos
linspace : genera puntos discretos desde en el intervalo [ x0 xf] incluido los extremos. en este 20 puntos.
figure(gcf): sirve para displayar el figura en el primer plano.
plot(X,f(X)) : la funcion plot sirve para graficar. X dominio y f(X) dominio
grid : Crea mallado en la ventana de figura

para mayor ayuda:
digita en la ventana de comandos de matlab las siguientes sentencias.


>> help clear
>> help num2str
>> help dsovle
>> help inline
>> help linspace
>> help figure
>> help gcf
>> help plot
>> help grid

Ejemplo:

>> help inline
INLINE Construct INLINE object.
INLINE(EXPR) constructs an inline function object from the
MATLAB expression contained in the string EXPR. The input
arguments are automatically determined by searching EXPR
for variable names (see SYMVAR). If no variable exists, 'x'
is used.
INLINE(EXPR, ARG1, ARG2, ...) constructs an inline
function whose input arguments are specified by the
strings ARG1, ARG2, ... Multicharacter symbol names may
be used.
INLINE(EXPR, N), where N is a scalar, constructs an
inline function whose input arguments are 'x', 'P1',
'P2', ..., 'PN'.

Examples:
g = inline('t^2')
g = inline('sin(2*pi*f + theta)')
g = inline('sin(2*pi*f + theta)', 'f', 'theta')
g = inline('x^P1', 1)

See also symvar.
Overloaded methods:
sym/inline
laurpoly/inline

Reference page in Help browser
doc inline

PARA MAYOR INFORMACION DE LA FUNCION inline Y VER EJEMPLOS DE ESTA FUNCION. HACER CLIC en doc inline que esta al final de la ayuda de este funcion inline. Es facil de ver, porque está de color azul y está subrayado.




Saludos.
JOSE JEREMIAS CABALLERO
ASESOR DE PROYECTOS CON MATLAB
PROFESOR DE METODOS NUMERICOS CON MATLAB
PROGRAMADOR EN MATLAB
[email protected]