author  hoelzl 
Tue, 06 Nov 2012 19:18:35 +0100  
changeset 50021  d96a3f468203 
parent 50003  8c213922ed49 
child 50038  8e32c9254535 
permissions  rwrr 
42146
5b52c6a9c627
split Product_Measure into Binary_Product_Measure and Finite_Product_Measure
hoelzl
parents:
42067
diff
changeset

1 
(* Title: HOL/Probability/Finite_Product_Measure.thy 
42067  2 
Author: Johannes Hölzl, TU München 
3 
*) 

4 

42146
5b52c6a9c627
split Product_Measure into Binary_Product_Measure and Finite_Product_Measure
hoelzl
parents:
42067
diff
changeset

5 
header {*Finite product measures*} 
42067  6 

42146
5b52c6a9c627
split Product_Measure into Binary_Product_Measure and Finite_Product_Measure
hoelzl
parents:
42067
diff
changeset

7 
theory Finite_Product_Measure 
5b52c6a9c627
split Product_Measure into Binary_Product_Measure and Finite_Product_Measure
hoelzl
parents:
42067
diff
changeset

8 
imports Binary_Product_Measure 
35833  9 
begin 
10 

47694  11 
lemma split_const: "(\<lambda>(i, j). c) = (\<lambda>_. c)" 
12 
by auto 

13 

40859  14 
abbreviation 
15 
"Pi\<^isub>E A B \<equiv> Pi A B \<inter> extensional A" 

39094  16 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

17 
syntax 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

18 
"_PiE" :: "[pttrn, 'a set, 'b set] => ('a => 'b) set" ("(3PIE _:_./ _)" 10) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

19 

3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

20 
syntax (xsymbols) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

21 
"_PiE" :: "[pttrn, 'a set, 'b set] => ('a => 'b) set" ("(3\<Pi>\<^isub>E _\<in>_./ _)" 10) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

22 

3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

23 
syntax (HTML output) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

24 
"_PiE" :: "[pttrn, 'a set, 'b set] => ('a => 'b) set" ("(3\<Pi>\<^isub>E _\<in>_./ _)" 10) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

25 

3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

26 
translations 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

27 
"PIE x:A. B" == "CONST Pi\<^isub>E A (%x. B)" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

28 

40859  29 
abbreviation 
30 
funcset_extensional :: "['a set, 'b set] => ('a => 'b) set" 

31 
(infixr ">\<^isub>E" 60) where 

32 
"A >\<^isub>E B \<equiv> Pi\<^isub>E A (%_. B)" 

33 

34 
notation (xsymbols) 

35 
funcset_extensional (infixr "\<rightarrow>\<^isub>E" 60) 

36 

37 
lemma extensional_insert[intro, simp]: 

38 
assumes "a \<in> extensional (insert i I)" 

39 
shows "a(i := b) \<in> extensional (insert i I)" 

40 
using assms unfolding extensional_def by auto 

41 

42 
lemma extensional_Int[simp]: 

43 
"extensional I \<inter> extensional I' = extensional (I \<inter> I')" 

44 
unfolding extensional_def by auto 

38656  45 

35833  46 
definition 
49780  47 
"merge I J = (\<lambda>(x, y) i. if i \<in> I then x i else if i \<in> J then y i else undefined)" 
40859  48 

49 
lemma merge_apply[simp]: 

49780  50 
"I \<inter> J = {} \<Longrightarrow> i \<in> I \<Longrightarrow> merge I J (x, y) i = x i" 
51 
"I \<inter> J = {} \<Longrightarrow> i \<in> J \<Longrightarrow> merge I J (x, y) i = y i" 

52 
"J \<inter> I = {} \<Longrightarrow> i \<in> I \<Longrightarrow> merge I J (x, y) i = x i" 

53 
"J \<inter> I = {} \<Longrightarrow> i \<in> J \<Longrightarrow> merge I J (x, y) i = y i" 

54 
"i \<notin> I \<Longrightarrow> i \<notin> J \<Longrightarrow> merge I J (x, y) i = undefined" 

40859  55 
unfolding merge_def by auto 
56 

57 
lemma merge_commute: 

49780  58 
"I \<inter> J = {} \<Longrightarrow> merge I J (x, y) = merge J I (y, x)" 
50003  59 
by (force simp: merge_def) 
40859  60 

61 
lemma Pi_cancel_merge_range[simp]: 

49780  62 
"I \<inter> J = {} \<Longrightarrow> x \<in> Pi I (merge I J (A, B)) \<longleftrightarrow> x \<in> Pi I A" 
63 
"I \<inter> J = {} \<Longrightarrow> x \<in> Pi I (merge J I (B, A)) \<longleftrightarrow> x \<in> Pi I A" 

64 
"J \<inter> I = {} \<Longrightarrow> x \<in> Pi I (merge I J (A, B)) \<longleftrightarrow> x \<in> Pi I A" 

65 
"J \<inter> I = {} \<Longrightarrow> x \<in> Pi I (merge J I (B, A)) \<longleftrightarrow> x \<in> Pi I A" 

40859  66 
by (auto simp: Pi_def) 
67 

68 
lemma Pi_cancel_merge[simp]: 

49780  69 
"I \<inter> J = {} \<Longrightarrow> merge I J (x, y) \<in> Pi I B \<longleftrightarrow> x \<in> Pi I B" 
70 
"J \<inter> I = {} \<Longrightarrow> merge I J (x, y) \<in> Pi I B \<longleftrightarrow> x \<in> Pi I B" 

71 
"I \<inter> J = {} \<Longrightarrow> merge I J (x, y) \<in> Pi J B \<longleftrightarrow> y \<in> Pi J B" 

72 
"J \<inter> I = {} \<Longrightarrow> merge I J (x, y) \<in> Pi J B \<longleftrightarrow> y \<in> Pi J B" 

40859  73 
by (auto simp: Pi_def) 
74 

49780  75 
lemma extensional_merge[simp]: "merge I J (x, y) \<in> extensional (I \<union> J)" 
40859  76 
by (auto simp: extensional_def) 
77 

78 
lemma restrict_Pi_cancel: "restrict x I \<in> Pi I A \<longleftrightarrow> x \<in> Pi I A" 

79 
by (auto simp: restrict_def Pi_def) 

80 

81 
lemma restrict_merge[simp]: 

49780  82 
"I \<inter> J = {} \<Longrightarrow> restrict (merge I J (x, y)) I = restrict x I" 
83 
"I \<inter> J = {} \<Longrightarrow> restrict (merge I J (x, y)) J = restrict y J" 

84 
"J \<inter> I = {} \<Longrightarrow> restrict (merge I J (x, y)) I = restrict x I" 

85 
"J \<inter> I = {} \<Longrightarrow> restrict (merge I J (x, y)) J = restrict y J" 

47694  86 
by (auto simp: restrict_def) 
40859  87 

88 
lemma extensional_insert_undefined[intro, simp]: 

89 
assumes "a \<in> extensional (insert i I)" 

90 
shows "a(i := undefined) \<in> extensional I" 

91 
using assms unfolding extensional_def by auto 

92 

93 
lemma extensional_insert_cancel[intro, simp]: 

94 
assumes "a \<in> extensional I" 

95 
shows "a \<in> extensional (insert i I)" 

96 
using assms unfolding extensional_def by auto 

97 

49780  98 
lemma merge_singleton[simp]: "i \<notin> I \<Longrightarrow> merge I {i} (x,y) = restrict (x(i := y i)) (insert i I)" 
41095  99 
unfolding merge_def by (auto simp: fun_eq_iff) 
100 

101 
lemma Pi_Int: "Pi I E \<inter> Pi I F = (\<Pi> i\<in>I. E i \<inter> F i)" 

102 
by auto 

103 

40859  104 
lemma PiE_Int: "(Pi\<^isub>E I A) \<inter> (Pi\<^isub>E I B) = Pi\<^isub>E I (\<lambda>x. A x \<inter> B x)" 
105 
by auto 

106 

107 
lemma Pi_cancel_fupd_range[simp]: "i \<notin> I \<Longrightarrow> x \<in> Pi I (B(i := b)) \<longleftrightarrow> x \<in> Pi I B" 

108 
by (auto simp: Pi_def) 

109 

110 
lemma Pi_split_insert_domain[simp]: "x \<in> Pi (insert i I) X \<longleftrightarrow> x \<in> Pi I X \<and> x i \<in> X i" 

111 
by (auto simp: Pi_def) 

39088
ca17017c10e6
Measurable on product space is equiv. to measurable components
hoelzl
parents:
39082
diff
changeset

112 

40859  113 
lemma Pi_split_domain[simp]: "x \<in> Pi (I \<union> J) X \<longleftrightarrow> x \<in> Pi I X \<and> x \<in> Pi J X" 
114 
by (auto simp: Pi_def) 

115 

116 
lemma Pi_cancel_fupd[simp]: "i \<notin> I \<Longrightarrow> x(i := a) \<in> Pi I B \<longleftrightarrow> x \<in> Pi I B" 

117 
by (auto simp: Pi_def) 

118 

41095  119 
lemma restrict_vimage: 
120 
assumes "I \<inter> J = {}" 

49780  121 
shows "(\<lambda>x. (restrict x I, restrict x J)) ` (Pi\<^isub>E I E \<times> Pi\<^isub>E J F) = Pi (I \<union> J) (merge I J (E, F))" 
41095  122 
using assms by (auto simp: restrict_Pi_cancel) 
123 

124 
lemma merge_vimage: 

125 
assumes "I \<inter> J = {}" 

49780  126 
shows "merge I J ` Pi\<^isub>E (I \<union> J) E = Pi I E \<times> Pi J E" 
41095  127 
using assms by (auto simp: restrict_Pi_cancel) 
128 

129 
lemma restrict_fupd[simp]: "i \<notin> I \<Longrightarrow> restrict (f (i := x)) I = restrict f I" 

47694  130 
by (auto simp: restrict_def) 
41095  131 

132 
lemma merge_restrict[simp]: 

49780  133 
"merge I J (restrict x I, y) = merge I J (x, y)" 
134 
"merge I J (x, restrict y J) = merge I J (x, y)" 

47694  135 
unfolding merge_def by auto 
41095  136 

137 
lemma merge_x_x_eq_restrict[simp]: 

49780  138 
"merge I J (x, x) = restrict x (I \<union> J)" 
47694  139 
unfolding merge_def by auto 
41095  140 

141 
lemma Pi_fupd_iff: "i \<in> I \<Longrightarrow> f \<in> Pi I (B(i := A)) \<longleftrightarrow> f \<in> Pi (I  {i}) B \<and> f i \<in> A" 

142 
apply auto 

143 
apply (drule_tac x=x in Pi_mem) 

144 
apply (simp_all split: split_if_asm) 

145 
apply (drule_tac x=i in Pi_mem) 

146 
apply (auto dest!: Pi_mem) 

147 
done 

148 

149 
lemma Pi_UN: 

150 
fixes A :: "nat \<Rightarrow> 'i \<Rightarrow> 'a set" 

151 
assumes "finite I" and mono: "\<And>i n m. i \<in> I \<Longrightarrow> n \<le> m \<Longrightarrow> A n i \<subseteq> A m i" 

152 
shows "(\<Union>n. Pi I (A n)) = (\<Pi> i\<in>I. \<Union>n. A n i)" 

153 
proof (intro set_eqI iffI) 

154 
fix f assume "f \<in> (\<Pi> i\<in>I. \<Union>n. A n i)" 

155 
then have "\<forall>i\<in>I. \<exists>n. f i \<in> A n i" by auto 

156 
from bchoice[OF this] obtain n where n: "\<And>i. i \<in> I \<Longrightarrow> f i \<in> (A (n i) i)" by auto 

157 
obtain k where k: "\<And>i. i \<in> I \<Longrightarrow> n i \<le> k" 

158 
using `finite I` finite_nat_set_iff_bounded_le[of "n`I"] by auto 

