离散时间信号系统的频域分析实验报告.pdfVIP

离散时间信号系统的频域分析

实验报告

LT

三、实验过程及实验结果

clearall;

w=-4*pi:8*pi/511:4*pi;

num=[21];den=[1-0.6];

h=freqz(num,den,w);

subplot(2,1,1)

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

title(实部)

xlabel(omega^pi);

ylabel(振幅);

subplot(2,1,2)

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

title(虚部)

xlabel(omega^pi);

ylabel(振幅);

figure;

subplot(2,1,1)

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

title(幅度谱)

xlabel(omega^pi);

ylabel(振幅);

subplot(2,1,2)

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

title(相位谱)

xlabel(omega^pi);

ylabel(以弧度为单位的相位);

h=[123456789];

w=0:pi/511:pi;

h=freqz(h,1,w);

subplot(2,1,1)

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

title(实部)

xlabel(omega^pi);

ylabel(振幅);

subplot(2,1,2)

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

title(虚部)

xlabel(omega^pi);

ylabel(振幅);

figure;

subplot(2,1,1)

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

title(幅度谱)

xlabel(omega^pi);

ylabel(振幅);

subplot(2,1,2)

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

title(相位谱)

xlabel(omega^pi);

ylabel(以弧度为单位的相位

);

clearall;

w=-pi:8*pi/511:1*pi;

num=[123456789];den=[1];

num2=[zeros(1,6),num];

h=freqz(num2,den,w);

h2=freqz(num,den,w);

h3=h2.*exp(-6*j*w);

subplot(3,2,1)

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

title(幅度谱)

xlabel(omega^pi);

ylabel(振幅);

subplot(3,2,2)

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

title(相位谱)

xlabel(omega^pi);

ylabel(以弧度为单位的相位);

subplot(3,2,3)

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

title(幅度谱)

xlabel(omega^pi);

ylabel(振幅);

subplot(3,2,4)

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

title(相位谱)

xlabel(omega^pi);

ylabel(以弧度为单位的相位);

subplot(3,2,5)

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

title(幅度谱)

xlabel(omega^pi);

ylabel(振幅);

subplot(3,2,6)

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

title(相位谱)

xlabel(omega^pi);

ylabel(以弧度为单位的相位);

四、实验小结

通过本次试验熟练了使用MATLAB软件的方式和技

巧，掌握了离散时间信号与系统的频域分析方法和离散

时间信号傅里叶变换与傅里叶逆变换的基本方法。系统

的频域分析法，将通过傅里叶变换，将信号分解成多个

正弦函数的和（或积分），得到信号的频谱，然后求系统

对各个正弦分量的响应得到响应的频谱，最后通过傅里

叶反变换，得到响应再加以分析。即将信号分解成一个

个的基信号，然后研究系统对于基信号的响应，再将这

些所有的基信号的响应叠加，便是系统对于一个完整的

复杂信号的响应。通过频域分析系统在物理上更为直观，

我们比较容易通过频域看出，系统与信号的特征。