EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
已知两个函数:正态分布x1和矩形分布x2,利用FFT和ifft求两个函数卷积x,并显示出图像,代码如下:" j! v6 `/ ~1 A$ k! S, F
fs = 100; ts = 1/fs;%时域抽样间隔 t = [-1:ts:1]; x1 = normpdf(t, 0, 0.07); x2 = unifpdf(t,-0.2,0.2);
% p2 e C) x" M- M, t%绘制时域波形 subplot(311) plot(t,x1); title('正态分布密度函数时域图像'); subplot(312) plot(t,x2); title('矩形分布密度函数时域图像'); - D8 j) r7 Q& ~$ T6 ?9 Q' B
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('正态分布与矩形分布联合密度函数时域图像'); * }, R; e5 z: v* C0 B- ^: |
结果无法显示x,请求高手指教
/ `" g6 n: @' i& D3 q' M |