赞
踩
利用matlap来求卷积:
代码如下:
clear all
close all
clc
dt=0.01;
t=0:dt:2;
N=length(t);%matlap只能做离散运算
f1=zeros(1,N);%先对f1进行取0
n=floor(1/dt); %1s内取样点
f1(1:n)=1;%取样点都取1
y1=conv(f1,f1)dt; %conv 卷积函数
t1=((1:2N-1)-1)*dt;%自己写时间轴
figure
plot(t1,y1)
f2=exp(-5*t);
u=ones(1,N);
y2=conv(f2,u)*dt;
figure
plot (t1,y2); (1)
%plot(t,y2(1:N)); (2)
打印结果如下
需要注意的是:
在第二部分求卷积的时候
如果我们采用的是(1)中打印程序结果则变成
而采用(2)中打印程序 结果变成
原因在于我们所选的时间跨度有所不同,对于理论计算卷积应为无限长的,而matlap仿真时不可能做到无穷时间。那么对于打印时候时间跨度的选择应该格外小心。
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。