[理学]DSP 课件.ppt

MPYA *AR2+ ；B=x2（1-x2/42（1-x2/72）） ST B,*AR2 ；（d_temp）=x2（1-x2/42（1-x2/72）） ‖LD *AR4,B ；B=1 MASR *AR2-,*AR3+,B,A ；A=1-x2/20（1-x2/42（1-x2/72）） MPYA *AR2+ ；B=x2（1-x2/20（1-x2/42（1-x2/72））） ST B,*AR2 ；（d_temp）=B=x2（1-x2/20（1-x2/42（1-x2/72））） ‖LD *AR4,B ；B=1 MASR *AR2-,*AR3+,B,A；A=1-x2/6（1-x2/20（1-x2/42（1-x2/72））） MPYA d_x ；B=x（1-x2/6（1-x2/20（1-x2/42（1-x2/72）））） STH B,d_sinx ；sin（theta） RET .end 2. 用z变换计算余弦值 下面采用z变换和反z变换来求cos（ x ），其的方法为： 根据z变换定义，序列x(n) 的变换公式为： 反z变换为x(n)=￡-1[X(z)]可用留数法、部分分式展开法和长除法求得。 则cos(ω0n)u(n)的z变换为： H（z）=￡[cos(n)u(n)]=　￡[ejω0n+e-jω0nu(n)] 其中C＝-cosω0，A=2cosω0T，B=-1。ω0为余弦输出信号的频率，T为离散余弦序列的采样频率。 如果以该函数设计一离散时间系统，为则其单位冲击响应就是余弦输出信号。此时的输出序列Y（k）为H（z）的反z变换。 Y（k）=￡-1[H（z）]=AY[k-1]+BY[k-2]+X[k]+CX[k-1] 当k=-1时　Y（k）=Y（-1）=AY[-2]+BY[-3]+X[-1]+CX[-2]=0 当k=0时　 Y（k）=Y（0）=AY[-1]+BY[-2]+X[0]+CX[-1]=0+0+1+0=1 当k=1时　 Y（k）=Y（1）=AY[0]+BY[-1]+X[1]+CX[0]=1+0+0+C=A+C 当k=2时　 Y（k）=Y（2）=AY[1]+BY[0]+X[2]+CX[1]=AY[1]+BY[0] 当k=3时　 Y（k）=Y（3）=AY[2]+BY[1]+X[3]+CX[2]=AY[2]+BY[1] 当k=n时　 Y（k）=Y（n）=AY[n-1]+BY[n-2] 　　在k＞2以后，Y（k）能用Y[k-1]和Y[k-2]算出，这是一个递归的差分方程。如果按第七章实验二的方式产生余弦信号，对应的初始化程序为： 初始化y[1]和y[2]： 　SSBX　 FRCT　 　　 ；置FRCT=1，准备进行小数乘法运算 　ST #INIT_A,AA ；将常数A装入变量AA 　ST #INIT_B,BB ；将常数B装入变量BB 　ST #INIT_C,CC ；将常数C装入变量CC PSHD CC ；将变量CC压入堆栈 POPD Y2 ；初始化Y2=CC LD 　AA, A　 　；装AA到A累加器 ADD CC, A ； A累加器=AA+CC STH A,Y2 ；Y2=Y[1]=AA+CC LD AA,T ；装AA到T寄存器 MPY Y2,A 　 ；Y2乘系数A，结果Y[1]*AA放入A累加器 ADD BB, A ； A累加器=Y[2]=Y[1]*AA+BB*Y[0] STH A,Y1 　；将A累加器中Y[2]的高16位存入变量Y1=Y[2] 以后的递推过程由中断服务程序完成Y[3]到Y[n]运算，相应的程序片段为： LD BB,T ；将系数B装入T寄存器 MPY Y2,A ；Y2乘系数B，结果BB*Y[1]放入A累加器 LTD Y1 ；将Y1=Y[2]装入T，同时复制到Y2，Y[2]退化为Y[1] MAC AA,A ；完成新余弦数据的计算，A累加器中为 ；Y1*AA+Y2*BB或Y[3]＝AA*Y[2]+BB*Y[1] STH A,1,Y1 ；将新数据存入Y1，因所有系数都除过2， ；所以在保存结果时左移一位，恢复数据正常大小。 STH A,1,Y0 ；将新正弦数据存入Y0 第四节