Skip to content
Permalink
Branch: master
Find file Copy path
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
35 lines (29 sloc) 830 Bytes
function [ ] = Richtungsfeld_zeichnen(f, xmin, xmax, ymin, ymax, hx, hy)
% Beispiel Parameter
% f = @(x, y) x.^2 + 0.1 * y;
% xmin = -2
% xmax = 2
% ymin = -1
% ymax = 2
% hx = 0.1
% hy = 0.1
%
% Richtungsfeld_zeichnen(@(x,y) x .^2 .* y .^2, -2, 2, -1, 2, 0.1, 0.1)
%
% bsp 2
% Richtungsfeld_zeichnen(@(x,y) sin(x) .* y./2, 0,6.282, 0, 4, 0.15, 0.15)
% bsp script 7.1
% Richtungsfeld_zeichnen(@(x,y) -0.5 .* y .* x.^2, 0,3,0, 4, 1, 1)
x = linspace(xmin, xmax, (xmax - xmin) / hx);
y = linspace(ymin, ymax, (ymax - ymin) / hy);
[x, y] = meshgrid(x, y);
vx = ones(size(x));
vy = f(x, y);
subplot(2,1,1)
quiver(x, y, vx, vy);
title('Slope - Original')
subplot(2,1,2)
len = sqrt(vx .^2 + vy.^2);
quiver(x, y, vx./len, vy./len);
title('Slope - normalised')
end
You can’t perform that action at this time.