untitled

Size: px

Start display at page:

Download "untitled"

缺榴毕
5 years ago
Views:

1 DSP idsp 1

2 DSP DSP5416 codec 4 CPU ARM MP3 2

3 DSP JP1 DSP JTAG 6.JP2 CPLD JTAG S1 4 ON S1 MP/MC CLKMD1 CLKMD2 CLKMD3 S1 J1 J2 HPI16 HPI_EN HPI 3

4 DSP

5 DSP

6 DSP

7 DSP

8 DSP

9 DSP

10 DSP

11 DSP

12 DSP

13 DSP DSP DSP DSP C DSP CUNSI-DSP DSP

14 DSP alldsp4\simulator\sub DSP program/simulate 14

15 DSP

16 DSP

17 DSP X ( K) = N 1 N = 0 x( n) e j(2π / N ) nk X ( K) = N 1 N = 0 x[ n] W N nk ( n mn k ln W + )( + ) N = W nk N 17

18 DSP X ( K) = N 1 n= 0 x[ n] W nk N + N 1 n= 0 x[ n] W nk N N / 2 1 n= 0 x[2n] W 2nk N + N / 2 1 n= 0 x[2n + 1] W (2n+ 1) k N W k N W = = k N j(2π ) / N 2 j(2π ) /( N / 2) [ e ] = e WN / 2 X ( K) = N / 2 1 N / 2 1 nk k x1 [ n] WN / 2 + WN x 2[ n] n= 0 n= 0 k = X K) + W X ( ) 1( N 2 K X ( ) X ( ) x [ ] x [ ] 1 K 2 K 1 n 2 n W nk N / 2 18

19 DSP

20 DSP

21 DSP

22 DSP

23 DSP

24 DSP () FIR FIR short oflag = fir (DATA *x, DATA *h, DATA *r, DATA **dbuffer,ushort nh, ushort nx) ( fir.asm) x[nx] nx h[nh] nx nh=3: H = b0, b1, b2 k ( k LSB 0) k log 2 (nh) TMS320C54X DSP Reference Set r[nx] nx Dbuffer[nh ] 0 r h Nx Nh Oflag oflag = 1 32 oflag = 0 32 h FIR x d(nx=1) nh r[j] = = k 0 h[ k] x[ j k] 0 <=j <=nx 1 MATLAB FIR WAVE alldsp4\simulator\ hpassfir 24

25 DSP DSP 3 ViewGraphTime/Frequence x r Display Type Dual Time ;Interleaved Data Sources No ; Start Address-upper display x Start Address-upper display r Acquistition Buffer Size 256 Display Data Size 256DSP Data Type 16-bit signed integer OK 25

26 DSP FIR FIR FIR program/simulator/lpassfir 26

27 DSP IIR (1) IIR (3) DSP IIR (4) IIR 16KHz 2KHz IIR 1KHZ 3KHZ 1KHZ 3KHZ x[k] y[k] z y N [] k = ai x[ k i] + bi y[ k i] H M i= 0 i= 0 ( z) N=M ( z) Y X ( z) M k i= 0 = = N 1 b k = 0 a z k k z k 1 a0 + a1z + K + a z H ( z) = 1 1+ b z + K + b 1 N n N N z H ( z) = C N z z j= 1 z p j j z j p j N N Z b j 0 FIR FIR IIR 27

28 DSP

29 DSP y(n)= b0*x(n) + b1*x(n-1) + b2*x(n-2) - a1*y(n-1) - a2*y(n-2) 29

30 DSP

31 DSP c54x simulator 31

32 DSP (1) (2) Matlab (3) LMS (4) DSP (5) LMS matlab DSP CCS d(k) y(k) e(k) 1 1 x(k) y(k) y(k)d(k) e(k)e(k) x(k) e(k) 32

33 DSP H(Z) d(n) x(n) + e(n) _ W(Z) y(n) 2 FIR IIR IIR FIR FIR Transversal StructureSymmetric Transversal Structure Lattice Structure FIR 3 x(n) z -1 x(n-1) z -1 x(n-2) z -1 x(n-n-1) b0 b1 b2 b(n-1) w0(n) w1(n) w2(n) wn-1(n) y(n) 3 x(n) w(n) : w(n)={w(0),w(1),,w(n-1)} y(n) : y(n)=x(n)* w(n) ) y( n) = W T ( n) X ( n) = N 1 i= 0 w ( n) x( n i) i 33

