讲义二：卷积与微分方程的数值法求解.pdfVIP

讲义二：卷积与微分方程的数值法求解 一、 从离散卷积和到连续卷积 序列f 1(k)和f 2(k) 的离散卷积定义式为 ∞ f k ∗f k f i f k =−i ( ) ( ) ∑ () ( ) 1 2 1 2 i=−∞ 用来计算离散卷积的函数为： f=conv(f1,f2) f1,f2 为参与卷积运算的两个序列，f 为卷积的 结果，长度为 length(f1)+length(f2)-1 。 [f,r]=deconv(f1,f2) 解卷运算，使 f1=conv(f,f2)+r 成立 EX 错误！文档中没有指定样式的文字。-1 x k sin k ,0 ≤k ≤10 ， 1 ( ) ( ) k ，计算离散卷积和 。 x k 0.8 ,0 =≤k ≤15 y k x k ∗x k 2 ( ) ( ) 1 ( ) 2 ( ) %程序 5_1 计算离散卷积和 k1=0:10; %x1 的变量取值范围 x1=sin(k1); %构建 x1 序列 k2=0:15; %x2 的变量取值范围 x2=0.8.^k2; %构建 x2 序列 y=conv(x1,x2); %计算卷积结果 %显示卷积结果 subplot(3,1,1);stem(k1,x1);title(x_1(k)); subplot(3,1,2);stem(k2,x2);title(x_2(k)); k=0:length(y)-1; subplot(3,1,3);stem(k,y);title(y(k)); 下面讨论连续卷积的计算： 连续时间函数f 1(t) 和f 2 (t) 的卷积定义为： ∞ f t f t =∗ f t f τ f t =−τ dτ () () () ( ) ( ) 1 2 ∫−∞ 1 2 由于计算机实际处理的数据必须满足：1、离散存储；2 、有限数 据量。连续信号的处理必须首先经过数值化的过程，以离散的形式被 分析、保存和处理。用数值方法计算卷积需要将卷积积分看作信号的 分段求和来实现，这样会得到一定的精确度要求下的卷积。 ∞ ∞ f t f t =∗f t f τ f t =−τd τ lim f k =Δ f t −k Δ Δ () () () ( ) ( ) ∑ ( ) ( ) 1 2 ∫−∞ 1 2 Δ→0 1 2 k =−∞ 如果我们只求当t n=Δ(n 为整数) 时f (t) 的值f (nΔ) ,则得： ∞ ∞ f n Δ ≈ f k Δ f n Δ−k Δ Δ=Δ f k Δ f [(n −k )Δ] ( ) ∑ ( ) (