source: proiecte/pmake3d/make3d_original/Make3dSingleImageStanford_version0.1/third_party/opt/yalmip/dev/tests-mbg/sdpvar/sdpvar_hankel.m @ 37

Last change on this file since 37 was 37, checked in by (none), 14 years ago

Added original make3d
File size: 1.9 KB
Line 
1function sdpvar_hankel
2
3c = sdpvar(4,1);
4assign(c,[2 4 6 8]');
5ok = all(all(double(hankel(c,c))-hankel(double(c),double(c)) == 0));
6mbg_asserttrue(ok)
7
8ok = all(all(double(hankel(c))-hankel(double(c)) == 0));
9mbg_asserttrue(ok)
10
11r = sdpvar(4,1);
12assign(r,-[2 4 6 8]');
13ok = all(all(double(hankel(c,r))-hankel(double(c),double(r)) == 0));
14mbg_asserttrue(ok)
15
16c = sdpvar(4,1);
17assign(c,[2 4 6 8]');
18r = sdpvar(2,1);
19assign(r,-[2 4]');
20ok = all(all(double(hankel(c,r))-hankel(double(c),double(r)) == 0));
21mbg_asserttrue(ok)
22
23ok = all(all(double(hankel(r,c'))-hankel(double(r),double(c')) == 0));
24mbg_asserttrue(ok)
25
26c = sdpvar(2,3);
27r = sdpvar(4,1);
28assign(c,[1 2 3;4 5 6]);
29assign(r,[1 2 3 4]');
30ok = all(all(double(hankel(c,r))-hankel(double(c),double(r))==0));
31mbg_asserttrue(ok)
32
33c = sdpvar(2,3);
34r = sdpvar(4,1);
35assign(c,[1 2 3;4 5 6]);
36assign(r,[1 2 3 4]');
37ok = all(all(double(hankel(r,c))-hankel(double(r),double(c))==0));
38mbg_asserttrue(ok)
39
40c = sdpvar(4,1,'full','complex');
41assign(c,[2 4 6 8]'+sqrt(-1)*[5 4 3 2]');
42ok = all(all(double(hankel(c,c))-hankel(double(c),double(c)) == 0));
43mbg_asserttrue(ok)
44
45ok = all(all(double(hankel(c))-hankel(double(c)) == 0));
46mbg_asserttrue(ok)
47
48c = sdpvar(4,1,'full','complex');
49assign(c,[2 4 6 8]'+sqrt(-1)*[5 4 3 2]');
50r = sdpvar(4,1);
51assign(r,-[2 4 6 8]');
52ok = all(all(double(hankel(c,r))-hankel(double(c),double(r)) == 0));
53mbg_asserttrue(ok)
54
55r = sdpvar(2,1);
56assign(r,-[2 4]'+sqrt(-1));
57ok = all(all(double(hankel(c,r))-hankel(double(c),double(r)) == 0));
58mbg_asserttrue(ok)
59
60ok = all(all(double(hankel(r,c'))-hankel(double(r),double(c')) == 0));
61mbg_asserttrue(ok)
62
63c = sdpvar(2,3);
64r = sdpvar(4,1);
65assign(c,[1 2 3;4 5 6]);
66assign(r,[1 2 3 4]'*sqrt(-1));
67ok = all(all(double(hankel(c,r))-hankel(double(c),double(r))==0));
68mbg_asserttrue(ok)
69
70c = sdpvar(2,3);
71r = sdpvar(4,1);
72assign(c,[1 2 3;4 5 6]*sqrt(-1));
73assign(r,[1 2 3 4]');
74ok = all(all(double(hankel(r,c))-hankel(double(r),double(c))==0));
75mbg_asserttrue(ok)
76
77
78
79
80
81
82
83
84
85
86
87
Note: See TracBrowser for help on using the repository browser.