File:Level grad.png
Level_grad.png (358 × 336 pixels, file size: 70 KB, MIME type: image/png)
File information
Structured data
Captions
Transferred from en.wikipedia to Commons by Maksim.
The original description page was here. All following user names refer to en.wikipedia.
Done with Matlab.
Public domainPublic domainfalsefalse |
I, the copyright holder of this work, release this work into the public domain. This applies worldwide. In some countries this may not be legally possible; if so: I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law. |
File:Level grad.svg is a vector version of this file. It should be used in place of this PNG file when not inferior.
File:Level grad.png → File:Level grad.svg
For more information, see Help:SVG.
|
.
Source code:
Lx1=-1.7; Lx2=2.3; Ly1=-2; Ly2=2; % box dimensions N=60; % split the box into N x N grid [X, Y]=meshgrid(Lx1:1/N:Lx2, Ly1:1/N:Ly2); % the grid f=inline('-((y+1).^4/25+(x-1).^4/10+x.^2+y.^2-1)'); % draw the level sets of f fx=inline('-2/5*(x-1).^3-2*x'); fy=inline('-4/25*(y+1).^3-2*y'); % partial deriv Z=f(X, Y); % the function value figure(1); clf; hold on; axis equal; axis off; % pop up a figure h=0.5; % spacing between heights v=[-20:h:0.8 0.85]; % the heights [c,h] = contour(X, Y, Z, v, 'b'); % the level sets at those heights x0=0.1333; y0=-0.0666; % coordinates of the top of the hill delta=0.01; % descend from the top of the hill with this step size Angles=linspace(0, 2*pi, 20); % will draw 19 descent curves with Angles(i) for i=1:length(Angles) x=x0+0.1*cos(Angles(i)); y=y0+0.1*sin(Angles(i)); % starting point Curve_x=[x]; Curve_y=[y]; % will hold a descent curve following the gradient % descend from the hill for j=1:500 x=x-delta*fx(x); y=y-delta*fy(y); Curve_x=[Curve_x x]; Curve_y=[Curve_y y]; % append the updated values if max(abs(x), abs(y)) > 5 % stop when going beyond the picture frame break; end end plot(Curve_x, Curve_y, 'r') % plot the curve of steepest descent end axis([Lx1 Lx2 Ly1 Ly2]); % the picture frame saveas(gcf, 'level_grad.eps', 'psc2') % save as color postscript. Use gimp to convert % to png. Does anybody know how to reduce aliasing (gimp helps, but only a bit)?
date/time | username | edit summary |
---|---|---|
19:47, 30 January 2006 | en:User:Varuna | (get rid of ugly tag warning) |
18:54, 23 December 2004 | en:User:Oleg Alexandrov | (added source code) |
18:41, 23 December 2004 | en:User:Oleg Alexandrov | (I made this. {{PD}}) |
Original upload log edit
Legend: (cur) = this is the current file, (del) = delete this old version, (rev) = revert to this old version.
Click on date to download the file or see the image uploaded on that date.
- (del) (cur) 18:41, 23 December 2004 . . en:User:Oleg_Alexandrov Oleg Alexandrov ( en:User_talk:Oleg_Alexandrov Talk) . . 358x336 (71326 bytes) (I made this. { { PD } })
File history
Click on a date/time to view the file as it appeared at that time.
Date/Time | Thumbnail | Dimensions | User | Comment | |
---|---|---|---|---|---|
current | 15:31, 19 March 2006 | 358 × 336 (70 KB) | Maksim (talk | contribs) | La bildo estas kopiita de wikipedia:en. La originala priskribo estas: Done with Matlab. {{PD-self}} Source code: <pre> Lx1=-1.7; Lx2=2.3; Ly1=-2; Ly2=2; % box dimensions N=60; % split the box into N x N grid [X, Y]=meshgrid(Lx1:1/N:Lx2, Ly1:1/N:Ly2); % |
You cannot overwrite this file.
File usage on Commons
There are no pages that use this file.
File usage on other wikis
The following other wikis use this file:
- Usage on es.wikipedia.org
- Usage on fr.wikipedia.org
- Usage on ja.wikipedia.org
- Usage on pl.wikipedia.org
- Usage on pt.wikibooks.org
- Usage on pt.wikiversity.org
- Usage on zh.wikipedia.org