femwind/hexa_coupling_test.m

   1 n=10;
   2 la=linspace(0,3,n);
   3 ld=linspace(0,3,n);
   4 a=10.^la;
   5 d=10.^ld;
   6 for i=1:n
   7     for j=1:n
   8         C=hexa_coupling([1,1,a(i)],[1,1,d(j)]);
   9         cv(i,j)=C(1,5);
  10         ch(i,j)=C(1,3);
  11         laa(i,j)=la(i);
  12         ldd(i,j)=ld(j);
  13         aa(i,j)=a(i);
  14         dd(i,j)=d(j);
  15     end
  16 end
  17 figure(1)
  18 mesh(laa,ldd,cv)
  19 hold on
  20 contour(laa,ldd,cv,[-1:0.1:1])
  21 hold off
  22 grid on
  23 title('Vertical coupling')
  24 xlabel('log vertical coefficient')
  25 ylabel('log vertical size')
  26
  27 figure(2)
  28 mesh(laa,ldd,ch)
  29 hold on
  30 contour(laa,ldd,ch,[-1:0.1:1])
  31 hold off
  32 grid on
  33 title('Horizontal coupling')
  34 xlabel('log vertical coefficient')
  35 ylabel('log vertical size')