matlab里sconv原理_matlab中的conv和sconv的区别 貌似我也看到conv里面是两个连续信号 所以为什么conv不能用来求连续信号卷积...

展开全部

sconv是别人写的函数吧？

参考资料636f70793231313335323631343130323136353331333363396438：网页链接

function  [f,k]=sconv(f1,f2,k1,k2,p)

%计算连续信号卷积积分f(t)=f1(t)*f2(t)

% f:  卷积积分f(t)对应的非零样值向量

% k： f(t)的对应时间向量

% f1: f1(t)非零样值向量

% f2: f2(t)的非零样值向量

% k1: f1(t)的对应时间向量

% k2:  序列f2(t)的对应时间向量

% p：  取样时间间隔

f=conv(f1,f2);          %计算序列f1与f2的卷积和f

f=f*p;

k0=k1(1)+k2(1);         %计算序列f非零样值的起点位置

k3=length(f1)+length(f2)-2;      %计算卷积和f的非零样值的宽度

k=k0:p:k3*p;                                   %确定卷积和f非零样值的时间向量

subplot(2,2,1)

plot(k1,f1)          %在子图1绘f1(t)时域波形图

title('f1(t)')

xlabel('t')

ylabel('f1(t)')

subplot(2,2,2)

plot(k2,f2)          %在子图2绘f2(t)时波形图

title('f2(t)')

xlabel('t')

ylabel('f2(t)')

subplot(2,2,3)

plot(k,f);                                      %画卷积f(t)的时域波形

h=get(gca,'position');

h(3)=2.5*h(3);

set(gca,'position',h)                      %将第三个子图的横坐标范围扩为原来的2.5倍

title('f(t)=f1(t)*f2(t)')

xlabel('t')

ylabel('f(t)')