fig3_2.m

数字信号处理Matlab演示文件,其中各个文件加放置了不同的matlab子文件

%	Figure 3.2
%	
%	BM Mar 98

clear
clf
colordef(1,'black')

rx = -3:0.1:0;
ix = -4:0.1:4.0;
[rrx,iix] = meshgrid(rx,ix);
xx = rrx + j*iix;
den = [ 1 2 5 ];
z = zeros(size(xx));
nn = max(size(rx));
for jj = 1:nn
	z(:,jj) = polyval(den,xx(:,jj));
end
h = surf(rx,ix,20*log10(1./abs(z)));hold
plot3([-1 -1 ],[ -2 -2 ],[-30 10 ])
plot3([-1 -1 ],[ 2 2 ], [ -30 10 ]) 

xlabel('real')
ylabel('imag')
zlabel('magnitude (dB)')
view([50,40]);
fprintf(1,'Figure 3.2: press return to rotate\n')
pause

rotvec = 50:5:90;
ss = max(size(rotvec));
step = 30/(ss-1);
hvec = 40:-step:10;
for ii = 1:length(rotvec)
	set(gca,'View',[rotvec(ii) hvec(ii)] )
	drawnow
	pause(2);
	end