starfun, starset, and other functions on NS types are now polymorphic;
many similar theorems have been generalized and merged;
(star_n X) replaces (Abs_star(starrel `` {X}));
many proofs have been simplified with the transfer tactic.
(* Title : Star.thy
Author : Jacques D. Fleuriot
Copyright : 1998 University of Cambridge
Conversion to Isar and new proofs by Lawrence C Paulson, 2003/4
*)
header{*Star-Transforms in Non-Standard Analysis*}
theory Star
imports NSA
begin
(* nonstandard extension of sets *)
syntax starset :: "'a set => 'a star set" ("*s* _" [80] 80)
translations "starset" == "Iset_of"
syntax starfun :: "('a => 'b) => 'a star => 'b star" ("*f* _" [80] 80)
translations "starfun" == "Ifun_of"
constdefs
(* internal sets *)
starset_n :: "(nat => 'a set) => 'a star set" ("*sn* _" [80] 80)
"*sn* As == Iset (star_n As)"
InternalSets :: "'a star set set"
"InternalSets == {X. \<exists>As. X = *sn* As}"
(* nonstandard extension of function *)
is_starext :: "['a star => 'a star, 'a => 'a] => bool"
"is_starext F f == (\<forall>x y. \<exists>X \<in> Rep_star(x). \<exists>Y \<in> Rep_star(y).
((y = (F x)) = ({n. Y n = f(X n)} : FreeUltrafilterNat)))"
(* internal functions *)
starfun_n :: "(nat => ('a => 'b)) => 'a star => 'b star"
("*fn* _" [80] 80)
"*fn* F == Ifun (star_n F)"
InternalFuns :: "('a star => 'b star) set"
"InternalFuns == {X. \<exists>F. X = *fn* F}"
(*--------------------------------------------------------
Preamble - Pulling "EX" over "ALL"
---------------------------------------------------------*)
(* This proof does not need AC and was suggested by the
referee for the JCM Paper: let f(x) be least y such
that Q(x,y)
*)
lemma no_choice: "\<forall>x. \<exists>y. Q x y ==> \<exists>(f :: nat => nat). \<forall>x. Q x (f x)"
apply (rule_tac x = "%x. LEAST y. Q x y" in exI)
apply (blast intro: LeastI)
done
subsection{*Properties of the Star-transform Applied to Sets of Reals*}
lemma STAR_UNIV_set [simp]: "*s*(UNIV::'a set) = (UNIV::'a star set)"
by (transfer UNIV_def, rule refl)
lemma STAR_empty_set [simp]: "*s* {} = {}"
by (transfer empty_def, rule refl)
lemma STAR_Un: "*s* (A Un B) = *s* A Un *s* B"
by (transfer Un_def, rule refl)
lemma STAR_Int: "*s* (A Int B) = *s* A Int *s* B"
by (transfer Int_def, rule refl)
lemma STAR_Compl: "*s* -A = -( *s* A)"
by (transfer Compl_def, rule refl)
lemma STAR_mem_Compl: "!!x. x \<notin> *s* F ==> x : *s* (- F)"
by (transfer Compl_def, simp)
lemma STAR_diff: "*s* (A - B) = *s* A - *s* B"
by (transfer set_diff_def, rule refl)
lemma STAR_subset: "A <= B ==> *s* A <= *s* B"
by (transfer subset_def, simp)
lemma STAR_mem: "a \<in> A ==> star_of a : *s* A"
by transfer
lemma STAR_mem_iff: "(star_of x \<in> *s* A) = (x \<in> A)"
by (transfer, rule refl)
lemma STAR_star_of_image_subset: "star_of ` A <= *s* A"
by (auto simp add: STAR_mem)
lemma STAR_hypreal_of_real_Int: "*s* X Int Reals = hypreal_of_real ` X"
by (auto simp add: SReal_def STAR_mem_iff)
lemma lemma_not_hyprealA: "x \<notin> hypreal_of_real ` A ==> \<forall>y \<in> A. x \<noteq> hypreal_of_real y"
by auto
lemma lemma_Compl_eq: "- {n. X n = xa} = {n. X n \<noteq> xa}"
by auto
lemma STAR_real_seq_to_hypreal:
"\<forall>n. (X n) \<notin> M ==> star_n X \<notin> *s* M"
apply (unfold Iset_of_def star_of_def)
apply (simp add: Iset_star_n)
done
lemma STAR_insert [simp]: "*s* (insert x A) = insert (star_of x) ( *s* A)"
by (transfer insert_def Un_def, rule refl)
lemma STAR_singleton: "*s* {x} = {star_of x}"
by simp
lemma STAR_not_mem: "x \<notin> F ==> star_of x \<notin> *s* F"
by transfer
lemma STAR_subset_closed: "[| x : *s* A; A <= B |] ==> x : *s* B"
by (blast dest: STAR_subset)
text{*Nonstandard extension of a set (defined using a constant
sequence) as a special case of an internal set*}
lemma starset_n_starset: "\<forall>n. (As n = A) ==> *sn* As = *s* A"
apply (drule expand_fun_eq [THEN iffD2])
apply (simp add: starset_n_def Iset_of_def star_of_def)
done
(*----------------------------------------------------------------*)
(* Theorems about nonstandard extensions of functions *)
(*----------------------------------------------------------------*)
(*----------------------------------------------------------------*)
(* Nonstandard extension of a function (defined using a *)
(* constant sequence) as a special case of an internal function *)
(*----------------------------------------------------------------*)
lemma starfun_n_starfun: "\<forall>n. (F n = f) ==> *fn* F = *f* f"
apply (drule expand_fun_eq [THEN iffD2])
apply (simp add: starfun_n_def Ifun_of_def star_of_def)
done
(*
Prove that abs for hypreal is a nonstandard extension of abs for real w/o
use of congruence property (proved after this for general
nonstandard extensions of real valued functions).
Proof now Uses the ultrafilter tactic!
*)
lemma hrabs_is_starext_rabs: "is_starext abs abs"
apply (simp add: is_starext_def, safe)
apply (rule_tac z = x in eq_Abs_star)
apply (rule_tac z = y in eq_Abs_star, auto)
apply (rule bexI, rule_tac [2] lemma_starrel_refl)
apply (rule bexI, rule_tac [2] lemma_starrel_refl)
apply (fold star_n_def)
apply (unfold star_abs_def Ifun_of_def star_of_def)
apply (simp add: Ifun_star_n star_n_eq_iff)
done
lemma Rep_star_FreeUltrafilterNat:
"[| X \<in> Rep_star z; Y \<in> Rep_star z |]
==> {n. X n = Y n} : FreeUltrafilterNat"
apply (rule_tac z = z in eq_Abs_star)
apply (auto, ultra)
done
text{*Nonstandard extension of functions*}
lemma starfun:
"( *f* f) (star_n X) = star_n (%n. f (X n))"
by (simp add: Ifun_of_def star_of_def Ifun_star_n)
lemma starfun_if_eq:
"w \<noteq> hypreal_of_real x
==> ( *f* (\<lambda>z. if z = x then a else g z)) w = ( *f* g) w"
apply (cases w)
apply (simp add: star_of_def starfun star_n_eq_iff, ultra)
done
(*-------------------------------------------
multiplication: ( *f) x ( *g) = *(f x g)
------------------------------------------*)
lemma starfun_mult: "!!x. ( *f* f) x * ( *f* g) x = ( *f* (%x. f x * g x)) x"
by (transfer, rule refl)
declare starfun_mult [symmetric, simp]
(*---------------------------------------
addition: ( *f) + ( *g) = *(f + g)
---------------------------------------*)
lemma starfun_add: "!!x. ( *f* f) x + ( *f* g) x = ( *f* (%x. f x + g x)) x"
by (transfer, rule refl)
declare starfun_add [symmetric, simp]
(*--------------------------------------------
subtraction: ( *f) + -( *g) = *(f + -g)
-------------------------------------------*)
lemma starfun_minus: "!!x. - ( *f* f) x = ( *f* (%x. - f x)) x"
by (transfer, rule refl)
declare starfun_minus [symmetric, simp]
(*FIXME: delete*)
lemma starfun_add_minus: "!!x. ( *f* f) x + -( *f* g) x = ( *f* (%x. f x + -g x)) x"
by (transfer, rule refl)
declare starfun_add_minus [symmetric, simp]
lemma starfun_diff: "!!x. ( *f* f) x - ( *f* g) x = ( *f* (%x. f x - g x)) x"
by (transfer, rule refl)
declare starfun_diff [symmetric, simp]
(*--------------------------------------
composition: ( *f) o ( *g) = *(f o g)
---------------------------------------*)
lemma starfun_o2: "(%x. ( *f* f) (( *f* g) x)) = *f* (%x. f (g x))"
by (transfer, rule refl)
lemma starfun_o: "( *f* f) o ( *f* g) = ( *f* (f o g))"
by (transfer o_def, rule refl)
text{*NS extension of constant function*}
lemma starfun_const_fun [simp]: "!!x. ( *f* (%x. k)) x = star_of k"
by (transfer, rule refl)
text{*the NS extension of the identity function*}
lemma starfun_Id [simp]: "!!x. ( *f* (%x. x)) x = x"
by (transfer, rule refl)
(* this is trivial, given starfun_Id *)
lemma starfun_Idfun_approx:
"x @= hypreal_of_real a ==> ( *f* (%x. x)) x @= hypreal_of_real a"
by (simp only: starfun_Id)
text{*The Star-function is a (nonstandard) extension of the function*}
lemma is_starext_starfun: "is_starext ( *f* f) f"
apply (simp add: is_starext_def, auto)
apply (rule_tac x = x in star_cases)
apply (rule_tac x = y in star_cases)
apply (auto intro!: bexI [OF _ Rep_star_star_n]
simp add: starfun star_n_eq_iff)
done
text{*Any nonstandard extension is in fact the Star-function*}
lemma is_starfun_starext: "is_starext F f ==> F = *f* f"
apply (simp add: is_starext_def)
apply (rule ext)
apply (rule_tac x = x in star_cases)
apply (drule_tac x = x in spec)
apply (drule_tac x = "( *f* f) x" in spec)
apply (auto dest!: FreeUltrafilterNat_Compl_mem simp add: starfun, ultra)
done
lemma is_starext_starfun_iff: "(is_starext F f) = (F = *f* f)"
by (blast intro: is_starfun_starext is_starext_starfun)
text{*extented function has same solution as its standard
version for real arguments. i.e they are the same
for all real arguments*}
lemma starfun_eq [simp]: "( *f* f) (star_of a) = star_of (f a)"
by (transfer, rule refl)
lemma starfun_approx: "( *f* f) (star_of a) @= hypreal_of_real (f a)"
by simp
(* useful for NS definition of derivatives *)
lemma starfun_lambda_cancel:
"!!x'. ( *f* (%h. f (x + h))) x' = ( *f* f) (star_of x + x')"
by (transfer, rule refl)
lemma starfun_lambda_cancel2:
"( *f* (%h. f(g(x + h)))) x' = ( *f* (f o g)) (star_of x + x')"
by (unfold o_def, rule starfun_lambda_cancel)
lemma starfun_mult_HFinite_approx: "[| ( *f* f) x @= l; ( *f* g) x @= m;
l: HFinite; m: HFinite
|] ==> ( *f* (%x. f x * g x)) x @= l * m"
apply (drule (3) approx_mult_HFinite)
apply (auto intro: approx_HFinite [OF _ approx_sym])
done
lemma starfun_add_approx: "[| ( *f* f) x @= l; ( *f* g) x @= m
|] ==> ( *f* (%x. f x + g x)) x @= l + m"
by (auto intro: approx_add)
text{*Examples: hrabs is nonstandard extension of rabs
inverse is nonstandard extension of inverse*}
(* can be proved easily using theorem "starfun" and *)
(* properties of ultrafilter as for inverse below we *)
(* use the theorem we proved above instead *)
lemma starfun_rabs_hrabs: "*f* abs = abs"
by (rule hrabs_is_starext_rabs [THEN is_starext_starfun_iff [THEN iffD1], symmetric])
lemma starfun_inverse_inverse [simp]: "( *f* inverse) x = inverse(x)"
by (unfold star_inverse_def, rule refl)
lemma starfun_inverse: "!!x. inverse (( *f* f) x) = ( *f* (%x. inverse (f x))) x"
by (transfer, rule refl)
declare starfun_inverse [symmetric, simp]
lemma starfun_divide: "!!x. ( *f* f) x / ( *f* g) x = ( *f* (%x. f x / g x)) x"
by (transfer, rule refl)
declare starfun_divide [symmetric, simp]
lemma starfun_inverse2: "!!x. inverse (( *f* f) x) = ( *f* (%x. inverse (f x))) x"
by (transfer, rule refl)
text{*General lemma/theorem needed for proofs in elementary
topology of the reals*}
lemma starfun_mem_starset:
"!!x. ( *f* f) x : *s* A ==> x : *s* {x. f x \<in> A}"
by (transfer, simp)
text{*Alternative definition for hrabs with rabs function
applied entrywise to equivalence class representative.
This is easily proved using starfun and ns extension thm*}
lemma hypreal_hrabs:
"abs (star_n X) = star_n (%n. abs (X n))"
by (simp only: starfun_rabs_hrabs [symmetric] starfun)
text{*nonstandard extension of set through nonstandard extension
of rabs function i.e hrabs. A more general result should be
where we replace rabs by some arbitrary function f and hrabs
by its NS extenson. See second NS set extension below.*}
lemma STAR_rabs_add_minus:
"*s* {x. abs (x + - y) < r} =
{x. abs(x + -hypreal_of_real y) < hypreal_of_real r}"
by (transfer, rule refl)
lemma STAR_starfun_rabs_add_minus:
"*s* {x. abs (f x + - y) < r} =
{x. abs(( *f* f) x + -hypreal_of_real y) < hypreal_of_real r}"
by (transfer, rule refl)
text{*Another characterization of Infinitesimal and one of @= relation.
In this theory since @{text hypreal_hrabs} proved here. Maybe
move both theorems??*}
lemma Infinitesimal_FreeUltrafilterNat_iff2:
"(x \<in> Infinitesimal) =
(\<exists>X \<in> Rep_star(x).
\<forall>m. {n. abs(X n) < inverse(real(Suc m))}
\<in> FreeUltrafilterNat)"
apply (cases x)
apply (auto intro!: bexI lemma_starrel_refl
simp add: Infinitesimal_hypreal_of_nat_iff star_of_def
star_n_inverse star_n_abs star_n_less hypreal_of_nat_eq)
apply (drule_tac x = n in spec, ultra)
done
lemma approx_FreeUltrafilterNat_iff: "star_n X @= star_n Y =
(\<forall>m. {n. abs (X n + - Y n) <
inverse(real(Suc m))} : FreeUltrafilterNat)"
apply (subst approx_minus_iff)
apply (rule mem_infmal_iff [THEN subst])
apply (auto simp add: star_n_minus star_n_add Infinitesimal_FreeUltrafilterNat_iff2)
apply (drule_tac x = m in spec, ultra)
done
lemma inj_starfun: "inj starfun"
apply (rule inj_onI)
apply (rule ext, rule ccontr)
apply (drule_tac x = "star_n (%n. xa)" in fun_cong)
apply (auto simp add: starfun star_n_eq_iff)
done
ML
{*
val starset_n_def = thm"starset_n_def";
val InternalSets_def = thm"InternalSets_def";
val is_starext_def = thm"is_starext_def";
val starfun_n_def = thm"starfun_n_def";
val InternalFuns_def = thm"InternalFuns_def";
val no_choice = thm "no_choice";
val STAR_UNIV_set = thm "STAR_UNIV_set";
val STAR_empty_set = thm "STAR_empty_set";
val STAR_Un = thm "STAR_Un";
val STAR_Int = thm "STAR_Int";
val STAR_Compl = thm "STAR_Compl";
val STAR_mem_Compl = thm "STAR_mem_Compl";
val STAR_diff = thm "STAR_diff";
val STAR_subset = thm "STAR_subset";
val STAR_mem = thm "STAR_mem";
val STAR_star_of_image_subset = thm "STAR_star_of_image_subset";
val STAR_hypreal_of_real_Int = thm "STAR_hypreal_of_real_Int";
val STAR_real_seq_to_hypreal = thm "STAR_real_seq_to_hypreal";
val STAR_singleton = thm "STAR_singleton";
val STAR_not_mem = thm "STAR_not_mem";
val STAR_subset_closed = thm "STAR_subset_closed";
val starset_n_starset = thm "starset_n_starset";
val starfun_n_starfun = thm "starfun_n_starfun";
val hrabs_is_starext_rabs = thm "hrabs_is_starext_rabs";
val Rep_star_FreeUltrafilterNat = thm "Rep_star_FreeUltrafilterNat";
val starfun = thm "starfun";
val starfun_mult = thm "starfun_mult";
val starfun_add = thm "starfun_add";
val starfun_minus = thm "starfun_minus";
val starfun_add_minus = thm "starfun_add_minus";
val starfun_diff = thm "starfun_diff";
val starfun_o2 = thm "starfun_o2";
val starfun_o = thm "starfun_o";
val starfun_const_fun = thm "starfun_const_fun";
val starfun_Idfun_approx = thm "starfun_Idfun_approx";
val starfun_Id = thm "starfun_Id";
val is_starext_starfun = thm "is_starext_starfun";
val is_starfun_starext = thm "is_starfun_starext";
val is_starext_starfun_iff = thm "is_starext_starfun_iff";
val starfun_eq = thm "starfun_eq";
val starfun_approx = thm "starfun_approx";
val starfun_lambda_cancel = thm "starfun_lambda_cancel";
val starfun_lambda_cancel2 = thm "starfun_lambda_cancel2";
val starfun_mult_HFinite_approx = thm "starfun_mult_HFinite_approx";
val starfun_add_approx = thm "starfun_add_approx";
val starfun_rabs_hrabs = thm "starfun_rabs_hrabs";
val starfun_inverse_inverse = thm "starfun_inverse_inverse";
val starfun_inverse = thm "starfun_inverse";
val starfun_divide = thm "starfun_divide";
val starfun_inverse2 = thm "starfun_inverse2";
val starfun_mem_starset = thm "starfun_mem_starset";
val hypreal_hrabs = thm "hypreal_hrabs";
val STAR_rabs_add_minus = thm "STAR_rabs_add_minus";
val STAR_starfun_rabs_add_minus = thm "STAR_starfun_rabs_add_minus";
val Infinitesimal_FreeUltrafilterNat_iff2 = thm "Infinitesimal_FreeUltrafilterNat_iff2";
val approx_FreeUltrafilterNat_iff = thm "approx_FreeUltrafilterNat_iff";
val inj_starfun = thm "inj_starfun";
*}
end