juanji.m

matlab实习卷积运算。很高兴与大家分享此代码。

while (1)
  sign=input('请输入要计算的卷积类型，离散为0，连续为1：');
  if (sign==0)                                                 %离散情况
  xstart=input('请输入信号的开始时间：');
  x=input('请输入离散信号的序列（以向量形式输入）：');
  hstart=input('请输入响应的开始时间：');
  h=input('请输入系统响应的序列（以向量形式输入）：');
  y=conv(x,h)
  n=(xstart+hstart):(length(y)+xstart+hstart-1)                %调整响应的起始和结束时间
  stem(n,y);                                                   %画出图形
 
  else                                                         %连续情况
  syms t
  a=input('请输入连续信号的开始时间：');
  b=input('请输入连续信号的结束时间：');
  t=a:0.001:b;                                                 %用离散信号近似，步长为0.001
  x=input('请输入信号函数x(t)=');
  x=x+t-t;                                 %处理信号是常数的情况（否则x就成为一个普通变量了）
  c=input('请输入系统响应的开始时间：');
  d=input('请输入系统响应的结束时间：');
  t=c:0.001:d;
  h=input('请输入系统响应h(t)=');
  h=h+t-t;
  y=conv(x,h);
  t=(a+c):0.001:(b+d);                                         %调整响应的起始和结束时间
  plot(t,y);                                                   %画出图形
  end
end