信号与信息处理综合实验实验报告二矩阵乘法.docVIP

信号与信息处理综合实验实验报告 实验二 矩阵乘法  一、实验要求： 完成矩阵A*B+C*D 其中： A, B, C和D矩阵都是六维矩阵 要求有汇编语言和C语言两种算法编写 要求有子程序的调用? 二、源代码 汇编部分： #includedefBF535.h .section data1; .BYTE2 aArray[36]=a.dat; .BYTE2 bArray[36]=b.dat; .BYTE2 cArray[36]=c.dat; .BYTE2 dArray[36]=d.dat; .BYTE2 eArray[36]; .section program ; .global _main ; .global _xiangcheng; _main : R0 = SYSCFG ;//系统配置寄存器 bitset (R0, 0x1) ;//位设置命令，把R0中0x1位置设置为1 SYSCFG = R0 ;//SYSCFG0x1位置为周期计数器使能，1使能64位自由运行周期计数器 M1=12;//利用M1来移动指针使其指向下一列数据 P1=6;//循环次数 I0.H = aArray ;//把a矩阵的地址给I0 I0.L = aArray ; I1.H = bArray ;//把b矩阵的地址给I1 I1.L = bArray ; I2.H = eArray ;//把e矩阵的地址给I2 I2.L = eArray ; CALL _xiangcheng; //调用子程序，完成aArray*bArray I0.H = cArray ;//把c矩阵的地址给I0 I0.L = cArray ; I1.H = dArray ;//把d矩阵的地址给I1 I1.L = dArray ; I2.H = eArray ;//把e矩阵的地址给I2 I2.L = eArray ; CALL _xiangcheng; //调用子程序，完成aArray*bArray +cArray*dArray idle;//空闲指令：通过设置SEQSTAT寄存器的idle_req位使Blackfin处理器进入空闲状态 _main.end: //矩阵相乘子函数 // _xiangcheng: //需要保存外层循环信息 loop lp1 LC0=P1;//定义循环计数器LC0为循环lp1计数，次数为6（for i=1:6） loop_begin lp1;//开始循环lp1 [--sp]=I1;/* save loop lp1 */ [--sp]=I0; [--sp]=LC0;//LC0寄存器保留剩余循环次数 [--sp]=LT0;//LT0寄存器保持循环中第一条语句的前地址 [--sp]=LB0;//LB0寄存器保持循环中最后一条语句的前地址 loop lp2 LC1=P1;//定义循环计数器LC0为循环lp2计数，次数为6（for j=1:6） loop_begin lp2;//开始循环lp2 [--sp]=I0;/* save loop lp2 */ [--sp]=I1; loop lp3 LC0=P1;//定义循环计数器LC0为循环lp3计数，次数为6（for k=1:6） loop_begin lp3; R1.L=W[I0];//把I0指针指向的数据a(i,k)给R1 I0+=2;//移动I0指针，使其指向矩阵该行的下一列元素a(i,k+1) R2.L=W[I1];//把I1指针指向的数据b(k,j)给R2 I1+=M1; //移动指针I1，使得指针后移六个数据，指向该列的下一行元素b(k+1,j) R0=R1.L*R2.L(IS);//做矩阵对应相乘的元素相乘运算a(i,k)*b(k,j）//IS：无符号整数 R3=R3+R0;//累加元素相乘结果，得到结果矩阵的元素c(i,j)=c(i,j)+a(i,k)*b(k,j) loop_end lp3; R4.L=W[I2];//取出结果矩阵的数据以便在第二次调用子函数的时候完成两次矩阵结果做和的运算 R3=R3+R4;//将本次运算结果累加结果矩阵中对应的元素 W[I2++]=R3.L;//保存运算结果 R3=0;//清0 R3，为下一次循环做准备 I1=[SP++];//输出