信号处理第一章实验报告.doc

2-1 R = 50; m = 0:R-1; s = 2*m.*(0.9.^m); % Generate the uncorrupted signal d = rand(R,1)-0.5; % Generate the random noise x1 = s+d; for n = 1:50; d = rand(R,1)-0.5; x = s + d; x1 = x1 + x; end x1 = x1/50; stem(m,d); xlabel(Time index n);ylabel(Amplitude); title(Noise); pause stem(m,x1); xlabel(Time index n);ylabel(Amplitude); title(Ensemble average); 2-2 (1) % Program 2_2 % Generation of complex exponential sequence % a = input(Type in real exponent = ); b = input(Type in imaginary exponent = ); c = a + b*i; K = input(Type in the gain constant = ); N = input (Type in length of sequence = ); n = 1:N; x = K*exp(c*n);%Generate the sequence stem(n,real(x));%Plot the real part xlabel(Time index n);ylabel(Amplitude); title(Real part); disp(PRESS RETURN for imaginary part); pause stem(n,imag(x));%Plot the imaginary part xlabel(Time index n);ylabel(Amplitude); title(Imaginary part); Type in real exponent = -1/12 Type in imaginary exponent = 3.14/6 Type in the gain constant = 1 Type in length of sequence = 40 PRESS RETURN for imaginary part (2) a = input(Type in real exponent = ); c = a; K = input(Type in the gain constant = ); N = input (Type in length of sequence = ); n = 1:N; x = K*c.^n;%Generate the sequence stem(n,real(x));%Plot the real part xlabel(Time index n);ylabel(Amplitude); title(Real part); disp(PRESS RETURN for imaginary part); Type in real exponent = 1.2 Type in the gain constant = 0.2 Type in length of sequence = 30 Type in real exponent = 0.9 Type in the gain constant = 20 Type in length of sequence = 30 (3) a = input(Type in w = ); c = a; K = 1.5; N = 40; n = 1:N; x = K*cos(c*n);%Generate the sequence stem(n,x);%Plot the real part xlabel(Time index n);ylabel(Amplitude); title(Real part); disp(PRESS RETURN for imaginary part); Type in w = 0 a = input(Type in w = ); c = a; K = 1.5; N = 40; n = 1:N; x = K*cos(c*n);%Generate the sequence stem(n,x);%Plot the real part xlabel(Time index n);ylabel(Amplitude); title(Real part); disp(PRESS RETURN fo