159 
have "f \<in> Pi I (A k)" 

160 
proof (intro Pi_I) 

161 
fix i assume "i \<in> I" 

162 
from mono[OF this, of "n i" k] k[OF this] n[OF this] 

163 
show "f i \<in> A k i" by auto 

164 
qed 

165 
then show "f \<in> (\<Union>n. Pi I (A n))" by auto 

166 
qed auto 

167 

168 
lemma PiE_cong: 

169 
assumes "\<And>i. i\<in>I \<Longrightarrow> A i = B i" 

170 
shows "Pi\<^isub>E I A = Pi\<^isub>E I B" 

171 
using assms by (auto intro!: Pi_cong) 

172 

173 
lemma restrict_upd[simp]: 

174 
"i \<notin> I \<Longrightarrow> (restrict f I)(i := y) = restrict (f(i := y)) (insert i I)" 

175 
by (auto simp: fun_eq_iff) 

176 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

177 
lemma Pi_eq_subset: 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

178 
assumes ne: "\<And>i. i \<in> I \<Longrightarrow> F i \<noteq> {}" "\<And>i. i \<in> I \<Longrightarrow> F' i \<noteq> {}" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

179 
assumes eq: "Pi\<^isub>E I F = Pi\<^isub>E I F'" and "i \<in> I" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

180 
shows "F i \<subseteq> F' i" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

181 
proof 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

182 
fix x assume "x \<in> F i" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

183 
with ne have "\<forall>j. \<exists>y. ((j \<in> I \<longrightarrow> y \<in> F j \<and> (i = j \<longrightarrow> x = y)) \<and> (j \<notin> I \<longrightarrow> y = undefined))" by auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

184 
from choice[OF this] guess f .. note f = this 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

185 
then have "f \<in> Pi\<^isub>E I F" by (auto simp: extensional_def) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

186 
then have "f \<in> Pi\<^isub>E I F'" using assms by simp 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

187 
then show "x \<in> F' i" using f `i \<in> I` by auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

188 
qed 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

189 

3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

190 
lemma Pi_eq_iff_not_empty: 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

191 
assumes ne: "\<And>i. i \<in> I \<Longrightarrow> F i \<noteq> {}" "\<And>i. i \<in> I \<Longrightarrow> F' i \<noteq> {}" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

192 
shows "Pi\<^isub>E I F = Pi\<^isub>E I F' \<longleftrightarrow> (\<forall>i\<in>I. F i = F' i)" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

193 
proof (intro iffI ballI) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

194 
fix i assume eq: "Pi\<^isub>E I F = Pi\<^isub>E I F'" and i: "i \<in> I" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

195 
show "F i = F' i" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

