资源描述:
《哈工程 数字信号处理实验 实验一》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实验一基本信号一:实验原理使用MATLAB内部向量程序来产生信号。用MATLAB的stem指令绘出离散时间信号.用MATLAB的stem指令会出离散时间信号。依据MATLAB的编址约定,标号n=0必须对应nn(1);必须给指定向量的第一个参数以得到正确的n轴。二:实验内容2.正弦信号X【n】=Acos(ωn+ρ)使用MATLAB的向量功能求解此问题,将向量赋予余弦或正弦函数,再利用一个函数调用。在每种指定区间上展开并标注水平轴n轴。使用stem指令显示每个序列。A.X₁【n】=sin(πn/17)0≤n≤25B.X₂【n】=sin(πn/1
2、7)-15≤n≤25C.X₃【n】=sin(3πn﹢π/2)-10≤n≤10D.X₄【n】=sin(πn/√23)0≤n≤503.指数信号衰减的指数信号是数字信号处理中的基本信号,因为它是线性常系数方程的解。A.利用functions研究下面的MATLAB函数,看他如何产生离线指数信号。然后是用函数在区间n=0,1,2,3,...,20上绘出指数信号x【n】=(0.9)ⁿ。B.指数信号序列aⁿu【n】须在有限区间上求和。这个和以下面闭合时表示:A.指数序列在信号处理中常常出现的一个原因是,时移并不改变其信号特征。Y【n】=ay【n-1】,1
3、≤n≤L-1B.产生指数信号另外的方法是使用查分方程给出的递归表示式。当输入x【n】是一个冲击信号的时候,信号y【n】=aⁿu【n】是下面查分方程的解:y【n】-a【n-1】=x【n】,初始条件y【-1】=0三.实验程序2.正弦函数A.n=0:0.01:25y=sin(pi*n/17)plot(n,y)ylabel(‘y=sin(pi*n/17)’)gridgtext(‘n’)B.n=-15:0.01:25y=sin(pi*n/17)plot(n,y)ylabel(‘y=sin(pi*n/17)’)gridgtext(‘n’)C.n=-10
4、:0.01:10y=sin(3*pi*n+pi/2)plot(n,y)ylabel(‘y=sin(3*pi*n+pi/2)’)gridgtext(‘n’)D.n=0:0.01:50y=cos(pi*n/(23^0.5))plot(n,y)ylabel(‘cos(pi*n/(23^0.5)’)gridgtext(‘n’)3.指数信号Afunctiony=genexp(b,n0,L)%GENEXPgenerateanexponentialsignal:b^n%usage:Y=genexp(B,N0,L)%Binputscalargivingra
5、tiobetweenterns%N0startingindex(integer)%Llengthofgeberatedsignal%YoutputsignalY(1:L)if(L<=0)error('GENEXP:lengthnotpositive')endnn=n0+[1:L]'-1;%---vectorofindicesy=b.^nn;endn=0:9;x1=genexp(0.9,0,20)stem(n,x1,'b')sum(x1(:))B.functiony=signal(a,n0,L)if(L<=0)error('SIGNAL:le
6、ngthnotpositive')endnn=n0+[1:L]'-1y=(1-a.^nn)/(1-a)endu1=genexp(0.9,0,21)ss(1)=0;fori=1:19ss(i+1)=u1(i)+ss(i);endss(19)nn=[1:21]'-1;a=0.9y=(1-a.^nn)/(1-a)C.formatcompact,subplot(111)n=0:20;m=1:21;x1=genexp(0.9,0,21)subplot(211)stem(n,x1,'b')sum(x1(:))grid,title('Test1_3_3_
7、1')x2=0.9*genexp(0.9,0,21)subplot(212)stem(m,x2,'b')grid,title('Test1_3_3_2')四.结果分析1.正弦函数A.BCD3.A.B.x1=1.00000.90000.81000.72900.65610.59050.53140.47830.43050.38740.34870.31380.28240.25420.22880.20590.18530.16680.15010.13510.1216ans=8.9058用题目所给公式所求结果为8.9058,结果一致。C.