EEQ1rot元程序的一些重要组件.docVIP

syms x y h a1 a2 a3 a4 a5 b1 b2 b3 b4 b5 f1 -1/4-y/2+3*y^2/4; f2 -1/4+x/2+3*x^2/4; f3 -1/4+y/2+3*y^2/4; f4 -1/4-x/2+3*x^2/4; f5 2-3*x^2/2-3*y^2/2; M [2*diff f1,x /h,2*diff f1,y /h;2*diff f2,x /h,2*diff f2,y /h;2*diff f3,x /h,2*diff f3,y /h;2*diff f4,x /h,2*diff f4,y /h;2*diff f5,x /h,2*diff f5,y /h]; Ke zeros 5,5 ; Ke sym Ke ; f [f1,f2,f3,f4,f5]; for i 1:5 for j 1:5 bj dot M i,: ,M j,: +f i *f j * h^2 /4; jf int int bj,x,-1,1 ,y,-1,1 ; Ke i,j jf; end end Ke simple Ke Ke [101/240*h1*h2+3/4/h2*h1+3/4/h1*h2, 3/80*h1*h2-3/4/h2*h1-3/4/h1*h2, -59/240*h1*h2+3/4/h2*h1+3/4/h1*h2, 3/80*h1*h2-3/4/h2*h1-3/4/h1*h2] [ 3/80*h1*h2-3/4/h2*h1-3/4/h1*h2, 101/240*h1*h2+3/4/h2*h1+3/4/h1*h2, 3/80*h1*h2-3/4/h2*h1-3/4/h1*h2, -59/240*h1*h2+3/4/h2*h1+3/4/h1*h2] [-59/240*h1*h2+3/4/h2*h1+3/4/h1*h2, 3/80*h1*h2-3/4/h2*h1-3/4/h1*h2, 101/240*h1*h2+3/4/h2*h1+3/4/h1*h2, 3/80*h1*h2-3/4/h2*h1-3/4/h1*h2] [ 3/80*h1*h2-3/4/h2*h1-3/4/h1*h2, -59/240*h1*h2+3/4/h2*h1+3/4/h1*h2, 3/80*h1*h2-3/4/h2*h1-3/4/h1*h2, 101/240*h1*h2+3/4/h2*h1+3/4/h1*h2] function K Q1 n N n* 2* n+1 ; h1 1/n;h2 1/n; K sparse N,N ; lnm bianhao n ; for i 1:n^2 x lnm :,i ; a1 x 1 ;a2 x 2 ;a3 x 3 ;a4 x 4 ; K a1,a1 K a1,a1 +101/240*h1*h2+3/4/h2*h1+3/4/h1*h2; K a1,a2 K a1,a2 +3/80*h1*h2-3/4/h2*h1-3/4/h1*h2; K a1,a3 K a1,a3 -59/240*h1*h2+3/4/h2*h1+3/4/h1*h2; K a1,a4 K a1,a4 +3/80*h1*h2-3/4/h2*h1-3/4/h1*h2; K a2,a1 K a2,a1 +3/80*h1*h2-3/4/h2*h1-3/4/h1*h2; K a2,a2 K a2,a2 +101/240*h1*h2+3/4/h2*h1+3/4/h1*h2; K a2,a3 K a2,a3 +3/80*h1*h2-3/4/h2*h1-3/4/h1*h2; K a2,a4 K a2,a4 -59/240*h1*h2+3/4/h2*h1+3/4/h1*h2; K a3,a1 K a3,a1 -59/240*h1*h2+3/4/h2*h1+3/4/h1*h2; K a3,a2 K a3,a2 +3/80*h1*h2-3/4/h2*h1-3/4/h1*h2; K a3,a3 K a3,a3 +101/240*h1*h2+3/4/h2*h1+3/4/h1*h2; K a3,a4 K a3,a4 +3/80*h1*h2-3/4/h2*h1-3/4/h1*h2; K a4,a1 K a4,a1 +3/80*h1*h2-3/4/h2*h1-3/4/h1*h2; K a4,a2 K a4,a2 -59/240*h1*h2+3/4/h2*h1+3/4/h1*h2; K a4,a3 K a4,a3 +3/80*h1*h2-3/4/h2*h1-3/4/h1*h2; K a4,a4 K a4,a4 +101/240*h1*h2+3/4/h2*h1+3/4/h1*h2; end function ln