实验3 FFT算法的应用
1、2N点实数序列
N=64。用一个64点的复数FFT程序,一次算出,并绘出 的图形。
原理:
①把x[n] (0
分析
定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析
求系统
的零、极点和幅度频率响应和相位响应。
运行结果:
M代码:
clear;clf;
num=[0.0528 0.0797 0.1295 0.1295 0.797 0.0528];
den=[1 -1.8107 2.4947 -1.8801 0.9537 -0.2336];
N=128;
w=0:pi/N:pi;
[z,p,k]=tf2zp(num,den);
disp('Áãµã');disp(z);
disp('¼«µã');disp(p);
disp('ÔöÒæ');disp(k);
subplot(2,2,1);
zplane(num,den);
xlabel('Á㼫µãͼ');
H=freqz(num,den,w);
subplot(2,2,3);
plot(w/pi,abs(H));grid;
xlabel('·ùƵÏìÓ¦');
H=freqz(num,den,w);
subplot(2,2,4);
plot(w/pi,angle(H));grid;
xlabel('ÏàƵÏìÓ¦');
附加题:时域的乘积对应于频域的圆周卷积
运行结果:
M代码:
clear;clf;
N=64;
n=0:N-1;
x1=randn([1,N]);
subplot(4,1,1);
stem(x1);
xlabel('x1[n]');
x2=randn([1,N]);
subplot(4,1,2);
stem(x2);
xlabel('x2[n]');
y=x1.*x2;
subplot(4,1,3)
stem(abs(fft(y)));
xlabel('|fft(x1[n].*x2[n])|');
X1=fft(x1);
X2=fft(x2);
Y=cconv(X1,X2,N)/N;
subplot(4,1,4);
stem(abs(Y));
xlabel('|cconv(X1,X2,N)/N|');
_1383204326.unknown