%-- 13.10.2010 14:01 --%
%-- 6.10.2010 14:51 --%
h = [  -4, 1, 1, 1, 0, 0, 1, 0, 0;
1, -4, 1, 0, 1, 0, 0, 1, 0;
1, 1, -4, 0, 0, 1, 0, 0, 1;
1, 0, 0, -4, 1, 1, 1, 0, 0;
0, 1, 0, 1, -4, 1, 0, 1, 0;
0, 0, 1, 1, 1, -4, 0, 0, 1;
1, 0, 0, 1, 0, 0, -4, 1, 1;
0, 1, 0, 0, 1, 0, 1, -4, 1;
0, 0, 1, 0, 0, 1, 1, 1, -4;]
f =[1;10;10;1;10;10;1;10;10]
f =[1;1;1;10;10;10;10;10;10]
g = h*f
mat = vekmat(3,3,g)
mat = veknamat(3,3,g)
hc = [2,0,0;0,2,0;0,0,2]
hr = hc
Hk = kron(hc,hr)
f = [1,2,3;5,7,11;13,17,19]
g = hc*f*hr'
f = [
0,0,0,0;
0,50,50,50,0;
0,0,0,0]
f = [
0,0,0,0;
0,50,50,0;
0,0,0,0]
hc = hr = [
1,1,1,1;
1,1,-1,-1;
1.414,-1.414,0,0;
0,0,1.414,-1.414]
hr = [
1,1,1,1;
1,1,-1,-1;
1.414,-1.414,0,0;
0,0,1.414,-1.414]
hc = hr
g = hc * f * hr'
H = kron(hr,hc);
vf = matnavek(4,4,f); % vektor F
vg = h*f; % vektor G
vg = H*f; % vektor G
vg = H*vf; % vektor G
vg = h*vf; % vektor G
vg = H*vf; % vektor G
g = veknamat(4,4,vg);
g
f
hr = [
1,1,1,1;
1,-1,1,-1;
1,1,-1,-1;
1,-1,-1,1]
hc = hr
g = hc * f * hr';
g
u1 = [1;1;1;1];
u2 = hc(:2); % druhý sloupec, všechny øádky
u3 = hc(:3);
u4 = hc(:4);
v1 = u1;
v2 = u2;
v3 = u3;
v4 = u4;
u1 = [1;1;1;1];
u2 = hc(:,2); % druhý sloupec, všechny øádky
u3 = hc(:,3);
u4 = hc(:,4);
v1 = u1;
v2 = u2;
v3 = u3;
v4 = u4;
u1 = [1;1;1;1];
u2 = hc(:,2); % druhý sloupec, všechny øádky
u3 = hc(:,3);
u4 = hc(:,4);
v1 = u1;
v2 = u2;
v3 = u3;
v4 = u4;
eo11 = 0*u1*v1';
eo12 = 0*u1*v2';
eo13 = 0*u1*v3';
eo14 = 0*u1*v4';
eo21 = 0*u2*v1';
eo22 = 50*u2*v2';
eo23 = 50*u2*v3';
eo24 = 0*u2*v4';
eo31 = 0*u2*v1';
eo32 = 50*u2*v2';
eo33 = 50*u2*v3';
eo34 = 0*u2*v4';
eo41 = 0*u1*v1';
eo42 = 0*u1*v2';
eo43 = 0*u1*v3';
eo44 = 0*u1*v4';
deo = oe11 +oe12 +oe13 +oe14 + oe21 +oe22 +oe23 +oe24 + oe31 +oe32 +oe33 +oe34+oe41 +oe42 +oe43 +oe44;
deo = eo11 +eo12 +eo13 +eo14 + eo21 +eo22 +eo23 +eo24 + eo31 +eo32 +eo33 +eo34+eo41 +eo42 +eo43 +eo44;
deo