34 DSP LMS Least Mean SquareLMS y(n) d(n)lms e(n)=d(n)y(n) MSEMean Square Error d(n) y(n) e(n) wi(n) LMS ε = E[ e 2 ( n)] = E[( d( n) y( n)) y(n) ε = E[ d 2 ( n)] + W T 2 ] ( n) RW ( n) 2W T ( n) P T R = E[ X ( n) X ( n)] = NN P = E[ d( n) X ( n)] N1 W W=R-1P R-1 P LMS W(n+1) W(n) W ( n + 1) = W ( n) u ( n) u (n) n LMS (n) E[ e 2 ( n)] 2 E[ e ( n)] ( n) = = 2E[ e( n) X ( n)] W ( n) W W (0) W ( n + 1) = W ( n) u ( n) W(n+1) W (n) (n) E[ ] 1 E{ e( n) X ( N)} = K K 1 i= 0 e( n i) X ( n i) 34

35 DSP K 2e( n) X ( n) 2E{ e( n) X ( n)} Widrow-Hoff LMS W ( n + 1) = W ( n) + 2ue( n) X ( n) LMS Widrow-Hoff LMS W ( n + 1) = W ( n) + 2ue( n)[ X ( n) + X ( n N + 1)] T X(n) [ x ( n), x( n 1),..., x( n N / 2 + 1)] X(n-N+1) [ x ( n N + 1), x( n N + 2),..., x( n N / 2)] Matlab Matlab T clear all; clf; signal = sin(2*pi*0.055*[0:1000-1]'); % noise=randn(1,128); % nfilt=fir1(11,0.4); % 11 FIR fnoise=filter(nfilt,1,noise); % d=signal.'+fnoise; % w0 = nfilt.' -0.01; % mu = 0.015; % s=initse(w0,mu); % adaptse [y,e,s] = adaptse(noise,d,s); plot(0:127,e(1:128)); DSP DSP c c #include <math.h> #include "tms320.h" #include "dsplib.h" 35

36 DSP #define NX 128 #define NH 12 #define STEP 278 #pragma DATA_SECTION (h,".coefs") #pragma DATA_SECTION (dbuffer,".dbuffer") float x1[nx] = { , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , }; float des1[nx] ={0, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , 36

37 DSP , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , }; DATA x[nx]; DATA dbuffer[nh]; DATA h[nh]; DATA r[nx]; DATA des[nx]; DATA *dp = dbuffer; short i; main() { //tranform data from float to fixed fltoq15(x1,x,nx); fltoq15(des1,des,nx); // clear for (i=0;i<nh;i++) h[i] =0; // clear coeff buffer (optional) for (i=0;i<nx;i++) r[i] =0; // clear output buffer (optional) for (i=0; i<nh; i++) dbuffer[i] = 0; // clear delay buffer (a must) // compute dlms(x,h,r,&dp,des,step, NH, NX); } return; c NX 128 NH 12 STEP h dbuffer.cofes.dbuffer x1 des1 r[nx] x1 des1 matlab fltoq15(x1,x,nx) fltoq15(des1,des,nx) hrdbuffer 0 dlms DSP dlms 37

38 DSP ;**************************************************************** ; Function: dlms ; Version: 1.00 ; Description: dlms fir filter ; ; Copyright Texas instruments Inc, 1998 ; ; Revision History: ; 1.00, R. Piedra, 8/31/98 - Original release. Started from code by ; Padma Mallela. ;****************************************************************.mmregs ; Far-mode adjustment.if far_mode offset.set 1 ; far mode uses one extra location for ret addr ll.else offset.set 0.endif.asg (0), mu_error ; 2*mu*error(i) variable.asg (1), save_st1 ; stack description.asg (2), save_st0.asg (3), save_ar6.asg (4), save_ar1.asg (5), ret_addr ; x in A.asg (6 + offset), arg_h.asg (7 + offset), arg_y.asg (8 + offset), arg_d.asg (9 + offset), arg_des.asg (10 + offset), arg_rate.asg (11 + offset), arg_nh.asg (12 + offset), arg_n ; register usage ; AR0 : circ addr. index.asg ar1, ar_count.asg ar2, ar_d.asg ar3, ar_h.asg ar4, ar_des.asg ar5, ar_x.asg ar6, ar_y 38

39 DSP ;*****************************************************************************.def _dlms _dlms: ; Preserve registers ; pshm ar1 pshm ar6 PSHM ST0 ; 1 cycle PSHM ST1 ; 1 cycle RSBX OVA ; 1 cycle RSBX OVB ; 1 cycle ; Preserve space for local variables ; frame -1 ; mu_error local var ; Set math and overflow modes ; ssbx CUNSIm ; sign extension on ssbx frct ; fract on ld #0,ASM ; clear ASM (used by st mpy) st #0, *sp(mu_error) ; clear error OJO : make argument ; ; Get arguments ; stlm a, ar_x ; pointer to x mvdk *sp(arg_h),*(ar_h) ; pointer to h mvdk *sp(arg_y),*(ar_y) ; pointer to y mvdk *sp(arg_d),*(ar_d) ; de-referencing mvdk *ar_d,*(ar_d) ; pointer to d mvdk *sp(arg_des),*(ar_des) ; for circ addr mvdk *sp(arg_n), *(ar_count) ; sample counter mar *ar_count- ; ar_count= nsamples-1 ld *sp(arg_nh),a ; a = nh ; Initialize registers ; stlm a,bk stm #1,AR0 sub #03,a _start: ; bk = circular buffer size = nh 39

40 DSP stl a, *sp(arg_nh) ; Loop through remaining (n-1) samples ; next_sample: mvdk *sp(arg_nh),brc ld *sp(mu_error), t ; t = 2 * mu * error (i) sub b,b ; b = 0 ; Loop through (nh-2) taps: compute filter and adapt at = time ; rptbd lms_end -1 mpy *ar_d+0%,a ; a = mu_error * x(1) ; point to x(2) lms *ar_h, *ar_d ; b = b + h(1) * x(2) ; a = a + h(1) lms_beg st a, *ar_h+0%; ; update h(1) ; point to h(2) mpy *ar_d+0%,a ; a = mu_error * x(2) ; point to x(0) lms *ar_h, *ar_d ; b = b + h(2) * x(0) ; a = a + h(2) lms_end ; Perform 1st tap outside to avoid that the oldest ; sample is over written before updating the h(0) coef ; st a, *ar_h+0% ; update last coef mpy *ar_d+0%,a ; error * x(0) mvdd *ar_x+, *ar_d ; new sample --> dbuffer x(0) lms *ar_h, *ar_d+0% ; b = b * h(0) * x(1) ; Write output sample ; sth b, *ar_y+ ; Compute new error ; st a,*ar_h+0% sub *ar_des+,a ; bh = error(i)= des(i)-y(i) 40

41 DSP banzd next_sample, *ar_countmpya *sp(arg_rate) ; b = step_size * error(i) sth b, *sp(mu_error) ; save as local variable ; Return ; _end.asg ar_h, ar_temp ; ar_h not used any more mvdk *sp(arg_d), *(ar_temp) mvkd *(ar_d), *ar_temp ; update new ar_d ; Return overflow flag ; ld 0,a xc 1,AOV ld #1,a ; Adjust stack pointer and restore registers ; frame +1 ; adjust for local vars POPM ST1 ; 1 cycle POPM ST0 ; 1 cycle popm ar6 popm ar1.if far_mode fretd.else retd.endif nop nop -c -l rts.lib -stack 0x200 -heap 0x200 MEMORY { PAGE 0: INT_PM_DRAM: origin = 3800h, length = 80h EXT_PM_RAM: origin = 3880h, length = 0700h PAGE 1: INT_DM_SCRATCH_PAD_DRAM: origin = 060h, length = 20h 41

42 DSP INT_DM_1: origin = 0080h, length = 1000h INT_DM_2: origin = 1080h, length = 380h EXT_DM_RAM: origin = 1400h, length = 02000h } SECTIONS {.text : {} > EXT_PM_RAM PAGE 0.cinit : {} > EXT_PM_RAM PAGE 0.switch : {} > EXT_PM_RAM PAGE 0.dbuffer: {} > INT_DM_1 PAGE 1.stack : {} > INT_DM_1 PAGE 1.bss : {} > INT_DM_1 PAGE 1.const : {} > INT_DM_1 PAGE 1.sysmem: {} > INT_DM_1 PAGE 1.cio : {} > INT_DM_1 PAGE 1.coefs : {} > INT_DM_2 PAGE 1 } 42

