clear;
data=xlsread('14-5.xls');
X=data(:,3);
Y=data(:,2);
%Z = peaks(X);
%[X,Y] = meshgrid(data(:,1),data(:,2));
%Z = sin(X) + cos(Y);
figure;
hold on;
%surfl(Z)
%colormap(pink) % change color map
%shading interp % interpolate colors across lines and faces
%surf(Y,X,Z)
plot(Y,X,'.');
plot(0:360,max(abs(X))*sin((0:360)*pi/180));
%patternCustom(Y,Z,Y);
%set(gca, 'XDir','reverse')
set(gca,'xtick',[0 90 180 270 360]);
grid on;
xlim([0 360])
%ylim([-0.02 0.02])
title ('Partial Discharge');
xlabel ('Phase angle ^o');
ylabel ('Amplitude');
legend ('');
hold off;