% PHD_PvD Investigate P & Q fields using Wait's formulations % for the following figures, use data generated from... load 'PHD_PvD_data'; filetype = '-dtiff'; points = 301; limit = 30; X = 0 : limit/(points-1) : limit; % --------------------------- plot argand diagrams ----------------- fig=1; figure(fig); D=0; T=1; Z=0; mode=3; parms = ['(0 : ', num2str(limit/(points-1)) ':', num2str(limit), ... ',', num2str(D), ',', num2str(T), ',',num2str(Z), ',',num2str(mode), ')' ]; cla; plot (sommerfeld(X,D,T,Z,mode),'k.-','linewidth',2,'markersize',12); set(gca,'LineWidth',1.5,'FontSize',12,'TickDir','out') xlabel('Real part of integrand','fontsize',12) ylabel('Imaginary part of integrand','fontsize',12) title(['Sommerfeld', parms],'fontsize',12) print(filetype, '', ['PvD-' num2str(fig)] ) fig=2; figure(fig); cla; plot (X,abs(sommerfeld(X,D,T,Z,mode)),'k-','linewidth',2,'markersize',12); set(gca,'LineWidth',1.5,'FontSize',12,'TickDir','out') xlabel('variable of integration','fontsize',12) ylabel('abs(integrand)','fontsize',12) title(['Sommerfeld', parms],'fontsize',12) print(filetype, '', ['PvD-' num2str(fig)] ) break % --------------------------- plot argand diagrams ----------------- fig=3; figure(fig); D=2; T=6; Z=0; mode=3; parms = ['(0 : ', num2str(limit/(points-1)) ':', num2str(limit), ... ',', num2str(D), ',', num2str(T), ',',num2str(Z), ',',num2str(mode), ')' ]; cla; plot (sommerfeld(X,D,T,Z,mode),'k.-','linewidth',2,'markersize',12); set(gca,'LineWidth',1.5,'FontSize',12,'TickDir','out') xlabel('Real part of integrand','fontsize',12) ylabel('Imaginary part of integrand','fontsize',12) title(['Sommerfeld', parms],'fontsize',12) print(filetype, '', ['PvD-' num2str(fig)] ) fig=4; figure(fig); cla; plot (X,abs(sommerfeld(X,D,T,Z,mode)),'k-','linewidth',2,'markersize',12); set(gca,'LineWidth',1.5,'FontSize',12,'TickDir','out') xlabel('variable of integration','fontsize',12) ylabel('abs(integrand)','fontsize',12) title(['Sommerfeld', parms],'fontsize',12) print(filetype, '', ['PvD-' num2str(fig)] ) % --------------------------- plot abs(Q) v D for different T ----------------- points = 201; limit = 5; X=linspace(0,limit,points); % tweak a couple of X values to show troughs X(find(X==1.425)) = sqrt(2); X(find(X==2.625)) = 2.62255; % ------------------------ T = 0 to 8 ------------------------ % ------------------------ abs(Q) ---------------------------- fig=5; figure(fig); cla; semilogy(X,abs(Qdata1),'k-','linewidth',1.5); axis([0 5 1E-8 1]); parms = ['D=[x], T=0:8, Z=',num2str(Z), ', mode=',num2str(mode) ]; set(gca,'LineWidth',1.5,'FontSize',12,'TickDir','out') xlabel('D: Normalised Offset','fontsize',12) ylabel('|Q| Normalised Vertical Field','fontsize',12) title(['PQintegral: ', parms],'fontsize',12) print(filetype, '', ['PvD-' num2str(fig)] ) % ------------------------ phase(Q) -------------------------- fig = fig+1; figure(fig); jcount = 9; start = 0; cla; plotphase1(X,points, jcount, start, Qdata1) axis([0 5 -1080 0]); set(gca, 'YTick', linspace(-1080,0,7)); parms = ['D=[x], T=0:8, Z=',num2str(Z), ', mode=',num2str(mode) ]; set(gca,'LineWidth',1.5,'FontSize',12,'TickDir','out') xlabel('D: Normalised Offset','fontsize',12) ylabel('angle(Q) Normalised Vertical Field','fontsize',12) title(['PQintegral: ', parms],'fontsize',12) print(filetype, '', ['PvD-' num2str(fig)] ) % accidentally deleted data for figs 7-12.