实验参考程序 数字信号处理.docVIP

实验一参考程序： 1 产生10点的单位抽样序列δ(n)； %function unit_pulse(N) % unit_pulse.m N=10； x=zeros(1,N); x(1)=1; n=0:N-1; figure(1); stem(n,x);xlabel(单位抽样序列) axis([-1 20 0 1.1]) 产生10点同时移位3位的单位抽样序列δ(n-3)； %function shift_unit_pulse(N,k) % shift_unit_pulse.m N=10; k=3; x=zeros(1,N); x(k+1)=1; n=0:N-1; figure(2); stem(n,x);xlabel(移位3位的单位抽样序列) axis([-1 20 0 1.1]) 或function [x, n]=i shift_unit_pulse (n0,ns,nf) n=[0:9]； x=[(n-3)==0] 产生任意序列 f(n)=8δ(n)+7δ(n-1)+6δ(n-2) +5δ(n-3)+ 4δ(n-4)+7δ(n-5)； %function arbitrary_pulse(N) % arbitrary_pulse.m N=10 x=zeros(1,N); x(1)=8;x(2)=7;x(3)=6;x(4)=5;x(5)=4;x(6)=7; n=0:N-1; figure(3); stem(n,x);xlabel(任意序列f(n)) axis([-1 20 0 9]) 产生N＝点的单位阶跃序列 %function unit_step(N) % unit_step.m N=32; x=ones(1,N); n=0:N-1; figure(4); stem(n,x);xlabel(单位阶跃序列) axis([-1 32 0 1.1]) 产生斜率为3，n0=4,点数为20点的斜坡序列g(n)=B(n-n0) %function slope(N,k,B) % slope.m N=20;k=4;B=3; x=[zeros(1,k) ones(1,N-k)]; for i=1:N x(i)=B*x(i)*(i-k); end n=0:N-1; figure(5); stem(n,x);xlabel(斜坡序列) axis([-1 10 0 90]) 产生幅度A=3，频率f＝100，初始相位＝1.2，点数为32点的正弦序列。 %function sine(N,A,f,fai) % sine.m A=3;f=100;fai=1.2;N=32; n=0:N-1; x=A*sin(2*pi*f*(n/N)+fai); figure(6); stem(n,x);xlabel(正弦序列) axis([-1 32 -3.2 3.2]) 产生幅度A=3，角频率ω＝314，点数为32点的复正弦序列。 %function complex_sine(N,A,w) % complex_sine.m A=3; w=314;N=32 n=0:N-1; x=A*exp(j*w*n); figure(7); stem(n,x);xlabel(复正弦序列) axis([-1 32 -3.2 3.2]) 产生幅度A=3，a＝0.7，点数为32点的实指数序列。 %function real_exponent(N,A,a) % real_exponent.m a=0.7;A=3; n=0:N-1; x=A*a.^n; figure(8); stem(n,x);xlabel(实指数序列) axis([-1 32 0 3.2]) 实验二 参考程序： 求序列的谱分析。160长度更准确，整周期采样估计更准确。 clear all close all clc %序列周期为r(2*pi/w)=r(2/0.225)=80/9 %选定序列的长度为80点可准确估计 N=160; n=0:N-1; x=cos(0.225*pi*n);%生成信号 y=fft(x); %求信号的频谱 subplot(2,1,1) %画图 plot(n,x) xlabel(n) ylabel(x(n)) subplot(2,1,2) plot(n/N*2,abs(y)) xlabel(nomorlized frequency) ylabel(Y(exp(j\omega))) xlim([0,1]) figure; subplot(2,1,1) stem(n,x) xlabel(n) ylabel(x(n)) subplot(2,1,2) stem(n,abs(y)) xlabel(k) ylabel(Y(k))) N=256; %点数不同，估计准确度不同 n=0:N-1; x=c