43 DSP

44 DSP c54x simulator program/simulator/lms 44

45 DSP DSP 54X DSP C c54x simulator oflag = short convol (DATA *x, DATA *h, DATA *r, ushort nh, ushort nr) ( convol.asm) x[nr+ nh-1] nr+nh-1 h[nh] r[nr] nr nh oflag nh nr oflag = 1 32 oflag = 0 32 r nr (nr = 1) r[j] = nh k = 0 h[ k] x[ j k] 0 <=j <=nr 1 CCS2.2Project/Open program\ Simulatort\ convol\ test.pjt test.pjt Source Project\ Rebuild All Debug Go Main, tesst.out 45

46 DSP File\Data\Load 500Hz1000z x,y 256 View/Graph/Time/Frequency xyr,

47 DSP Hz 47

48 DSP

49 DSP

50 DSP DSP 54X DSP C c54x simulator short oflag = acorr_raw (DATA *x, DATA *r, ushort nx, ushort nr) short oflag = acorr_bias (DATA *x, DATA *r, ushort nx, ushort nr) short oflag = acorr_unbias (DATA *x, DATA *r, ushort nx, ushort nr) ( araw.asm, abias.asm, aubias.asm) short oflag = acorr (DATA *x, DATA *r, ushort nx, ushort nr, type) x [nx] r [nr] nx nr oflag type rawbias unbias nx nr nr <= nx r x oflag = 1 32 oflag = 0 32 nr 0 nr r[j] = nx j 1 k = 0 r[j] = 1/nx x[ j + k] x[ k] nx j 1 k = 0 x[ j + k] x[ k] 0 <=j <=nr 0 <=j <=nr r[j] = 1/(nx-abs(j)) nx j 1 k = 0 x[ j + k] x[ k] 0 <=j <=nr 50

