From 68b4d9871f52077dc58edcf4efbcea582ef5bf8b Mon Sep 17 00:00:00 2001 From: Jean-Paul Dax <karlmarx79@hotmail.com> Date: Mon, 3 Apr 2017 17:27:28 +0200 Subject: [PATCH] correction date --- README.md | 2 ++ cours.tex | 44 ++++++++++++------------- tpFourier/mydata.txt | Bin 0 -> 72498 bytes tpFourier/tpFourier.tex | 70 ++++++++++++++++++++++++++++++++++++++++ 4 files changed, 94 insertions(+), 22 deletions(-) create mode 100644 tpFourier/mydata.txt create mode 100644 tpFourier/tpFourier.tex diff --git a/README.md b/README.md index 05e8af0..9fda9b0 100644 --- a/README.md +++ b/README.md @@ -1,3 +1,5 @@ # Mathématiques pour deuxième année ITI Ce projet contient la tentative de polycopié du cours de Mathématiques pour la filière ITI de hepia. + +Un bonus de 0.1 points sur la note de l'examen sera obtenu pour tout pull request réussi. \ No newline at end of file diff --git a/cours.tex b/cours.tex index 72a33b9..3afa941 100644 --- a/cours.tex +++ b/cours.tex @@ -2457,7 +2457,7 @@ b_j&\equiv A_j\sin(\phi_j), \end{align} on obtient \begin{equation} - f(t)=\sum_{j=0}^\infty \left(a_j\cos(j\omega t)+b_j\sin(j\omega t)\right). \label{eq_decomp_sincos} + f(t)=\sum_{j=0}^\infty \left(a_j\sin(j\omega t)+b_j\cos(j\omega t)\right). \label{eq_decomp_sincos} \end{equation} On a donc transformé une équation où on devait déterminer une amplitude et une phase, ce qui est très compliqué, en une autre équation où on doit déterminer uniquement deux amplitude. Par ailleurs, comme $\cos$ et $\sin$ sont indépendants, @@ -2469,7 +2469,7 @@ que les $a_j\cos(j\omega t)$ et $b_j\sin(j\omega t)$ approximent au mieux la fon On va donc considérer les fonctions d'erreur suivantes \begin{equation} - E^s_j=\int_0^T(f(t)-b_j\sin(j\omega t))^2\dd t,\quad E^c_j=\int_0^T(f(t)-a_j\cos(j\omega t))^2\dd t. + E^s_j=\int_0^T(f(t)-a_j\sin(j\omega t))^2\dd t,\quad E^c_j=\int_0^T(f(t)-b_j\cos(j\omega t))^2\dd t. \end{equation} Puis on va déterminer $a_j,b_j$ tels que $E_j^s$ et $E_j^c$ sont minimales. Pour ce faire on va utiliser les dérivées et déterminer nos coefficients en résolvant les équations @@ -2479,22 +2479,22 @@ et déterminer nos coefficients en résolvant les équations \end{align} Pour l'équation \eqref{eq_deriv_aj}, on a \begin{align} - \dDeriv{E^c_j}{a_j}&=\dDeriv{\int_0^T(f(t)-a_j\cos(j\omega t))^2\dd t}{a_j},\nonumber\\ - &=\underbrace{\dDeriv{(\int_0^Tf^2(t)\dd t)}{a_j}}_{=0}+\dDeriv{(a_j^2\int_0^T(\cos^2(j\omega t)\dd t))}{a_j}-\dDeriv{(2a_j\int_0^T(f(t)\cos(j\omega t)\dd t))}{a_j},\nonumber\\ - &=2a_j\int_0^T\cos^2(j\omega t)\dd t-2\int_0^Tf(t)\cos(j\omega t)\dd t,\nonumber\\ - &=2a_j\frac{T}{2}-2\int_0^T\cos(j\omega t)f(t)\dd t. + \dDeriv{E^c_j}{b_j}&=\dDeriv{\int_0^T(f(t)-b_j\cos(j\omega t))^2\dd t}{b_j},\nonumber\\ + &=\underbrace{\dDeriv{(\int_0^Tf^2(t)\dd t)}{b_j}}_{=0}+\dDeriv{(b_j^2\int_0^T(\cos^2(j\omega t)\dd t))}{b_j}-\dDeriv{(2b_j\int_0^T(f(t)\cos(j\omega t)\dd t))}{b_j},\nonumber\\ + &=2b_j\int_0^T\cos^2(j\omega t)\dd t-2\int_0^Tf(t)\cos(j\omega t)\dd t,\nonumber\\ + &=2b_j\frac{T}{2}-2\int_0^T\cos(j\omega t)f(t)\dd t. \end{align} Finalement on obtient \begin{equation} - a_j=\frac{2}{T}\int_0^T\cos(j\omega t)f(t)\dd t. + b_j=\frac{2}{T}\int_0^T\cos(j\omega t)f(t)\dd t. \end{equation} -Pour $b_j$ on a de façon similaire +Pour $a_j$ on a de façon similaire \begin{equation} - b_j=\frac{2}{T}\int_0^T\sin(j\omega t)f(t)\dd t. + a_j=\frac{2}{T}\int_0^T\sin(j\omega t)f(t)\dd t. \end{equation} En particulier si $j=0$, on a \begin{equation} -b_0=0,\quad a_0=\frac{2}{T}\int_0^T f(t)\dd t. +a_0=0,\quad b_0=\frac{2}{T}\int_0^T f(t)\dd t. \end{equation} On constate que $a_0/2$ correspond à la valeur moyenne de $f(t)$ dans $[0,T]$. Cela permet d'approximer des fonctions dons la valeur moyenne n'est pas nulle (les sinus et cosinus ont toujours @@ -2518,8 +2518,8 @@ Cela est dû à la propriété d'othorgonalité des fonctions sinus/cosinus. En multipliant l'équation \eqref{eq_decomp_sincos} par $\frac{2}{T}\sin(k \omega t)$ et en intégrant entre $0$ et $T$, on obtient \begin{align} -\frac{2}{T}\int_0^T f(t)\sin(k\omega t)\dd t&=\frac{2}{T}\sum_{j=0}^\infty \left(a_j\underbrace{\int_0^T\cos(j\omega t)\sin(k\omega t)\dd t}_{=0}+b_j\underbrace{\int_0^T\sin(j\omega t)\sin(k \omega t)}_{=\frac{T}{2}\delta_{jk}}\right),\nonumber\\ -\frac{2}{T}\int_0^T f(t)\sin(k\omega t)\dd t&=\sum_{j=0}^\infty b_j \delta_{jk}=b_k, +\frac{2}{T}\int_0^T f(t)\sin(k\omega t)\dd t&=\frac{2}{T}\sum_{j=0}^\infty \left(b_j\underbrace{\int_0^T\cos(j\omega t)\sin(k\omega t)\dd t}_{=0}+a_j\underbrace{\int_0^T\sin(j\omega t)\sin(k \omega t)\dd t}_{=\frac{T}{2}\delta_{jk}}\right),\nonumber\\ +\frac{2}{T}\int_0^T f(t)\sin(k\omega t)\dd t&=\sum_{j=0}^\infty a_j \delta_{jk}=a_k, \end{align} où $\delta_{jk}$ est le ``delta de Kronecker'', dont la définition est \begin{equation} @@ -2532,8 +2532,8 @@ où $\delta_{jk}$ est le ``delta de Kronecker'', dont la définition est En multipliant l'équation \eqref{eq_decomp_sincos} par $\frac{2}{T}\cos(k \omega t)$ et en intégrant entre $0$ et $T$, on obtient \begin{align} -\frac{2}{T}\int_0^T f(t)\cos(k\omega t)\dd t&=\frac{2}{T}\sum_{j=0}^\infty \left(a_j\underbrace{\int_0^T\cos(j\omega t)\sin(k\omega t)\dd t}_{=\frac{T}{2}\delta_{jk}}+b_j\underbrace{\int_0^T\sin(j\omega t)\sin(k \omega t)}_{=0}\right),\nonumber\\ -\frac{2}{T}\int_0^T f(t)\cos(k\omega t)\dd t&=\sum_{j=0}^\infty a_j \delta_{jk}=a_k. +\frac{2}{T}\int_0^T f(t)\cos(k\omega t)\dd t&=\frac{2}{T}\sum_{j=0}^\infty \left(a_j\underbrace{\int_0^T\cos(j\omega t)\sin(k\omega t)\dd t}_{=0}+b_j\underbrace{\int_0^T\cos(j\omega t)\cos(k \omega t)\dd t}_{=\frac{T}{2}\delta_{jk}}\right),\nonumber\\ +\frac{2}{T}\int_0^T f(t)\cos(k\omega t)\dd t&=\sum_{j=0}^\infty b_j \delta_{jk}=b_k. \end{align} \subsubsection{Les séries de Fourier en notations complexes} @@ -2543,9 +2543,9 @@ Cette écriture nous fait penser qu'il pourrait être possible de réécrire cet de nouveaux coefficients $c_n$, \begin{equation} c_n=\left\{\begin{array}{ll} - \frac{a_n+ib_n}{2}, & $\mbox{ si }$n<0\\ - \frac{a_0}{2}, & $\mbox{ si }$n=0\\ - \frac{a_n-ib_n}{2}, & $\mbox{ si }$n>0 + \frac{b_n+ia_n}{2}, & $\mbox{ si }$n<0\\ + \frac{b_0}{2}, & $\mbox{ si }$n=0\\ + \frac{b_n-ia_n}{2}, & $\mbox{ si }$n>0 \end{array}\right. \end{equation} Avec cette notation, on peut réécrire l'équation \eqref{eq_decomp_sincos} (exercice) comme @@ -2594,7 +2594,7 @@ ainsi que la notation $\omega_j=j\omega$, on peut réécrire cette équation &=\frac{1}{2\pi}\sum_{j=-\infty}^\infty (\Delta \omega_j)\fh(\omega_j) e^{i\omega_j t}. \end{align} Maintenant pour passer dans le cas où la fonction n'est pas périodique (la période est infinie), nous -devons prendre la limite $\omega_j\rightarrow 0$ dans l'équation précédente, et on voit apparaître une somme de Riemann +devons prendre la limite $\Delta \omega_j\rightarrow 0$ dans l'équation précédente, et on voit apparaître une somme de Riemann \begin{align} f(t)&=\frac{1}{2\pi}\sum_{j=-\infty}^\infty \lim\limits_{\Delta \omega_j\rightarrow 0}\Delta \omega_j\fh(\omega_j) e^{i\omega_j t},\nonumber\\ &=\frac{1}{2\pi}\int_{-\infty}^\infty \fh(\omega) e^{i\omega t}\dd\omega. @@ -2830,9 +2830,9 @@ Montrons à présent que la transformée inverse discrète de la transformée de discrète donne bien la suite de départ \begin{align} f[n]&=\frac{1}{N}\sum_{k=0}^{N-1} \fh[k] e^{\frac{2\pi i k n}{N}},\nonumber\\ - &=\frac{1}{N}\sum_{k=0}^{N-1} \sum_{m=0}^{N-1} f[m] e^{-\frac{2\pi k m}{N}} e^{\frac{2\pi i k n}{N}},\nonumber\\ - &=\frac{1}{N}\sum_{k=0}^{N-1} \sum_{m=0}^{N-1} f[m] e^{\frac{2\pi k (n-m)}{N}},\nonumber\\ - &=\frac{1}{N}\sum_{m=0}^{N-1} f[m] \sum_{k=0}^{N-1} e^{\frac{2\pi k (n-m)}{N}},\nonumber\\ + &=\frac{1}{N}\sum_{k=0}^{N-1} \sum_{m=0}^{N-1} f[m] e^{-\frac{2\pi i k m}{N}} e^{\frac{2\pi i k n}{N}},\nonumber\\ + &=\frac{1}{N}\sum_{k=0}^{N-1} \sum_{m=0}^{N-1} f[m] e^{\frac{2\pi i k (n-m)}{N}},\nonumber\\ + &=\frac{1}{N}\sum_{m=0}^{N-1} f[m] \sum_{k=0}^{N-1} e^{\frac{2\pi i k (n-m)}{N}},\nonumber\\ &=\frac{1}{N}\sum_{m=0}^{N-1} f[m] N \delta_{nm},\nonumber\\ &=f[n]. \end{align} @@ -2867,7 +2867,7 @@ La transformée de Fourier discrète étant un échantillonage de la transformé à temps discret, toutes les propriétés discutées pour la transformée de Fourier à temps discret restent valides. En particulier la transformée de Fourier discrète est périodique, de période $N$ \begin{equation} - f[n]=f[n+N]. + \fh[k]=\fh[k+N]. \end{equation} A démontrer en exercice. diff --git a/tpFourier/mydata.txt b/tpFourier/mydata.txt new file mode 100644 index 0000000000000000000000000000000000000000..137ee761c50bd13aaf37d06f5a6c78077a1bcfa5 GIT binary patch literal 72498 zcmeI5J(lFK4TK%fDK?+C{rE5OKn~{q7w}gADRz$?Ikd4+FGPW00R%y6RLZYDYWe5i zT)*-wwK??jyS_%Mlb`M9>u-IvOD)=Wd6k;7S|Gk&HPzeX?KMWJt9g0lGOf|y{Hn7z zZ_n?MP5a;S>a}zi^ZQEn`+L07yjG{<)qD0%Gx~<a>QpvLVgHuCdnx(bUM)8%Z{Phj z=ad}N2im!I<)mxf@6Sn=_V4R8)08faQI@)Om@$+6ba@)tdii^6FxS&Qmu5yWn%0{k z`14i3txJ6ct1;$=Fuz(YK6r-H=W3o`y-8zL|1}9r21(;|d*kuSV;aWNoolhyn0U=1 zh#3-hns+KAzf)d4(n~GY3qbX0Lq_#{jVqg@82)wDJ^cE*dU0QJXSObB-d@E_Z73OQ z=zY9OCqT;Gms4(jo9*rIASV3h&sTrdUd<6;k&;0t7t;L7W;~HK%TTwBqrP%8*EW>5 zzLf)So|g|hV`bkiFR$tD|Hi?|6uQH(#RjfesnhN>2e+K<KhiiXSg?kOzJ@mk)Mc+3 zNz6;GfWmjOSL*w#j5XGduiaodp$r445^km3#|y(l1|y3{AzB(SUp)x9@w$(VbxiXP zlM6JKhe7okHV_$8Yu;W{C@sjqpr)r+@Q}`ilZip=3EP>lv)Vv#C!s<M`Genec!=Yt z{vdhDH0z39`_)z$(!R{g+Y@`Y#gOz&1@n%f(X3m^l{73lR85u3BO7^p&CBzLMZe8< zwO1^fJstUMNhTn_Uxrgd*sh^5A)L&6hXZj+-i_Rmi4--VOUArRM?5Z5_l-86A9AW@ z{BU6Es}la)Dfvwtda?k!R(25J19$NE#Md6m7Qfxa=)0%;7bSFBY?um-bemIA&M59B zW)OJ5`(y%L-OLpZV|n^WnH!JvsKAh;;;l!b?)1n6TLzLpj)Bi;aNrru+t>-BGxj&E zFz{e?pykd3?QsM-;5g#AjYvFh+b%|UT>KUf$WJ`H#Se03h0S6BKk-JNZ}E2G6NKL4 z^Tf0CoOnL*&#eE5|2<!n*WTyLo?hsGrWZMDO-5qqZ8Cz)n!60-%-y%SPOl$xedh-J zadN{+C3JdHiQM6d2kj^E&}5}4<)bG}y=66~`pIevf|Uhna<bq_np%lY(!^~m$AZks zu?C0dWDtJK$xjLw1Sf@`jL)4X<DUhBka$}l&U(V+IO_@MEz5|Av;3UJ3ftcnE3Ini zJ*yf3`PzpOd+oDW^mP*HStrl3ij&W>N|Ut+)A}sJ&x#GGXT@gJYk5};ye;p$7{K0B z3_O*C#s01w=w{a$$WNVd$`Y{nC`)z~qp9*#F{s&94~g=l>e*G2<$?DqX;oE>|5Q~f zzqC4^^6S)TjQ!MU>d?zNex0)J6oH6O5qL^QtT-j3?p!4tEuI3C5$}3bSUvS9mhSqM z@O!IYpK{kS$zCR9p%=#l!6}Ya=1f}x2gZ!Pi#1><#j-c1x4DaB9Nvp9FL;yEM{AU? zK>7Bf6&vT+D?BlrS3$sB>ArZA;!tBXdt(rHO}5snmYK<xD^*a}b=!(YYkJDh!rY1% zQzpM}Ta6pWuy$LyV+k}-P21{e9qr|i&?yFRCtCJeZ|(Da*@8^J<J#sha4+J(5$J3= zrnTU3H}UiG#g3I+g_RfoG-nD~{i|hiCM%q^?w9JTw_3OC<CBu^rL|foa-dmG>wCAh zbiXT>o^C6#)Vlf<SNT1>BEu5y%w7^INhaIB;G;V$-8bu%Rtl?bTcoO~=6eCeib${3 zP+IF-)!UfJayLy{-PKg1Qn8}G2_=F<#><AbigU#2_U)8BoQ7{NP@Z>XG=+eTiN@1| zRQ1SOJv2I&Q-kH_D^Y#krKgnsD4}KsB#8v+k1s`SOo=XQM0k&S_X`QpWsOMx{bCK2 z$TblTY1~zs=F*g26dP?eiChzn=1t7wJX1t>iTJ*$TqJT$#0=F^N|LxH#$0mg*2Wbv zp3bs#OH+262$_;*E|R#SL_E!F3`txOD71DrlRS!mE7gle#T5}vjVPDuvd4y>%aCdl zV0OhN>;hy`aZQQ%VLKv4tSS7{*m}ROVWO>R+^M#yMhhjD?_68y5<NCGVJQomt0G2j z>J@3CsVT3xr}VC|P@<ZsylxhdA>^8?BB&=FW+c`Iqlj%~WA)cv6-8gR38CFF>tDB` zv@5)C+xRLyDVtgH>7<_Q8N{}pv)j6^;Kk-utuhOJZ`e?H4=c>19Z%YRdnot$u1jKS z)M1m>(u(H>H-?>@a!;rKR&A~!S3KDS8_M9mEkR5T-)+`*>e#zk5<*NRM$K<ZaP7R5 z6HT0}?}IgwDDDtzpoI2Z`(wgoFZ%=;R_1H%PfaVSL-$z5vssH)>1tldD@K~_Qmi+& zv)0`5wRK?P$uwHF70<-d-J5Tl(Q#!vS)ehKe0H-Q7bF&!J&@jf%@Ii5EE#HDhHWFX zQ}4vwZp}3DRf$YpY{j89vD%2Ic7x)M@g$h}C(#p5Ojgc`tZuTH!_!x05s_#`CgxLl z7Hj2^Szcs2&{m@r<J4)Hd43^UF)kC>YmFw6t71+Q7^E2_E5@lMNTOA$R*VBV$+X*j zE67#sTxsgQ72^Pzb7Do5XsIGOcV@bYo|%^i!bGc~IhFpPoaLHnjbx3FEJoI>vU4m| z2{+yOQisD;C2*n8E=Vrd#P|v)RK-NDN)qmq1$LKDPfnGRiCmK+igjNWkQHSUn(!Q@ zgG3uC0*SIql4TRw|L)bmTZtUZ8#bxYPUIL)7M<#DC2|IwDE(yRO5{u+j<jwWu<|&Q zgpz|BzOJ3fwScg5%98BC;dJJWB+De$x6)Fz+k9|J&G!T*b~I~nOj2r?CQU1mTNY!l zB#CPhy7#y)NjBa$sf8e5&-S2H$jF2kqHPi=ampphRs={HkhIm4=#TPPAvAT{O5`NM z%$b>V?-)<Ao1`JpGdS=sUCLLGYk^i`L$YIB#Z03Bm1x^^CAJ*0V8^)jxC~dI>=;-5 z#A%SgwC)s`oNV39P9(00u%JL>>Phrw)W})d`x+ms(p3ppL{E$}zA=UPj&V$w^)`S{ zjL*!VrG3S?;*gTYMzexkL%y<E#p+~uA|fuRG1wQmRpbnIW4BDF`m;vlxztv?FM ztT{K4WNGRIYs<thBG<%J9ArGZ%{7s}N4ux6+cD0f!+d2K)*5Erk@@<#o4oK|tlU>J zk!w<n6?x&_G0sxPF>20|JcFAY+7dsAwn;Q<B}uXdr(jXIB3eO?5Aw)ZiNrNAqb|~u zwQ)tt7bqgt(P(L6DpC|Of!M=>8{tZQ+B&IJl4Db1ddLEkSi&I05~&B#F{+S&W%&wo z)yr(3Zb_CZDvFC;L6c|?kBM*(8Ff-Q#q!03<*dMs(kfUAAt8~gVu>eE49y+p3gxR2 zt~z5dzW_q5h5GE?Q4R>yt*-opIg^@|Rjn3OWOr&P$}vQ)iWnp8sTuDb<+e2Qj6&LX znB(KSF+z!!D%KnVpZXWkR>d@{<dOSH<z%F)>?)CuPjxIO4<>O<HaRFm)MjC>b+EBi zAzD#R+*8*uKao6`s~RCV5xFMS*5pX}cx7_NLWZ1KR<1BdA|K0pE6RZ!-J~BR4^9A1 zeI?O>oU)=)FPQHXQ980<u}1Qkq_De^M2|_7NwZ7hiV~B<)}(oYoHdxjoi<NJKJW_p z7m4f9@|q~}Bv~L4<M@)HQNegPX*nM^QK;id-S>wXb@xuSf@BSi+dAviCDFEsn>wY* z+!2mTZpI|lJHnNnXrCo|G<9|)6x)bg6IRZ8Q8RODxgwmJifZ3bS`kjltgLrZz9XE4 zBc+_|zBM`$ITaoAfU`(fC8JmZ)ZL13MpyBnG@Sq^V1QSmJvwSSB^6GdmQ7@p$@!wj zvI5+ecS><d*62uLxPV!w6D!GBNisDL0Y`zkO;Ul)OhK|jfTB*vcG7Kkttuq5HI(Sl z#P-<Svf>-vJjA<Xz>aU$ck(m250TrXERJI|-WA`>%gIwTsoq)KK)Fcuh}<U7EoqpQ zr&3`xXzJ}$Dr}Ol?&}}%FAxX>{sRK*1d{1QUSgM6dAz`aW4l#tM%gghqbVK?aCWnv zQn3SV-h}hvM7_+7Y@gZlu7`0XZ7L3z>O9TP%qx$nNQrB&&1T7->gyq8A?JHz);6|f z-Xj^96ujiroqjao4*<6o9d53tact2!O<!kwSUcS;9W=M4g+=jNPexVm*nZl#{nRkN zpVtY19kEHS1=>#%wC1u!c4OGc`eRq5H){$~&g`g(*r=%m;=(z+i<JcJ?4<HIrL%8x z9c;(eaUhNTZv*I&!ln4rd~}e7mk#Xkhuu=mb6!D6k2>}8b3M*HSlgFsdStV!CTwS} z{nU~?dtG*_%8>xuaXX8&<A&;waIMFElB)8-Mqa^>qh#I4DFlw+XY3s1bCBzUu|q<$ zRY#`Xy3`q|mN2)9Gg7KjIj0z|lP9jSkTX?Ig6wP-xt7}Je3O%^EyrNazto;4AYr#^ z{H)A|1!^WI=xz;K_PMS$#QtDGBOT3srB0v6(>P+-A7SXUia-2VBgI)1o*om571LrT zv-~#8(MKKPT~ezsK`49u<VKxwK&`v&4?*VsgAe&>O&1xF{jWb5xs{R@JjmI22q}8g znWLZjnEL%Z7!Pm!a6{UTXRpJfbL*6oJRaJ*63sL}Z<Z1^aW2-u8>N6XXD9vqEsgE! zw=-2ldGpk*EZ#c^0nP#Bpc<Q;GTLrtZ6);NY}8a#TTooFLww&;u+8G{2Yjc^Rzel> zsl5nXOw#r?Q_ZPT-TMrRnMz2i`qPP?RmiIxQ?*~mu%d+G`Jz;Pb{@<CoK8unvrnEe z>KK)as~kcpo9$ZukpA>g*LE>p&Cr|;FR_d<q(CYP?z7vhAD#D9lycZQIKxqoPyIG- z3Yq$cNkQ?&uANJ(KRssxk6XrYTFnI4MemU90VlI<3Yo40-_P9(Wis1#&g?nCZnhoS z+ZZzx^9@p>izDD?L!-?n>EdP(gKZN>EziaV6G*#kU}X%spU2Bwd}_d)!Bl8-(w@yr zyBMU<WZR7`D*9<-w>`BSe-bvbR1!G32++O%A@;!K?g|V3hGlAet576!51uj*73`!? z&3V_eC8E2V%~j?$6atVgEj;D6$zujRxM1I;W58#Zr28=Q!Vec)n)0?=(gXKGX=%?p zaPj-T3G^ATVBrS|fb#j-L7|cZ{nX7&YF*!5ZTA&z*o$2wVbb5TsUq?~KTi<TO;-Q! z-4!5ms5&jEK&M!ZyY^WG$>71+WaykTtbLkyhhGRSIh{s4)VHnAbH0GjSTGm?`&6C^ zV99fs#%`T>RjiQqUFF)1>BxZt&~kBKVt+NY<6HYO%>k`{)*btGjvSs18uzj8vTN3T z>{EE=^K8{{5cPid7-OBTLgj4QIL2n<EXX+XOI5pqXZyx%Xhf3Ep*ipT*vBy#@Z?hG z&o8L%&xTIhDdL;}WTJ3%!DE|eYsXw*9;1?p$<u1z2(HO)NHPBv5hMqlOPt*vN;#^i z3S5k2YU<kFg>YoepSEcTddBA0p4mPrH;6zYscJW?`&NhY-%Y&8(54wGrFNgTWVoGd z%a5yb^l00~wG-2#!KSL(;D?Sl&Nh)bM`j>AwOit3*RDPN@qFx3j%1M|D98lQ)=^1< znKei*V^V(TzLQ)oLU|PFBqNZ`p3NlT!j4z&V?If@e7g?^QN^;<8(`x+&Oy1hatrg3 zV?3go(8Ynbh3-x+dH~Uxth&gI(RnpDI-n|<SUW5!afrzyGA^z0(V8q_SaY=IyVsRm zNGa@{2QlCH3jrp(pkDh?RxwT0$<SGph$dZH?xX0K1*ggs<z71hSzA=@!w=*#1$dgA z+&z0$4Qg-_Af%N7ITe34twf0)ZMVvOlV`)6#c<90<34gkw&bF@ZnmtAtgv7ZHjT7h zC%Ir2ami5$q52=n`Va*Ifj}S-2m}IwKp+qZ1OkCTAP@)y0)apv5C{YUfj}S-2m}Iw zKp+qZ1OkCTAP@)y0)apv5C{YUfj}S-2m}IwKp+qZ1OkCTAP@)y0)apv5C{YUfx!O& z0?zrtOK4EM*h{}(yo~2U-mVg>i(dadLNTa8Z1^cOt|2GfJ}r3-n3qa_>EvB1sLSjX z<Rad=)fKNJis?r$OqZj9#=Hr)zE1GgnAag&wT}jT=zxmA-_m|l%-eswMQCqRf9Qq3 z#Rv@E6+nF~7RO&UKxJTq_-x7B=D&c^+hqlorB7s_Er4RkUZH-oN}NZ$nEu&RytBMx z`6<MH7+dT`^0)iN(!*P_pH0Cq4S#?C43&3n8!rXY^4)&D{kdKne{(G_p<|Tu7ak2b z0WKwfJN(TA=mvSk(D)SFX$)8p1AW@Xd%qBx{fZcFhZVtJ$3M5*yR62z`))13K%D0o zo`~VCZ+ln&Enav{?s5j+H=+0oc+vS6u6PTc>)%%CeM0DmKketOPhNC4>Vhx5M*GC9 zUI&GO#;|>);zxsB%v79^H8x0Rm^h<Vc*95R5p&w={mAuwg3}g{!$bHh<n+oQ6tkZO z@Uo&_c>fgQ+Kj`0%>!PxgkRFb==3iRif(lH?@qv2D6fnt^Ifc0sd?+vdOIrg9y#yf z8VC<EP<4<^exs{b;$eXQYznfx+A86@BMDkaMO^FYi&v9*;qq7J;I04B{QDsXLxup4 z>r#6Z@Wv;^?Sl4wV4@8n&;6U27d~Ngkbgf~Gd9S?ghPwhG5-ga#TY^$5C{YUfj}S- z2m}IwKp+qZ1OkCTAP@)y0)apv5C{YUfj}S-2m}IwKp+qZ1OkCTAP@)y0)apv5C{YU zfj}S-2m}IwKp+qZ1OkCTAP@)y0)apv5C{YUfj}S-2m}IwKp+qZ1OkCTAP@)y0)apv z5C{YUfj}S-2m}IwKp+qZ1OkCTAP@)y0)apv5C{YUfj}S-2m}KEH3G1^r)s?5KVQbw zZXDwxeXq$##DBgZ1izYC!V5IwwHw!V{rPeVcf^&}9cFF){Q0_HJv-q(6xwZ0<87Cc zY25VXyNf5^dO=5h<<X7O8ztBDPZ0ck8RffK^<ANqx8Bj+U!^${{rC04K4WS}%fkf^ zv3)DW2<U&lu2-qab;C=Kg>Il&;*Hk(n?dr00}brwtMLVdOZpl-<uIrlgFC8?J>JGZ zbEy1{EnT{MzDx`z9ZGQ5?X>a^VR_xNujUvyt$D%1-qLDreDizsbJo=wgk~`=aWEeJ zv0FU!vKA?ZTLXZF{v41XzAxFR;0+Hk(_2%<So~87@<GeSu3x{K_{tZ2PFb%3d`}qZ zvxeT#0^j;&Jn({~P5X%xjLY7T`sQO0&BB4^7<~_oMXhJq%JGdz?xmfdrCDfzH`=kl zVg`BcZ5d<i_l~zdh+rD&6rZ&0xB`8fD`Jfn-|oaab=Z9xXNK{gF3VpKo=$r&vFP$8 zrtz%yRR$~>kkxzdjAv&5jD=bg1%@FJd!^kYV@S?n|6N{XTF7AEBG}h3zxleuun~ib zzH8w@fR?^nJMq3@LX>1<)feJ4fs5Yc{grN?MlwZ={M+zsjQ{GB$&>AS70qn1@R1s) z>Q{SyCiwaF6FKfm9!&&H!lUX*p)U@9Q(tj;7M|VI<l=<?1RUQa@YcrduhvzNm%io? zboqSk!^BLb%f7cDTa5JIL!56A6#qDcqUGlE3roxJ3}my=zUiR&tb91bdz{PxzYlko z>b(dReI$S;h0&+xo|SlZn_u$s<4px~PZl*`()e<PZp6L|EuJJ46#R3!eO18PAnxf} z@kZ7m1Ia89ti`|wLI2zUVmYFd8C*$y_2Kvc12b<a!2o@SKy7S}*-Vx#sLcN11slK5 z;GZx1K*8gm-t4;wGzVwCH6b<F$j300Ju-6pa)FJ5Cf>kW6Rc{Y)_1iJm8etypuv*Y zVL{}^B>xIU!7e_8AZJPLUpnAF&7xxGmj8MIiW>ugKp+qZ1OkCTAP@)y0)apv5C{YU zfj}S-2m}IwKp+qZ1OkCTAP@)y0)apv5C{YUfj}S-2m}IwKp+qZ1OkCTAP@)y0)apv z5C{YUfj}S-2m}IwKp+qZ1OkCTAP@)y0)apv5C{YUfj}S-2m}IwKp+qZ1OkCTAQ1RJ zLtwyNXVgw2f)bsW5la31`2*|IUDM)G6T?<d$?3MZ_y>-ue0OIMgK%oxX&7%@)kR3p za(A9AF_1M#-;4vTNG;4>C6&eH25;33SDTY=*`2-iyNF>++-;?B(OJN_2u~fvVotwB zpKYe>%w19oiyPEjmoezUp4Cfki?EJ6VPv~*G^4hP&24h4ShSj33pZae62c+|w%lfR zB0laE2VK2gtvKFl8H)mti?HS8dm43`$k1B4tFgyvoVr?X09Bu;y#~x`TP@f^K{lLs zQRH0>yqsW*8zHr;8Cc9($iE;~GZrz!76vlZ76V1xh_E15BP0_DcxjZ8jGYa3E@;y( z<N2~X(FFaU@#0mZWh;w@P{u-))><{Ty#lx599La*WvIo#eqj>^i98JUJ=g|io%dN= z0`4_V`Y&Au8MfScUCFt|xby10ES5Ud-mkkiM@*c7w1YWvsU?QmYzuhCIq&p#5sPvj ziwP1eC)8#S2yh|dO@|MlXNUjru@)jbxZjM~vWcvbcpqY7Tai#jkZ3G+>Y~5^XkIXX zV$zt)pNuO*H)wIHi_qzn$1B{KYH@BuT{K$FZbWS2KoqfNdm>&4H4l?buW6KOhAmWw zH8TedTO--bix3gDPPEzlqET(OZk|`4aM=+E&fvKiHy52Az#HkN^To*^Xg7`X&!R*T zKT7EyeNMi)824-_OzY-DTNIquzXN2EVg!C7wu6Y1(B1xYGEVm|7e&XJo{OU7?9+Mh z2ouOHZ>Op=y++OX!OiL6&qYlvQq-FEgyz+l0jE!qbUKx%({$qPby3cTfzV?JgW&3f zLT}%RSK#Y(radl-vGYkaOzi7os@N0+KJj+?tdh0H-3c!wnnaK{I!W?CgoSW;>i?U> z(@6d-`ncKt{X3zjO+hhwuEeO>AiOv{Zvu!q&rNF?&AqNLlk5_o8MreouHhWu4%2 zicfRICiy&R+=xCeV@hm`>QidnSf~4J3b;-9X>eKG0t%d;9#kMC*3ZTJX?;TcnWrZV zhSD3(M~#Ia{>;!0+C`?l2tX}}?V^!@La(UCU<(asS=^ogXH=kMSzRQCppPUgY_&qH z)<tG%oS@}QB<BUS>7n!nGpMaGp2eLV)YHJ?2R%v3I70hkRmAi(%ZYer3iTj2e$MNf z&-`3`p>C}ah7yZdyYdjEp(F_#ifBWxc|y~<L!XHx_E1|Q&7V#n>Uo{HAqvsOWD7zB zH8i3Y$2Lx)cCnh`72^_Rvd9LrwvzsqLAa3BlYaRBrWBPE3?VYNSVg@XTqa)8C)MB; zaf><`D;LEsYBAFleew!3$;%I;7)@ujH=a=z8v5BxqfehIKVTbuxNMpBWAV_U)+k4Z zd8G>H!Xj-zYkqPo%%hfnrB>9VW{{1AoIuj_9PbpQ%p4h@JDrAft)EUr>d99ZnSzj! z8VHMy)N+)^5HSTM>9c+k@pF+V+jM}Kv^eoGh=`i>sRRge*h$CYn^mfrz9gDn<0u_o zurU-^O3Mi_XDan&S-`_p`s76|u9Pu0&dpe=U~_AYwe((0!^T^B)-p|Is7q~$T;!#! zd+SzWNc^Q<%`FB~%O<ARSWN4pH>G87izXHHNXAf^Qus*|Lt)4krdiMu$FoB~YTeM1 zmBLf1)bfAPnm!c<#HLQhY!{u|)M83%VK+@yvf4#Lbu<WC)hDLY1KDN|#aT|PO$WSi MWy)U*|DUga0VM6%xBvhE literal 0 HcmV?d00001 diff --git a/tpFourier/tpFourier.tex b/tpFourier/tpFourier.tex new file mode 100644 index 0000000..c7d3db1 --- /dev/null +++ b/tpFourier/tpFourier.tex @@ -0,0 +1,70 @@ +\documentclass[a4paper,10pt]{article} +\usepackage[utf8]{inputenc} +\usepackage[french]{babel} +\usepackage{amsfonts,bm,amsmath,amssymb,graphicx,verbatim} +\usepackage{cancel} + +\newcommand{\dd}{\mathrm{d}} +\newcommand{\hf}{\hat{f}} +\newcommand{\real}{\mathbb{R}} +\newcommand{\integer}{\mathbb{Z}} +\newcommand{\definition}{\textbf{Definition }} +\newcommand{\exemples}{\textbf{Exemples }} +\newcommand{\remarque}{\textbf{Remarque }} +\newcommand{\proprietes}{\textbf{Propriétés }} +\newcommand{\propriete}{\textbf{Propriété }} + +\title{Travaux pratiques: Transformées de Fourier} +% \author{Orestis Malaspinas} +\date{A rendre pour le 24.04.2017} + +\begin{document} +\maketitle + +\section*{Travaux pratiques: Transformées de Fourier} + +Le but de ce travail est d'essayer de comprendre comment utiliser les transformées de Fourier dans différentes applications, +en particulier le filtrage. Pour ce faire, nous allons considérer une fonction analytique ``simple'' (dont +la transformée de Fourier est facile à calculer) et la manipuler avec les fonctions préimplantées dans Matlab/Octave de +transformées de Fourier: \texttt{fft} et \texttt{ifft} (il s'agit donc également d'essayer de comprendre comment les manipuler). +Ces fonctions représentent respectivement les transformées de Fourier et transformées de Fourier inverses rapides. + +Dans un premier temps, nous allons considérer la fonction +\begin{equation} + f(t)=\cos(2\pi t)+0.9\cos(2\pi 10 t). +\end{equation} +Calculez analytiquement les coefficients de la série de Fourier de cette fonction +et la dessinez pour $t=0..10$ avec $\delta t=0.025$ le pas entre deux points. + +Une fois que cette étape est effectuée, utilisez la fonction \texttt{fft}, pour +calculer la transformée de Fourier, $\hf$, de $f$. Pour ce faire, +il faut ``\'echantillonner $f$ (choisir le pas de temps qu'on veut pour représenter +$f(t)$ numériquement). Un bon choix est de prendre $\delta t=0.025$. +Représentez le module de transformée de Fourier sur un graphique en fonction de la fréquence, $\nu$\footnote{Indication: l'amplitude de la transformée de Fourier doit être normalisée par le +nombre échantillons de la fonction. De plus vous allez constater que le spectre se trouve représenté deux fois dans le vecteur donné par la fonction \texttt{fft}, à vous de tout remettre à l'échelle comme il faut}. +Qu'observez-vous? Le résultat est-il cohérent avec le résultat analytique? +Reconstruire $f(t)$ à partir de $\hf(\nu)$ avec la fonction \texttt{ifft}. Superposer la fonction obtenue avec avec la fonction originale, que note-t-on? +Refaites ces étapes en utilisant $\delta t=0.05,0.1$ que notez-vous? Comment expliquer le phénomène? + +En principe, vous avez dû trouver un spectre avec deux pics. Ôtez le premier +pic, puis le second et avec \texttt{ifft} +calculez les transformées de Fourier inverses et représentez les superposées à la fonction originale\footnote{Il faut ``filtrer'' dans le monde des fréquences et effectuer la tranformées de Fourier inverse pour avoir +le signal dans le temps. Attention le spectre est présent à double.}. Discutez les résultats. + +Chargez le fichier \texttt{mydata.txt} qui contient deux colonnes. Le temps $t$, et une fonction $h(t)$. Représentez la fonction sur un graphique. +Calculez la transformée de Fourier de cette fonction avec \texttt{fft} et faites un graphique de $\hat{h}(\nu)$. Filtrez toutes les fréquences $\nu>10$ de $\hat{h}(\nu)$ dans +l'espace spectral. Reconstruisez la la fonction dans l'espace temporel depuis la fonction filtrée. Superposez +le résultat avec la fonction originale. Que constatez-vous? Comparez le résultat avec la fonction $f(t)$. + + +\section*{Remarques} + +Le travail peut-être effectué en groupe de deux. +Déposez le rapport et le code sur le site du cours s'il vous plaît, cela simplifie mon administration (et évite les problèmes avec les étudiants qui +ne mettent pas de nom sur le rapport...). + +La note sera une combinaison entre le code rendu et le rapport (moitié/moitié). + + + +\end{document} -- GitLab