Boncelet Homepagehttp://www.ece.udel.edu/~boncelet/Home page for Charles Boncelet.enSun, 09 Apr 2017 02:51:56 GMTNikola (getnikola.com)http://blogs.law.harvard.edu/tech/rssSSIS Notebook Tutorialhttp://www.ece.udel.edu/~boncelet/blog/ssis-notebook.htmlCharles Boncelet<div><p>SSIS is a steganography technique developed by Lisa Marvel (my graduate
student and scientist with the US Army Research Laboratory), Charles Retter (also ARL), and myself in the late 1990's and early 2000's. We published several papers describing the technique. We also were awarded a US Patent. Those papers and the patent have been cited over a thousand times.</p>
<p>The notebook, <a href="http://www.ece.udel.edu/~boncelet/blog/research/SSIS.html">[html]</a> and
<a href="http://www.ece.udel.edu/~boncelet/blog/research/SSIS.ipynb">[ipynb]</a>, presents the basics of SSIS.</p>
<p>References: <a href="http://www.ece.udel.edu/~boncelet/papers/SSIS99.pdf">[TrIP99]</a> <a href="http://www.ece.udel.edu/~boncelet/papers/HidingICIP98.pdf">[ICIP98]</a></p></div>http://www.ece.udel.edu/~boncelet/blog/ssis-notebook.htmlSun, 09 Apr 2017 02:29:05 GMTBlock Arithmetic Coding Tutorialhttp://www.ece.udel.edu/~boncelet/blog/block-arithmetic-coding-tutorial.htmlCharles Boncelet<div><h2>Block Arithmetic Coding Tutorial</h2>
<p>I've written two Jupyter notebooks describing my Block Arithmetic Coding (BAC).
BAC is a variable to fixed encoder (parse the input into variable length
strings and encode each with a fixed length output string). I consider it
the best work I've done, but hardly anyone read the paper. Hopefully, some
people will at least read my tutorial. </p>
<p>The first notebook goes over the basics of BAC and shows how to compute
it's coding efficiency. <a href="http://www.ece.udel.edu/~boncelet/research/BAC_Basics.ipynb">[ipynb]</a>
<a href="http://www.ece.udel.edu/~boncelet/research/BAC_Basics.html">[html]</a></p>
<p>The second notebook provides a simple encoder and decoder for binary inputs.
<a href="http://www.ece.udel.edu/~boncelet/research/BAC_Coder_Decoder.ipynb">[ipynb]</a>
<a href="http://www.ece.udel.edu/~boncelet/research/BAC_Coder_Decoder.html">[html]</a></p></div>http://www.ece.udel.edu/~boncelet/blog/block-arithmetic-coding-tutorial.htmlFri, 31 Mar 2017 21:03:06 GMTMy Book is Out and Website Updatehttp://www.ece.udel.edu/~boncelet/blog/my-book-is-out-and-website-update.htmlCharles Boncelet<div><h2>My Book is Out and Website Update</h2>
<p>My book, <strong>Probability, Statistics, and Random Signals</strong>, came out in 2016. I
was busy doing lots of things last spring, summer, fall, and this winter
(visiting China, doing a sabbatical in Idaho, and taking a study abroad
trip to London) and I'm finally getting around to updating my website. </p>
<p>It is available through my publisher, <em>Oxford University
Press</em>, and through Amazon.com.</p>
<p>My book is designed for a junior or senior year class, but many reviewers
suggested they may use it in their graduate classes. It has
lots of carefully drawn illustrations, computational
examples using Matlab, Python, and R, and many worked out examples. </p>
<p>The first nine chapters cover discrete and continuous probability.
The remaining five chapters cover statistics, hypothesis testing, random
vectors and linear regression, and random processes. </p>
<p>Contact me if you are interested in using the text in your
class or your studies. The Solutions Manual has solutions for every homework
problem (writing solutions to every problem took me much longer than I anticipated!).</p>
<p><a href="http://tinyurl.com/kxff7xm">Oxford University Press link to <em>Probability, Statistics, and Random Signals</em></a></p>
<p><a href="http://a.co/7QCFBP9" title="Probabilty, Statistics, and Random Signals">Amazon link to <em>Probability, Statistics, and Random Signals</em></a></p>
<p><img src="http://www.ece.udel.edu/~boncelet/images/boncelet_book_cover.jpg" alt="Book Cover" style="width:200px;height:240px;"></p></div>http://www.ece.udel.edu/~boncelet/blog/my-book-is-out-and-website-update.htmlFri, 31 Mar 2017 20:56:07 GMTSolar Insolationhttp://www.ece.udel.edu/~boncelet/blog/solar-insolation.htmlCharles Boncelet<div><p><em>I read Steve Goddard's blog, <em>Real Science</em>. He is doing a wonderful job exposing the fraudulent manipulation of the temperature record by the climate alarmist crowd (e.g., NASA). </em></p>
<p><em>But lately, he has gotten off-track talking about climate physics. I sent the following comment to one of his posts:</em></p>
<p>There’s way too much imprecision in this discussion–from all sides.</p>
<p>I used the solar calculator at http://www.esrl.noaa.gov/gmd/grad/solcalc/ to record the sun’s elevation throughout the day. I chose Jun 21, 2014, the summer solstice and recorded the sun’s elevation hourly for three locations: the North Pole, Anchorage Alaska, and a point on the equator.</p>
<p>At the north pole, the sun’s elevation is about 23.5 degrees all day. The sun at Anchorage reaches an elevation of about 52 degrees. The sun at the equator reaches an elevation of 66.6 degrees.</p>
<p>Solar insolation is proportional to the cos(90-elevation angle) for angles above the horizon. Integrating over the 24 hours, I get the following: the north pole sees 9.5 hours direct sun equivalent, Anchorage sees 8.7 hours, and the equator 7.0 hours.</p>
<p>So Steve/Tony is right: at the summer solstice, the Arctic gets more solar insolation than the tropics.</p>
<p>However, this statement is misleading in many ways. I’ve neglected the effect of the lower sun’s rays at the pole being absorbed more by the atmosphere than at the equator. The peak solar elevation is much higher at the equator than the pole.</p>
<p>Furthermore, the summer solstice is the most favorable day of the year for the north pole. That’s the day it gets most solar insolation and the day the equator gets least (though the equator doesn’t vary much). (BTW, the south pole gets no insolation that day.)</p>
<p>Finally, a day’s solar insolation affects that location’s change in temperature far more than it affects it’s absolute temperature. The tropics are warm all year around, but the pole is very cold coming out of winter.</p>
<p>Blaming the cold weather at the poles on the lack of greenhouse gases is wrong. (BTW, the standard assumption is that CO2 is distributed uniformly.) The poles are cold because they get less sunlight over the year.</p>
<p>Water is immensely important to the planet’s temperature distribution. The oceans transfer heat from the tropics poleward, clouds and thunderstorms cool the tropics, and clouds help retain heat at night. The atmosphere also contributes with convection and winds. It’s complicated and oversimplifying is misleading at best.</p></div>climatepythonweatherhttp://www.ece.udel.edu/~boncelet/blog/solar-insolation.htmlWed, 19 Nov 2014 00:20:42 GMTFun with FFT IPython Notebookhttp://www.ece.udel.edu/~boncelet/blog/fun-with-fft-ipython-notebook.htmlCharles Boncelet<p>Another notebook from my DSP class: Fun with FFT. It's part of a lecture I gave on Halloween. It really needs to be cleaned up, but I didn't want to wait any longer before posting it. See my Ipython page <a href="http://www.ece.udel.edu/~boncelet/ipython.html">IPython</a>.</p>http://www.ece.udel.edu/~boncelet/blog/fun-with-fft-ipython-notebook.htmlSat, 08 Nov 2014 18:10:10 GMTFiltering IPython Notebookhttp://www.ece.udel.edu/~boncelet/blog/filtering-ipython-notebook.htmlCharles Boncelet<p>I've added a simple IPython notebook illustrating using Python for signal processing. I use a few of the functions in scipy.signal. See my Ipython page <a href="http://www.ece.udel.edu/~boncelet/ipython.html">IPython</a>.</p>http://www.ece.udel.edu/~boncelet/blog/filtering-ipython-notebook.htmlSat, 08 Nov 2014 17:39:48 GMTBloom Energy Subsidyhttp://www.ece.udel.edu/~boncelet/blog/bloom-energy-subsidy.htmlCharles Boncelet<div><p>I sent the following letter to the News Journal:</p>
<p>A recent report in the News Journal indicates that Delmarva Power customers are subsidizing Bloom Energy in excess of $3 million per month. That is about $40 million per year.</p>
<p>I have a couple of questions for the governor and legislature:</p>
<p>1) If Bloom actually creates 1000 jobs (likely an optimistic scenario), the subsidy works out to $40,000 per job per year. Are these jobs actually worth $40,000 each to the state of Delaware? If not, why are we paying this subsidy?</p>
<p>2) If these jobs are important to the State of Delaware, why is it fair that only the Delmarva Power customers have to pay the subsidy, why not all the citizens in Delaware?</p></div>Delawarehttp://www.ece.udel.edu/~boncelet/blog/bloom-energy-subsidy.htmlMon, 01 Sep 2014 21:11:27 GMTNew IPython Notebooks: Getting Started and Linear Regressionhttp://www.ece.udel.edu/~boncelet/blog/new-ipython-notebooks-getting-started-and-linear-regression.htmlCharles Boncelet<div><p>I added two more IPython Notebooks. </p>
<p>The first is a <em>Getting Started</em> notebook with lots of advice on where to get python (anaconda or canopy), how to run IPython remotely (wakari.io and sagemathcloud.com), and where to find tutorials. </p>
<p>The second is an example showing four different ways to compute <em>linear regression</em> estimates with Python. Bottom line: for most problems, use statsmodels with the QR option.</p>
<p>See my <a href="http://www.ece.udel.edu/~boncelet/ipython.html">IPython</a> page for more information and to get the files. </p></div>http://www.ece.udel.edu/~boncelet/blog/new-ipython-notebooks-getting-started-and-linear-regression.htmlThu, 28 Aug 2014 03:38:20 GMTSampling Theorem -- Aliasinghttp://www.ece.udel.edu/~boncelet/blog/sampling-theorem-aliasing.htmlCharles Boncelet<div><p>The sampling theorem consists of two parts: the first is sampling and aliasing and the second is reconstruction. Taken together, the sampling theorem is the
most important concept in digital signal processing: it is why we can
use digital computers to analyze continuous time signals.</p>
<p>I wrote up a quick IPython script demonstrating aliasing. Basically, we have
a series of plots showing two continuous time sinusoids and the resulting
samples. The digital frequencies are the same.</p>
<p>In later work, we will demonstrate reconstruction.</p>
<p>Here are the files: <a href="http://www.ece.udel.edu/~boncelet/SamplingTheoremAliasing.ipynb">IPynb</a>, <a href="http://www.ece.udel.edu/~boncelet/SamplingTheoremAliasing.html">HTML</a></p></div>dspipythonsignal processinghttp://www.ece.udel.edu/~boncelet/blog/sampling-theorem-aliasing.htmlMon, 18 Aug 2014 02:28:12 GMTRyan Howard and clutch hittinghttp://www.ece.udel.edu/~boncelet/blog/ryan-howard-and-clutch-hitting.htmlCharles Boncelet<div><p>It's been a frustrating couple of years for this Phillies fan. The Phillies
management (mostly the GM Amaro) has made many mistakes, but none bigger than
signing Ryan Howard to a huge contract. </p>
<p>He can't hit for average, can't run, can't throw, and can't field (in fairness,
he's okay at catching throws in the dirt, but has no range and doesn't
field enough batted balls). </p>
<p>But, people will say, he hits for power and gets RBI's. Let's look at
the numbers:</p>
<p>As of today (17 August 2014), Howard is in a three-way tie for 14th in the NL
with 18 homeruns through 117 games. He's on pace to hit
25 = (162/117)*18 homeruns on the season. </p>
<p>Is that enough to compensate for all his other failings? Not in my book.</p>
<p>How about RBI's? Howard is currently 3rd in the NL with 77 RBI's. That's
pretty good, right? Doesn't it point to his "clutch hitting"?</p>
<p>No it doesn't. Clutch hitting is mostly a fiction. Players get RBI's
because they bat with runners on base, not because of some mystical ability
to hit better with runners one base. (If they could hit better with runners
on base, why don't they hit better without runners on base? Are they just
lazy?)</p>
<p>Howard leads the NL with 373 baserunners on base during his at bats. That's
about 0.74 baserunners per at bat. It's an amazing testament to the first
three hitters in the Phillies lineup (generally Revere, Rollins, and Utley).</p>
<p>Those runners score at a 16% rate during or because of Howard's at bat. That's
currently 81st in the NL. In other words, 80 players in the NL knock in runs
at a greater rate than Howard.</p>
<p>Howard's OPS is .676, good enough for 59th best in the NL. (MLB.com only lists
71 batters, so 59th is not very good.)</p>
<p>But he gets RBI's, they say.</p>
<p>(All data is from MLB.com and baseball-reference.com before today's game.)</p></div>baseballphillieshttp://www.ece.udel.edu/~boncelet/blog/ryan-howard-and-clutch-hitting.htmlSun, 17 Aug 2014 19:31:56 GMT