DSPs原理及应用 CCS实验报告

DSPs原理及应用——实验报告实验一第17章利用CCS开发DSP程序程序V olume.c#include <stdio.h>#include "volume.h"/* Global declarations */int inp_buffer[BUFSIZE]; /* processing data buffers */int out_buffer[BUFSIZE];int gain = MINGAIN; /* volume control variable */ unsigned int processingLoad = BASELOAD; /* processing routine load value */ struct PARMS str ={2934,9432,213,9432,&str};/* Functions */extern void load(unsigned int loadValue);static int processing(int *input, int *output);static void dataIO(void);/** ======== main ========*/void main(){int *input = &inp_buffer[0];int *output = &out_buffer[0];puts("volume example started\n");/* loop forever */while(TRUE){/** Read input data using a probe-point connected to a host file.* Write output data to a graph connected through a probe-point.*/dataIO();#ifdef FILEIOputs("begin processing") /* deliberate syntax error */#endif/* apply gain */processing(input, output);}}/** ======== processing ========** FUNCTION: apply signal processing transform to input signal.** PARAMETERS: address of input and output buffers.** RETURN V ALUE: TRUE.*/static int processing(int *input, int *output){int size = BUFSIZE;while(size--){*output++ = *input++ * gain;}/* additional processing load */load(processingLoad);return(TRUE);}/** ======== dataIO ========** FUNCTION: read input signal and write processed output signal.** PARAMETERS: none.** RETURN V ALUE: none.*/static void dataIO(){/* do data I/O */return;}输入波形输出波形此时Gain值为1.输入波形输出波形Gain值为5Gain值为10时，输入输出波形如下：实验二第18章6701核心系统的定浮点运算程序的设计与分析1.定点小数乘法product.asm 代码如下：.def start.data.def _c_int00_c_int00multiplicand .word 5*32768/10multiplier .word 5*32768/10.bss result,4,1,1.textstart: MVKL multiplicand,A11MVKH multiplicand,A11MVKL multiplier,B11MVKh multiplier,B11LDW *A11,A12LDW *B11,B12NOP 5MPYU A12,B12,A11MPYHLU A12,B12,B11MPYH A12,B12,B12MPYLHU A12,B12,A12NOPADD B11,A12,B11SHR B11,16,A12SHL B11,16,B11ADD A11,B11,A11ADD A12,B12,A12MVKL result,B11MVKH result,B11STW A11,*B11++STW A12,*B11NOP 5LOOP: B LOOPNOP 5.endproduct.cmd 代码如下：MEMORY{ IPRAM:origin=0x00008000,len=0x01000 IDRAM: origin=0x80000000,len=0x10000 } SECTIONS{ .text : > IPRAM.data : > IDRAM.bss : > IDRAM }反汇编窗口打开memory windows根据.cmd文件的内容做设置此时存储器窗口显示的图如下Q值改为15 Q值改为30 2．64位整数加减法程序加法程序：.global_add64;.text_add64:ADDU.L1 A4,A6,A1:A0;||ADD.L2 B4,B6,B7;|| ADD.S2X A8,4,B0;ADD.L2X A1,B7.B7;STW.D1 A0,*A8;|| STW.D2 B7,*B0;||B.S2 B3NOP 5.endPC改为0x00008020，加数和被加数A4、A6改为0x12345678,0xf2345678减法程序：.global _sub64.text_sub64:SUBU.L1 A4,A6, A1:A0 || SUB.L2 B4,B6,B7|| SUB.S2X A8,4,B0SUB.L2X A1,B7,B7STW .D1 A0,* A8|| STW .D2 B7,* B0|| B.S2 B3NOP 5.endPC改为0x000080003. 浮点乘除法运算程序浮点运算程序如下：#include "math.h"#pragma DATA_SECTION(data1,"data_buf1") float data1;#pragma DATA_SECTION(data2,"data_buf2") float data2;#pragma DATA_SECTION(data3,"data_buf3") float data3;#pragma DATA_SECTION(data4,"data_buf4") float data4;#pragma DATA_SECTION(data5,"data_buf5") float data5;void main(void){data1=1.1;data2=2.84;data3=data1 * data2;data4=data1/data2;data5=sqrt(1.1);}.cmd文件如下：MEMORY{IPRAM : origin=0x0,len=0x10000IDRAM : origin=0x10000000,len=0xF0000 }SECTIONS{.text >IPRAMdata_buf1>IDRAMdata_buf2>IDRAMdata_buf3>IDRAMdata_buf4>IDRAMdata_buf5>IDRAM.stack >IDRAMcinit >IDRAM.far >IDRAM}Memory window 设置如下图实验三6701核心系统的FFT运算程序FFT主程序#define PI 3.14159265358979#define M 10#define N 1024#include<math.h>void radix2(int * xr,int * xi,short * wr,short * wi);void ChangeOrder(int * xr,int * xi);int xr[N],xi[N];int x[2* N];int xm[N];int xlr[N],xli[N];short w1[N],w2[N];int y[2* N];int ym[N];int radix;int nx1;int i;main(){double delta;radix=2;nx1=N;delta=2* PI/nx1;for(i=0;i<nx1/2;i++){w1[i]=32767* (cos(i* delta));w2[i]=32767* (sin(i* delta));}for(i=0;i<nx1;i++){xr[i]=xlr[i]=(int)((cos(PI* i/5.0))* 0x80);xi[i]=xli[i]=1;x[2* i]=xlr[i];x[2* i+1]=xli[i];}for(i=0;i<nx1;i++){xm[i]=(int)sqrt((xr[i]* xr[i])+(xi[i]* xi[i]));}ChangeOrder(xlr,xli);radix2(xlr,xli,w1,w2);for(i=0;i<nx1;i++){int sum=0;sum=(int)sqrt((xlr[i]* xlr[i])+(xli[i]* xli[i]));ym[i]=sum;y[2* i]=(int)sqrt((xlr[i])* (xlr[i]));y[2* i+1]=(int)sqrt((xli[i])* (xli[i]));}}/*基2FFT算法程序*/void radix2(int xr[],int xi[],short wr[],short wi[]){short L,B,J,P,k;int rPartKB,iPartKB;for(L=1;L<=M;L++){B=(int)(pow(2,(L-1))+0.5);for(J=0;J<=B-1;J++){P=J* ((int)(pow(2,(M-L))+0.5));for(k=J;k<=N-1;k+=(int)(pow(2,L)+0.5)){rPartKB=xr[k+B]* wr[P]+xi[k+B]* wi[P];iPartKB=xi[k+B]* wr[P]-xr[k+B]* wi[P];rPartKB>>=15;iPartKB>>=15;xr[k+B]=xr[k]-rPartKB;xi[k+B]=xi[k]-iPartKB;xr[k]=xr[k]+rPartKB;xi[k]=xi[k]+iPartKB;}}}}/*倒序子程序*/void ChangeOrder(int xr[],int xi[]){int LH,N1,I,J,K;short T;LH=N/2;J=LH;N1=N-2;for(I=1;I<=N1;I++){if(I<J){T=xr[I];xr[I]=xr[J];xr[J]=T;T=xi[I];xi[I]=xi[J];xi[J]=T;}K=LH;while(J>=K){J=J-K;K=(int)(K/2+0.5);}J=J+K;}}fft.cmd:MEMORY{VECTORS:origin=0x0,len=0x400IPRAM:origin=0x400,len=0xf000IDRAM:origin=0x80000000,len=0x10000 }SECTIONS{.vec>VECTORS.text>IPRAM.data>IDRAM.bss>IDRAM.cinit>IDRAM.comst>IDRAM.far>IDRAMxlr变量,发现无法找到x1r变量。