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

信号与信息处理综合实验实验报告 实验二矩阵乘法 一、实验要求: 完成矩阵A*B+C*D其中： A, B, C和D矩阵都是六维矩阵 要求有汇编语言和C语言两种算法编写 要求有子程序的调用 源代码 汇编部分： #include,,defBF535.h .section datal; .BYTE2 aArray[36]=a.dat; .BYTE2 bArray[36]=,,b.dat; .BYTE2 cArray[36]=nc.dat; .BYTE2 dArray[36]=d.datH; .BYTE2 eArray[36]; .section program; .global _main; .global _xiangcheng; main : RO = SYSCFG ;//系统配置寄存器 bitset (RO, 0x1);//位设置命令，把R0中0x1位置设置为1 SYSCFG = RO ;//SYSCFGOxl位置为周期计数器使能，1使能64位自由运行周期计数器 Ml=12;〃利用Ml来移动指针使其指向下一列数据 Pl=6;//循环次数 I0.H = aArray ;〃把a矩阵的地址给10 10 丄=aArray; H = bArray;〃把b矩阵的地址给11 11 丄二 bArray; H = eArray;〃把e矩阵的地址给12 I2.L = eArray; CALL _xiangcheng; // 调用子程序，完成 a Array* bArray H = cArray;〃把c矩阵的地址给10 10 丄=cArray; H = dArray ;//把d矩阵的地址给11 L = dArray; H = eArray;〃把e矩阵的地址给12 12丄=eArray ; CALL _xiangcheng; // 调用子程序，完成 a Array* bArray +cArray*dArray idle;//空闲指令：通过设置SEQSTAT寄存器的idle_req位使Blackfin处理器进入空闲状态 _main.e nd: 〃矩阵相乘子函数〃 _xia ngche ng: 〃需要保存外层循环信息 loop Ipl LCO=P1;//定义循环计数器LCO为循坏Ipl计数，次数为6 (for i=l:6) loop_begin Ipl;//开始循环 Ipl [-sp]=ll;/* save loop Ipl */ [--sp]=IO; [-sp]=LCO;//LCO寄存器保留剩余循环次数 [-sp]=LTO;//LTO寄存器保持循环中第一条语句的前地址 [-sp]=LBO;//LBO寄存器保持循环中最后一条语句的前地址 loop Ip2 LC1=P1;//定义循环计数器LCO为循环Ip2计数，次数为6 (forj=l:6) loop_begin Ip2;〃开始循环 Ip2 [-sp]=IO;/* save loop Ip2 */ [~sp]=ll； loop Ip3 LCO二Pl;〃定义循环计数器LCO为循环Ip3计数,次数为6(for k=l:6) loop_begin Ip3; R1.L=W[IO];//把10指针指向的数据a(i,k)给R1 10+=2;//移动10指针，使其指向矩阵该行的下一列元素a(i,k+l) R2.L=W[I1];//把11指针指向的数据b(k,j)给R2 I1+=M1; 〃移动指针11,使得指针后移六个数据，指向该列的下一行元素 b(k+lj) RO=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 Ip3; R4.L=W[I2];//取出结果矩阵的数据以便在第二次调用子函数的时候完成两次矩阵 结果做和的运算 R3二R3+R4;//将本次运算结果累加结果矩阵中对应的元素 W[I2++]=R3.L;//保存运算结果 R3=0;〃清0R3,为下一次循环做准备 H=[SP++]；//输出 Ip2 参数 10 二[SP++]; 11+=2;//移动指针II,此次运算设置为乘以第二个矩阵的下一列 loop_end Ip2; LB0=[SP++];//输出 Ip2 参数 LT0=[SP++]; LC0=[SP++]; IO=[SP++]; ll=[SP++]; 〃为此次乘法循环设置为第一个矩阵的下一行元素 10+二Ml;〃移动指针10,使得指针后移六位数据，指向该列的下一行元素 loop_end Ipl; RTS; _xiangche ng.end: C语言部分： #include stdio.h #include std