EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
已知两个函数:正态分布x1和矩形分布x2,利用FFT和ifft求两个函数卷积x,并显示出图像,代码如下:
U% _0 C: p( ?* G2 J# H7 Pfs = 100; ts = 1/fs;%时域抽样间隔 t = [-1:ts:1]; x1 = normpdf(t, 0, 0.07); x2 = unifpdf(t,-0.2,0.2); ! @4 M7 V3 f- ^
%绘制时域波形 subplot(311) plot(t,x1); title('正态分布密度函数时域图像'); subplot(312) plot(t,x2); title('矩形分布密度函数时域图像');
8 W& X: N' K2 G9 _6 \N = 512;%所做的FFT点数,2的次幂能实现快速算法 X1 = fft(x1,N);%求得x(n)的DFT结果 X2 = fft(x2,N);%求得x(n)的DFT结果 X1 = X1/fs; X2 = X2/fs; x=ifft(X1.*X2); subplot(313) plot(t,x); title('正态分布与矩形分布联合密度函数时域图像');
: D; K/ D7 _) @/ L" f) ^/ G' Q结果无法显示x,请求高手指教
( N# y, d$ i! D5 L6 S |