-
Notifications
You must be signed in to change notification settings - Fork 0
/
visual.m
35 lines (26 loc) · 757 Bytes
/
visual.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
function visual(f,x,x0,xex)
% Postprocessing: Plot contour lines for the objective function f and the
% solution path of a particular sequence of iterates.
%
% input: f - function
% x - solution path of iterates
% x0 - initial guess
% xex - exact solution
% Specify the window in parameter space that we want to plot.
x1 = [-1.5:0.01:1.5];
x2 = [-1.5:0.01:1.5];
% Plot the contours of the objective function
Z=zeros(length(x2),length(x1));
for i=1:length(x1)
for j=1:length(x2)
Z(j,i) = f([x1(i);x2(j)]);
end
end
contour(x1,x2,Z,50);
hold on
% Then plot the initial guess, the exact solution and the iterates
plot(x0(1),x0(2),'md');
plot(xex(1),xex(2),'b*');
plot(x(1,:),x(2,:),'r-');
hold off
end