function?mul_result=mul（ab）????????%8bit‘s??multiply????mul[]?is?the?result
%?syms?a0?a1?a2?a3?a4?a5?a6?a7;
%?a=[a0?a1?a2?a3?a4?a5?a6?a7];
%b=[1?1?1?0?1?0?0?0];%所乘的系數要倒置，因為matlab是默認高位在左邊，比如23的二進制數是0?0?0??1?0?1?1?1，矩陣相乘時要寫為b=[1?1?1?0?1?0?0?0]；
mul（1）=b（1）*a（1）+b（2）*a（8）+b（3）*a（7）+b（4）*a（6）+b（5）*a（5）+b（6）*a（4）+b（6）*a（8）+b（7）*a（3）+b（7）*a（7）+b（7）*a（8）+b（8）*a（2）+b（8）*a（6）+b（8）*a（7）+b（8）*a（8）;
mul（2）=b（1）*a（2）+b（2）*a（1）+b（3）*a（8）+b（4）*a（7）+b（5）*a（6）+b（6）*a（5）+b（7）*a（4）+b（7）*a（8）+b（8）*a（3）+b（8）*a（7）+b（8）*a（8）;
mul（3）=b（1）*a（3）+b（2）*a（2）+b（2）*a（8）+b（3）*a（1）+b（3）*a（7）+b（4）*a（6）+b（4）*a（8）+b（5）*a（5）+b（5）*a（7）+b（6）*a（4）+b（6）*a（6）+b（6）*a（8）+b（7）*a（3）+b（7）*a（5）+b（7）*a（7）+b（7）*a（8）+b（8）*a（2）+b（8）*a（4）+b（8）*a（6）+b（8）*a（7）;
mul（4）=b（1）*a（4）+b（2）*a（3）+b（2）*a（8）+b（3）*a（2）+b（3）*a（7）+b（3）*a（8）+b（4）*a（1）+b（4）*a（6）+b（4）*a（7）+b（5）*a（5）+b（5）*a（6）+b（5）*a（8）+b（6）*a（4）+b（6）*a（5）+b（6）*a（7）+b（6）*a（8）+b（7）*a（3）+b（7）*a（4）+b（7）*a（6）