196 
using Pi_eq_subset[of I F F', OF ne eq i] 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

197 
using Pi_eq_subset[of I F' F, OF ne(2,1) eq[symmetric] i] 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

198 
by auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

199 
qed auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

200 

3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

201 
lemma Pi_eq_empty_iff: 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

202 
"Pi\<^isub>E I F = {} \<longleftrightarrow> (\<exists>i\<in>I. F i = {})" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

203 
proof 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

204 
assume "Pi\<^isub>E I F = {}" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

205 
show "\<exists>i\<in>I. F i = {}" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

206 
proof (rule ccontr) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

207 
assume "\<not> ?thesis" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

208 
then have "\<forall>i. \<exists>y. (i \<in> I \<longrightarrow> y \<in> F i) \<and> (i \<notin> I \<longrightarrow> y = undefined)" by auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

209 
from choice[OF this] guess f .. 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

210 
then have "f \<in> Pi\<^isub>E I F" by (auto simp: extensional_def) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

211 
with `Pi\<^isub>E I F = {}` show False by auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

212 
qed 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

213 
qed auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

214 

3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

215 
lemma Pi_eq_iff: 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

216 
"Pi\<^isub>E I F = Pi\<^isub>E I F' \<longleftrightarrow> (\<forall>i\<in>I. F i = F' i) \<or> ((\<exists>i\<in>I. F i = {}) \<and> (\<exists>i\<in>I. F' i = {}))" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

217 
proof (intro iffI disjCI) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

218 
assume eq[simp]: "Pi\<^isub>E I F = Pi\<^isub>E I F'" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

219 
assume "\<not> ((\<exists>i\<in>I. F i = {}) \<and> (\<exists>i\<in>I. F' i = {}))" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

220 
then have "(\<forall>i\<in>I. F i \<noteq> {}) \<and> (\<forall>i\<in>I. F' i \<noteq> {})" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

221 
using Pi_eq_empty_iff[of I F] Pi_eq_empty_iff[of I F'] by auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

222 
with Pi_eq_iff_not_empty[of I F F'] show "\<forall>i\<in>I. F i = F' i" by auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

223 
next 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

224 
assume "(\<forall>i\<in>I. F i = F' i) \<or> (\<exists>i\<in>I. F i = {}) \<and> (\<exists>i\<in>I. F' i = {})" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

225 
then show "Pi\<^isub>E I F = Pi\<^isub>E I F'" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

226 
using Pi_eq_empty_iff[of I F] Pi_eq_empty_iff[of I F'] by auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

227 
qed 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

228 

40859  229 
section "Finite product spaces" 
230 

231 
section "Products" 

232 

47694  233 
definition prod_emb where 
234 
"prod_emb I M K X = (\<lambda>x. restrict x K) ` X \<inter> (PIE i:I. space (M i))" 

235 

236 
lemma prod_emb_iff: 

237 
"f \<in> prod_emb I M K X \<longleftrightarrow> f \<in> extensional I \<and> (restrict f K \<in> X) \<and> (\<forall>i\<in>I. f i \<in> space (M i))" 

238 
unfolding prod_emb_def by auto 

40859  239 

47694  240 
lemma 
241 
shows prod_emb_empty[simp]: "prod_emb M L K {} = {}" 

242 
and prod_emb_Un[simp]: "prod_emb M L K (A \<union> B) = prod_emb M L K A \<union> prod_emb M L K B" 

243 
and prod_emb_Int: "prod_emb M L K (A \<inter> B) = prod_emb M L K A \<inter> prod_emb M L K B" 

244 
and prod_emb_UN[simp]: "prod_emb M L K (\<Union>i\<in>I. F i) = (\<Union>i\<in>I. prod_emb M L K (F i))" 

245 
and prod_emb_INT[simp]: "I \<noteq> {} \<Longrightarrow> prod_emb M L K (\<Inter>i\<in>I. F i) = (\<Inter>i\<in>I. prod_emb M L K (F i))" 

246 
and prod_emb_Diff[simp]: "prod_emb M L K (A  B) = prod_emb M L K A  prod_emb M L K B" 

247 
by (auto simp: prod_emb_def) 

40859  248 

47694  249 
lemma prod_emb_PiE: "J \<subseteq> I \<Longrightarrow> (\<And>i. i \<in> J \<Longrightarrow> E i \<subseteq> space (M i)) \<Longrightarrow> 
250 
prod_emb I M J (\<Pi>\<^isub>E i\<in>J. E i) = (\<Pi>\<^isub>E i\<in>I. if i \<in> J then E i else space (M i))" 

251 
by (force simp: prod_emb_def Pi_iff split_if_mem2) 

252 

253 
lemma prod_emb_PiE_same_index[simp]: "(\<And>i. i \<in> I \<Longrightarrow> E i \<subseteq> space (M i)) \<Longrightarrow> prod_emb I M I (Pi\<^isub>E I E) = Pi\<^isub>E I E" 

254 
by (auto simp: prod_emb_def Pi_iff) 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

255 

47694  256 
definition PiM :: "'i set \<Rightarrow> ('i \<Rightarrow> 'a measure) \<Rightarrow> ('i \<Rightarrow> 'a) measure" where 
257 
"PiM I M = extend_measure (\<Pi>\<^isub>E i\<in>I. space (M i)) 

258 
{(J, X). (J \<noteq> {} \<or> I = {}) \<and> finite J \<and> J \<subseteq> I \<and> X \<in> (\<Pi> j\<in>J. sets (M j))} 

259 
(\<lambda>(J, X). prod_emb I M J (\<Pi>\<^isub>E j\<in>J. X j)) 

260 
(\<lambda>(J, X). \<Prod>j\<in>J \<union> {i\<in>I. emeasure (M i) (space (M i)) \<noteq> 1}. if j \<in> J then emeasure (M j) (X j) else emeasure (M j) (space (M j)))" 

261 

262 
definition prod_algebra :: "'i set \<Rightarrow> ('i \<Rightarrow> 'a measure) \<Rightarrow> ('i \<Rightarrow> 'a) set set" where 

263 
"prod_algebra I M = (\<lambda>(J, X). prod_emb I M J (\<Pi>\<^isub>E j\<in>J. X j)) ` 

264 
{(J, X). (J \<noteq> {} \<or> I = {}) \<and> finite J \<and> J \<subseteq> I \<and> X \<in> (\<Pi> j\<in>J. sets (M j))}" 

265 

266 
abbreviation 

267 
"Pi\<^isub>M I M \<equiv> PiM I M" 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

268 

40859  269 
syntax 
47694  270 
"_PiM" :: "pttrn \<Rightarrow> 'i set \<Rightarrow> 'a measure \<Rightarrow> ('i => 'a) measure" ("(3PIM _:_./ _)" 10) 
40859  271 

272 
syntax (xsymbols) 

47694  273 
"_PiM" :: "pttrn \<Rightarrow> 'i set \<Rightarrow> 'a measure \<Rightarrow> ('i => 'a) measure" ("(3\<Pi>\<^isub>M _\<in>_./ _)" 10) 
40859  274 

275 
syntax (HTML output) 

47694  276 
"_PiM" :: "pttrn \<Rightarrow> 'i set \<Rightarrow> 'a measure \<Rightarrow> ('i => 'a) measure" ("(3\<Pi>\<^isub>M _\<in>_./ _)" 10) 
40859  277 

278 
translations 

47694  279 
"PIM x:I. M" == "CONST PiM I (%x. M)" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

280 

47694  281 
lemma prod_algebra_sets_into_space: 
282 
"prod_algebra I M \<subseteq> Pow (\<Pi>\<^isub>E i\<in>I. space (M i))" 

283 
using assms by (auto simp: prod_emb_def prod_algebra_def) 

40859  284 

47694  285 
lemma prod_algebra_eq_finite: 
286 
assumes I: "finite I" 

287 
shows "prod_algebra I M = {(\<Pi>\<^isub>E i\<in>I. X i) X. X \<in> (\<Pi> j\<in>I. sets (M j))}" (is "?L = ?R") 

288 
proof (intro iffI set_eqI) 

289 
fix A assume "A \<in> ?L" 

290 
then obtain J E where J: "J \<noteq> {} \<or> I = {}" "finite J" "J \<subseteq> I" "\<forall>i\<in>J. E i \<in> sets (M i)" 

291 
and A: "A = prod_emb I M J (PIE j:J. E j)" 

292 
by (auto simp: prod_algebra_def) 

293 
let ?A = "\<Pi>\<^isub>E i\<in>I. if i \<in> J then E i else space (M i)" 

294 
have A: "A = ?A" 

295 
unfolding A using J by (intro prod_emb_PiE sets_into_space) auto 

296 
show "A \<in> ?R" unfolding A using J top 

297 
by (intro CollectI exI[of _ "\<lambda>i. if i \<in> J then E i else space (M i)"]) simp 

298 
next 

299 
fix A assume "A \<in> ?R" 

300 
then obtain X where "A = (\<Pi>\<^isub>E i\<in>I. X i)" and X: "X \<in> (\<Pi> j\<in>I. sets (M j))" by auto 

301 
then have A: "A = prod_emb I M I (\<Pi>\<^isub>E i\<in>I. X i)" 

302 
using sets_into_space by (force simp: prod_emb_def Pi_iff) 

303 
from X I show "A \<in> ?L" unfolding A 

304 
by (auto simp: prod_algebra_def) 

305 
qed 

41095  306 

47694  307 
lemma prod_algebraI: 
308 
"finite J \<Longrightarrow> (J \<noteq> {} \<or> I = {}) \<Longrightarrow> J \<subseteq> I \<Longrightarrow> (\<And>i. i \<in> J \<Longrightarrow> E i \<in> sets (M i)) 

309 
\<Longrightarrow> prod_emb I M J (PIE j:J. E j) \<in> prod_algebra I M" 

310 
by (auto simp: prod_algebra_def Pi_iff) 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

311 

47694  312 
lemma prod_algebraE: 
313 
assumes A: "A \<in> prod_algebra I M" 

314 
obtains J E where "A = prod_emb I M J (PIE j:J. E j)" 

315 
"finite J" "J \<noteq> {} \<or> I = {}" "J \<subseteq> I" "\<And>i. i \<in> J \<Longrightarrow> E i \<in> sets (M i)" 

316 
using A by (auto simp: prod_algebra_def) 

42988  317 

47694  318 
lemma prod_algebraE_all: 
319 
assumes A: "A \<in> prod_algebra I M" 

320 
obtains E where "A = Pi\<^isub>E I E" "E \<in> (\<Pi> i\<in>I. sets (M i))" 

321 
proof  

322 
from A obtain E J where A: "A = prod_emb I M J (Pi\<^isub>E J E)" 

323 
and J: "J \<subseteq> I" and E: "E \<in> (\<Pi> i\<in>J. sets (M i))" 

324 
by (auto simp: prod_algebra_def) 

325 
from E have "\<And>i. i \<in> J \<Longrightarrow> E i \<subseteq> space (M i)" 

326 
using sets_into_space by auto 

327 
then have "A = (\<Pi>\<^isub>E i\<in>I. if i\<in>J then E i else space (M i))" 

328 
using A J by (auto simp: prod_emb_PiE) 

329 
moreover then have "(\<lambda>i. if i\<in>J then E i else space (M i)) \<in> (\<Pi> i\<in>I. sets (M i))" 

330 
using top E by auto 

331 
ultimately show ?thesis using that by auto 

332 
qed 

40859  333 

47694  334 
lemma Int_stable_prod_algebra: "Int_stable (prod_algebra I M)" 
335 
proof (unfold Int_stable_def, safe) 

336 
fix A assume "A \<in> prod_algebra I M" 

337 
from prod_algebraE[OF this] guess J E . note A = this 

338 
fix B assume "B \<in> prod_algebra I M" 

339 
from prod_algebraE[OF this] guess K F . note B = this 

340 
have "A \<inter> B = prod_emb I M (J \<union> K) (\<Pi>\<^isub>E i\<in>J \<union> K. (if i \<in> J then E i else space (M i)) \<inter> 

341 
(if i \<in> K then F i else space (M i)))" 

342 
unfolding A B using A(2,3,4) A(5)[THEN sets_into_space] B(2,3,4) B(5)[THEN sets_into_space] 

343 
apply (subst (1 2 3) prod_emb_PiE) 

344 
apply (simp_all add: subset_eq PiE_Int) 

345 
apply blast 

346 
apply (intro PiE_cong) 

347 
apply auto 

348 
done 

349 
also have "\<dots> \<in> prod_algebra I M" 

350 
using A B by (auto intro!: prod_algebraI) 

351 
finally show "A \<inter> B \<in> prod_algebra I M" . 

352 
qed 

353 

354 
lemma prod_algebra_mono: 

355 
assumes space: "\<And>i. i \<in> I \<Longrightarrow> space (E i) = space (F i)" 

356 
assumes sets: "\<And>i. i \<in> I \<Longrightarrow> sets (E i) \<subseteq> sets (F i)" 

357 
shows "prod_algebra I E \<subseteq> prod_algebra I F" 

358 
proof 

359 
fix A assume "A \<in> prod_algebra I E" 

360 
then obtain J G where J: "J \<noteq> {} \<or> I = {}" "finite J" "J \<subseteq> I" 

361 
and A: "A = prod_emb I E J (\<Pi>\<^isub>E i\<in>J. G i)" 

362 
and G: "\<And>i. i \<in> J \<Longrightarrow> G i \<in> sets (E i)" 

363 
by (auto simp: prod_algebra_def) 

364 
moreover 

365 
from space have "(\<Pi>\<^isub>E i\<in>I. space (E i)) = (\<Pi>\<^isub>E i\<in>I. space (F i))" 

366 
by (rule PiE_cong) 

367 
with A have "A = prod_emb I F J (\<Pi>\<^isub>E i\<in>J. G i)" 

368 
by (simp add: prod_emb_def) 

369 
moreover 

370 
from sets G J have "\<And>i. i \<in> J \<Longrightarrow> G i \<in> sets (F i)" 

371 
by auto 

372 
ultimately show "A \<in> prod_algebra I F" 

373 
apply (simp add: prod_algebra_def image_iff) 

374 
apply (intro exI[of _ J] exI[of _ G] conjI) 

375 
apply auto 

376 
done 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

377 
qed 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

378 

47694  379 
lemma space_PiM: "space (\<Pi>\<^isub>M i\<in>I. M i) = (\<Pi>\<^isub>E i\<in>I. space (M i))" 
380 
using prod_algebra_sets_into_space unfolding PiM_def prod_algebra_def by (intro space_extend_measure) simp 

381 

382 
lemma sets_PiM: "sets (\<Pi>\<^isub>M i\<in>I. M i) = sigma_sets (\<Pi>\<^isub>E i\<in>I. space (M i)) (prod_algebra I M)" 

383 
using prod_algebra_sets_into_space unfolding PiM_def prod_algebra_def by (intro sets_extend_measure) simp 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

384 

47694  385 
lemma sets_PiM_single: "sets (PiM I M) = 
386 
sigma_sets (\<Pi>\<^isub>E i\<in>I. space (M i)) {{f\<in>\<Pi>\<^isub>E i\<in>I. space (M i). f i \<in> A}  i A. i \<in> I \<and> A \<in> sets (M i)}" 

387 
(is "_ = sigma_sets ?\<Omega> ?R") 

388 
unfolding sets_PiM 

389 
proof (rule sigma_sets_eqI) 

390 
interpret R: sigma_algebra ?\<Omega> "sigma_sets ?\<Omega> ?R" by (rule sigma_algebra_sigma_sets) auto 

391 
fix A assume "A \<in> prod_algebra I M" 

392 
from prod_algebraE[OF this] guess J X . note X = this 

393 
show "A \<in> sigma_sets ?\<Omega> ?R" 

394 
proof cases 

395 
assume "I = {}" 

396 
with X have "A = {\<lambda>x. undefined}" by (auto simp: prod_emb_def) 

397 
with `I = {}` show ?thesis by (auto intro!: sigma_sets_top) 

398 
next 

399 
assume "I \<noteq> {}" 

400 
with X have "A = (\<Inter>j\<in>J. {f\<in>(\<Pi>\<^isub>E i\<in>I. space (M i)). f j \<in> X j})" 

401 
using sets_into_space[OF X(5)] 

402 
by (auto simp: prod_emb_PiE[OF _ sets_into_space] Pi_iff split: split_if_asm) blast 

403 
also have "\<dots> \<in> sigma_sets ?\<Omega> ?R" 

404 
using X `I \<noteq> {}` by (intro R.finite_INT sigma_sets.Basic) auto 

405 
finally show "A \<in> sigma_sets ?\<Omega> ?R" . 

406 
qed 

407 
next 

408 
fix A assume "A \<in> ?R" 

409 
then obtain i B where A: "A = {f\<in>\<Pi>\<^isub>E i\<in>I. space (M i). f i \<in> B}" "i \<in> I" "B \<in> sets (M i)" 

410 
by auto 

411 
then have "A = prod_emb I M {i} (\<Pi>\<^isub>E i\<in>{i}. B)" 

412 
using sets_into_space[OF A(3)] 

413 
apply (subst prod_emb_PiE) 

414 
apply (auto simp: Pi_iff split: split_if_asm) 

415 
apply blast 

416 
done 

417 
also have "\<dots> \<in> sigma_sets ?\<Omega> (prod_algebra I M)" 

418 
using A by (intro sigma_sets.Basic prod_algebraI) auto 

419 
finally show "A \<in> sigma_sets ?\<Omega> (prod_algebra I M)" . 

420 
qed 

421 

422 
lemma sets_PiM_I: 

423 
assumes "finite J" "J \<subseteq> I" "\<forall>i\<in>J. E i \<in> sets (M i)" 

424 
shows "prod_emb I M J (PIE j:J. E j) \<in> sets (PIM i:I. M i)" 

425 
proof cases 

426 
assume "J = {}" 

427 
then have "prod_emb I M J (PIE j:J. E j) = (PIE j:I. space (M j))" 

428 
by (auto simp: prod_emb_def) 

429 
then show ?thesis 

430 
by (auto simp add: sets_PiM intro!: sigma_sets_top) 

431 
next 

432 
assume "J \<noteq> {}" with assms show ?thesis 

50003  433 
by (force simp add: sets_PiM prod_algebra_def) 
40859  434 
qed 
435 

47694  436 
lemma measurable_PiM: 
437 
assumes space: "f \<in> space N \<rightarrow> (\<Pi>\<^isub>E i\<in>I. space (M i))" 

438 
assumes sets: "\<And>X J. J \<noteq> {} \<or> I = {} \<Longrightarrow> finite J \<Longrightarrow> J \<subseteq> I \<Longrightarrow> (\<And>i. i \<in> J \<Longrightarrow> X i \<in> sets (M i)) \<Longrightarrow> 

439 
f ` prod_emb I M J (Pi\<^isub>E J X) \<inter> space N \<in> sets N" 

440 
shows "f \<in> measurable N (PiM I M)" 

441 
using sets_PiM prod_algebra_sets_into_space space 

442 
proof (rule measurable_sigma_sets) 

443 
fix A assume "A \<in> prod_algebra I M" 

444 
from prod_algebraE[OF this] guess J X . 

445 
with sets[of J X] show "f ` A \<inter> space N \<in> sets N" by auto 

446 
qed 

447 

448 
lemma measurable_PiM_Collect: 

449 
assumes space: "f \<in> space N \<rightarrow> (\<Pi>\<^isub>E i\<in>I. space (M i))" 

450 
assumes sets: "\<And>X J. J \<noteq> {} \<or> I = {} \<Longrightarrow> finite J \<Longrightarrow> J \<subseteq> I \<Longrightarrow> (\<And>i. i \<in> J \<Longrightarrow> X i \<in> sets (M i)) \<Longrightarrow> 

451 
{\<omega>\<in>space N. \<forall>i\<in>J. f \<omega> i \<in> X i} \<in> sets N" 

452 
shows "f \<in> measurable N (PiM I M)" 

453 
using sets_PiM prod_algebra_sets_into_space space 

454 
proof (rule measurable_sigma_sets) 

455 
fix A assume "A \<in> prod_algebra I M" 

456 
from prod_algebraE[OF this] guess J X . note X = this 

457 
have "f ` A \<inter> space N = {\<omega> \<in> space N. \<forall>i\<in>J. f \<omega> i \<in> X i}" 

458 
using sets_into_space[OF X(5)] X(2) space unfolding X(1) 

459 
by (subst prod_emb_PiE) (auto simp: Pi_iff split: split_if_asm) 

460 
also have "\<dots> \<in> sets N" using X(3,2,4,5) by (rule sets) 

461 
finally show "f ` A \<inter> space N \<in> sets N" . 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

462 
qed 
41095  463 

47694  464 
lemma measurable_PiM_single: 
465 
assumes space: "f \<in> space N \<rightarrow> (\<Pi>\<^isub>E i\<in>I. space (M i))" 

466 
assumes sets: "\<And>A i. i \<in> I \<Longrightarrow> A \<in> sets (M i) \<Longrightarrow> {\<omega> \<in> space N. f \<omega> i \<in> A} \<in> sets N" 

467 
shows "f \<in> measurable N (PiM I M)" 

468 
using sets_PiM_single 

469 
proof (rule measurable_sigma_sets) 

470 
fix A assume "A \<in> {{f \<in> \<Pi>\<^isub>E i\<in>I. space (M i). f i \<in> A} i A. i \<in> I \<and> A \<in> sets (M i)}" 

471 
then obtain B i where "A = {f \<in> \<Pi>\<^isub>E i\<in>I. space (M i). f i \<in> B}" and B: "i \<in> I" "B \<in> sets (M i)" 

472 
by auto 

473 
with space have "f ` A \<inter> space N = {\<omega> \<in> space N. f \<omega> i \<in> B}" by auto 

474 
also have "\<dots> \<in> sets N" using B by (rule sets) 

475 
finally show "f ` A \<inter> space N \<in> sets N" . 

476 
qed (auto simp: space) 

40859  477 

50003  478 
lemma sets_PiM_I_finite[measurable]: 
47694  479 
assumes "finite I" and sets: "(\<And>i. i \<in> I \<Longrightarrow> E i \<in> sets (M i))" 
480 
shows "(PIE j:I. E j) \<in> sets (PIM i:I. M i)" 

481 
using sets_PiM_I[of I I E M] sets_into_space[OF sets] `finite I` sets by auto 

482 

50021
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

483 
lemma measurable_component_singleton: 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

484 
assumes "i \<in> I" shows "(\<lambda>x. x i) \<in> measurable (Pi\<^isub>M I M) (M i)" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

485 
proof (unfold measurable_def, intro CollectI conjI ballI) 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

486 
fix A assume "A \<in> sets (M i)" 
47694  487 
then have "(\<lambda>x. x i) ` A \<inter> space (Pi\<^isub>M I M) = prod_emb I M {i} (\<Pi>\<^isub>E j\<in>{i}. A)" 
488 
using sets_into_space `i \<in> I` 

489 
by (fastforce dest: Pi_mem simp: prod_emb_def space_PiM split: split_if_asm) 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

490 
then show "(\<lambda>x. x i) ` A \<inter> space (Pi\<^isub>M I M) \<in> sets (Pi\<^isub>M I M)" 
47694  491 
using `A \<in> sets (M i)` `i \<in> I` by (auto intro!: sets_PiM_I) 
492 
qed (insert `i \<in> I`, auto simp: space_PiM) 

493 

50021
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

494 
lemma measurable_component_singleton'[measurable_app]: 
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

495 
assumes f: "f \<in> measurable N (Pi\<^isub>M I M)" 
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

496 
assumes i: "i \<in> I" 
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

497 
shows "(\<lambda>x. (f x) i) \<in> measurable N (M i)" 
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

498 
using measurable_compose[OF f measurable_component_singleton, OF i] . 
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

499 

d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

500 
lemma measurable_nat_case[measurable (raw)]: 
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

501 
assumes [measurable (raw)]: "i = 0 \<Longrightarrow> f \<in> measurable M N" 
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

502 
"\<And>j. i = Suc j \<Longrightarrow> (\<lambda>x. g x j) \<in> measurable M N" 
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

503 
shows "(\<lambda>x. nat_case (f x) (g x) i) \<in> measurable M N" 
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

504 
by (cases i) simp_all 
d96a3f468203
add support for function application to measurability prover
hoelzl
parents:
50003
diff
changeset

505 

50003  506 
lemma measurable_add_dim[measurable]: 
49776  507 
"(\<lambda>(f, y). f(i := y)) \<in> measurable (Pi\<^isub>M I M \<Otimes>\<^isub>M M i) (Pi\<^isub>M (insert i I) M)" 
47694  508 
(is "?f \<in> measurable ?P ?I") 
509 
proof (rule measurable_PiM_single) 

510 
fix j A assume j: "j \<in> insert i I" and A: "A \<in> sets (M j)" 

511 
have "{\<omega> \<in> space ?P. (\<lambda>(f, x). fun_upd f i x) \<omega> j \<in> A} = 

512 
(if j = i then space (Pi\<^isub>M I M) \<times> A else ((\<lambda>x. x j) \<circ> fst) ` A \<inter> space ?P)" 

513 
using sets_into_space[OF A] by (auto simp add: space_pair_measure space_PiM) 

514 
also have "\<dots> \<in> sets ?P" 

515 
using A j 

516 
by (auto intro!: measurable_sets[OF measurable_comp, OF _ measurable_component_singleton]) 

517 
finally show "{\<omega> \<in> space ?P. prod_case (\<lambda>f. fun_upd f i) \<omega> j \<in> A} \<in> sets ?P" . 

518 
qed (auto simp: space_pair_measure space_PiM) 

41661  519 

50003  520 
lemma measurable_component_update: 
521 
"x \<in> space (Pi\<^isub>M I M) \<Longrightarrow> i \<notin> I \<Longrightarrow> (\<lambda>v. x(i := v)) \<in> measurable (M i) (Pi\<^isub>M (insert i I) M)" 

522 
by simp 

523 

524 
lemma measurable_merge[measurable]: 

49780  525 
"merge I J \<in> measurable (Pi\<^isub>M I M \<Otimes>\<^isub>M Pi\<^isub>M J M) (Pi\<^isub>M (I \<union> J) M)" 
47694  526 
(is "?f \<in> measurable ?P ?U") 
527 
proof (rule measurable_PiM_single) 

528 
fix i A assume A: "A \<in> sets (M i)" "i \<in> I \<union> J" 

49780  529 
then have "{\<omega> \<in> space ?P. merge I J \<omega> i \<in> A} = 
47694  530 
(if i \<in> I then ((\<lambda>x. x i) \<circ> fst) ` A \<inter> space ?P else ((\<lambda>x. x i) \<circ> snd) ` A \<inter> space ?P)" 
49776  531 
by (auto simp: merge_def) 
47694  532 
also have "\<dots> \<in> sets ?P" 
533 
using A 

534 
by (auto intro!: measurable_sets[OF measurable_comp, OF _ measurable_component_singleton]) 

49780  535 
finally show "{\<omega> \<in> space ?P. merge I J \<omega> i \<in> A} \<in> sets ?P" . 
49776  536 
qed (auto simp: space_pair_measure space_PiM Pi_iff merge_def extensional_def) 
42988  537 

50003  538 
lemma measurable_restrict[measurable (raw)]: 
47694  539 
assumes X: "\<And>i. i \<in> I \<Longrightarrow> X i \<in> measurable N (M i)" 
540 
shows "(\<lambda>x. \<lambda>i\<in>I. X i x) \<in> measurable N (Pi\<^isub>M I M)" 

541 
proof (rule measurable_PiM_single) 

542 
fix A i assume A: "i \<in> I" "A \<in> sets (M i)" 

543 
then have "{\<omega> \<in> space N. (\<lambda>i\<in>I. X i \<omega>) i \<in> A} = X i ` A \<inter> space N" 

544 
by auto 

545 
then show "{\<omega> \<in> space N. (\<lambda>i\<in>I. X i \<omega>) i \<in> A} \<in> sets N" 

546 
using A X by (auto intro!: measurable_sets) 

547 
qed (insert X, auto dest: measurable_space) 

548 

50003  549 
lemma sets_in_Pi_aux: 
550 
"finite I \<Longrightarrow> (\<And>j. j \<in> I \<Longrightarrow> {x\<in>space (M j). x \<in> F j} \<in> sets (M j)) \<Longrightarrow> 

551 
{x\<in>space (PiM I M). x \<in> Pi I F} \<in> sets (PiM I M)" 

552 
by (simp add: subset_eq Pi_iff) 

553 

554 
lemma sets_in_Pi[measurable (raw)]: 

555 
"finite I \<Longrightarrow> f \<in> measurable N (PiM I M) \<Longrightarrow> 

556 
(\<And>j. j \<in> I \<Longrightarrow> {x\<in>space (M j). x \<in> F j} \<in> sets (M j)) \<Longrightarrow> 

557 
Sigma_Algebra.pred N (\<lambda>x. f x \<in> Pi I F)" 

558 
unfolding pred_def 

559 
by (rule measurable_sets_Collect[of f N "PiM I M", OF _ sets_in_Pi_aux]) auto 

560 

561 
lemma sets_in_extensional_aux: 

562 
"{x\<in>space (PiM I M). x \<in> extensional I} \<in> sets (PiM I M)" 

563 
proof  

564 
have "{x\<in>space (PiM I M). x \<in> extensional I} = space (PiM I M)" 

565 
by (auto simp add: extensional_def space_PiM) 

566 
then show ?thesis by simp 

567 
qed 

568 

569 
lemma sets_in_extensional[measurable (raw)]: 

570 
"f \<in> measurable N (PiM I M) \<Longrightarrow> Sigma_Algebra.pred N (\<lambda>x. f x \<in> extensional I)" 

571 
unfolding pred_def 

572 
by (rule measurable_sets_Collect[of f N "PiM I M", OF _ sets_in_extensional_aux]) auto 

573 

47694  574 
locale product_sigma_finite = 
575 
fixes M :: "'i \<Rightarrow> 'a measure" 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

576 
assumes sigma_finite_measures: "\<And>i. sigma_finite_measure (M i)" 
40859  577 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

578 
sublocale product_sigma_finite \<subseteq> M: sigma_finite_measure "M i" for i 
40859  579 
by (rule sigma_finite_measures) 
580 

47694  581 
locale finite_product_sigma_finite = product_sigma_finite M for M :: "'i \<Rightarrow> 'a measure" + 
582 
fixes I :: "'i set" 

583 
assumes finite_index: "finite I" 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

584 

40859  585 
lemma (in finite_product_sigma_finite) sigma_finite_pairs: 
586 
"\<exists>F::'i \<Rightarrow> nat \<Rightarrow> 'a set. 

587 
(\<forall>i\<in>I. range (F i) \<subseteq> sets (M i)) \<and> 

47694  588 
(\<forall>k. \<forall>i\<in>I. emeasure (M i) (F i k) \<noteq> \<infinity>) \<and> incseq (\<lambda>k. \<Pi>\<^isub>E i\<in>I. F i k) \<and> 
589 
(\<Union>k. \<Pi>\<^isub>E i\<in>I. F i k) = space (PiM I M)" 

40859  590 
proof  
47694  591 
have "\<forall>i::'i. \<exists>F::nat \<Rightarrow> 'a set. range F \<subseteq> sets (M i) \<and> incseq F \<and> (\<Union>i. F i) = space (M i) \<and> (\<forall>k. emeasure (M i) (F k) \<noteq> \<infinity>)" 
592 
using M.sigma_finite_incseq by metis 

40859  593 
from choice[OF this] guess F :: "'i \<Rightarrow> nat \<Rightarrow> 'a set" .. 
47694  594 
then have F: "\<And>i. range (F i) \<subseteq> sets (M i)" "\<And>i. incseq (F i)" "\<And>i. (\<Union>j. F i j) = space (M i)" "\<And>i k. emeasure (M i) (F i k) \<noteq> \<infinity>" 
40859  595 
by auto 
596 
let ?F = "\<lambda>k. \<Pi>\<^isub>E i\<in>I. F i k" 

47694  597 
note space_PiM[simp] 
40859  598 
show ?thesis 
41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

599 
proof (intro exI[of _ F] conjI allI incseq_SucI set_eqI iffI ballI) 
40859  600 
fix i show "range (F i) \<subseteq> sets (M i)" by fact 
601 
next 

47694  602 
fix i k show "emeasure (M i) (F i k) \<noteq> \<infinity>" by fact 
40859  603 
next 
47694  604 
fix A assume "A \<in> (\<Union>i. ?F i)" then show "A \<in> space (PiM I M)" 
605 
using `\<And>i. range (F i) \<subseteq> sets (M i)` sets_into_space 

606 
by auto blast 

40859  607 
next 
47694  608 
fix f assume "f \<in> space (PiM I M)" 
41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

609 
with Pi_UN[OF finite_index, of "\<lambda>k i. F i k"] F 
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

610 
show "f \<in> (\<Union>i. ?F i)" by (auto simp: incseq_def) 
40859  611 
next 
612 
fix i show "?F i \<subseteq> ?F (Suc i)" 

41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

613 
using `\<And>i. incseq (F i)`[THEN incseq_SucD] by auto 
40859  614 
qed 
615 
qed 

616 

49780  617 
lemma 
618 
shows space_PiM_empty: "space (Pi\<^isub>M {} M) = {\<lambda>k. undefined}" 

619 
and sets_PiM_empty: "sets (Pi\<^isub>M {} M) = { {}, {\<lambda>k. undefined} }" 

620 
by (simp_all add: space_PiM sets_PiM_single image_constant sigma_sets_empty_eq) 

621 

622 
lemma emeasure_PiM_empty[simp]: "emeasure (PiM {} M) {\<lambda>_. undefined} = 1" 

623 
proof  

624 
let ?\<mu> = "\<lambda>A. if A = {} then 0 else (1::ereal)" 

625 
have "emeasure (Pi\<^isub>M {} M) (prod_emb {} M {} (\<Pi>\<^isub>E i\<in>{}. {})) = 1" 

626 
proof (subst emeasure_extend_measure_Pair[OF PiM_def]) 

627 
show "positive (PiM {} M) ?\<mu>" 

628 
by (auto simp: positive_def) 

629 
show "countably_additive (PiM {} M) ?\<mu>" 

630 
by (rule countably_additiveI_finite) 

631 
(auto simp: additive_def positive_def sets_PiM_empty space_PiM_empty intro!: ) 

632 
qed (auto simp: prod_emb_def) 

633 
also have "(prod_emb {} M {} (\<Pi>\<^isub>E i\<in>{}. {})) = {\<lambda>_. undefined}" 

634 
by (auto simp: prod_emb_def) 

635 
finally show ?thesis 

636 
by simp 

637 
qed 

638 

639 
lemma PiM_empty: "PiM {} M = count_space {\<lambda>_. undefined}" 

640 
by (rule measure_eqI) (auto simp add: sets_PiM_empty one_ereal_def) 

641 

49776  642 
lemma (in product_sigma_finite) emeasure_PiM: 
643 
"finite I \<Longrightarrow> (\<And>i. i\<in>I \<Longrightarrow> A i \<in> sets (M i)) \<Longrightarrow> emeasure (PiM I M) (Pi\<^isub>E I A) = (\<Prod>i\<in>I. emeasure (M i) (A i))" 

644 
proof (induct I arbitrary: A rule: finite_induct) 

40859  645 
case (insert i I) 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

646 
interpret finite_product_sigma_finite M I by default fact 
40859  647 
have "finite (insert i I)" using `finite I` by auto 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

648 
interpret I': finite_product_sigma_finite M "insert i I" by default fact 
41661  649 
let ?h = "(\<lambda>(f, y). f(i := y))" 
47694  650 

651 
let ?P = "distr (Pi\<^isub>M I M \<Otimes>\<^isub>M M i) (Pi\<^isub>M (insert i I) M) ?h" 

652 
let ?\<mu> = "emeasure ?P" 

653 
let ?I = "{j \<in> insert i I. emeasure (M j) (space (M j)) \<noteq> 1}" 

654 
let ?f = "\<lambda>J E j. if j \<in> J then emeasure (M j) (E j) else emeasure (M j) (space (M j))" 

655 

49776  656 
have "emeasure (Pi\<^isub>M (insert i I) M) (prod_emb (insert i I) M (insert i I) (Pi\<^isub>E (insert i I) A)) = 
657 
(\<Prod>i\<in>insert i I. emeasure (M i) (A i))" 

658 
proof (subst emeasure_extend_measure_Pair[OF PiM_def]) 

659 
fix J E assume "(J \<noteq> {} \<or> insert i I = {}) \<and> finite J \<and> J \<subseteq> insert i I \<and> E \<in> (\<Pi> j\<in>J. sets (M j))" 

660 
then have J: "J \<noteq> {}" "finite J" "J \<subseteq> insert i I" and E: "\<forall>j\<in>J. E j \<in> sets (M j)" by auto 

661 
let ?p = "prod_emb (insert i I) M J (Pi\<^isub>E J E)" 

662 
let ?p' = "prod_emb I M (J  {i}) (\<Pi>\<^isub>E j\<in>J{i}. E j)" 

663 
have "?\<mu> ?p = 

664 
emeasure (Pi\<^isub>M I M \<Otimes>\<^isub>M (M i)) (?h ` ?p \<inter> space (Pi\<^isub>M I M \<Otimes>\<^isub>M M i))" 

665 
by (intro emeasure_distr measurable_add_dim sets_PiM_I) fact+ 

666 
also have "?h ` ?p \<inter> space (Pi\<^isub>M I M \<Otimes>\<^isub>M M i) = ?p' \<times> (if i \<in> J then E i else space (M i))" 

667 
using J E[rule_format, THEN sets_into_space] 

668 
by (force simp: space_pair_measure space_PiM Pi_iff prod_emb_iff split: split_if_asm) 

669 
also have "emeasure (Pi\<^isub>M I M \<Otimes>\<^isub>M (M i)) (?p' \<times> (if i \<in> J then E i else space (M i))) = 

670 
emeasure (Pi\<^isub>M I M) ?p' * emeasure (M i) (if i \<in> J then (E i) else space (M i))" 

671 
using J E by (intro M.emeasure_pair_measure_Times sets_PiM_I) auto 

672 
also have "?p' = (\<Pi>\<^isub>E j\<in>I. if j \<in> J{i} then E j else space (M j))" 

673 
using J E[rule_format, THEN sets_into_space] 

674 
by (auto simp: prod_emb_iff Pi_iff split: split_if_asm) blast+ 

675 
also have "emeasure (Pi\<^isub>M I M) (\<Pi>\<^isub>E j\<in>I. if j \<in> J{i} then E j else space (M j)) = 

676 
(\<Prod> j\<in>I. if j \<in> J{i} then emeasure (M j) (E j) else emeasure (M j) (space (M j)))" 

677 
using E by (subst insert) (auto intro!: setprod_cong) 

678 
also have "(\<Prod>j\<in>I. if j \<in> J  {i} then emeasure (M j) (E j) else emeasure (M j) (space (M j))) * 

679 
emeasure (M i) (if i \<in> J then E i else space (M i)) = (\<Prod>j\<in>insert i I. ?f J E j)" 

680 
using insert by (auto simp: mult_commute intro!: arg_cong2[where f="op *"] setprod_cong) 

681 
also have "\<dots> = (\<Prod>j\<in>J \<union> ?I. ?f J E j)" 

682 
using insert(1,2) J E by (intro setprod_mono_one_right) auto 

683 
finally show "?\<mu> ?p = \<dots>" . 

47694  684 

49776  685 
show "prod_emb (insert i I) M J (Pi\<^isub>E J E) \<in> Pow (\<Pi>\<^isub>E i\<in>insert i I. space (M i))" 
686 
using J E[rule_format, THEN sets_into_space] by (auto simp: prod_emb_iff) 

687 
next 

688 
show "positive (sets (Pi\<^isub>M (insert i I) M)) ?\<mu>" "countably_additive (sets (Pi\<^isub>M (insert i I) M)) ?\<mu>" 

689 
using emeasure_positive[of ?P] emeasure_countably_additive[of ?P] by simp_all 

690 
next 

691 
show "(insert i I \<noteq> {} \<or> insert i I = {}) \<and> finite (insert i I) \<and> 

692 
insert i I \<subseteq> insert i I \<and> A \<in> (\<Pi> j\<in>insert i I. sets (M j))" 

693 
using insert by auto 

694 
qed (auto intro!: setprod_cong) 

695 
with insert show ?case 

696 
by (subst (asm) prod_emb_PiE_same_index) (auto intro!: sets_into_space) 

50003  697 
qed simp 
47694  698 

49776  699 
lemma (in product_sigma_finite) sigma_finite: 
700 
assumes "finite I" 

701 
shows "sigma_finite_measure (PiM I M)" 

702 
proof  

703 
interpret finite_product_sigma_finite M I by default fact 

704 

705 
from sigma_finite_pairs guess F :: "'i \<Rightarrow> nat \<Rightarrow> 'a set" .. 

706 
then have F: "\<And>j. j \<in> I \<Longrightarrow> range (F j) \<subseteq> sets (M j)" 

707 
"incseq (\<lambda>k. \<Pi>\<^isub>E j \<in> I. F j k)" 

708 
"(\<Union>k. \<Pi>\<^isub>E j \<in> I. F j k) = space (Pi\<^isub>M I M)" 

709 
"\<And>k. \<And>j. j \<in> I \<Longrightarrow> emeasure (M j) (F j k) \<noteq> \<infinity>" 

47694  710 
by blast+ 
49776  711 
let ?F = "\<lambda>k. \<Pi>\<^isub>E j \<in> I. F j k" 
47694  712 

49776  713 
show ?thesis 
47694  714 
proof (unfold_locales, intro exI[of _ ?F] conjI allI) 
49776  715 
show "range ?F \<subseteq> sets (Pi\<^isub>M I M)" using F(1) `finite I` by auto 
47694  716 
next 
49776  717 
from F(3) show "(\<Union>i. ?F i) = space (Pi\<^isub>M I M)" by simp 
47694  718 
next 
719 
fix j 

49776  720 
from F `finite I` setprod_PInf[of I, OF emeasure_nonneg, of M] 
721 
show "emeasure (\<Pi>\<^isub>M i\<in>I. M i) (?F j) \<noteq> \<infinity>" 

722 
by (subst emeasure_PiM) auto 

40859  723 
qed 
724 
qed 

725 

47694  726 
sublocale finite_product_sigma_finite \<subseteq> sigma_finite_measure "Pi\<^isub>M I M" 
727 
using sigma_finite[OF finite_index] . 

40859  728 

729 
lemma (in finite_product_sigma_finite) measure_times: 

47694  730 
"(\<And>i. i \<in> I \<Longrightarrow> A i \<in> sets (M i)) \<Longrightarrow> emeasure (Pi\<^isub>M I M) (Pi\<^isub>E I A) = (\<Prod>i\<in>I. emeasure (M i) (A i))" 
731 
using emeasure_PiM[OF finite_index] by auto 

41096  732 

40859  733 
lemma (in product_sigma_finite) positive_integral_empty: 
41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

734 
assumes pos: "0 \<le> f (\<lambda>k. undefined)" 
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

735 
shows "integral\<^isup>P (Pi\<^isub>M {} M) f = f (\<lambda>k. undefined)" 
40859  736 
proof  
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

737 
interpret finite_product_sigma_finite M "{}" by default (fact finite.emptyI) 
47694  738 
have "\<And>A. emeasure (Pi\<^isub>M {} M) (Pi\<^isub>E {} A) = 1" 
40859  739 
using assms by (subst measure_times) auto 
740 
then show ?thesis 

47694  741 
unfolding positive_integral_def simple_function_def simple_integral_def[abs_def] 
742 
proof (simp add: space_PiM_empty sets_PiM_empty, intro antisym) 

41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

743 
show "f (\<lambda>k. undefined) \<le> (SUP f:{g. g \<le> max 0 \<circ> f}. f (\<lambda>k. undefined))" 
44928
7ef6505bde7f
renamed Complete_Lattices lemmas, removed legacy names
hoelzl
parents:
44890
diff
changeset

744 
by (intro SUP_upper) (auto simp: le_fun_def split: split_max) 
41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

745 
show "(SUP f:{g. g \<le> max 0 \<circ> f}. f (\<lambda>k. undefined)) \<le> f (\<lambda>k. undefined)" using pos 
44928
7ef6505bde7f
renamed Complete_Lattices lemmas, removed legacy names
hoelzl
parents:
44890
diff
changeset

746 
by (intro SUP_least) (auto simp: le_fun_def simp: max_def split: split_if_asm) 
40859  747 
qed 
748 
qed 

749 

47694  750 
lemma (in product_sigma_finite) distr_merge: 
40859  751 
assumes IJ[simp]: "I \<inter> J = {}" and fin: "finite I" "finite J" 
49780  752 
shows "distr (Pi\<^isub>M I M \<Otimes>\<^isub>M Pi\<^isub>M J M) (Pi\<^isub>M (I \<union> J) M) (merge I J) = Pi\<^isub>M (I \<union> J) M" 
47694  753 
(is "?D = ?P") 
40859  754 
proof  
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

755 
interpret I: finite_product_sigma_finite M I by default fact 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

756 
interpret J: finite_product_sigma_finite M J by default fact 
40859  757 
have "finite (I \<union> J)" using fin by auto 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

758 
interpret IJ: finite_product_sigma_finite M "I \<union> J" by default fact 
47694  759 
interpret P: pair_sigma_finite "Pi\<^isub>M I M" "Pi\<^isub>M J M" by default 
49780  760 
let ?g = "merge I J" 
47694  761 

41661  762 
from IJ.sigma_finite_pairs obtain F where 
763 
F: "\<And>i. i\<in> I \<union> J \<Longrightarrow> range (F i) \<subseteq> sets (M i)" 

41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

764 
"incseq (\<lambda>k. \<Pi>\<^isub>E i\<in>I \<union> J. F i k)" 
47694  765 
"(\<Union>k. \<Pi>\<^isub>E i\<in>I \<union> J. F i k) = space ?P" 
766 
"\<And>k. \<forall>i\<in>I\<union>J. emeasure (M i) (F i k) \<noteq> \<infinity>" 

41661  767 
by auto 
768 
let ?F = "\<lambda>k. \<Pi>\<^isub>E i\<in>I \<union> J. F i k" 

47694  769 

770 
show ?thesis 

771 
proof (rule measure_eqI_generator_eq[symmetric]) 

772 
show "Int_stable (prod_algebra (I \<union> J) M)" 

773 
by (rule Int_stable_prod_algebra) 

774 
show "prod_algebra (I \<union> J) M \<subseteq> Pow (\<Pi>\<^isub>E i \<in> I \<union> J. space (M i))" 

775 
by (rule prod_algebra_sets_into_space) 

776 
show "sets ?P = sigma_sets (\<Pi>\<^isub>E i\<in>I \<union> J. space (M i)) (prod_algebra (I \<union> J) M)" 

777 
by (rule sets_PiM) 

778 
then show "sets ?D = sigma_sets (\<Pi>\<^isub>E i\<in>I \<union> J. space (M i)) (prod_algebra (I \<union> J) M)" 

779 
by simp 

780 

781 
show "range ?F \<subseteq> prod_algebra (I \<union> J) M" using F 

782 
using fin by (auto simp: prod_algebra_eq_finite) 

783 
show "(\<Union>i. \<Pi>\<^isub>E ia\<in>I \<union> J. F ia i) = (\<Pi>\<^isub>E i\<in>I \<union> J. space (M i))" 

784 
using F(3) by (simp add: space_PiM) 

41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

785 
next 
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

786 
fix k 
47694  787 
from F `finite I` setprod_PInf[of "I \<union> J", OF emeasure_nonneg, of M] 
788 
show "emeasure ?P (?F k) \<noteq> \<infinity>" by (subst IJ.measure_times) auto 

41661  789 
next 
47694  790 
fix A assume A: "A \<in> prod_algebra (I \<union> J) M" 
50003  791 
with fin obtain F where A_eq: "A = (Pi\<^isub>E (I \<union> J) F)" and F: "\<forall>i\<in>J. F i \<in> sets (M i)" "\<forall>i\<in>I. F i \<in> sets (M i)" 
47694  792 
by (auto simp add: prod_algebra_eq_finite) 
793 
let ?B = "Pi\<^isub>M I M \<Otimes>\<^isub>M Pi\<^isub>M J M" 

794 
let ?X = "?g ` A \<inter> space ?B" 

795 
have "Pi\<^isub>E I F \<subseteq> space (Pi\<^isub>M I M)" "Pi\<^isub>E J F \<subseteq> space (Pi\<^isub>M J M)" 

50003  796 
using F[rule_format, THEN sets_into_space] by (force simp: space_PiM)+ 
47694  797 
then have X: "?X = (Pi\<^isub>E I F \<times> Pi\<^isub>E J F)" 
798 
unfolding A_eq by (subst merge_vimage) (auto simp: space_pair_measure space_PiM) 

799 
have "emeasure ?D A = emeasure ?B ?X" 

800 
using A by (intro emeasure_distr measurable_merge) (auto simp: sets_PiM) 

801 
also have "emeasure ?B ?X = (\<Prod>i\<in>I. emeasure (M i) (F i)) * (\<Prod>i\<in>J. emeasure (M i) (F i))" 

50003  802 
using `finite J` `finite I` F unfolding X 
49776  803 
by (simp add: J.emeasure_pair_measure_Times I.measure_times J.measure_times Pi_iff) 
47694  804 
also have "\<dots> = (\<Prod>i\<in>I \<union> J. emeasure (M i) (F i))" 
41661  805 
using `finite J` `finite I` `I \<inter> J = {}` by (simp add: setprod_Un_one) 
47694  806 
also have "\<dots> = emeasure ?P (Pi\<^isub>E (I \<union> J) F)" 
41661  807 
using `finite J` `finite I` F unfolding A 
808 
by (intro IJ.measure_times[symmetric]) auto 

47694  809 
finally show "emeasure ?P A = emeasure ?D A" using A_eq by simp 
810 
qed 

41661  811 
qed 
41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

812 

bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

813 
lemma (in product_sigma_finite) product_positive_integral_fold: 
47694  814 
assumes IJ: "I \<inter> J = {}" "finite I" "finite J" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

815 
and f: "f \<in> borel_measurable (Pi\<^isub>M (I \<union> J) M)" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

816 
shows "integral\<^isup>P (Pi\<^isub>M (I \<union> J) M) f = 
49780  817 
(\<integral>\<^isup>+ x. (\<integral>\<^isup>+ y. f (merge I J (x, y)) \<partial>(Pi\<^isub>M J M)) \<partial>(Pi\<^isub>M I M))" 
41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

818 
proof  
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

819 
interpret I: finite_product_sigma_finite M I by default fact 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

820 
interpret J: finite_product_sigma_finite M J by default fact 
41831  821 
interpret P: pair_sigma_finite "Pi\<^isub>M I M" "Pi\<^isub>M J M" by default 
49780  822 
have P_borel: "(\<lambda>x. f (merge I J x)) \<in> borel_measurable (Pi\<^isub>M I M \<Otimes>\<^isub>M Pi\<^isub>M J M)" 
49776  823 
using measurable_comp[OF measurable_merge f] by (simp add: comp_def) 
41661  824 
show ?thesis 
47694  825 
apply (subst distr_merge[OF IJ, symmetric]) 
49776  826 
apply (subst positive_integral_distr[OF measurable_merge f]) 
49999
dfb63b9b8908
for the product measure it is enough if only one measure is sigmafinite
hoelzl
parents:
49784
diff
changeset

827 
apply (subst J.positive_integral_fst_measurable(2)[symmetric, OF P_borel]) 
47694  828 
apply simp 
829 
done 

40859  830 
qed 
831 

47694  832 
lemma (in product_sigma_finite) distr_singleton: 
833 
"distr (Pi\<^isub>M {i} M) (M i) (\<lambda>x. x i) = M i" (is "?D = _") 

834 
proof (intro measure_eqI[symmetric]) 

41831  835 
interpret I: finite_product_sigma_finite M "{i}" by default simp 
47694  836 
fix A assume A: "A \<in> sets (M i)" 
837 
moreover then have "(\<lambda>x. x i) ` A \<inter> space (Pi\<^isub>M {i} M) = (\<Pi>\<^isub>E i\<in>{i}. A)" 

838 
using sets_into_space by (auto simp: space_PiM) 

839 
ultimately show "emeasure (M i) A = emeasure ?D A" 

840 
using A I.measure_times[of "\<lambda>_. A"] 

841 
by (simp add: emeasure_distr measurable_component_singleton) 

842 
qed simp 

41831  843 

41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

844 
lemma (in product_sigma_finite) product_positive_integral_singleton: 
40859  845 
assumes f: "f \<in> borel_measurable (M i)" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

846 
shows "integral\<^isup>P (Pi\<^isub>M {i} M) (\<lambda>x. f (x i)) = integral\<^isup>P (M i) f" 
40859  847 
proof  
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

848 
interpret I: finite_product_sigma_finite M "{i}" by default simp 
47694  849 
from f show ?thesis 
850 
apply (subst distr_singleton[symmetric]) 

851 
apply (subst positive_integral_distr[OF measurable_component_singleton]) 

852 
apply simp_all 

853 
done 

40859  854 
qed 
855 

41096  856 
lemma (in product_sigma_finite) product_positive_integral_insert: 
49780  857 
assumes I[simp]: "finite I" "i \<notin> I" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

858 
and f: "f \<in> borel_measurable (Pi\<^isub>M (insert i I) M)" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

859 
shows "integral\<^isup>P (Pi\<^isub>M (insert i I) M) f = (\<integral>\<^isup>+ x. (\<integral>\<^isup>+ y. f (x(i := y)) \<partial>(M i)) \<partial>(Pi\<^isub>M I M))" 
41096  860 
proof  
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

861 
interpret I: finite_product_sigma_finite M I by default auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

862 
interpret i: finite_product_sigma_finite M "{i}" by default auto 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

863 
have IJ: "I \<inter> {i} = {}" and insert: "I \<union> {i} = insert i I" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

864 
using f by auto 
41096  865 
show ?thesis 
49780  866 
unfolding product_positive_integral_fold[OF IJ, unfolded insert, OF I(1) i.finite_index f] 
867 
proof (rule positive_integral_cong, subst product_positive_integral_singleton[symmetric]) 

47694  868 
fix x assume x: "x \<in> space (Pi\<^isub>M I M)" 
49780  869 
let ?f = "\<lambda>y. f (x(i := y))" 
870 
show "?f \<in> borel_measurable (M i)" 

47694  871 
using measurable_comp[OF measurable_component_update f, OF x `i \<notin> I`] 
872 
unfolding comp_def . 

49780  873 
show "(\<integral>\<^isup>+ y. f (merge I {i} (x, y)) \<partial>Pi\<^isub>M {i} M) = (\<integral>\<^isup>+ y. f (x(i := y i)) \<partial>Pi\<^isub>M {i} M)" 
874 
using x 

875 
by (auto intro!: positive_integral_cong arg_cong[where f=f] 

876 
simp add: space_PiM extensional_def) 

41096  877 
qed 
878 
qed 

879 

880 
lemma (in product_sigma_finite) product_positive_integral_setprod: 

43920  881 
fixes f :: "'i \<Rightarrow> 'a \<Rightarrow> ereal" 
41096  882 
assumes "finite I" and borel: "\<And>i. i \<in> I \<Longrightarrow> f i \<in> borel_measurable (M i)" 
41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

883 
and pos: "\<And>i x. i \<in> I \<Longrightarrow> 0 \<le> f i x" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

884 
shows "(\<integral>\<^isup>+ x. (\<Prod>i\<in>I. f i (x i)) \<partial>Pi\<^isub>M I M) = (\<Prod>i\<in>I. integral\<^isup>P (M i) (f i))" 
41096  885 
using assms proof induct 
886 
case (insert i I) 

887 
note `finite I`[intro, simp] 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

888 
interpret I: finite_product_sigma_finite M I by default auto 
41096  889 
have *: "\<And>x y. (\<Prod>j\<in>I. f j (if j = i then y else x j)) = (\<Prod>j\<in>I. f j (x j))" 
890 
using insert by (auto intro!: setprod_cong) 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

891 
have prod: "\<And>J. J \<subseteq> insert i I \<Longrightarrow> (\<lambda>x. (\<Prod>i\<in>J. f i (x i))) \<in> borel_measurable (Pi\<^isub>M J M)" 
41096  892 
using sets_into_space insert 
47694  893 
by (intro borel_measurable_ereal_setprod 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

894 
measurable_comp[OF measurable_component_singleton, unfolded comp_def]) 
41096  895 
auto 
41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

896 
then show ?case 
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

897 
apply (simp add: product_positive_integral_insert[OF insert(1,2) prod]) 
47694  898 
apply (simp add: insert(2) * pos borel setprod_ereal_pos positive_integral_multc) 
899 
apply (subst positive_integral_cmult) 

900 
apply (auto simp add: pos borel insert(2) setprod_ereal_pos positive_integral_positive) 

41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

901 
done 
47694  902 
qed (simp add: space_PiM) 
41096  903 

41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

904 
lemma (in product_sigma_finite) product_integral_singleton: 
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

905 
assumes f: "f \<in> borel_measurable (M i)" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

906 
shows "(\<integral>x. f (x i) \<partial>Pi\<^isub>M {i} M) = integral\<^isup>L (M i) f" 
41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

907 
proof  
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

908 
interpret I: finite_product_sigma_finite M "{i}" by default simp 
43920  909 
have *: "(\<lambda>x. ereal (f x)) \<in> borel_measurable (M i)" 
910 
"(\<lambda>x. ereal ( f x)) \<in> borel_measurable (M i)" 

41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

911 
using assms by auto 
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

912 
show ?thesis 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

913 
unfolding lebesgue_integral_def *[THEN product_positive_integral_singleton] .. 
41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

914 
qed 
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

915 

bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

916 
lemma (in product_sigma_finite) product_integral_fold: 
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

917 
assumes IJ[simp]: "I \<inter> J = {}" and fin: "finite I" "finite J" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

918 
and f: "integrable (Pi\<^isub>M (I \<union> J) M) f" 
49780  919 
shows "integral\<^isup>L (Pi\<^isub>M (I \<union> J) M) f = (\<integral>x. (\<integral>y. f (merge I J (x, y)) \<partial>Pi\<^isub>M J M) \<partial>Pi\<^isub>M I M)" 
41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

920 
proof  
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

921 
interpret I: finite_product_sigma_finite M I by default fact 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

922 
interpret J: finite_product_sigma_finite M J by default fact 
41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

923 
have "finite (I \<union> J)" using fin by auto 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

924 
interpret IJ: finite_product_sigma_finite M "I \<union> J" by default fact 
47694  925 
interpret P: pair_sigma_finite "Pi\<^isub>M I M" "Pi\<^isub>M J M" by default 
49780  926 
let ?M = "merge I J" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

927 
let ?f = "\<lambda>x. f (?M x)" 
47694  928 
from f have f_borel: "f \<in> borel_measurable (Pi\<^isub>M (I \<union> J) M)" 
929 
by auto 

49780  930 
have P_borel: "(\<lambda>x. f (merge I J x)) \<in> borel_measurable (Pi\<^isub>M I M \<Otimes>\<^isub>M Pi\<^isub>M J M)" 
49776  931 
using measurable_comp[OF measurable_merge f_borel] by (simp add: comp_def) 
47694  932 
have f_int: "integrable (Pi\<^isub>M I M \<Otimes>\<^isub>M Pi\<^isub>M J M) ?f" 
49776  933 
by (rule integrable_distr[OF measurable_merge]) (simp add: distr_merge[OF IJ fin] f) 
41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

934 
show ?thesis 
47694  935 
apply (subst distr_merge[symmetric, OF IJ fin]) 
49776  936 
apply (subst integral_distr[OF measurable_merge f_borel]) 
47694  937 
apply (subst P.integrable_fst_measurable(2)[symmetric, OF f_int]) 
938 
apply simp 

939 
done 

41026
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

940 
qed 
bea75746dc9d
folding on arbitrary Lebesgue integrable functions
hoelzl
parents:
41023
diff
changeset

941 

49776  942 
lemma (in product_sigma_finite) 
943 
assumes IJ: "I \<inter> J = {}" "finite I" "finite J" and A: "A \<in> sets (Pi\<^isub>M (I \<union> J) M)" 

944 
shows emeasure_fold_integral: 

49780  945 
"emeasure (Pi\<^isub>M (I \<union> J) M) A = (\<integral>\<^isup>+x. emeasure (Pi\<^isub>M J M) ((\<lambda>y. merge I J (x, y)) ` A \<inter> space (Pi\<^isub>M J M)) \<partial>Pi\<^isub>M I M)" (is ?I) 
49776  946 
and emeasure_fold_measurable: 
49780  947 
"(\<lambda>x. emeasure (Pi\<^isub>M J M) ((\<lambda>y. merge I J (x, y)) ` A \<inter> space (Pi\<^isub>M J M))) \<in> borel_measurable (Pi\<^isub>M I M)" (is ?B) 
49776  948 
proof  
949 
interpret I: finite_product_sigma_finite M I by default fact 

950 
interpret J: finite_product_sigma_finite M J by default fact 

951 
interpret IJ: pair_sigma_finite "Pi\<^isub>M I M" "Pi\<^isub>M J M" .. 

49780  952 
have merge: "merge I J ` A \<inter> space (Pi\<^isub>M I M \<Otimes>\<^isub>M Pi\<^isub>M J M) \<in> sets (Pi\<^isub>M I M \<Otimes>\<^isub>M Pi\<^isub>M J M)" 
49776  953 
by (intro measurable_sets[OF _ A] measurable_merge assms) 
954 

955 
show ?I 

956 
apply (subst distr_merge[symmetric, OF IJ]) 

957 
apply (subst emeasure_distr[OF measurable_merge A]) 

958 
apply (subst J.emeasure_pair_measure_alt[OF merge]) 

959 
apply (auto intro!: positive_integral_cong arg_cong2[where f=emeasure] simp: space_pair_measure) 

960 
done 

961 

962 
show ?B 

963 
using IJ.measurable_emeasure_Pair1[OF merge] 

964 
by (simp add: vimage_compose[symmetric] comp_def space_pair_measure cong: measurable_cong) 

965 
qed 

966 

41096  967 
lemma (in product_sigma_finite) product_integral_insert: 
47694  968 
assumes I: "finite I" "i \<notin> I" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

969 
and f: "integrable (Pi\<^isub>M (insert i I) M) f" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

970 
shows "integral\<^isup>L (Pi\<^isub>M (insert i I) M) f = (\<integral>x. (\<integral>y. f (x(i:=y)) \<partial>M i) \<partial>Pi\<^isub>M I M)" 
41096  971 
proof  
47694  972 
have "integral\<^isup>L (Pi\<^isub>M (insert i I) M) f = integral\<^isup>L (Pi\<^isub>M (I \<union> {i}) M) f" 
973 
by simp 

49780  974 
also have "\<dots> = (\<integral>x. (\<integral>y. f (merge I {i} (x,y)) \<partial>Pi\<^isub>M {i} M) \<partial>Pi\<^isub>M I M)" 
47694  975 
using f I by (intro product_integral_fold) auto 
976 
also have "\<dots> = (\<integral>x. (\<integral>y. f (x(i := y)) \<partial>M i) \<partial>Pi\<^isub>M I M)" 

977 
proof (rule integral_cong, subst product_integral_singleton[symmetric]) 

978 
fix x assume x: "x \<in> space (Pi\<^isub>M I M)" 

979 
have f_borel: "f \<in> borel_measurable (Pi\<^isub>M (insert i I) M)" 

980 
using f by auto 

981 
show "(\<lambda>y. f (x(i := y))) \<in> borel_measurable (M i)" 

982 
using measurable_comp[OF measurable_component_update f_borel, OF x `i \<notin> I`] 

983 
unfolding comp_def . 

49780  984 
from x I show "(\<integral> y. f (merge I {i} (x,y)) \<partial>Pi\<^isub>M {i} M) = (\<integral> xa. f (x(i := xa i)) \<partial>Pi\<^isub>M {i} M)" 
47694  985 
by (auto intro!: integral_cong arg_cong[where f=f] simp: merge_def space_PiM extensional_def) 
41096  986 
qed 
47694  987 
finally show ?thesis . 
41096  988 
qed 
989 

990 
lemma (in product_sigma_finite) product_integrable_setprod: 

991 
fixes f :: "'i \<Rightarrow> 'a \<Rightarrow> real" 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

992 
assumes [simp]: "finite I" and integrable: "\<And>i. i \<in> I \<Longrightarrow> integrable (M i) (f i)" 
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

993 
shows "integrable (Pi\<^isub>M I M) (\<lambda>x. (\<Prod>i\<in>I. f i (x i)))" (is "integrable _ ?f") 
41096  994 
proof  
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

995 
interpret finite_product_sigma_finite M I by default fact 
41096  996 
have f: "\<And>i. i \<in> I \<Longrightarrow> f i \<in> borel_measurable (M i)" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

997 
using integrable unfolding integrable_def by auto 
47694  998 
have borel: "?f \<in> borel_measurable (Pi\<^isub>M I M)" 
999 
using measurable_comp[OF measurable_component_singleton[of _ I M] f] by (auto simp: comp_def) 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

1000 
moreover have "integrable (Pi\<^isub>M I M) (\<lambda>x. \<bar>\<Prod>i\<in>I. f i (x i)\<bar>)" 
41096  1001 
proof (unfold integrable_def, intro conjI) 
47694  1002 
show "(\<lambda>x. abs (?f x)) \<in> borel_measurable (Pi\<^isub>M I M)" 
41096  1003 
using borel by auto 
47694  1004 
have "(\<integral>\<^isup>+x. ereal (abs (?f x)) \<partial>Pi\<^isub>M I M) = (\<integral>\<^isup>+x. (\<Prod>i\<in>I. ereal (abs (f i (x i)))) \<partial>Pi\<^isub>M I M)" 
43920  1005 
by (simp add: setprod_ereal abs_setprod) 
1006 
also have "\<dots> = (\<Prod>i\<in>I. (\<integral>\<^isup>+x. ereal (abs (f i x)) \<partial>M i))" 

41096  1007 
using f by (subst product_positive_integral_setprod) auto 
41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

1008 
also have "\<dots> < \<infinity>" 
47694  1009 
using integrable[THEN integrable_abs] 
1010 
by (simp add: setprod_PInf integrable_def positive_integral_positive) 

1011 
finally show "(\<integral>\<^isup>+x. ereal (abs (?f x)) \<partial>(Pi\<^isub>M I M)) \<noteq> \<infinity>" by auto 

1012 
have "(\<integral>\<^isup>+x. ereal ( abs (?f x)) \<partial>(Pi\<^isub>M I M)) = (\<integral>\<^isup>+x. 0 \<partial>(Pi\<^isub>M I M))" 

41981
cdf7693bbe08
reworked Probability theory: measures are not type restricted to positive extended reals
hoelzl
parents:
41831
diff
changeset

1013 
by (intro positive_integral_cong_pos) auto 
47694  1014 
then show "(\<integral>\<^isup>+x. ereal ( abs (?f x)) \<partial>(Pi\<^isub>M I M)) \<noteq> \<infinity>" by simp 
41096  1015 
qed 
1016 
ultimately show ?thesis 

1017 
by (rule integrable_abs_iff[THEN iffD1]) 

1018 
qed 

1019 

1020 
lemma (in product_sigma_finite) product_integral_setprod: 

1021 
fixes f :: "'i \<Rightarrow> 'a \<Rightarrow> real" 

49780  1022 
assumes "finite I" and integrable: "\<And>i. i \<in> I \<Longrightarrow> integrable (M i) (f i)" 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

1023 
shows "(\<integral>x. (\<Prod>i\<in>I. f i (x i)) \<partial>Pi\<^isub>M I M) = (\<Prod>i\<in>I. integral\<^isup>L (M i) (f i))" 
49780  1024 
using assms proof induct 
1025 
case empty 

1026 
interpret finite_measure "Pi\<^isub>M {} M" 

1027 
by rule (simp add: space_PiM) 

1028 
show ?case by (simp add: space_PiM measure_def) 

41096  1029 
next 
1030 
case (insert i I) 

1031 
then have iI: "finite (insert i I)" by auto 

1032 
then have prod: "\<And>J. J \<subseteq> insert i I \<Longrightarrow> 

41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

1033 
integrable (Pi\<^isub>M J M) (\<lambda>x. (\<Prod>i\<in>J. f i (x i)))" 
49780  1034 
by (intro product_integrable_setprod insert(4)) (auto intro: finite_subset) 
41689
3e39b0e730d6
the measure valuation is again part of the measure_space type, instead of an explicit parameter to the locale;
hoelzl
parents:
41661
diff
changeset

1035 
interpret I: finite_product_sigma_finite M I by default fact 
41096  1036 
have *: "\<And>x y. (\<Prod>j\<in>I. f j (if j = i then y else x j)) = (\<Prod>j\<in>I. f j (x j))" 
1037 
using `i \<notin> I` by (auto intro!: setprod_cong) 

1038 
show ?case 

49780  1039 
unfolding product_integral_insert[OF insert(1,2) prod[OF subset_refl]] 
47694  1040 
by (simp add: * insert integral_multc integral_cmult[OF prod] subset_insertI) 
41096  1041 
qed 
1042 

49776  1043 
lemma sets_Collect_single: 
1044 
"i \<in> I \<Longrightarrow> A \<in> sets (M i) \<Longrightarrow> { x \<in> space (Pi\<^isub>M I M). x i \<in> A } \<in> sets (Pi\<^isub>M I M)" 

50003  1045 
by simp 
49776  1046 

1047 
lemma sigma_prod_algebra_sigma_eq_infinite: 

1048 
fixes E :: "'i \<Rightarrow> 'a set set" 

49779
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1049 
assumes S_union: "\<And>i. i \<in> I \<Longrightarrow> (\<Union>j. S i j) = space (M i)" 
49776  1050 
and S_in_E: "\<And>i. i \<in> I \<Longrightarrow> range (S i) \<subseteq> E i" 
1051 
assumes E_closed: "\<And>i. i \<in> I \<Longrightarrow> E i \<subseteq> Pow (space (M i))" 

1052 
and E_generates: "\<And>i. i \<in> I \<Longrightarrow> sets (M i) = sigma_sets (space (M i)) (E i)" 

1053 
defines "P == {{f\<in>\<Pi>\<^isub>E i\<in>I. space (M i). f i \<in> A}  i A. i \<in> I \<and> A \<in> E i}" 

1054 
shows "sets (PiM I M) = sigma_sets (space (PiM I M)) P" 

1055 
proof 

1056 
let ?P = "sigma (space (Pi\<^isub>M I M)) P" 

1057 
have P_closed: "P \<subseteq> Pow (space (Pi\<^isub>M I M))" 

1058 
using E_closed by (auto simp: space_PiM P_def Pi_iff subset_eq) 

1059 
then have space_P: "space ?P = (\<Pi>\<^isub>E i\<in>I. space (M i))" 

1060 
by (simp add: space_PiM) 

1061 
have "sets (PiM I M) = 

1062 
sigma_sets (space ?P) {{f \<in> \<Pi>\<^isub>E i\<in>I. space (M i). f i \<in> A} i A. i \<in> I \<and> A \<in> sets (M i)}" 

1063 
using sets_PiM_single[of I M] by (simp add: space_P) 

1064 
also have "\<dots> \<subseteq> sets (sigma (space (PiM I M)) P)" 

1065 
proof (safe intro!: sigma_sets_subset) 

1066 
fix i A assume "i \<in> I" and A: "A \<in> sets (M i)" 

1067 
then have "(\<lambda>x. x i) \<in> measurable ?P (sigma (space (M i)) (E i))" 

1068 
apply (subst measurable_iff_measure_of) 

1069 
apply (simp_all add: P_closed) 

1070 
using E_closed 

1071 
apply (force simp: subset_eq space_PiM) 

1072 
apply (force simp: subset_eq space_PiM) 

1073 
apply (auto simp: P_def intro!: sigma_sets.Basic exI[of _ i]) 

1074 
apply (rule_tac x=Aa in exI) 

1075 
apply (auto simp: space_PiM) 

1076 
done 

1077 
from measurable_sets[OF this, of A] A `i \<in> I` E_closed 

1078 
have "(\<lambda>x. x i) ` A \<inter> space ?P \<in> sets ?P" 

1079 
by (simp add: E_generates) 

1080 
also have "(\<lambda>x. x i) ` A \<inter> space ?P = {f \<in> \<Pi>\<^isub>E i\<in>I. space (M i). f i \<in> A}" 

1081 
using P_closed by (auto simp: space_PiM) 

1082 
finally show "\<dots> \<in> sets ?P" . 

1083 
qed 

1084 
finally show "sets (PiM I M) \<subseteq> sigma_sets (space (PiM I M)) P" 

1085 
by (simp add: P_closed) 

1086 
show "sigma_sets (space (PiM I M)) P \<subseteq> sets (PiM I M)" 

1087 
unfolding P_def space_PiM[symmetric] 

1088 
by (intro sigma_sets_subset) (auto simp: E_generates sets_Collect_single) 

1089 
qed 

1090 

49779
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1091 
lemma bchoice_iff: "(\<forall>a\<in>A. \<exists>b. P a b) \<longleftrightarrow> (\<exists>f. \<forall>a\<in>A. P a (f a))" 
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1092 
by metis 
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1093 

47694  1094 
lemma sigma_prod_algebra_sigma_eq: 
49779
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1095 
fixes E :: "'i \<Rightarrow> 'a set set" and S :: "'i \<Rightarrow> nat \<Rightarrow> 'a set" 
47694  1096 
assumes "finite I" 
49779
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1097 
assumes S_union: "\<And>i. i \<in> I \<Longrightarrow> (\<Union>j. S i j) = space (M i)" 
47694  1098 
and S_in_E: "\<And>i. i \<in> I \<Longrightarrow> range (S i) \<subseteq> E i" 
1099 
assumes E_closed: "\<And>i. i \<in> I \<Longrightarrow> E i \<subseteq> Pow (space (M i))" 

1100 
and E_generates: "\<And>i. i \<in> I \<Longrightarrow> sets (M i) = sigma_sets (space (M i)) (E i)" 

1101 
defines "P == { Pi\<^isub>E I F  F. \<forall>i\<in>I. F i \<in> E i }" 

1102 
shows "sets (PiM I M) = sigma_sets (space (PiM I M)) P" 

1103 
proof 

1104 
let ?P = "sigma (space (Pi\<^isub>M I M)) P" 

49779
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1105 
from `finite I`[THEN ex_bij_betw_finite_nat] guess T .. 
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1106 
then have T: "\<And>i. i \<in> I \<Longrightarrow> T i < card I" "\<And>i. i\<in>I \<Longrightarrow> the_inv_into I T (T i) = i" 
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1107 
by (auto simp add: bij_betw_def set_eq_iff image_iff the_inv_into_f_f) 
47694  1108 
have P_closed: "P \<subseteq> Pow (space (Pi\<^isub>M I M))" 
1109 
using E_closed by (auto simp: space_PiM P_def Pi_iff subset_eq) 

1110 
then have space_P: "space ?P = (\<Pi>\<^isub>E i\<in>I. space (M i))" 

1111 
by (simp add: space_PiM) 

1112 
have "sets (PiM I M) = 

1113 
sigma_sets (space ?P) {{f \<in> \<Pi>\<^isub>E i\<in>I. space (M i). f i \<in> A} i A. i \<in> I \<and> A \<in> sets (M i)}" 

1114 
using sets_PiM_single[of I M] by (simp add: space_P) 

1115 
also have "\<dots> \<subseteq> sets (sigma (space (PiM I M)) P)" 

1116 
proof (safe intro!: sigma_sets_subset) 

1117 
fix i A assume "i \<in> I" and A: "A \<in> sets (M i)" 

1118 
have "(\<lambda>x. x i) \<in> measurable ?P (sigma (space (M i)) (E i))" 

1119 
proof (subst measurable_iff_measure_of) 

1120 
show "E i \<subseteq> Pow (space (M i))" using `i \<in> I` by fact 

1121 
from space_P `i \<in> I` show "(\<lambda>x. x i) \<in> space ?P \<rightarrow> space (M i)" 

1122 
by (auto simp: Pi_iff) 

1123 
show "\<forall>A\<in>E i. (\<lambda>x. x i) ` A \<inter> space ?P \<in> sets ?P" 

1124 
proof 

1125 
fix A assume A: "A \<in> E i" 

1126 
then have "(\<lambda>x. x i) ` A \<inter> space ?P = (\<Pi>\<^isub>E j\<in>I. if i = j then A else space (M j))" 

1127 
using E_closed `i \<in> I` by (auto simp: space_P Pi_iff subset_eq split: split_if_asm) 

1128 
also have "\<dots> = (\<Pi>\<^isub>E j\<in>I. \<Union>n. if i = j then A else S j n)" 

1129 
by (intro PiE_cong) (simp add: S_union) 

49779
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1130 
also have "\<dots> = (\<Union>xs\<in>{xs. length xs = card I}. \<Pi>\<^isub>E j\<in>I. if i = j then A else S j (xs ! T j))" 
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1131 
using T 
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1132 
apply (auto simp: Pi_iff bchoice_iff) 
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1133 
apply (rule_tac x="map (\<lambda>n. f (the_inv_into I T n)) [0..<card I]" in exI) 
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1134 
apply (auto simp: bij_betw_def) 
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1135 
done 
47694  1136 
also have "\<dots> \<in> sets ?P" 
1137 
proof (safe intro!: countable_UN) 

49779
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1138 
fix xs show "(\<Pi>\<^isub>E j\<in>I. if i = j then A else S j (xs ! T j)) \<in> sets ?P" 
47694  1139 
using A S_in_E 
1140 
by (simp add: P_closed) 

49779
1484b4b82855
remove incseq assumption from sigma_prod_algebra_sigma_eq
hoelzl
parents:
49776
diff
changeset

1141 
(auto simp: P_def subset_eq intro!: exI[of _ "\<lambda>j. if i = j then A else S j (xs ! T j)"]) 
47694  1142 
qed 
1143 
finally show "(\<lambda>x. x i) ` A \<inter> space ?P \<in> sets ?P" 

1144 
using P_closed by simp 

1145 
qed 

1146 
qed 

1147 
from measurable_sets[OF this, of A] A `i \<in> I` E_closed 

1148 
have "(\<lambda>x. x i) ` A \<inter> space ?P \<in> sets ?P" 

1149 
by (simp add: E_generates) 

1150 
also have "(\<lambda>x. x i) ` A \<inter> space ?P = {f \<in> \<Pi>\<^isub>E i\<in>I. space (M i). f i \<in> A}" 

1151 
using P_closed by (auto simp: space_PiM) 

1152 
finally show "\<dots> \<in> sets ?P" . 

1153 
qed 

1154 
finally show "sets (PiM I M) \<subseteq> sigma_sets (space (PiM I M)) P" 

1155 
by (simp add: P_closed) 

1156 
show "sigma_sets (space (PiM I M)) P \<subseteq> sets (PiM I M)" 

1157 
using `finite I` 

50003  1158 
by (auto intro!: sigma_sets_subset sets_PiM_I_finite simp: E_generates P_def) 
47694  1159 
qed 
1160 

1161 
end 