% 
% ГсбцЮмбфб
%
% ГЯнефбй чсЮуз фщн ухнбсфЮуещн: plot, hist, polar, mesh, surf
%

%%%%%%%%%%%%%%%%%%%%%%%
% X-Y гсбцЮмбфб: PLOT
%%%%%%%%%%%%%%%%%%%%%%%
y = '0.25 + 0.75*cos(x.^3)';	% псйумьт ухнЬсфзузт ме блцбсйимзфйкь
angle = 0: 1: 90;					% псйумьт редЯпх фймюн гщнйюн
x = pi * angle / 180;			% псйумьт фпх бгнюуфпх x кбй мефбфспрЮ уе бкфЯнйб
z = eval( y);						% екфЯмзуз фзт ухнЬсфзузт y
figure;								% дзмйпхсгЯб кенпэ рбсбиэспх учедЯбузт
plot( angle, z);					% учедЯбуз Ч-Х гсбцЮмбфпт фзт y ухнбсфЮуей фзт гщнЯбт
drawnow;

%%%%%%%%%%%%%%%%%%%%%%%
% ЙуфпгсЬммбфб: HIST
%%%%%%%%%%%%%%%%%%%%%%%
y = rand( 1, 100);							% рбсбгщгЮ 100 фхчбЯщн бсйимюн
[ n, x] = hist( y );							% хрплпгйумьт йуфпгсЬммбфпт фймюн
figure;											% дзмйпхсгЯб кенпэ рбсбиэспх учедЯбузт
bar( x, n);										% учедЯбуз гсбцЮмбфпт уфзлюн
colormap( cool);								% кбипсйумьт фпх чсюмбфпт брейкьнйузт
% хрплпгйумьт фзт иЭузт рбспхуЯбузт фщн рпупуфюн уфпн кбмвЬ учедЯбузт
S = sum( n);
for k=1:10
	pct = 100*n( k)/S;
	text( x( k)-.02, n( k)-.5, [ num2str( pct) '%' ]);	% кеЯменп уфпн кбмвЬ учедЯбузт
end
drawnow;

%%%%%%%%%%%%%%%%%%%%%%%
% ГсбцЮмбфб рплйкюн ухнфефбгмЭнщн: POLAR
%%%%%%%%%%%%%%%%%%%%%%%
a = 60*pi/180;
theta = (pi-a)/2: a/60: (pi+a)/2;
rho = ones( size( theta));
figure;											% дзмйпхсгЯб кенпэ рбсбиэспх учедЯбузт
subplot(121);									% псйумьт рпллбрлпэ гсбцЮмбфпт кбй ерйлпгЮ иЭузт
polar( theta, rho, 'r');					% гсЬцзмб рплйкюн ухнфефбгмЭнщн
hold on;											% енфплЮ учедЯбузт уфзн Ядйб иЭуз
polar( [ (pi-a)/2 (pi-a)/2 ], [ 0 1]);	% гсЬцзмб рплйкюн ухнфефбгмЭнщн
polar( [ (pi+a)/2 (pi+a)/2 ], [ 0 1]);	% гсЬцзмб рплйкюн ухнфефбгмЭнщн
hold off;										% бкэсщуз енфплЮт учедЯбузт уфзн Ядйб иЭуз

a = 60*pi/180;
theta = (pi-a)/2: a/60: (pi+a)/2;
rho = ones( size( theta));
rho1 = rho * sin( (pi-a)/2) ./ sin( theta);
subplot(121);									% псйумьт рпллбрлпэ гсбцЮмбфпт кбй ерйлпгЮ иЭузт
polar( theta, rho, 'r');
hold on;
polar( theta, rho1);
hold off;
drawnow;

%%%%%%%%%%%%%%%%%%%%%%%
% ГсбцЮмбфб фсйюн дйбуфЬуещн: MESH, SURF
%%%%%%%%%%%%%%%%%%%%%%%
x = -13000:520:13000;
y = x;
[ X, Y] = meshgrid( x, y);					% дзмйпхсгЯб рлЭгмбфпт дэп дйбуфЬуещн
a=1; b=2; p=1000;
z = (-X.^2/a^2 + Y.^2/b^2) / (2*p);		% дзмйпхсгЯб ухнЬсфзузт рбсбдеЯгмбфпт
figure;											% дзмйпхсгЯб кенпэ рбсбиэспх учедЯбузт
subplot(121);
mesh( z);										% учедЯбуз 3D гсбцЮмбфпт рлЭгмбфпт
subplot(122);
surf( z);										% учедЯбуз 3D гсбцЮмбфпт ерйцЬнейбт
shading interp;								% ерйлпгЮ фэрпх чсщмбфйкЮт брьдпузт
