Borrar filtros
Borrar filtros

Error: The size of X must match the size of U or the number of columns of U.

6 visualizaciones (últimos 30 días)
What's the problem with this code? when I run it the following Error occurs:
Error using quiver
The size of X must match the size of U or the number of columns of U.
syms x y;
f=input('enter function: ','s');
f = symfun(eval(f), [x y]);
u=-diff(f,x);
v=-diff(f,y);
[X,Y]=meshgrid(-10:.5:10);
quiver(X,Y,u,v);

Respuesta aceptada

Stephen23
Stephen23 el 14 de Mayo de 2017
Editada: Stephen23 el 14 de Mayo de 2017
There is no need to make the code slower and more complicated by using syms and eval. Exactly as I wrote in my answer to your related question, simple numeric calculations will do everything you need:
% str = input('enter a function in x and y: ','s');
str = 'log(x^2+y^2)';
fun = str2func(['@(x,y)',str]);
[X,Y] = meshgrid(-10:.51:10);
[DX,DY] = gradient(arrayfun(fun,X,Y));
quiver(X,Y,DX,DY,5);
  3 comentarios
Star Strider
Star Strider el 14 de Mayo de 2017
This was a previous attempt by OP (in another Question) that apparently did not work.
See How to use gradient (link) for context.
OP is using another approach here, so it is best to follow that approach.
geometry geometry
geometry geometry el 15 de Mayo de 2017
Thanks; sorry to make a similar post.
after you changed the name of that post I thought someone deleted that post! hence I made this post.

Iniciar sesión para comentar.

Más respuestas (1)

Star Strider
Star Strider el 14 de Mayo de 2017
Remember that ‘u’ and ‘v’ are also functions of (x,y), so use them as ‘u(x,y)’ and ‘v(x,y)’. With one change (changing the step in the meshgrid vector to 0.51 to prevent ‘division by zero’ errors), this works:
syms x y;
% f=input('enter function: ','s');
f = 'log(x^2+y^2)';
f = symfun(eval(f), [x y]);
u=-diff(f,x);
v=-diff(f,y);
[X,Y]=meshgrid(-10:.51:10);
quiver(X,Y,u(X,Y),v(X,Y),5);

Categorías

Más información sobre Numerical Integration and Differential Equations en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by