matlab微机编程,数字信号处理基于计算机的方法第三章MATLAB编程题

第三章MATLAB程序

刘伶

创建一M文件，三个程序均是函数unit3的调用function unit3(num,den,k)

w = 0:pi/(k-1):pi;

h = freqz(num, den, w);

subplot(2,2,1)

plot(w/pi,real(h));grid

title('实部')

xlabel('\omega/\pi'); ylabel('Amplitude')

subplot(2,2,2)

plot(w/pi,imag(h));grid

title('虚部')

xlabel('\omega/\pi'); ylabel('Amplitude')

subplot(2,2,3)

plot(w/pi,abs(h));grid

title('幅度 ')

xlabel('\omega/\pi'); ylabel('Magnitude')

subplot(2,2,4)

plot(w/pi,angle(h));grid

title('相位')

xlabel('\omega/\pi'); ylabel('Phase, radians')

习题3.1

六幅图分别对应r=[0.3 0.7]，theta=[0 pi/2 pi]，r T*theta的值>> num1=[1 0 0];

>>for r=[0.3 0.7]

for theta=[0 pi/2 pi];

figure;

den1=[1 r*theta r^2*theta];

unit3(num1,den1,50)

end

end