$2Project/Open program\ Simulatort\ acorr\ test.$ $pjt Source Project\ Rebuild All Debug Go Main, test.$

51 DSP CCS2.2Project/Open program\ Simulatort\ acorr\ test.pjt test.pjt Source Project\ Rebuild All Debug Go Main, test.out 2File\Data\Load 500Hz x, 256 View/Graph/Time/Frequency x,r,

52 DSP

53 DSP DSP 54X DSP C CCS short oflag = corr_raw (DATA *x, DATA *y, DATA *r, ushort nx, ushort ny) short oflag = corr_bias (DATA *x, DATA *y, DATA *r, ushort nx, ushort ny) short oflag = corr_unbias (DATA *x, DATA *y, DATA *r, ushort nx, ushort ny) ( craw.asm, cbias.asm, cubias.asm) short oflag = corr(data *x, DATA *y, DATA *r, ushort nx, ushort nytype) type rawbias unbias x [nx] 1 nx y [ny] r[nx+ny-1] nx ny oflag 2 ny nx+ny oflag = 1 32 oflag = 0 32 x y r r[j] = nr j 1 k= 0 r[j] = 1/nr x[ j + k]* y[ k] nr j 1 k = 0 r[j] = 1/(nx abs(j)) x[ j + k]* y[ k] nr j 1 k = 0 x[ j + k]* y[ k] 0 <=j <=nr = nx+ny-1 0 <=j <=nr = nx+ny-1 0 <=j <=nr = nx+ny-1 53

