实验参考程序.doc

实验一参考程序： 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＝10点的单位阶跃序列 function unit_step(N) % unit_step.m N=10; 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]) 观察并分析采用不同频率时，对函数的频谱影响。 　　　　（a）：以，对其进行采样得到x1(n)。 　　　　（b）：以，对其进行采样得到x2(n) 　 （c）：以，对其进行采样得到x3(n) a=218.2; b=50*pi; fs=1000; %采样频率1000hz T=1/fs; Tp=50*0.001; %观察时间50微秒 M=Tp*fs; %采样点数 n=0:M-1 y=a*exp(-b*n*T).*sin(b*n*T ) %函数表达式 subplot(3,2,1) stem(n,y) xlabel(n);ylabel(xa(nT));title(fs=1000HZ); axis([0,M,-2,1.2*max(abs(y))]) yk=T*fft(y,M) %M点FFT K=0:M-1; fk=K/Tp; subplot(3,2,2) plot(fk,abs(yk)) xlabel(f(Hz));ylabel(幅度);