EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
已知两个函数:正态分布x1和矩形分布x2,利用FFT和ifft求两个函数卷积x,并显示出图像,代码如下:5 R& r- }; i" t+ T; _
fs = 100; ts = 1/fs;%时域抽样间隔 t = [-1:ts:1]; x1 = normpdf(t, 0, 0.07); x2 = unifpdf(t,-0.2,0.2); 7 q, b3 T2 f$ r1 M
%绘制时域波形 subplot(311) plot(t,x1); title('正态分布密度函数时域图像'); subplot(312) plot(t,x2); title('矩形分布密度函数时域图像');
6 t$ h, {( D% E! W& D3 t( WN = 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('正态分布与矩形分布联合密度函数时域图像'); 7 F: _' @1 D0 i! U7 M# m2 \( }; x- s
结果无法显示x,请求高手指教
9 Z% E# {9 R$ u$ F5 @ |