54 DSP CCS2.2Project/Open program\ Simulatort\ corr\ test.pjt test.pjt Source Project\ Rebuild All Debug Go Main, tesst.out 2File\Data\Load 500Hz xy,

55 DSP y 3 View/Graph/Time/Frequency x,y, X 1y 2 r r 512 x,y 55

56 DSP r 56

57 DSP

58 DSP

59 DSP

60 DSP DSP McBSP0McBSP2 AIC23 MP3 LED3 60

61 DSP DSP DSP CUNSI-DSP J6 J7, TLC320AIC23 TI AD DA 96Kbps, IC manual DSP Mcbsp0 Mcbsp2 AIC23 program /yuyinfirlpass 61

62 DSP YUV411 DSP CUNSI-DSP DSP, f(x)=255-f(x) 1 program\yuv1 test.pjt,ccs Debug go main 2 DSP picture ok.dat, OK 62

63 DSP

64 DSP CCS 1 lines per Display

65 DSP

66 DSP YUV411 DSP CUNSI-DSP DSP f1 0 f f(x)<100 fx= 255 f(x)>100 1 program\yuv2 test.pjt,ccs, Debug go main 2 DSP picture ok.dat, OK 66

67 DSP

68 DSP CCS 1 lines per Display

69 DSP

70 DSP DSP CUNSI-DSP DSP program\yuv6 test.pjt,ccs Debug go main 6 DSP picture flowers.dat, OK 70

71 DSP

72 DSP CCS 1 lines per Display

73 DSP

$com DSP CUNSI-DSP DSP S1S2 S1S2 9 program\yuv5 test.$

74 DSP DSP CUNSI-DSP DSP S1S2 S1S2 9 program\yuv5 test.pjt,ccs Debug go main 10 DSP picture flower.dat, OK 74

75 DSP

76 DSP CCS 1 lines per Display

77 DSP

78 DSP DSP5416 DSP CUNSI-DSP DSP 0 s1 s2 s3 255 s1 s2 s1s2t1t2 13 program/yuv4 test.pjt,ccs Debug go main 14 DSP picture 78

79 DSP flower.dat, OK 15 79

80 DSP CCS 1 lines per Display

81 DSP

82 DSP dsp DSP CUNSI-DSP DSP 1 f 1,0 f 2 82

83 DSP DSP Hz 1000Hz 83

84 DSP Graphical Display program\ Graphical Display 84

85 DSP

86 DSP

87 DSP (PAM) dsp (PAM) DSP CUNSI-DSP DSP STt ms(t) LPF 87

88 DSP STt 88

89 DSP Pω k m(t) Q Q15 =A*SINxA SINx DSP Graphical Display 89

90 DSP Graphical Display 90

91 DSP

92 DSP

93 DSP (PAM) (PAM) dsp DSP CUNSI-DSP DSP PAM PAM =A*SINxA SINx APAM *SINx A*SINx*SINx=A*1-COS2X/2=A/2-A*COS2X/2 AA*COS2X/2A DSP A*SINx SINx FIR FIR DSP A A A A*SINx PAM 93

94 DSP A*SINx*SINx program\ Graphical Display 94

95 DSP

96 DSP

97 DSP

98 DSP (QPSK) dsp (QPSK) DSP CUNSI-DSP DSP QPSK QPSK - QPSK QPSK 98

99 DSP QPSK QPSK 4 90 DSP QPSK 99

100 DSP program\ Graphical Display 100

101 DSP

102 DSP

103 DSP dsp CCS DSP CUNSI-DSP E J41 J27 P1 P2 P2 DSP 103

104 DSP OK 1 104

C++ 程式設計

C++ 程式設計 C C 料, 數, - 列串理列 main 數串列什 pointer) 數, 數, 數數省不不, 數 (1) 數, 不數 * 料 * 數 int *int_ptr; char *ch_ptr; float *float_ptr; double *double_ptr; 數 (2) int i=3; int *ptr; ptr=&i; 1000 1012 ptr 數, 數 1004