第一次实验：.doc

第一次实验： %实验一 program exp2_1 k=-2:5; x=[3,2,-1,0,2,-3,0,-1]; stem(k,x); xlabel(k); %实验一 program exp2_2 A=1;a=50;w0=50*pi; t=[0:0.002:1]; %fs=500Hz x=A*exp(-a*t).*sin(w0*t); stem(t,x); %实验一 program exp2_3 N=100;k=0:N; x=10*sin(0.02*pi*k)+randn(1,N+1); plot(k,x); xlabel(k);ylabel(x[k]); %实验一 program exp3 x=[1,2,1,1,0,-3]; kx=0:5; h=[1,-1,1]; kh=-1:1; y=conv(x,h); k=kx(1)+kh(1):kx(end)+kh(end); stem(k,y); xlabel(k);ylabel(y); %实验一 program exp4-1 b=[1 2 3 4 5 4 3 2 1]; a=[25]; w=linspace(-pi,pi,100); H=freqz(b,a,w); subplot(2,1,1) plot(w/pi,abs(H)) subplot(2,1,2) plot(w/pi,angle(H)) %实验一 program exp4-2 N=10; k=0:N; w0=0.02*pi; x=10*sin(w0*k)+randn(1,N+1); a=[1 2 3 4 5 4 3 2 1]; b=[25]; y=filter(a,b,x); subplot(2,1,1); stem(k,x); subplot(2,1,2); plot(k,y); %实验一 program exp4-3 N=10; k=0:N; w0=0.4*pi; x=10*sin(w0*k)+randn(1,N+1); a=[1 2 3 4 5 4 3 2 1]; b=[25]; y=filter(a,b,x); subplot(2,1,1); stem(k,x); subplot(2,1,2); plot(k,y); %实验一 program exp5 N=20; W1=0.1*pi; W2=0.4*pi; W3=0.7*pi; A=[2*cos(2*W1) 2*cos(W1);2*cos(2*W2) 2*cos(W2);2*cos(2*W3) 2*cos(W3)]; B=[0;1;0]; x=A\B; y=[x(1) x(2) 0 x(2) x(1)] k=0:N; x1=cos(W1*k); x2=cos(W2*k); x3=cos(W3*k); a=1; b=y; y=filter(b,a,x1+x2+x3); plot(k,y,r,k,x2,k--,k,x1+x2+x3,k: ); legend(y[k],x2[k],x1+x2+x3) y = -0.4667 0.3580 0 0.3580 -0.4667 第二次实验： %实验二 program exp1 N=10; k=-N:N ; x=cos(2*k*pi/N); b=x; a=[zeros(1,N) 1]; w=linspace(-pi,pi,512); H=freqz(b,a,w); plot(w/pi,abs(H)); hold on; M=length(x); X=fftshift(fft(x,M)); k1=(1-M)/M:2/M:(M-1)/M; stem(k1,abs(X),r); %实验二 program exp2-1 function fn=circonvt(x1,x2,N) if(length(x1)N|length(x2)N) error; end x1=[x1,zeros(1,N-length(x1))]; x2=[x2,zeros(1,N-length(x2))]; m=0:N-1; x=zeros(N,N); for n=0:N-1 x(:,n+1)=x2(mod((n-m),N)+1); end fn=x1*x; %实验二 program exp2-2 x=[1 -2 3 1]; h=[3 0 1 -1 2 1]; N=6; k=0:N-1; y=circonvt(x,h,N); Y=fft(y,N); subplot(2,1,1); stem(k*2*pi/N,Y); X=fft(x,N); H=fft(h,N); F=X.*H; subplot(2,1,2) stem(k*2*pi/N,F); ylabel(F); %实验二 p