EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
已知两个函数:正态分布x1和矩形分布x2,利用FFT和ifft求两个函数卷积x,并显示出图像,代码如下:6 G9 M5 m, s0 e6 ]+ l3 A
fs = 100; ts = 1/fs;%时域抽样间隔 t = [-1:ts:1]; x1 = normpdf(t, 0, 0.07); x2 = unifpdf(t,-0.2,0.2);
3 I( k2 o1 \, p! ^. A* ^%绘制时域波形 subplot(311) plot(t,x1); title('正态分布密度函数时域图像'); subplot(312) plot(t,x2); title('矩形分布密度函数时域图像'); ! I2 {) Z: m2 S4 @4 m* P' t4 u
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('正态分布与矩形分布联合密度函数时域图像'); 0 e; l' s& c/ B' D
结果无法显示x,请求高手指教 , Q7 z$ J6 k* |2 K
|