%read data from simlation

[f,a,p] = textread('C:\Documents and Settings\Bruce Land\My Documents\WebSites\ece576\DDA\AnalogSimNiosII\SimData.txt',...
    'freq=%f, amp=%f, phase=%f')

figure(1), clf;

subplot(2,1,1)
semilogy(log10(f),a,'bx')
line([log10(485) log10(485)], [100,1],'color','red')
title('Second order system')
set(gca,'xticklabels',num2str(10.^(str2num(get(gca,'xticklabels'))),'%3.0f'));

subplot(2,1,2)
plot(log10(f),p,'bx')
line([log10(485) log10(485)], [100,1],'color','red')
set(gca,'xticklabels',num2str(10.^(str2num(get(gca,'xticklabels'))),'%3.0f'));
xlabel('log10(frequency)')