09物电.docVIP

MatLab与信息处理实验三 西华师大物电学院 学生信息：班级：2009级10班 学号：2009341047 姓名：钟一先 注意: 1,请按要求作题,填写求解的MatLab语句和输出的答案,并保存文件. 2,仔细填写 “学生信息”. 3,把此word文档的文件名更名为你的 “系别班级姓名”.上传到教师上相应的班级的文件夹中。例： “MatLab%%%.doc” ——- “0901敖兵.doc” 已知各离散序列的波形如图所示，试用Matlab求下列卷积和，并绘出卷积和序列的时域波形。 （1） （2） (3) 答：（1）n1=-1:1; x1=[1 2 1]; n2=-2:2; x2=ones(1,5); n=(n1(1)+n2(1)):(n1(1)+n2(1)+length(n1)+length(n2)-2) n = -3 -2 -1 0 1 2 3 f1=conv(x1,x2) f1 = 1 3 4 4 4 3 1 stem(n,f1,filled) title(f1=f1k*f2k) xlabel(n) (2) n3=0:2; x3=[3 2 1]; n4=0:3; x4=[1 -1 1 -1]; n=(n3(1)+n4(1)):(n3(1)+n4(1)+length(n3)+length(n4)-2) n = 0 1 2 3 4 5 f2=conv(x3,x4) f2 = 3 -1 2 -2 -1 -1 stem(n,f2,filled) title(f2=f3k*f4k) xlabel(n) (3)编辑函数 edit jxl function[x,n]=jxl(x1,x2,n1,n2) %两离散序列相加 n=min(min(n1),min(n2)):max(max(n1),max(n2)); s1=zeros(1,length(n)); s2=s1; s1(find((n=min(n1))(n=max(n1))==1))=x1; s2(find((n=min(n2))(n=max(n2))==1))=x2; x=s1+s2; axis([(min(min(n1),min(n2))-1),(max(max(n1),max(n2))+1),(min(x)-0.5),(max(x)+0.5)]) n2=-3:3; x2=[0 1 1 1 1 1 0]; n1=-2:3; x1=[0 0 3 2 1 0]; x3=-x1; [x,n]=jxl(x1,x2,n1,n2); x3=x; n3=n; n2=-2:3; x2=[0 0 3 2 1 0]; n=(n2(1)+n3(1)):(n2(1)+n3(1)+length(n2)+length(n3)-2) n = -5 -4 -3 -2 -1 0 1 2 3 4 5 6 f3=conv(x2,x3) f3 = 0 0 0 3 5 15 18 16 7 2 0 0 stem(n,f3,filled) title(f3=(f2k-f1k)*f3k) xlabel(n) 已知各连续信号的波形如下图所示，试用解析方法求下列卷积，并用Matlab绘出卷积信号波形，将其与解析计算结果进行比较 (2) 答：编辑函数edit gggfconv function[f,t]=gggfconv(f1,f2,t1,t2); %计算连续信号卷积积分实用函数 d=input( 请输入采样时间间隔d); f=conv(f1,f2); f=f*d; ts=t1(1)+t2(1); l=length(f1)+length(f2)-2; t=ts:d:(ts+l*d); plot(t,f); axis([min(t),max(t),min(f)-min(f)*0.2,max(f)+max(f)*0.2]) title(f(t)=f1(t)*f2(t)); （1）t1=-1:0.01:3; t2=0:0.01:4; f1