我想利用矩阵来给定义多项式的系数。因为参数太多了,有30个,所以利用矩阵(别人教我的)具体的程序如下:function a=part1(p)
a1=p(1);a2=p(2);b1=p(3);b2=p(4);c1=p(5);c2=p(6);d1=p(7);d2=p(8);f1=p(9);f2=(10);j1=p(11);j2=p(12);a3=p(13);a4=p(14);m1=p(15);m2=p(16);n1=p(17);n2=p(18);g1=p(19);g2=p(20);h1=p(21);h2=p(22);b3=p(23);c3=p(24);q1=p(25);q2=p(26);k1=p(27);k2=p(28);e1=p(29);e2=p(30);
A1=h1*(-b1*(e2*g1-e1*g2)+c1*(e2*m1+g2*n1)+c2(n1*g1+e1*m1));
B1=h1*(b2*(e2*g1-e1*g2)+c1*(e2*m2+g2*n2)+c2*(n2*g1+e1*m2));
C1=(a1*h1-b1*q1)*(e2*g1-e1*g2)-((1-k1)*a3*h1-m1*q1)*(c1*e2+c2*e1)-((1-k2)*a4*h1-n1*q1)*(c1*g2+c2*g1)+j1*k1*b3*(e2*m1+g2*n1);D1=b2*q1*(e2*g1-e1*g2)+(c1*q1+j1*k1*b3)*(e2*m2+g2*n2)+c2*q1*(n2*g1+e1*m2);
F1=a1*q1*(e2*g1-e1*g2)-(c1*q1+j1*k1*b3)*((1-k1)*a3*e2)+(1-k2)*a4*g2-c2*q1*((1-k2)*a4*g1+(1-k1)*a3*e1)
;A2=h2*(-d2*(e2*g1-e1*g2)+f1*(e2*m2+g2*n2)+f2*(n2*g1+e1*m2));
B2=h2*(d1*(e2*g1-e1*g2)+f1*(e2*m1+g2*n1)+f2*(e1*m1+g1*n1));
C2=d1*q2*(e2*g1-e1*g2)+f1*q2*(e2*m1+g2*n1)+(f2*q2+j2*k2*c3)*(e1*m1+g1*n1);D2=(a2*h2-d2*q2)*(e2*g1-e1*g2)-((1-k1)*a3*h2-m2*q2)*(f1*e2+f2*e1)-((1-k2)*a4*h2-n2*q2)*(f1*g2+f2*g1)+j2*k2*c3*(n2*g1+e1*m2);
F2=a2*q2*(e2*g1-e1*g2)-f1*q2*((1-k1)*a3*e2+(1-k2)*a4*g2)-(f2*q2+j2*k2*c3)*((1-k2)*a4*g1+(1-k1)*a3*e1);
a=[A1,B1,C1,D1,F1,A2,B2,C2,D2,F2];
我的想法是:命令窗口输入一个矩阵,然后调用函数part1,就可以得到A1,B1等等的值,如果程序有问题请具体指出,我是一个初学者,太深奥的话我就不知道怎么改了,先谢谢大家的帮助,如果时间允许的话,最好给出可以复制过去就可以运行的程序!!



雷达卡



京公网安备 11010802022788号







