0001 N=4;
0002 theta=(0:0.01:pi)';
0003 stPinTaun=PwePinTaun(N,theta);
0004
0005
0006 scrsz = get(0,'ScreenSize');
0007 figAspectRatio=1;
0008 figHeight=scrsz(4)-150;
0009 figWidth=figAspectRatio*figHeight;
0010 figure('Name','PweScriptPlotPinTaun: pi_n(theta) and tau_n(theta)', ...
0011 'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0012 for n=1:N;
0013
0014 hax=subplot(N,2,2*n-1);
0015 polar([-theta(end:-1:1);theta],[(stPinTaun.taun(end:-1:1,n));(stPinTaun.taun(:,n))]);
0016 set(hax,'Outerposition',[0,1-n/N,0.5,1.0/N]);
0017 title(['\tau_n, n=' int2str(n) ]);
0018
0019 hax=subplot(N,2,2*n);
0020 polar([-theta(end:-1:1);theta],[(stPinTaun.pin(end:-1:1,n));(stPinTaun.pin(:,n))]);
0021 set(hax,'Outerposition',[0.5,1-n/N,0.5,1.0/N]);
0022 title(['\pi_n, n=' int2str(n) ]);
0023 end;
0024
0025
0026
0027
0028
0029
0030 scrsz = get(0,'ScreenSize');
0031 figAspectRatio=1;
0032 figHeight=scrsz(4)-150;
0033 figWidth=figAspectRatio*figHeight;
0034 figure('Name','PweScriptPlotPinTaun: abs(pi_n(theta)) and abs(tau_n(theta))', ...
0035 'Position',[(scrsz(3)-figWidth)/2 50 figWidth figHeight]);
0036 for n=1:N;
0037
0038 hax=subplot(N,2,2*n-1);
0039 polar([-theta(end:-1:1);theta],[abs(stPinTaun.taun(end:-1:1,n));abs(stPinTaun.taun(:,n))]);
0040 set(hax,'Outerposition',[0,1-n/N,0.5,1.0/N]);
0041 title(['|\tau_n|, n=' int2str(n) ]);
0042
0043 hax=subplot(N,2,2*n);
0044 polar([-theta(end:-1:1);theta],[abs(stPinTaun.pin(end:-1:1,n));abs(stPinTaun.pin(:,n))]);
0045 set(hax,'Outerposition',[0.5,1-n/N,0.5,1.0/N]);
0046 title(['|\pi_n|, n=' int2str(n) ]);
0047 end;
0048