DSP上机实验报告

（1）

>> A=[1 2 3 4];B=[3 4 5 6];

C=A+B; D=A-B; E=A.*B; F=A./B; G=A.^B;

subplot(2,4,1);stem(A,'.'); subplot(2,4,2);stem(B,'.'); subplot(2,4,3);stem(C,'.'); subplot(2,4,4);stem(D,'.'); subplot(2,4,5);stem(E,'.'); subplot(2,4,6);stem(F,'.'); subplot(2,4,7);stem(G,'.')

（2）

closeall;clearall;clc; n=0:15; x1n=0.8.^n;

x2n=exp((0.2+3*1i)*n);

x3n=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi); subplot(2,2,1);stem(n,x1n,'.');

subplot(2,2,2);stem(n,abs(x2n),'.'); subplot(2,2,3);stem(n,angle(x2n),'.'); subplot(2,2,4);stem(n,x3n,'.');

(6)>> b=[1,sqrt(2),1];a=[1,-0.67,0.9]; [h,w]=freqz(b,a);

am=20*log10(abs(h));

subplot(2,1,1); -1-，

plot(w,abs(h));xlabel('w');ylabel('幅频响应');title('系统响应') ph=angle(h); subplot(2,1,2);

plot(w,ph);xlabel('w');ylabel('相频响应'); （7）

>> a=[8 -2 -1 2 3];

b=[2 3 -1 -3]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c); xlabel('n');ylabel('幅度'); （8） >> N=50;

a=[1 -2]; b=[1 0.1 -0.06]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(a,b,x); stem(k,y); xlabel('n');ylabel('幅度 ');

实验二 信号的采样与重建

（5）

>> n1=0:1:15; p1=8; q1=2;

x=exp(-(n1-p1).^2/q1); n2=0:1:15; a=0.1; f2=0.0625;

y=(exp(-a*n2)).*sin(2*pi*f2*n2); N=length(x); n=0:N-1; n3=0:30; X=fft(x); Y=fft(y);

x32=[x zeros(1,16)]; y32=[y zeros(1,16)]; X32=fft(x32); Y32=fft(y32); z16=ifft(X.*Y);

z32=ifft(X32.*Y32); subplot(2,2,1); plot(n,z16,'-*'); xlabel('n'); ylabel('z(n)');

title('循环卷积结果'); subplot(2,2,2);

plot(n3,z32(1:2*N-1),'-o'); xlabel('n'); ylabel('z(n)');

title('线性卷积结果');

rm16=real(ifft(conj(X).*Y));

rm32_0=real(ifft(conj(X32).*Y32)); rm32=[rm32_0(N+2:2*N) rm32_0(1:N)]; m=n;

subplot(2,2,3); plot(m,rm16,'--'); xlabel('m'); ylabel('rm'); title('循环相关结果'); m=-(N-1):N-1; subplot(2,2,4); plot(m,rm32,'--o'); xlabel('m'); ylabel('rm');title('线性相关结果');