/**/
function alert4() //9812060000
{
document.write(''
+'Output may contain error, Please verify first. '
+'Program environment is MSIE 6.0, please use MSIE '
+''
);
} //function alert4() 9812060002
/**
Cauchy Inequality used frequently
create a command to build it any
where.
2009-06-26-10-52 here
2009-06-26-11-20 must call HelloCauchy(hcPar)
outside of
..
otherwise width
wrong.
9903101457 add w1 to HelloCauchy()
in holder02.htm
99,02,07,21,45 the use of w1
from 9902072131
to 9902072143
solve infinite loop problem
and no need change from uppercase
to lowercase.
/**/
function HelloCauchy(hcPar1) //9806261053
{
strCauchy1=''
+'<a name="Cauchy'
+hcPar1
+'"> '
+'Cauchy's Inequality'
+'
'
+'
'
+'
k=∞
\n
∑
\n
k=1
'
+'
a k b k
'
+'
≦
'
+'
[
'
+'
k=∞
\n
∑
\n
k=1
'
+'
a k2
'
+'
]
'
+'
1/2
\n
\n
'
+'
[
'
+'
k=∞
\n
∑
\n
k=1
'
+'
b k2
'
+'
]
'
+'
1/2
\n
\n
'
+'
'
+'
---Page 5 ---eqn.1.7 '
+'width of above equation'
+''
+' calling ID '+hcPar1
;
var iter0=0; //9902161139
var iterMax=5; //[=][][]
if(arguments.length==1)
document.write(strCauchy1);
else
{
argLen=arguments.length; //9806261457
var w0; //9806261502
var w2; //9806261336
var w1; //9902072131
var strCauchy2=strCauchy1;
for(w0=1;w0=0)
{ //9902161158 add this '{', '}'
strCauchy2= //9806261337
strCauchy2.replace(arguments[w0][0],arguments[w0][1]);
w1=arguments[w0][1].length; //9902072133
if(w1==0)w1=1; //9902072143
//iter0++;
} //9902161158 add this '{', '}'
} //if(arguments[w0].length==1) else
} //for(w0=0;w0.. otherwise width
wrong.
one test call code is next, but delete four '#'
99,02,07,21,45 the use of w1
from 9902072131
to 9902072143
solve infinite loop problem
and no need change from uppercase
to lowercase.
/**/
function HelloHilbert(hhPar1) //9904041826
{
var strHilbert1=''
+'<a name="Hilbert'
+hhPar1
+'"> '
+'Hilbert's Inequality'
+'
'
+'
'
+'
m=∞
∑
m=1
'
+'
n=∞
∑
n=1
'
+'
ambn
m+n
'
+'
<C
(
'
+'
m=∞
∑
m=1
'
+'
am2
'
+'
)
'
+'
1/2
'
+'
(
'
+'
n=∞
∑
n=1
'
+'
bn2
'
+'
)
'
+'
1/2
'
+'
'
+'
---Page 155 ---line 16 ---eqn.10.1 '
+'width of above equation'
+''
+' calling ID '+hhPar1
;
var iter0=0; //9902161139
var iterMax=5; //[=][][]
if(arguments.length==1)
document.write(strHilbert1);
else
{
argLen=arguments.length; //9806261457
var w0; //9806261502
var w2; //9806261336
var w1; //9902072131
var strHilbert2=strHilbert1;
for(w0=1;w0=0)
{ //9902161158 add this '{', '}'
strHilbert2= //9806261337
strHilbert2.replace(arguments[w0][0],arguments[w0][1]);
w1=arguments[w0][1].length; //9902072133
if(w1==0)w1=1; //9902072143
//iter0++;
} //9902161158 add this '{', '}'
} //if(arguments[w0].length==1) else
} //for(w0=0;w0
Inequality
Study 40th file
Update 2010-04-26
indexthis
program
DocA
Limit
XYGraph v2.3 - web page graph
☜☞
donate
get code
The Cauchy-Schwarz Master Class
J. Michael
Steele
★★★★★
This file is personal home work. No one
proofread. Cannot promise correctness.
If you suspect any view point wrong,
please ask a math expert near by.
Freeman 2009-06-19-10-46
Please use MSIE browser to read this file.
Did not test other browser. This file is
written under MSIE 6.0
<a name="docA001">Index beginIndex this file
2009-06-08-19-10 start
This file http://freeman2.com/tute0008.htm
is Freeman's reading notes. Although Freeman
always keep correct view point. But Freeman's
capability is limited, plus no one proofread
this file. Then you can still find wrong view
point. When read, please put question mark as
often as possible. If you suspect any view
point wrong, please ask a math expert near by.
<a name="docA002">,<a name="textbook">
This file is a note for reading inequality
book written by ProfessorJ. MichaelSteeleThe Cauchy-Schwarz Master Class★★★★★
Below use 'textbook' as abbreviation.
Freeman also read web pages online, and will
indicate the source URL at discussion point.
<a name="docA003">
Freeman study mechanical engineering.
Engineering mathematics do not teach
inequality. Above book is first inequality
book. First time read, it was very hard.
Although high school time learned
a*a + b*b >= 2*a*b
But this little knowledge do not help.
<a name="docA004">
This file follow textbook chapter section
order, but not continuous, Freeman skip
those uncertain sections/problems.
This file first function is to learn
inequality. Second function is to learn
how to use html code to write math
equations.
<a name="docA005">Index beginIndex this file
On 2009-01-27-10-08 Freeman accessed
the next page
http://www.sftw.umac.mo/~fstitl/10mmo/inequality.html
save as sftw.umac.mo_text_math_eqn_good.htm
Above page is the main reference for html
math equation.
In order to let reader to build html math
equation, previous file tute0007.htm page
end has math symbol and internal code.
This file third function is to display how
to draw curves in web page. The main engine
is XYGraph v2.3 - Technical Figures
Thank you for read Freeman's inequality page.
2009-06-08-19-47 stop
<a name="docA006">
2009-08-23-15-00 start
□ Exercise 1.14 solution
Exercise 1.14 problem (a) is solved by hint
Exercise 1.14 problem (b) is out of LiuHH's
reach. What is "cardinality of a set B⊂Z3" ?
Without clear understand the definition,
LiuHH is unable to solve problem (b)
LiuHH major in mechanical engineering.
LiuHH is mathematics admirer and outsider.
<a name="docA007">
To solve problems in
The Cauchy-Schwarz Master Class
LiuHH's goal is to solve 50% of the exercise
problems. In reality, solve 40% is good enough.
Please do not be surprise that LiuHH skip some
text section, skip some problems.
2009-08-23-15-12 stop
<a name="ch11a001">Index beginIndex this file
2010-04-21-16-27 start
■■Chapter 11: Hardy's Inequality
and the Flop
Textbook page 166
Hardy's Inequality work with
integration of a function f(x).
Chapter title say "Hardy's
Inequality and the Flop".
<a name="ch11a002">
What is Flop? //here do flop
Textbook page 176 top say:
"If we have a factor on both
side of an equation and if it
appears to a smaller power on
the 'right' than on the 'left'
then we can clear the factor
to the left to obtain a new--
and potentially useful--bound."
<a name="ch11a003">
Textbook target eqn.11.1 has
two terms [∫f(u)du]/x and f(x)
They are (finally) separated.
During proving process, the
step of integration by parts
generate [∫f(u)du]/x * f(x)
two terms are multiplied. Until
eqn.11.7 two are separatable.
This separation process is
called "Flop".
<a name="ch11a004">Index beginIndex this file
2010-04-21-16-47 here
■ Problem 11.1
(Hardy's Inequality)
Show that every integrable
function f:(0,T)→Real
satisfies the inequality
---page 166
---line 18
---eqn.11.1
width of above equation
<a name="ch11a006">
and show, moreover, that the
constant 4 cannot be replaced
with any smaller value.
2010-04-21-17-00 stop
<a name="ch11a007">
2010-04-21-17-01 start
Following is a reading notes
parallel to textbook page 166
to 168 with some change.
Online paper use Hardy's
Inequality which is slightly
generalized as following
<a name="ch11a009">Index beginIndex this file
2010-04-21-17-27 here
If p=2, eqn.BH001 become eqn.11.1
eqn.11.1 use p=2 (squared), then
if f(x)≦0, f(x)_squared≧0.
eqn.11.1 do not use absolute
sign for f(x).
<a name="ch11a010">
If choose eqn.BH001 to start,
if p=3, negative^3 is negative.
Some paper require that f(x)≧0
Some paper use absolute sign
for |f(x)|. To simplify equation
the following use eqn.BH001 and
require f(x)≧0.<a name="ch11a011">
The term {∫[u=0,x]f(u)du}/x is
an integration average of f(u)
over u=0 to u=x.
If f(u) is force, and if x is
distance, variable u upper
bound is u=x. Then u and x
have same physics dimension
distance, du/x is pure number.
{f(u)du}/x has same dimension
as f(u) force.
{∫[u=0,x]f(u)du}/x is average
force between distance=0 to x.
<a name="ch11a012">
eqn.11.1 tell us that the
integral of average_squared
is not too far from
integral of f(x)_squared
<a name="ch11a013">
eqn.11.1 greater than side has
coefficient 4. Whether this
number 4 is critical? The
following is stress testing
(textbook page 167 line 3)
<a name="ch11a014">Index beginIndex this file
Textbook set
f(x)=xα ---eqn.BH002
to test constant C. To set f(x)
to specific function xα, it is
for stress testing only. Target
at constant C=4. To prove eqn.11.1
f(x) is general. Start at here.
For stress testing, we have next
equation.
---page 167
---line 7
---eqn.11.2
width
<a name="ch11a016">
Now replace f(x) with xα get
x=T
∫
x=0
{
1
x
u=x
∫
u=0
uα
du
}
2
dx
≦ C
x=T
∫
x=0
x2α
dx
---page 167
---line 7
---eqn.BH003
width of above equation
<a name="ch11a017">
u integral is
∫[u=0,x]uαdu
=∫[u=0,x]d(uα+1)/(α+1)
=(xα+1-0α+1)/(α+1)
=(xα+1)/(α+1) ---eqn.BH004
//Power integration is simple
Integration average become
{∫[u=0,x]uαdu}/x
=(xα+1)/(α+1)/x
=(xα)/(α+1) ---eqn.BH005
<a name="ch11a018">
eqn.11.4 is {term} in eqn.11.2
{term} shall be squared and
integrated from x=0 to x=T.
The result is
∫[x=0,T]dx*[(xα)/(α+1)]2
={∫[x=0,T]x2αdx}/[(α+1)2]
={∫[x=0,T]d(x2α+1)/(2α+1)}/[(α+1)2]
={T2α+1}/(2α+1)/[(α+1)2] ---eqn.BH006
<a name="ch11a019">Index beginIndex this file
Similarly eqn.11.2 right side
integration has the value
∫[x=0,T]x2αdx
=∫[x=0,T]d(x2α+1)/(2α+1)
=(T2α+1-02α+1)/(2α+1)
=(T2α+1)/(2α+1) ---eqn.BH007
<a name="ch11a020">
Now put eqn.BH006 and eqn.BH007
back to eqn.11.2 get
T2α+1
(2α+1)(α+1)2
≦ C
T2α+1
(2α+1)
---page 167
---line 10
---eqn.BH008
width of above equation
<a name="ch11a021">
2010-04-21-18-29 here
Why it is called "stress testing"?
Because eqn.BH008 both side have
"(2α+1)" when it approach to zero
equation is stressed (divide by 0)
If α=-1/2, we can cancel non-zero
term T2α+1 in eqn.BH008. But we can
not cancel zero term (2α+1). We
require 2α+1>0, that is require
α>-1/2 ---eqn.BH009
<a name="ch11a022">
In the limit case α→-1/2
eqn.BH008 give us
lim[α→-1/2]C→1/[(-1/2)+1]2
(non-zero 2α+1 drop out)
[(-1/2)+1]2 has value 4
We get C=4 as the limit value.
<a name="ch11a023">
"limit value" means that C=4 is
critical or is best. C a little
bit higher than 4, α>-1/2 OK
If C=4, α = -1/2, then violate
the rule α>-1/2
2010-04-21-18-36 here
Above is stress testing for
eqn.11.1
<a name="ch11a024">Index beginIndex this file
Parallel to eqn.11.2, now see
eqn.BH001. (in eqn.BH001
set p=2 get eqn.11.2)
In eqn.BH001, replace red
(p/(p-1))p with C
and replace f(x) with xα
calculate as following.
---eqn.BH010
require p>1, f(x)≧0
width of above equation
<a name="ch11a026">
Less than side {term} can be
copied from eqn.BH005.
Instead of square eqn.BH005,
now take p-th power to
eqn.BH005 and integrate
from x=0 to x=T.
<a name="ch11a027">
The result is
∫[x=0,T]dx*[(xα)/(α+1)]p
={∫[x=0,T]xpαdx}/[(α+1)p]
={∫[x=0,T]d(xpα+1)/(pα+1)}/[(α+1)p]
={Tpα+1}/(pα+1)/[(α+1)p] ---eqn.BH011
<a name="ch11a028">
Similarly eqn.BH010 right side
integration has the value
∫[x=0,T]xpαdx
=∫[x=0,T]d(xpα+1)/(pα+1)
=(Tpα+1-0pα+1)/(pα+1)
=(Tpα+1)/(pα+1) ---eqn.BH012
<a name="ch11a029">Index beginIndex this file
Now put eqn.BH011 and eqn.BH012
back to eqn.BH010 get
Tpα+1
(pα+1)(α+1)p
≦ C
Tpα+1
(pα+1)
---eqn.BH013
width of above equation
<a name="ch11a030">
2010-04-21-19-29 here
If α=-1/p, we can cancel non-zero
term Tpα+1 in eqn.BH013.
But we can not cancel zero
term (pα+1). We require pα+1>0
that is require
α>-1/p ---eqn.BH014
<a name="ch11a031">
In the limit case α→-1/p
eqn.BH013 give us
lim[α→-1/p]C→1/[(-1/p)+1]p
(non-zero pα+1 drop out)
[(-1/p)+1]p has value (p/(p-1))p
We get C=(p/(p-1))p as the
limit value.
If p=2, (p/(p-1))p=4.
<a name="ch11a032">
Square in eqn.11.1 and
p-th power in eqn.BH001
are parallel
Coefficient C in eqn.11.1 and
coef. (p/(p-1))p in eqn.BH001
are parallel.
Parallel analysis is success.
<a name="ch11a033">
Above is special case
f(x)=xα ---eqn.BH002
for stress testing.
The following is general case
f(x) be arbitrary.
2010-04-21-19-48 stop
<a name="ch11a034">Index beginIndex this file
2010-04-21-20-42 start
LiuHH try to prove eqn.BH001
instead of eqn.11.1.
In eqn.BH001, if set p=2,
it is same as eqn.11.1.
Prove method still parallel
with textbook.
<a name="ch11a035">
We start from eqn.BH001 less
than side. Textbook suggest
use integration by parts.
//Previous notes: ch07a040
Calculus product rule tell us
that
d[u(x)*v(x)]/dx
= u(x)*d[v(x)]/dx
+ v(x)*d[u(x)]/dx ---eqn.BH015
<a name="ch11a036">
Let u'(x)=d[u(x)]/dx ---eqn.BH016
let v'(x)=d[v(x)]/dx ---eqn.BH017
then
d[u(x)*v(x)] ---eqn.BH018
=[u(x)*v'(x)+u'(x)*v(x)]dx
<a name="ch11a037">
Integration get
∫d[u(x)*v(x)]
=∫u(x)*v'(x)dx
+∫u'(x)*v(x)dx ---eqn.BH019
or
<a name="ch11a038">
∫u(x)*v'(x)dx
=-∫u'(x)*v(x)dx
+∫d[u(x)*v(x)] ---eqn.BH020
eqn.BH020 is common used form.
∫d[u(x)*v(x)] is boundary term
<a name="ch11a039">Index beginIndex this file
re-write eqn.BH001 less than
side as next (not in textbook)
here create coefficient '1/(-p+1)'; there create coefficient 'p'
x=T
∫
x=0
{
u=x
∫
u=0
f(u)du
}
p
dx
xp
=
x=T
∫
x=0
{
u=x
∫
u=0
f(u)du
}
p
d(x-p+1)
-p+1
---eqn.BH021
require p>1, f(x)≧0
red is u(x); blue is v'(x); purple is d[v(x)]
width
<a name="ch11a040">
The next equality sign is integration by parts
=-
x=T
∫
x=0
x-p+1
-p+1
d [
{
u=x
∫
u=0
f(u)du
}
p
]
+
x=T
|
x=0
{
u=x
∫
u=0
f(u)du
}
p
x-p+1
-p+1
---eqn.BH022
red bar '|' and right side is boundary term
width
<a name="insert41">
2010-04-22-10-08 insert start
Although the proof is for eqn.BH001, still record
textbook equation by set p=2, make sure compatibility.
In eqn.BH021 left side set p=2, get next eqn.BH023
<a name="insert42">
I
define
=
x=T
∫
x=0
{
u=x
∫
u=0
f(u)du
}
2
dx
x2
=
-
x=T
∫
x=0
{
u=x
∫
u=0
f(u)du
}
2
(
1
x
)
'
dx
---page 167
---line 20
---eqn.BH023
require p>1, f(x)≧0; red is u(x); blue is v'(x)dx
width
<a name="insert43">Index beginIndex this file Integration by parts process drop prime (differentiation)
from above blue term, and differentiate to above red term.
Result a factor '2', red term reduce power by one and factor
out a 'f(x)' appears in next line.
Below red from above red. Below blue from above blue.
<a name="insert44">
Please compare eqn.BH023 with eqn.BH020.
Begin integration by parts eqn.BH020 has positive sign.
But eqn.BH023 second line has a negative sign and begin.
When change from v'(x)dx to u'(x)dx, a negative sign
appears, which make red '2' (below) positive and make
boundary term (purple) negative.
<a name="insert45">f(x) from Leibniz's rule
I=
2
x=T
∫
x=0
{
u=x
∫
u=0
f(u)du
}
f(x)
1
x
dx
-
x=T
|
x=0
{
u=x
∫
u=0
f(u)du
}
2
1
x
---page 167
---line 24
---eqn.11.3
purple bar '|' and right side is boundary term (2nd line)
width
2010-04-22-11-16 insert stop
<a name="ch11a043">
2010-04-21-22-03 here
Upper side x=T has no problem.
Lower side x=0 is 0/0 which is
undetermined.
We need step back and see what
is the limit value of this 0/0
<a name="ch11a045">
We need L'Hopital's Rule
If g(a)=h(a)=0 and g'(a), h'(a)
exist, assume h'(a)≠0, then
lim[x→a]g(x)/h(x)
=g'(x)/h'(x) ---eqn.BH027
<a name="ch11a046">Index beginIndex this file
Now in eqn.BH026
g(x) is {∫[u=0,x]f(u)du}p
h(x) is xp-1
point a is a=0
g(0)=0 because integration domain
∫[u=0,0] is zero.
h(0)=0 because p-1>0, then 0p-1=0
<a name="ch11a047">
L'Hopital's Rule tell us to
differentiate both g(x) and h(x)
and evaluate g'(a)/h'(a)
( here []' represent multiple
differentiation)
until not see 0/0. The last
g'(a)/h'(a) is answer.
<a name="ch11a048">
Look at eqn.BH026,
{∫[u=0,x]f(u)du}p has power p,
differentiate p times to constant
xp-1 has power p-1,
differentiate p-1 times to
constant.
<a name="ch11a049">
The competition result is that
after differentiate p-1 times
{∫[u=0,x]f(u)du}p still has
power one, ∫[u=0,0] is zero.
xp-1 become constant (p-1)!
We get 0/(p-1)! which is zero.
<a name="ch11a050">
Lower boundary term eqn.BH025 is
zero !!
Upper boundary term eqn.BH024≦0
Because p>1, then -p+1<0.
eqn.BH024 integration≧0
<a name="ch11a051">Index beginIndex this file
eqn.BH022 is equality with a
negative boundary term.
What to do?
This is inequality class.
Drop negative boundary term.
eqn.BH022 change from equality
to inequality
<a name="ch11a052">
eqn.BH022 remaining non-boundary
term will be eqn.BH001 greater
than side.
But still need flop operation.
2010-04-21-22-46 stop
<a name="ch11a053">
2010-04-22-11-36 start
We start at eqn.BH021=eqn.BH022
Integration by parts tell us that
boundary terms are non-positive,
we dropped boundary term from
eqn.BH022.
then eqn.BH021=eqn.BH022
become eqn.BH021≦eqn.BH022<a name="ch11a054">
The whole inequality is next
x=T
∫
x=0
{
u=x
∫
u=0
f(u)du
}
p
dx
xp
≦━
x=T
∫
x=0
x-p+1
-p+1
d
{
u=x
∫
u=0
f(u)du
}
p
---eqn.BH028
require f(x)≧0 and p>1, two red become positive.
Blue term participate 'flop' operation.
width
2010-04-22-11-51 here
<a name="ch11a055">
2010-04-22-12-16 start
eqn.BH028 greater than side
has a differentiation to an
integration. Variable x located
at integration upper bound. We
need Leibniz's rule
eqn.BH028 become next
---eqn.BH029
Red term come from Leibniz's rule.
width
<a name="ch11a057">
red dx/dx come from Leibniz's rule, black dx come from
integration by parts. Simplify eqn.BH029 as following.
eqn.BH028
left side
≦
x=T
∫
x=0
p
p-1
{
1
x
u=x
∫
u=0
f(u)du
}
p-1
f(x)
dx
---eqn.BH030
width
<a name="ch11a058">
In eqn.BH030 set p=2, get textbook equation (11.4)
x=T
∫
x=0
{
1
x
u=x
∫
u=0
f(u)du
}
2
dx
≦2
x=T
∫
x=0
{
1
x
u=x
∫
u=0
f(u)du
}
f(x)
dx
---eqn.11.4
width of above equation
<a name="ch11a059">
2010-04-22-12-58 here
Let us pay attention to the
blue term {∫[u=0,x]f(u)du}/x
and red term f(x)
Blue term show up at both less
than side and greater than side.
<a name="ch11a060">
Red term show up at greater
than side only. Can we make
blue at one side and red at
the other side of inequality?
Flop operation is to separate
blue from red.
2010-04-22-13-15 stop
<a name="ch11a061">Index beginIndex this file
2010-04-22-14-31 start
To simplify work, textbook
define blue term as φ(x) and
define red term as ψ(x)
---eqn.BH031
width
<a name="ch11a064">
for p=2, eqn.BH031 become textbook eqn.11.6
x=T
∫
x=0
φ(x)2
dx
≦
C
x=T
∫
x=0
φ(x) ψ(x) dx
---page 168
---line 14
---eqn.11.6
width of above equation
<a name="ch11a065">
where C=2.
The critical feature of this
inequality (11.6) is that the
function φ(x) is raised to a
higher power on the left side
of the equation than on the
right. This is far from a
minor detail; it open up the
possibility of a maneuver
which has featured in thousands
of investigations.
<a name="ch11a066">Index beginIndex this file
Now apply Schwarz's inequality to
the right hand side of eqn.BH031.
Get the following [ignore p/(p-1)]
---eqn.BH032
width
<a name="ch11a068">
Whole equation is eqn.BH031 less than side plus
eqn.BH032 greater than side and coefficient p/(p-1)
x=T
∫
x=0
φ(x)p
dx
≦
p
p-1
{
x=T
∫
x=0
[φ(x)p-1]2
dx
}
1/2
{
x=T
∫
x=0
[ψ(x)]2
dx
}
1/2
---eqn.BH033
width
<a name="ch11a069">
If φ(x)≡0, whole equation is 0=0. We can assume φ(x)≠0.
For f(x)>0 and x>0. eqn.11.5A tell us that φ(x)>0
For φ(x)>0, we divide equation both side by φ(x)
inequality not change direction.
<a name="ch11a070">
2010-04-22-15-32
OOPS ! choked !! eqn.BH033 fits p=2 only !!
If p≠2, left side φ(x)p and right side
φ(x)2p-2 are different.
2010-04-22-15-36 stop
2010-04-22-16-52 start
Continue with the case p=2, which is textbook
problem 11.1. For p=2, eqn.BH033 become
<a name="ch11a071">Index beginIndex this file
x=T
∫
x=0
φ(x)2
dx
≦
C
{
x=T
∫
x=0
φ(x)2
dx
}
1/2
{
x=T
∫
x=0
ψ(x)2
dx
}
1/2
---page 168
---line 21
---eqn.11.7
width of above equation
<a name="ch11a072">
Red terms are common at both side.
Cancel a square root of them, only
less than side left a square root
of red term. We have
---page 168
---line 26
---eqn.11.8
width of above equation
<a name="ch11a074">
Refer to eqn.BH033,
p/(p-1) become C in eqn.11.8
with p=2 , C=2/(2-1)=2
eqn.11.8 take square, recover
φ(x) and ψ(x) we have target
equation eqn.11.1.
Problem 11.1 is solved.
2010-04-22-17-11 stop
<a name="ch11a075">
2010-04-22-19-25 start
Now try solve eqn.BH001
Get trouble at this step
[[
Now apply Schwarz's inequality to
the right hand side of eqn.BH031.
]]
<a name="ch11a076">Index beginIndex this file
Because Schwarz's inequality
use square root operation.
That is in
1/p + 1/q = 1 ---eqn.BH034
Schwarz's inequality insist
p=q=2.<a name="ch11a077">Exercise 9.7
Holder's Inequality for Integrals
eqn.BC001 can be used.
Re-write eqn.BC001 as following
but change p,q to m,n. Because
eqn.BH031 use p and can not be
mixed.
Red term is I1
Blue term is I2
---page 150
---line 16
---eqn.BC001
width of above equation
<a name="ch11a079">
where, as usual, 1<m<∞ and
m-1+n-1=1 ---eqn.BC002
Now apply Integral Holder to
the right hand side of eqn.BH031
Consider
eqn.BC001 f(x) = φ(x)p-1 ---eqn.BH035
eqn.BC001 g(x) = ψ(x) ---eqn.BH036
eqn.BC001 w(x) ≡ 1 ---eqn.BH037
get the following
<a name="ch11a081">Index beginIndex this file
eqn.BH031 less than side is
our concern
eqn.BH038 less than side is
intermediate step
eqn.BH038 greater than side
is our concern
Now put eqn.BH031 less than
side and eqn.BH038 greater
than side into one equation
as following
<a name="ch11a083">
Pay attention to φ(x) power.
Less than side φ(x) power is p
Greater than side φ(x) power is
(p-1)*m where m is variable.
We want
p=(p-1)*m ---eqn.BH040
so that eqn.BH039 two sides have
similar term which can be merged.
<a name="ch11a084">
Solve eqn.BH040 for m get
m=p/(p-1) ---eqn.BH041
Because
1/m + 1/n = 1 ---eqn.BH042
solve for n get
1/n = 1-1/m = 1-(p-1)/p = 1/p
get
n = p ---eqn.BH043
<a name="ch11a085">
Now put m and n into eqn.BH039
get
---eqn.BH045
width
<a name="ch11a087">
From eqn.BH045 to eqn.BH046
is FLOP operation.
What is flop?
Two red terms in eqn.BH045 are the same.
It is possible to merge two to one. Get next equation.
{
x=T
∫
x=0
[φ(x)p]
dx
}
1-[(p-1)/p]
≦
p
p-1
{
x=T
∫
x=0
[ψ(x)]p
dx
}
1/p
---eqn.BH046
width
<a name="ch11a088">
In eqn.BH046, red power has value 1/p
eqn.BH046 whole equation raise power p, get
<a name="ch11a089">
Recover φ(x) and ψ(x) from
eqn.11.5A and eqn.11.5B.
we find that
eqn.BH047 is the solution of
generalized problem eqn.BH001
2010-04-22-20-43 stop
<a name="ch11a090">Index beginIndex this file
2010-04-23-12-33 start
■ Problem 11.2
The Discrete Hardy's Inequality
Show that for any sequence of
nonnegative real numbers a1,a2,
...,an one has the inequality
---page 169
---line 13
---eqn.11.9
width
<a name="ch11a092">
eqn.11.9 in expanded form is next
(
a1
1
)
2
+
(
a1+a2
2
)
2
+
(
a1+a2+a3
3
)
2
+
(
.....
)
2
+
(
a1+a2+a3+.....+aN
N
)
2
≦
4
(a12+a22+...+aN2 )
---eqn.BH048
width of above equation
2010-04-23-12-59 stop
<a name="ch11a093">
2010-04-23-14-11 start
eqn.BH048 tell us that
Problem 11.2 is a two
sequence comparison problem.
First sequence is a1,a2,...,an
Second sequence is
gradually AMed first sequence.
(see eqn.BH048 less than side)
<a name="ch11a094">
eqn.BH048 compare the sum
of element squared result.
If without factor 4, second
sequence has greater value.
Second sequence square sum is
not greater than four times
of first sequence square sum.
<a name="ch11a095">Index beginIndex this filegradually AMed sequence is
not new. We saw gradually GMed
sequence eqn.2.15 before.
<a name="ch11a096">
Integral Hardy's inequality
start from target equation
eqn.11.1 less than side
using integration by parts
After we drop boundary term,
we get "preflop" inequality
eqn.11.4<a name="ch11a097">
Discrete Hardy's Inequality
start from target equation
eqn.11.9 less than side via
summation by parts, we get
the next eqn.11.10
---page 169
---line 20
---eqn.11.10
width of above equation
<a name="ch11a099">
Let us define new {b} sequence
Its element bn is defined as
bn=(a1+a2+...+an)/n ---eqn.BH049
//bn and Aneqn.BH072 are identical
{b} sequence is
gradually AMed {a} sequence.
// AM=Arithmetic Mean
---page 169
---line 13
---eqn.BH050
width
<a name="ch11a101">
Discrete "preflop" inequality eqn.11.10 can be written as
n=N
∑
n=1
bn*bn
≦ 2
n=N
∑
n=1
bn*an
---page 169
---line 20
---eqn.BH051
{b} sequence is gradually AMed {a} sequence.
width of above equation
2010-04-23-15-09 stop
<a name="ch11a102">
2010-04-23-16-29 start
Summation by parts is the
natural analog of integration
by parts, although it is a bit
less mechanical. Here, for
example, we must decide how to
represent 1/n2 as a difference.
<a name="ch11a103">
There are two choices.
First define sn as sum of 1/k2
from k=n to k=∞
Second define s~n as sum of 1/k2
from k=1 to k=n
Mathematics expression is next
<a name="ch11a105">Index beginIndex this file
Above two summation are very
different.
Red sn sum from k=n to k=∞
Blue s~n sum from k=1 to k=n
Red sn sum better? or
Blue s~n sum better?
Textbook suggest experimentation
and choose red sn sum.
<a name="ch11a106">
We define TN=eqn.11.9 less than
side, use sn to express TN as
following
---page 170
---line 7
---eqn.BH055
width
<a name="ch11a110">
Above equation has sn and sn+1. We unify them
to sn only. Shift index of negative (blue) term, get
TN=
n=N
∑
n=1
sn*
(a1+a2+...+an)2
-
N+1
∑
n=2
sn*
(a1+a2+...+an-1)2
---page 170
---L.8 borrow
---eqn.BH056
width of above equation
<a name="ch11a111">Index beginIndex this file
eqn.BH055 positive to eqn.BH056
positive (black) not change.
eqn.BH055 negative (blue) to
eqn.BH056 negative (red) changed.
<a name="ch11a112">
Verify.
Blue start from n=1, for n=1
sn+1*(a1+a2+...+an)2
become s2*(a1)2 ---eqn.BH057
Red start from n=2, for n=2
sn*(a1+a2+...+an-1)2
become s2*(a1)2 ---eqn.BH058
eqn.BH057 and eqn.BH058 are
the same.
<a name="ch11a113">
Blue end at n=N, for n=N
sn+1*(a1+a2+...+an)2
become sN+1*(a1+a2+...+aN)2 ---eqn.BH059
Red end at n=N+1, for n=N+1
sn*(a1+a2+...+an-1)2
become sN+1*(a1+a2+...+aN)2 ---eqn.BH060
eqn.BH059 and eqn.BH060 are
the same.
<a name="ch11a114">
Shift index of negative term
from eqn.BH055 to eqn.BH056
is correct.
<a name="ch11a115">
Now let us target at eqn.BH056
sn unified for positive and
negative terms. But their sum
range are different.
We keep common (between positive
and negative in eqn.BH056) and
isolate not paired terms.
2010-04-23-17-55 here
<a name="ch11a116">Index beginIndex this file
In our focus eqn.BH056
positive sum from n=1 to n=N
negative sum from n=2 to n=N+1
positive n=1 is single
negative n=N+1 is single
from n=2 to n=N are paired
<a name="ch11a117">
TN= ---eqn.BH061 begin
+s1*a12
+∑[n=2,N]sn*(a1+a2+...+an)2
-∑[n=2,N]sn*(a1+a2+...+an-1)2
-sN+1*(a1+a2+...+aN)2 ---eqn.BH061 end
<a name="ch11a118">
First two lines come from positive.
Last two lines come from negative.
Middle two lines ∑[n=2,N] has the
form c*c-d*d, which is (c+d)(c-d)
<a name="ch11a119">
TN= ---eqn.BH062 begin
+{ ∑[n=2,N]sn*
[(a1+a2+...+an)+(a1+a2+...+an-1)]
*[(a1+a2+...+an)-(a1+a2+...+an-1)]
}
+s1*a12
-sN+1*(a1+a2+...+aN)2 ---eqn.BH062 end
<a name="ch11a120">
TN= ---eqn.BH063 begin
+{ ∑[n=2,N]sn*
[2*(a1+a2+...+an-1) + an]
*[an]
}
+s1*a12
-sN+1*(a1+a2+...+aN)2 ---eqn.BH063 end
<a name="ch11a121">Index beginIndex this file
TN= ---eqn.BH064 begin
+{∑[n=2,N]sn*[2*an(a1+a2+...+an-1)
+ an*an]
}
+s1*a12
-sN+1*(a1+a2+...+aN)2 ---eqn.BH064 end
<a name="ch11a122">
Here, we insert two zeros.
TN= ---eqn.BH065 begin
+{∑[n=2,N]sn*[2*an(a1+a2+...+an-1)
+ an*an]
} +s1*a12
+(s1*a12-s1*a12) //red is zero
+∑{n=2,N}(sn*an2-sn*an2) //blue is zero
-sN+1*(a1+a2+...+aN)2 ---eqn.BH065 end
<a name="ch11a123">
TN= ---eqn.BH066 begin
+{∑[n=2,N]sn*[2*an(a1+a2+...+an-1)
+ an*an]
} +s1*a12
+(s1*a12-s1*a12)
+∑{n=2,N}(sn*an2-sn*an2)
-sN+1*(a1+a2+...+aN)2 ---eqn.BH066 end
<a name="ch11a124">
Purple s1*a12 contribute to change
from ∑[n=2,N] to ∑[n=1,N]
Green sn*an2 contribute to extend
from (a1+a2+...+an-1)
to (a1+a2+...+an)
After they move to right place
<a name="ch11a125">
TN= ---eqn.BH067 begin
+{∑[n=1,N]sn*2*an(a1+a2+...+an)}
-s1*a12
-∑{n=2,N}sn*an2
-sN+1*(a1+a2+...+aN)2 ---eqn.BH067 end
<a name="ch11a126">Index beginIndex this file
Up to here eqn.BH067 is equality.
We drop three negative terms.
TN become less than or equal to
the positive term.
// 5=7-2, drop '-2', get 5<7
Recover TNoriginal face,
we have
---page 170
---line 11
---eqn.11.11
width of above equation
<a name="ch11a128">
eqn.11.11 is very close to our
goal eqn.11.10.
If we want flop eqn.11.11, that
is to move (a1+a2+...+an) from
right side to left side. We need
the red sn change to 1/n, to meet
left side (a1+a2+...+an)/n
2010-04-23-19-10 stop
<a name="ch11a129">
2010-04-23-20-06 start
sn is defined in eqn.BH052
Next find sn and 1/n relation.
Telescope expansion 010203, 01 close to eqn.BH069
width of above equation
<a name="ch11a132">
eqn.BH069 red equality is
telescoping result.
Expand eqn.BH069 summation term,
except first 1/(n-1), all other
terms are cancelled by the next
twin terms. Result is 1/(n-1)
which is less than equal to 2/n
for n greater than one.
<a name="ch11a133">
The case n=1, s1 is
s1= 1/1 + s2---eqn.BH070
eqn.BH068 and eqn.BH069 tell
us that s2≦2/2, so
s1≦ 1/1 + 2/2 = 2
or
s1≦ 2/n = 2/1---eqn.BH071
For all n≧1, we have
---page 171
---line 3
---eqn.11.13
width of above equation
<a name="ch11a137">
Busy few hours, change from
sn to 2/n and get the 1/n.
Our goal is eqn.11.10
Compare eqn.11.10 with eqn.11.13
there is one point different.
eqn.11.10 require coefficient 2
eqn.11.13 give us coefficient 4.
Integral version eqn.11.4 use 2
2010-04-23-21-02 here
2010-04-23-21-17 start
If we flop eqn.11.13, we get 8
but target eqn.11.9 require 4.
<a name="ch11a138">
Go back to eqn.11.10 which
textbook reprint as eqn.11.14.
n=N
∑
n=1
{
1
n
(a1+a2+...+an)
}
2
≦ 2
n=N
∑
n=1
{
1
n
(a1+a2+...+an)
}
an
---page 171
---line 18
---eqn.11.14
same as eqn.11.10
width of above equation
<a name="ch11a139">
Above two failures both start
from eqn.11.9 less than side
and use summation by parts. The
next third method, textbook use
delta calculation. To simplify
expression, define
An=(a1+a2+...+an)/n ---eqn.BH072
//bneqn.BH049 and An are identical
Attention: An definition do not
contain summation (∑[n=1,n=N])
<a name="ch11a140">Index beginIndex this file
Define
Δn=An*An-2*An*an ---eqn.BH073
Δn do not contain summation.
Δn is eqn.11.14 less than
side minus greater than side
and no n-sum.
Less than side minus greater
than side, that is Δn≦0
Below, we verify Δn≦0
eqn.11.14 is simply summation
of Δn<a name="ch11a141">
Δndefinition contains an. First
thing to do is to eliminate an.
Express Δn in terms of An and
An-1<a name="ch11a142">
an=(a1+a2+...+an)
-(a1+a2+...+an-1) ---eqn.BH074
an=An*n - An-1*(n-1) ---eqn.BH075
Use eqn.BH075 to eliminate an in
eqn.BH073, get
Δn=An*An-2*An*[An*n - An-1*(n-1)]
Δn=(1-2n)*An*An+2*(n-1)*An*An-1 ---eqn.BH076
<a name="ch11a143">
We do not want the term An*An-1
We want An2+An-12Humble bound will give us
An*An-1≦(An2+An-12)/2 ---eqn.BH077
Then we have
Δn≦(1-2n)*An2+(n-1)*(An2+An-12) ---eqn.BH078
or
Δn≦(n-1)*An-12-n*An2 ---eqn.BH079
<a name="ch11a144">
Now we sum Δn from n=1 to n=N
n=N
∑
n=1
Δn
≦
n=N
∑
n=1
{(n-1)*An-12-n*An2}
The last inequality is one that telescopes beautifully.
(textbook page 172 line 11) ---eqn.BH080 begin
<a name="ch11a145">Index beginIndex this file
=(1-1)*A1-12-1*A12 //this line: n=1
+(2-1)*A2-12-2*A22 //this line: n=2
+(3-1)*A3-12-3*A32 //this line: n=3
+.....
+(N-1)*AN-12-N*AN2 //this line: n=N
---page 172
---line 13
---eqn.BH080 end
Telescope expansion 010203 ; 03 is right here.
width of above equation
<a name="ch11a146">
Red cancel red, blue cancel blue,
purple cancel purple, all the way
to the end. Only two terms left
are first and last
(1-1)*A1-12-N*AN2
But (1-1)=0, the net result is
∑[n=1,N]Δn≦-N*AN2 ---eqn.BH081
<a name="ch11a147">
∑[n=1,N]Δn is short writing of
eqn.11.14 which is same as
eqn.11.10
Delta calculation let us verified
eqn.11.10.
2010-04-23-22-28 here
<a name="ch11a148">
eqn.11.10 is same as eqn.BH051
It is easier to see from eqn.BH051.
(eqn.BH051 is verified with Delta
calculation)
Apply Cauchy's inequality to
eqn.BH051 right hand side
(eqn.11.10 right hand side)
we get
∑[n=1,N]bn2≦2*∑[n=1,N]bnan
≦2*√{∑[n=1,N]bnbn}*√{∑[n=1,N]anan} ---eqn.BH082
//bn and An are identical; an not
<a name="ch11a149">
eqn.BH082 first line (black≦red)
is eqn.11.10
eqn.BH082 red part is Cauchy's
inequality.
Head and tail of eqn.BH082
∑[n=1,N]bn2 ---eqn.BH083
≦2*√{∑[n=1,N]bn2}*√{∑[n=1,N]an2}
is equation just about flop.
<a name="ch11a150">
We move blue term to one side
of the equation, this separation
is flop. The result is
√{∑[n=1,N]bn2} ---eqn.BH084
≦2*√{∑[n=1,N]an2}
<a name="ch11a151">
bn is defined at eqn.BH049
Recover bn and square eqn.BH084
we get our goal eqn.11.9Problem 11.2 solved.
2010-04-23-22-50 stop
2010-04-24-21-40 done first proofread
2010-04-25-12-02 done second proofread
2010-04-25-12-40 done spelling check
<a name="ch11a152">
2010-04-25-23-42
"Update 2010-04-26" add three "∑{n=2,N}"
You can still find many small bugs.
No one proofread for me.
<a name="Copyright">Index beginIndex this file
2009-06-19-10-48
If you are interested in inequality,
suggest you buy the book
The Cauchy-Schwarz Master Class
written by Prof. J. Michael Steele
The Cauchy-Schwarz Master Class is
this web page's textbook.
To buy textbook, that is to show thanks
to Prof. J.M. Steele's great work.
and it is also respect copyright law.
Thank you. Freeman 2009-06-19
The Cauchy-Schwarz Master Class
J. MichaelSteele★★★★★
ISBN 978-0-521-54677-5
2009-06-19-10-56