Commit b0db00a7 authored by Adam Taheraly's avatar Adam Taheraly

je sais pas trop ce que c'est

parent 12be553b
This diff is collapsed.
This diff is collapsed.
% Created 2020-03-30 lun. 22:50
% Intended LaTeX compiler: pdflatex
\documentclass[11pt]{article}
\usepackage[utf8]{inputenc}
\usepackage[T1]{fontenc}
\usepackage{graphicx}
\usepackage{grffile}
\usepackage{longtable}
\usepackage{wrapfig}
\usepackage{rotating}
\usepackage[normalem]{ulem}
\usepackage{amsmath}
\usepackage{textcomp}
\usepackage{amssymb}
\usepackage{capt-of}
\usepackage{hyperref}
\usepackage[utf8]{inputenc}
\usepackage[T1]{fontenc}
\usepackage[francais, frenchb]{babel}
\usepackage{amsmath}
\usepackage{amssymb}
\usepackage{mathrsfs}
\author{Adam Taheraly\thanks{taheraly.adam@gmail.com}}
\date{Lundi 30 Mars}
\title{À propros du calcul de \(\pi\)}
\hypersetup{
pdfauthor={Adam Taheraly},
pdftitle={À propros du calcul de \(\pi\)},
pdfkeywords={MOOC-RR exo1},
pdfsubject={Exercice Mooc RR sur le calcul de \(\pi\)},
pdfcreator={Emacs 26.3 (Org mode 9.1.9)},
pdflang={Frenchb}}
\begin{document}
\maketitle
\tableofcontents
\section{En demandant à la lib maths}
\label{sec:orgc4eb4d4}
Mon ordinateur m'indique que \(\pi\) vaut \emph{approximativement}
\begin{verbatim}
pi
\end{verbatim}
\begin{verbatim}
[1] 3.141593
\end{verbatim}
\section{En utilisant la méthode des aiguilles de Buffon}
\label{sec:org3a95611}
Mais calculé avec la \textbf{méthode} des \href{https://fr.wikipedia.org/wiki/Aiguille\_de\_Buffon}{aiguilles de Buffon}, on obtiendrait
comme \textbf{approximation} :
\begin{verbatim}
set.seed(42)
N = 100000
x = runif(N)
theta = pi/2*runif(N)
2/(mean(x+sin(theta)>1))
\end{verbatim}
\begin{verbatim}
[1] 3.14327
\end{verbatim}
\begin{verbatim}
[1] 3.14327
\end{verbatim}
\section{Avec un argument "fréquentiel" de surface}
\label{sec:org757811b}
Sinon, une méthode plus simple à comprendre et ne faisant pas
intervenir d'appel à la fonction sinus se base sur le fait que
si \(X \sim U(0,1)\) et \(Y \sim
U(0,1)\) alors \(P[X^{2} + Y^{2} \le 1] = \pi/4\)
(voir \href{https://fr.wikipedia.org/wiki/M\%C3\%A9thode\_de\_Monte-Carlo\#D\%C3\%A9termination\_de\_la\_valeur\_de\_\%CF\%80}{méthode
de Monte Carlo sur Wikipédia}). Le code suivant illustre ce fait :
\begin{verbatim}
set.seed(42)
N = 1000
df = data.frame(X = runif(N), Y = runif(N))
df$Accept = (df$X**2 + df$Y**2 <=1)
library(ggplot2)
ggplot(df, aes(x=X,y=Y,color=Accept)) + geom_point(alpha=.2) + coord_fixed() + theme_bw()
\end{verbatim}
\begin{center}
\includegraphics[width=.9\linewidth]{/tmp/babel-A2ZKik/figureoqEQiY.png}
\end{center}
Il est alors aisé d'obtenir une approximation (pas terrible) de \(\pi\) en
comptant combien de fois, en moyenne, \(X^{2} + Y^{2}\)
est inférieur à 1 :
\begin{verbatim}
4*mean(df$Accept)
\end{verbatim}
\begin{verbatim}
[1] 3.156
\end{verbatim}
Emacs 26.3 (Org mode 9.1.9)
\end{document}
\ No newline at end of file
% Created 2020-03-30 lun. 22:45
% Intended LaTeX compiler: pdflatex
\documentclass[11pt]{article}
\usepackage[utf8]{inputenc}
\usepackage[T1]{fontenc}
\usepackage{graphicx}
\usepackage{grffile}
\usepackage{longtable}
\usepackage{wrapfig}
\usepackage{rotating}
\usepackage[normalem]{ulem}
\usepackage{amsmath}
\usepackage{textcomp}
\usepackage{amssymb}
\usepackage{capt-of}
\usepackage{hyperref}
\author{Adam Taheraly\thanks{taheraly.adam@gmail.com}}
\date{Lundi 30 Mars}
\title{À propros du calcul de \(\pi\)}
\hypersetup{
pdfauthor={Adam Taheraly},
pdftitle={À propros du calcul de \(\pi\)},
pdfkeywords={MOOC-RR exo1},
pdfsubject={Exercice Mooc RR sur le calcul de \(\pi\)},
pdfcreator={Emacs 26.3 (Org mode 9.1.9)},
pdflang={Frenchb}}
\begin{document}
\maketitle
\tableofcontents
\section{En demandant à la lib maths}
\label{sec:orgc5414c9}
Mon ordinateur m'indique que \(\pi\) vaut \emph{approximativement}
\begin{verbatim}
pi
\end{verbatim}
\begin{verbatim}
[1] 3.141593
\end{verbatim}
\section{En utilisant la méthode des aiguilles de Buffon}
\label{sec:orgdde632f}
Mais calculé avec la \textbf{méthode} des \href{https://fr.wikipedia.org/wiki/Aiguille\_de\_Buffon}{aiguilles de Buffon}, on obtiendrait
comme \textbf{approximation} :
\begin{verbatim}
set.seed(42)
N = 100000
x = runif(N)
theta = pi/2*runif(N)
2/(mean(x+sin(theta)>1))
\end{verbatim}
\begin{verbatim}
[1] 3.14327
\end{verbatim}
\begin{verbatim}
[1] 3.14327
\end{verbatim}
\section{Avec un argument "fréquentiel" de surface}
\label{sec:orgdef1a48}
Sinon, une méthode plus simple à comprendre et ne faisant pas
intervenir d'appel à la fonction sinus se base sur le fait que
si \(X \sim U(0,1)\) et \(Y \sim
U(0,1)\) alors \(P[X^{2} + Y^{2} \le 1] = \pi/4\)
(voir \href{https://fr.wikipedia.org/wiki/M\%C3\%A9thode\_de\_Monte-Carlo\#D\%C3\%A9termination\_de\_la\_valeur\_de\_\%CF\%80}{méthode
de Monte Carlo sur Wikipédia}). Le code suivant illustre ce fait :
\begin{verbatim}
set.seed(42)
N = 1000
df = data.frame(X = runif(N), Y = runif(N))
df$Accept = (df$X**2 + df$Y**2 <=1)
library(ggplot2)
ggplot(df, aes(x=X,y=Y,color=Accept)) + geom_point(alpha=.2) + coord_fixed() + theme_bw()
\end{verbatim}
\begin{center}
\includegraphics[width=.9\linewidth]{/tmp/babel-A2ZKik/figureGj86qm.png}
\end{center}
Il est alors aisé d'obtenir une approximation (pas terrible) de \(\pi\) en
comptant combien de fois, en moyenne, \(X^{2} + Y^{2}\)
est inférieur à 1 :
\begin{verbatim}
4*mean(df$Accept)
\end{verbatim}
\begin{verbatim}
[1] 3.156
\end{verbatim}
Emacs 26.3 (Org mode 9.1.9)
\end{document}
\ No newline at end of file
#+TITLE: Savoir faire un calcul simple soi-même
#+AUTHOR: Adam Taheraly
#+DATE: Mardi 31 Mars 2020
#+LANGUAGE: fr
# #+PROPERTY: header-args :eval never-export
#+HTML_HEAD: <link rel="stylesheet" type="text/css" href="http://www.pirilampo.org/styles/readtheorg/css/htmlize.css"/>
#+HTML_HEAD: <link rel="stylesheet" type="text/css" href="http://www.pirilampo.org/styles/readtheorg/css/readtheorg.css"/>
#+HTML_HEAD: <script src="https://ajax.googleapis.com/ajax/libs/jquery/2.1.3/jquery.min.js"></script>
#+HTML_HEAD: <script src="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.4/js/bootstrap.min.js"></script>
#+HTML_HEAD: <script type="text/javascript" src="http://www.pirilampo.org/styles/lib/js/jquery.stickytableheaders.js"></script>
#+HTML_HEAD: <script type="text/javascript" src="http://www.pirilampo.org/styles/readtheorg/js/readtheorg.js"></script>
* Chargement des librairies
#+begin_src R :results output :session *R* :exports both
library(prettyR)
#+end_src
* Chargement des données
#+begin_src R :results output :session *R* :exports results
df <- data.frame(x = c(14.0, 7.6, 11.2, 12.8, 12.5, 9.9, 14.9, 9.4, 16.9, 10.2, 14.9, 18.1, 7.3, 9.8, 10.9,12.2, 9.9, 2.9, 2.8, 15.4, 15.7, 9.7, 13.1, 13.2, 12.3, 11.7, 16.0, 12.4, 17.9, 12.2, 16.2, 18.7, 8.9, 11.9, 12.1, 14.6, 12.1, 4.7, 3.9, 16.9, 16.8, 11.3, 14.4, 15.7, 14.0, 13.6, 18.0, 13.6, 19.9, 13.7, 17.0, 20.5, 9.9, 12.5, 13.2, 16.1, 13.5, 6.3, 6.4, 17.6, 19.1, 12.8, 15.5, 16.3, 15.2, 14.6, 19.1, 14.4, 21.4, 15.1, 19.6, 21.7, 11.3, 15.0, 14.3, 16.8, 14.0, 6.8, 8.2, 19.9, 20.4, 14.6, 16.4, 18.7, 16.8, 15.8, 20.4, 15.8, 22.4, 16.2, 20.3, 23.4, 12.1, 15.5, 15.4, 18.4, 15.7, 10.2, 8.9, 21.0))
#+end_src
* Calcul
#+begin_src R :results output :session *R* :exports results
describe(df, num.desc = c("sd", "min", "mean", "median", "max"))
#+end_src
This diff is collapsed.
This diff is collapsed.
...@@ -60,6 +60,7 @@ Voici l'explication des colonnes donnée sur [[https://ns.sentiweb.fr/incidence/ ...@@ -60,6 +60,7 @@ Voici l'explication des colonnes donnée sur [[https://ns.sentiweb.fr/incidence/
| ~geo_insee~ | Code de la zone géographique concernée (Code INSEE) http://www.insee.fr/fr/methodes/nomenclatures/cog/ | | ~geo_insee~ | Code de la zone géographique concernée (Code INSEE) http://www.insee.fr/fr/methodes/nomenclatures/cog/ |
| ~geo_name~ | Libellé de la zone géographique (ce libellé peut être modifié sans préavis) | | ~geo_name~ | Libellé de la zone géographique (ce libellé peut être modifié sans préavis) |
** Téléchargement ** Téléchargement
La première ligne du fichier CSV est un commentaire, que nous ignorons en précisant ~skip=1~. La première ligne du fichier CSV est un commentaire, que nous ignorons en précisant ~skip=1~.
......
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment