t=(0:257)*10^9/fs;
aj=plot(t(1:180),GPRJointData(:,1:1000),'b');
hold on
am=plot(t(1:180),Ndata.Moisture(1:180,1:1000),'r');
h=[am(1),aj(1)];
legend(h,'Moisture','Joint')
xlabel('Time/ns')
ylabel('Amplitude/mv')
hold off
Sometime, we should use legend to seperate different classess. Thus, we should use properties of plot and legend.