EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
已知两个函数:正态分布x1和矩形分布x2,利用FFT和ifft求两个函数卷积x,并显示出图像,代码如下:" u$ W' N& J' F% m' v% P/ o
fs = 100; ts = 1/fs;%时域抽样间隔 t = [-1:ts:1]; x1 = normpdf(t, 0, 0.07); x2 = unifpdf(t,-0.2,0.2);
1 z, I$ Q/ x# Z7 B8 h7 I- `* L%绘制时域波形 subplot(311) plot(t,x1); title('正态分布密度函数时域图像'); subplot(312) plot(t,x2); title('矩形分布密度函数时域图像');
' i, T6 N; T: }! L" f+ K6 DN = 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('正态分布与矩形分布联合密度函数时域图像');
) C9 x7 O/ Z0 x) D2 W: ~1 X结果无法显示x,请求高手指教
4 [) ^8 Q& n+ q4 v: o: c |