实验一 信号 数字信号处理实验报告
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
七程序及图像
1.分析理想采样信号序列的特性。
1.理想采样信号序列
>> n=0:50;A=444.128;a=50*sqrt(2)*pi;
>> w0=50*sqrt(2)*pi;T=1/1000;
>> x=A*exp(-a*n*T).*sin(w0*n*T);
>> k=-25:25; W=(pi/12.5)*k; f=(1/25)*k*1000;
>> X=x*(exp(-j*pi/12.5)).^(n'*k);
>> magX=abs(X);angX=angle(X);
>> subplot(3,1,1);stem(x);title('x的时域')
>> subplot(3,1,2);stem(f,magX);title('x的幅度谱') >> subplot(3,1,3);stem(f,angX);title('x的相位谱')
>> T=1/300;
>> x=A*exp(-a*n*T).*sin(w0*n*T);
>> magX=abs(X);angX=angle(X);
>> subplot(3,1,1);stem(x);title('x的时域')
>> subplot(3,1,2);stem(f,magX);title('x的幅度谱') >> subplot(3,1,3);stem(f,angX);title('x的相位谱')
>> T=1/200;
>> x=A*exp(-a*n*T).*sin(w0*n*T);
>> X=x*(exp(-j*pi/12.5)).^(n'*k);
>> magX=abs(X);angX=angle(X);
>> subplot(3,1,1);stem(x);title('x的时域')
>> subplot(3,1,2);stem(f,magX);title('x的幅度谱') >> subplot(3,1,3);stem(f,angX);title('x的相位谱')
没有出现“混淆”现象，采样频率均大于信号频率的二倍。
2.离散信号，系统和系统响应的分析
（1）>> n=0:50;xb=[1 zeros(1,50)];k=-25:25;
>> Xb=xb*(exp(-j*pi/25)).^(n'*k);
>> magX=abs(Xb);angX=angle(Xb);
>> subplot(3,1,1);stem(xb);title('xb的时域')
>> subplot(3,1,2);stem(magX);title('xb的幅度谱')
>> subplot(3,1,3);stem(angX);title('xb的相位谱')
>> n=1:50; hb=zeros(1,50);
>> hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;k=-25:25; >> X=hb*(exp(-j*pi/25)).^(n'*k);
>> magX=abs(X);angX=angle(X);
>> subplot(3,1,1);stem(hb);title('hb的时域') >> subplot(3,1,2);stem(magX);title('hb的幅度谱') >> subplot(3,1,3);stem(angX);title('hb的相位谱')
>> y=conv(xb,hb);
>> subplot(1,1,1);stem(y);title('y的时域')
>> n=0:99;
>> k=-50:50;
>> X=y*(exp(-j*pi/25)).^(n'*k);
>> magX=abs(X);
>> angX=angle(X);
>> subplot(3,1,2);stem(magX);title('y的幅度谱') subplot(3,1,3);stem(angX);title('y的相位谱')
(2)>> n=0:50;xc=[ones(1,10) zeros(1,41)];
>> subplot(3,1,1);stem(n,xc);title('xc的时域') >> axis([0 50 0 1.2]);
>> k=-25:25;
>> X=xc*(exp(-j*pi/25)).^(n'*k);
>> magX=abs(X);angX=angle(X);
>> subplot(3,1,2);stem(magX);title('xc的幅度谱') >>subplot(3,1,3);stem(angX);title('xc的相位谱') >> ha=[ones(1,10) zeros(1,41)];
>> y=conv(xc,ha);subplot(3,1,1);
>>stem(y);title('线性卷积y的时域')
>> n=0:100;k=-50:50;
>> X=y*(exp(-j*pi/25)).^(n'*k);
>> magX=abs(X);angX=angle(X);
>> subplot(3,1,2);stem(magX);title('y的幅度谱') >>subplot(3,1,3);stem(angX);title('y的相位谱')
(3)
>> n=0:50;A=1;a=0.4;w=2.0734;T=1;
>> x=A*exp(-a*n*T).*sin(w*n*T);
>> k=-25:25;W=(pi/12.5)*k;
>> f=(1/25)*k*1000;
>> X=x*(exp(-j*pi/12.5)).^(n'*k);
>> magX=abs(X);angX=angle(X);
>> subplot(3,1,1);stem(x);title('x的时域')
>> subplot(3,1,2);stem(f,magX);title('x的幅度谱') >> subplot(3,1,3);stem(f,angX);title('x的相位谱')
>> ha=[ones(1,10) zeros(1,41)];
>> y1=conv(x,ha);
>> subplot(3,1,1);stem(y1);title('线性卷积y1的时域') >>n=0:100;k=-50:50;
>>X=y1*(exp(-j*pi/25)).^(n'*k);
>>magX=abs(X);angX=angle(X);
>>subplot(3,1,2);stem(magX);title('y1的幅度谱') >>subplot(3,1,3);stem(angX);title('y1的相位谱')
>> n=0:50;A=1;a=0.1;w=2.0734;T=1;
>>x=A*exp(-a*n*T).*sin(w*n*T);
>> ha=[ones(1,10) zeros(1,41)];
>>y2=conv(x,ha);
>>subplot(3,1,1);stem(y2);title('线性卷积y2的时域') >>n=0:100;k=-50:50;X=y2*(exp(-j*pi/25)).^(n'*k); >>magX=abs(X);angX=angle(X);
>>subplot(3,1,2);stem(magX);title('y2的幅度谱') >>subplot(3,1,3);stem(angX);title('y2的相位谱')
>> n=0:50;A=1;a=0.1;w=1.2516;T=1;
>>x=A*exp(-a*n*T).*sin(w*n*T);
>> ha=[ones(1,10) zeros(1,41)];
>>y3=conv(x,ha);
>>subplot(3,1,1);stem(y3);title('线性卷积y3的时域') >>n=0:100;k=-50:50;X=y3*(exp(-j*pi/25)).^(n'*k); >>magX=abs(X);angX=angle(X);
>>subplot(3,1,2);stem(magX);title('y3的幅度谱') >>subplot(3,1,3);stem(angX);title('y3的相位谱')