Fork me on GitHub

P4基函数

 1 July 29,2015
 2 P4
 3 syms p1 p2 p3
 4 phi(:,1) = 1/6*p1.*(4*p1-1).*(4*p1-2).*(4*p1-3);
 5 phi(:,2) = 1/6*p2.*(4*p2-1).*(4*p2-2).*(4*p2-3);
 6 phi(:,3) = 1/6*p3.*(4*p3-1).*(4*p3-2).*(4*p3-3);
 7 phi(:,4) = 32*p1.*p2.*p3.*(4*p1-1);
 8 phi(:,5) = 32*p1.*p2.*p3.*(4*p2-1);
 9 phi(:,6) = 32*p1.*p2.*p3.*(4*p3-1);
10 phi(:,7) = 8/3*p1.*p2.*(4*p1-1).*(4*p1-2);
11 phi(:,8) = 4*p1.*p2.*(4*p1-1).*(4*p2-1);
12 phi(:,9) = 8/3*p1.*p2.*(4*p2-1).*(4*p2-2);
13 phi(:,10) = 8/3*p1.*p3*(4*p1-1).*(4*p1-2);
14 phi(:,11) = 4*p1.*p3.*(4*p1-1).*(4*p3-1);
15 phi(:,12) = 8/3*p1.*p3.*(4*p3-1).*(4*p3-1);
16 phi(:,13) = 8/3*p2.*p3.*(4*p2-1).*(4*p2-2);
17 phi(:,14) = 4*p2.*p3.*(4*p2-1).*(4*p3-1);
18 phi(:,15) = 8/3*p2.*p3.*(4*p3-1).*(4*p3-2);
P4基函数

 1 July 29,2015
 2 P4
 3 syms p1 p2 p3
 4 phi(:,1) = 1/6*p1.*(4*p1-1).*(4*p1-2).*(4*p1-3);
 5 phi(:,2) = 1/6*p2.*(4*p2-1).*(4*p2-2).*(4*p2-3);
 6 phi(:,3) = 1/6*p3.*(4*p3-1).*(4*p3-2).*(4*p3-3);
 7 phi(:,4) = 8/3*p2.*p3.*(4*p2-1).*(4*p2-2);
 8 phi(:,5) = 8/3*p1.*p3.*(4*p3-1).*(4*p3-1);
 9 phi(:,6) = 8/3*p1.*p2.*(4*p1-1).*(4*p1-2);
10 phi(:,7) =  4*p2.*p3.*(4*p2-1).*(4*p3-1);
11 phi(:,8) =  4*p1.*p3.*(4*p1-1).*(4*p3-1);
12 phi(:,9) =  4*p1.*p2.*(4*p1-1).*(4*p2-1);
13 phi(:,10) = 8/3*p2.*p3.*(4*p3-1).*(4*p3-2);
14 phi(:,11) = 8/3*p1.*p3.*(4*p1-1).*(4*p1-2);
15 phi(:,12) = 8/3*p1.*p2.*(4*p2-1).*(4*p2-2);
16 phi(:,13) = 32*p1.*p2.*p3.*(4*p1-1);
17 phi(:,14) = 32*p1.*p2.*p3.*(4*p2-1);
18 phi(:,15) = 32*p1.*p2.*p3.*(4*p3-1);
P4元基函数

 

posted @ 2016-07-29 10:55  小奔奔  阅读(151)  评论(0编辑  收藏  举报