funmatinmmdpel.m

计算基本矩阵的matlab例程

		         
		   a=Fini(i1,j1);
		   b=Fini(i1,j2);
		   c=Fini(i2,j1);
		   d=Fini(i2,j2);
		
			switch im,
         case 1,
            F= [ (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2,                    -b*x2-d*y2;
					                     -a*x1-b*y1,                             a,                             b;
					                     -c*x1-d*y1,                             c,                             d];
         case 2,
            F= [ (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2,                    -b*x2-d*y2;
					                     -c*x1-d*y1,                             c,                             d;
					                     -a*x1-b*y1,                             a,                             b];
         case 3,
            F= [                    -a*x1-b*y1,                             a,                             b;
					  (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2,                    -b*x2-d*y2;
					                     -c*x1-d*y1,                             c,                             d];
         case 4,
            F= [                    -c*x1-d*y1,                             c,                             d;
					  (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2,                    -b*x2-d*y2;
					                     -a*x1-b*y1,                             a,                             b];
         case 5,
            F= [                    -a*x1-b*y1,                             a,                             b;
					                     -c*x1-d*y1,                             c,                             d;
					  (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2,                    -b*x2-d*y2];
         case 6,
            F= [                    -c*x1-d*y1,                             c,                             d;
					                     -a*x1-b*y1,                             a,                             b;
					  (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2,                    -b*x2-d*y2];
         case 7,
            F= [ (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2,                    -a*x2-c*y2;
					                     -a*x1-b*y1,                             b,                             a;
					                     -c*x1-d*y1,                             d,                             c];
         case 8,
            F= [ (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2,                    -a*x2-c*y2;
					                     -c*x1-d*y1,                             d,                             c;
					                     -a*x1-b*y1,                             b,                             a];
         case 9,
            F= [                    -a*x1-b*y1,                             b,                             a;
					  (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2,                    -a*x2-c*y2;
					                     -c*x1-d*y1,                             d,                             c];
         case 10,
            F= [                    -c*x1-d*y1,                             d,                             c;
					  (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2,                    -a*x2-c*y2;
					                     -a*x1-b*y1,                             b,                             a];
         case 11,
            F= [                    -a*x1-b*y1,                             b,                             a;
					                     -c*x1-d*y1,                             d,                             c;
					  (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2,                    -a*x2-c*y2];
         case 12,
            F= [                    -c*x1-d*y1,                             d,                             c;
					                     -a*x1-b*y1,                             b,                             a;
					  (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2,                    -a*x2-c*y2];
         case 13,
            F= [                    -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2;
					                              a,                    -a*x1-b*y1,                             b;
					                              c,                    -c*x1-d*y1,                             d];
         case 14,
            F= [                    -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2;
					                              c,                    -c*x1-d*y1,                             d;
					                              a,                    -a*x1-b*y1,                             b];
         case 15,
            F= [                             a,                    -a*x1-b*y1,                             b;
					                     -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2;
					                              c,                    -c*x1-d*y1,                             d];
         case 16,
            F= [                             c,                    -c*x1-d*y1,                             d;
					                     -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2;
					                              a,                    -a*x1-b*y1,                             b];
         case 17,
            F= [                             a,                    -a*x1-b*y1,                             b;
					                              c,                    -c*x1-d*y1,                             d;
					                     -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2];
         case 18,
            F= [                             c,                    -c*x1-d*y1,                             d;
					                              a,                    -a*x1-b*y1,                             b;
					                     -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -b*x2-d*y2];
         case 19,
            F= [                    -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2;
					                              b,                    -a*x1-b*y1,                             a;
					                              d,                    -c*x1-d*y1,                             c];
         case 20,
            F= [                    -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2;
					                              d,                    -c*x1-d*y1,                             c;
					                              b,                    -a*x1-b*y1,                             a];
         case 21,
            F= [                             b,                    -a*x1-b*y1,                             a;
					                     -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2;
					                              d,                    -c*x1-d*y1,                             c];
         case 22,
            F= [                             d,                    -c*x1-d*y1,                             c;
					                     -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2;
					                              b,                    -a*x1-b*y1,                             a];
         case 23,
            F= [                             b,                    -a*x1-b*y1,                             a;
					                              d,                    -c*x1-d*y1,                             c;
					                     -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2];
         case 24,
            F= [                             d,                    -c*x1-d*y1,                             c;
					                              b,                    -a*x1-b*y1,                             a;
					                     -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2,                    -a*x2-c*y2];
         case 25,
            F= [                    -a*x2-c*y2,                    -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2;
					                              a,                             b,                    -a*x1-b*y1;
					                              c,                             d,                    -c*x1-d*y1];
         case 26,
            F= [                    -a*x2-c*y2,                    -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2;
					                              c,                             d,                    -c*x1-d*y1;
					                              a,                             b,                    -a*x1-b*y1];
         case 27,
            F= [                             a,                             b,                    -a*x1-b*y1;
					                     -a*x2-c*y2,                    -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2;
					                              c,                             d,                    -c*x1-d*y1];
         case 28,
            F= [                             c,                             d,                    -c*x1-d*y1;
					                     -a*x2-c*y2,                    -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2;
					                              a,                             b,                    -a*x1-b*y1];
         case 29,
            F= [                             a,                             b,                    -a*x1-b*y1;
					                              c,                             d,                    -c*x1-d*y1;
					                     -a*x2-c*y2,                    -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2];
         case 30,
            F= [                             c,                             d,                    -c*x1-d*y1;
					                              a,                             b,                    -a*x1-b*y1;
					                     -a*x2-c*y2,                    -b*x2-d*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2];
         case 31,
            F= [                    -b*x2-d*y2,                    -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2;
					                              b,                             a,                    -a*x1-b*y1;
					                              d,                             c,                    -c*x1-d*y1];
         case 32,
            F= [                    -b*x2-d*y2,                    -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2;
					                              d,                             c,                    -c*x1-d*y1;
					                              b,                             a,                    -a*x1-b*y1];
         case 33,
            F= [                             b,                             a,                    -a*x1-b*y1;
					                     -b*x2-d*y2,                    -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2;
					                              d,                             c,                    -c*x1-d*y1];
         case 34,
            F= [                             d,                             c,                    -c*x1-d*y1;
					                     -b*x2-d*y2,                    -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2;
					                              b,                             a,                    -a*x1-b*y1];
         case 35,
            F= [                             b,                             a,                    -a*x1-b*y1;
					                              d,                             c,                    -c*x1-d*y1;
					                     -b*x2-d*y2,                    -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2];
         case 36,
            F= [                             d,                             c,                    -c*x1-d*y1;
					                              b,                             a,                    -a*x1-b*y1;
					                     -b*x2-d*y2,                    -a*x2-c*y2, (a*x1+b*y1)*x2+(c*x1+d*y1)*y2];
         case 37,
				disp('Error')
         end         
         Fini=F;    
         vmini=vm;
	   else
         F=Fini;   
         break;
	   end	      
	end;
end;