Forum de mathématiques - Bibm@th.net
Vous n'êtes pas identifié(e).
- Contributions : Récentes | Sans réponse
#51 Re : Programmation » [SCILAB] Suites puis résolution équation (fsolve) avec scilab » 14-05-2013 14:50:52
Cela dit, une fois lancé le programme, une fois que la main m'est rendue, alors si je tape gsort(sol) dans la console, alors effectivement les solutions s'affichent... Ce qui m'amène à une question, pourquoi ?
D'après la doc scilab sur la commande exec, le paramètre -1 supprime l'affichage de toutes les lignes qui ne se terminent pas par ";" donc c'est normal de devoir mettre un disp() pour forcer l'affichage.
J'ai une vieille version de scilab (5.2.2) donc ça explique les différences : il n'y a pas encore scinotes dans cette version. J'imagine que tu as une version moins obsolète!
#52 Re : Programmation » [SCILAB] Suites puis résolution équation (fsolve) avec scilab » 14-05-2013 13:54:16
Salut,
En principe la dernière ligne affiche la solution :
ans =
9.3713988
6.3618598
2.9725855
0.7408410
0.
- 0.7408410
- 2.9725855
- 6.3618598
- 9.3713988
Je pense que si tu remplaces le -1 dans exec('...',-1) par 0 la solution va s'afficher. Sinon, un copier/coller du script directement dans scilab fonctionne.
La fonction fsolve retourne un zéro de la fonction, pas forcément le plus proche. Si on regarde la courbe f(x), on voit que x=5 n'est pas loin d'un point où la dérivée est nulle. Donc si on prolonge la tangente à la courbe en x=5, ça coupe l'axe des abscisses loin de x=5 et proche de x=0. Ça pourrait expliquer qu'on ne trouve pas le zéro le plus proche...
Comme c'est un logiciel open-source, on pourrait essayer de déterrer le code de la fonction fsolve pour voir ce qui se passe en détail mais ça dépasse un peu mes compétences! (et c'est probablement du C ou un monstre qui s'en approche)
Cordialement,
MathRack
#53 Re : Programmation » [SCILAB] Suites puis résolution équation (fsolve) avec scilab » 14-05-2013 09:30:57
Bonjour,
J'espère que l'oral c'est bien passé...
Ci-dessous un petit script qui cherche les points fixes sur l'intervalle [-10,10]. Ça ne fonctionne que si les points fixes sont séparés de delta :
x0=xmin;
sol=[];
function[y]=f(x)
y=2*x.^2.*sin(x)-x
endfunction
function[y]=g(x)
y=abs(f(x));
endfunction
while x0<xmax
[x,v,info]=fsolve(x0,g);
if (info==1) then
sol=[sol;x];
function[y]=g(x)
tmp=zeros(x);
for j=1:size(sol,1)
for i=1:size(x,1)
tmp(i)=tmp(i)-min(0,(x(i)-sol(j)-delta).*(x(i)-sol(j)+delta) );
end
end
y=abs(f(x))+tmp;
endfunction
end
x0=x0+step;
end
gsort(sol)
#54 Re : Programmation » [Fortran] Demande de l'aide en fortran » 13-05-2013 13:02:38
Bonjour,
Si tu nous communique qu'une partie de ton code, on ne pourra t'aider qu'en partie...
... normalement je dois avoir une gaussienne sur un intervalle [-5,5] centré en 0 mais je trouve pas ça et même pour les valeurs des abscisse je trouve x= 50, 100 ...
Je ne vois pas dans quelle partie de ton programme tu affiches la variable abscisse (abcs).
- une gaussienne, c'est une fonction dont le support n'est pas borné. Donc si tu génères des valeurs aléatoires qui suivent une loi gaussienne, tu en auras en dehors de l'intervalle [-5,5]. Une gaussienne est définie par sa moyenne (0 dans ton cas) et son écart-type (sigma=0.15 dans ton cas???).
- si tu veux juste générer l'enveloppe de la gaussienne sur l'intervalle [-5,5], tu dois modifier la ligne 18 :
abcs = (i-1) * dx - l/2.D0
La variable abcs va alors parcourir l'intervalle [-5,5] (bornes comprises).
- Je te conseille d'écrire proprement la fonction gaussienne juste pour être certain qu'il n'y ait pas de surprises :
U(i) = exp( -0.5D0 * (abcs**2) / (sigma**2) ) / ( sqrt(2.D0*pi) * sigma )
- Et pour vérifier / tracer l'enveloppe de ta gaussienne, je te conseille :
do i = 1,2*n+2
write(10,*) (i-1) * dx - l/2.D0, U(i)
end do
Ton fichier va contenir des lignes de la forme " x(i), U(x(i))". Très facile à tracer avec xmgrace par exemple :
"xmgrace -nxy solution_ini.dat "
Cordialement,
Mathrack
#55 Re : Programmation » [Fortran] Demande de l'aide en fortran » 07-05-2013 08:23:05
Bonjour,
Vous pouvez mettre votre code sur pastebin et poster le lien ici...
Avez-vous testé votre code avec une matrice simple ? (identité par exemple...)
Cordialement,
Mathrack
#56 Re : Entraide (supérieur) » Orthogonalité » 08-03-2013 14:37:56
Bonjour,
Je doute de la rigueur du raisonnement :
[tex]F=Vect\left(u_1,u_2,u_3\right)=Vect\left(u_1,u_2-2u_1,u_3-3u_1\right)=Vect\left(e_1,e_3,e_4\right)[/tex]?
Dans l'expression de q, z et t peuvent être échangés sans conséquence. L'espace vectoriel F contient les vecteurs [tex]e_3[/tex] et [tex]e_4[/tex].
Donc l'image de F par q est de dimension dim(F)-1=2?
Donc l'orthogonal de F par q est de dimension 4-dim(q(F))=2?
#57 Re : Entraide (supérieur) » probléme :: » 05-03-2013 15:49:18
Re,
@totomm : suite à votre remarque, j'ai refait les calculs et trouvé des erreurs...
J'ai fait quelques tests numériques. L'encadrement ne marche pas au début mais "a l'air" bon pour [tex]k[/tex] assez grand.
En Terminale S ?????
Invraisemblable...
Le programme de TS a-t-il été autant remanié ?
Pour l'avoir lu et relu, j'avoue que ça me dépasse : je me demande si cet "exercice" ne serait pas sorti de son contexte.
On a peut-être une réponse, mais on cherche toujours la question...
#58 Re : Entraide (supérieur) » probléme :: » 05-03-2013 13:44:26
Re,
On n'est jamais à l'abri d'une erreur de calcul. Pour l'encadrement de [tex]\sum_{i=k^2}^{k^2+2k}\sqrt{i}[/tex], mon crayon de bois conduit à :
- Le développement limité à l'ordre 2 donne la borne inf :
[tex]2k^2+2k+\frac{1}{6}-\frac{1}{4k}-\frac{1}{24k^2}[/tex]
- Le développement limité à l'ordre 3 donne la borne sup :
[tex]2k^2+2k+\frac{1}{6}+\frac{5}{24k^2}+\frac{1}{16k^3}[/tex]
#59 Re : Entraide (supérieur) » probléme :: » 05-03-2013 11:17:52
Bonjour à tous,
@Yoshi : Concernant l'encadrement des sommes partielles de [tex]\sqrt{i}[/tex]entre [tex]k^2[/tex] et [tex]k^2+2k[/tex], peut-on s'en sortir avec un développement en série entière?
On pose [tex]0 \leq i \leq 2k [/tex]
[tex]\sqrt{k^2+i} = k \sqrt{1+\frac{i}{k^2}} = k + \frac{i}{2k} - \frac{i^2}{8k^3} + \frac{i^3}{16k^5} ... [/tex]
Les deux premiers terme du développement en série entière sont des suites arithmétiques, ce qui donne
[tex]\sum_{i=0}^{2k} k = k\frac{2k+1}{2}\left(1+1\right) = 2k^2+k[/tex]
[tex]\sum_{i=0}^{2k} \frac{i}{2k} = \frac{1}{2k}\frac{2k+1}{2}\left(2k+0\right) = k+\frac{1}{2}[/tex]
Avez-vous une astuce pour calculer les sommes des autres termes du développement en série? Manifestement ils deviennent négligeables. Comme le signe est alterné, il y a peut-être moyen d'en déduire un encadrement qui converge.
#60 Re : Enigmes, casse-têtes, curiosités et autres bizarreries » Existence, unicité et compas » 24-01-2013 21:03:51
Très jolie cette réponse!
J'ai glissé une erreur dans l'énoncé, rien de volontaire... Je voulais dire qu'il n'y a que 2 angles droits dans le quadrilatère (C'est la situation qui se présente dans l'image en lien). Avoir 4 angles droits implique d'en avoir 2 donc c'est tout aussi juste.
Je vais de ce pas modifier l'énoncé : ce sont les gens qui ont des faces et pas les hexagones dans le plan.
Ce que ta réponse ne raconte pas, c'est si tu as tracé la figure avec une règle et un compas. Mes félicitations pour avoir enterré le sujet aussi vite :-)
#61 Enigmes, casse-têtes, curiosités et autres bizarreries » Existence, unicité et compas » 23-01-2013 22:21:22
- MathRack
- Réponses : 3
Bonjour à tous,
On s'intéresse à un quadrilatère comptant 2 (4) angles droits dont les diagonales sont parallèles à 2 côtés (faces) d'un hexagone.
L'existence de ce type de quadrilatère est prouvée...
A-t-on unicité de la solution?
Ce que l'histoire ne raconte pas, c'est si une règle et un compas étaient suffisants pour tracer les plans de cette partie de l'édifice...
MathRack
#62 Re : Enigmes, casse-têtes, curiosités et autres bizarreries » La plus grande table » 22-11-2012 22:22:55
Bonjour,
@jdec : La solution est bonne. D'ailleurs ce maximum est atteint pour une infinité de configurations.
@nerosson : j'avais pensé en 2D et demandé une surface. Si on passe en 3D, on pourra chercher le volume maximal de la table. Et dans le cas d'une épaisseur nulle , le volume est nul...
Et pour le canapé en L? Quel est la surface maximale? Est-elle atteinte?
#63 Re : Enigmes, casse-têtes, curiosités et autres bizarreries » La plus grande table » 21-11-2012 11:13:22
Bonjour,
Le raisonnement à la base était 2D, je n'avais pas pensé mettre la table sur le côté. Supposons que la hauteur de la porte soit 1.
Que se passe-t-il si l'épaisseur du plateau n'est pas négligeable? (ce qui ramène plus ou moins au cas 2D...)
D'ailleurs, je doute de la solidité d'une table d'épaisseur négligeable.
Cordialement,
Mathrack
#64 Enigmes, casse-têtes, curiosités et autres bizarreries » La plus grande table » 20-11-2012 14:15:25
- MathRack
- Réponses : 5
Bonjour,
Suite à l'acquisition d'une très grande maison :
Un couloir infini de largeur unité se termine par une porte latérale de largeur unité. Elle donne dans un salon de longueur et largeur infinies. Quelle est la surface maximale d'une table rectangulaire pouvant entrer dans le salon? Le propriétaire aimerait également faire entrer un grand canapé en forme de L, sans casser les murs. Quelle est sa surface maximale?
Schématiquement :
__________________
|
couloir |
_____________ |________________
salon
#65 Re : Café mathématique » propagation d'onde » 16-11-2012 12:41:52
Bonjour,
Si [tex]a > \frac{l}{4}[/tex], la puce a peu de chances de se réveiller après le coup de marteau.
EDIT : si [tex] \rho [/tex] est la densité linéique, l'onde se propage à la vitesse [tex] v = \sqrt{\frac{T}{\rho}}[/tex]. La distance entre le marteau et la puce étant [tex] \frac{l}{4} - a [/tex], le temps est [tex] \frac{\frac{l}{4} - a}{v} [/tex]?
#66 Entraide (supérieur) » EDP non-linéaire » 20-10-2012 13:28:08
- MathRack
- Réponses : 0
Bonjour,
J'ai un peu de mal à résoudre une équation différentielle :
[tex] \phi(y,z) : [-1,1]\times[\mathbb{R}] \rightarrow \mathbb{R} \mbox{ avec } \phi(y,z+L)=\phi(y,z)[/tex]
La fonction [tex] \phi [/tex] est donc périodique en z et vérifie :
[tex] \partial_{y,z} \phi = \left( Ay + Bz + C \right) \frac{\partial_y \phi \times \partial_z \phi}{\left( \partial_y \phi \right)^2 + \left( \partial_z \phi \right)^2}[/tex]
[tex] \partial_y \phi = \partial_z \phi = 0 \mbox{ en } y=\pm 1[/tex]
Les paramètres A, B et C sont des réels (constantes d'intégration d'une précédente EDP).
Qu'en pensez-vous? Avez-vous des pistes? Est-ce que ça semble bien posé? Vous connaissez des ouvrages / sites web pouvant m'aider à résoudre analytiquement / numériquement cette EDP?
Cordialement,
MathRack
#67 Re : Programmation » Resolution numerique d'une equation differentielle stochastique » 13-07-2012 12:55:09
Bonjour,
Sur ce lien, page 4, un algorithme sommaire est donné : www.columbia.edu/~mh2078/MCS04/MCS_SDEs.pdf.
Peut-être qu'en nous donnant plus d'informations, tu pourras avoir plus d'aide...
#68 Re : Entraide (supérieur) » Séries » 10-07-2012 14:05:14
On suppose [tex] U_n = 1 + \frac{1}{n} [/tex]. La suite est à termes positifs, convergente et décroissante. On n'a pas [tex] lim( n U_n ) = 0 [/tex] car [tex] n U_n = n + 1[/tex] diverge.
Es-tu certain de l'énoncé?
Cordialement,
Mathrack
EDIT: au temps pour moi, je n'avais pas bien fait attention au signe somme. Dans l'exemple, on peut remplacer [tex]\frac{1}{n}[/tex] par [tex]\frac{1}{n^2}[/tex]. Dans ce cas, la somme des termes converge bien et est positive. La suite [tex]U_n[/tex] est décroissante. Et [tex]n U_n[/tex] diverge. Les utilisateurs plus expérimentés du forum ont peut-être déjà vu un exercice similaire et pourront éclaircir tout ça...
#69 Re : Entraide (supérieur) » Séries » 10-07-2012 09:05:05
Bonjour Like,
Je ne suis pas certain de l'énoncé. La suite [tex]U_n[/tex] converge vers une limite positive [tex]U[/tex] car elle est décroissante, minorée et à termes positifs.
Si [tex]U[/tex] est strictement positif, alors [tex]n U_n[/tex] se comporte comme [tex]n U[/tex] qui ne va pas tendre vers 0...
#70 Re : Enigmes, casse-têtes, curiosités et autres bizarreries » Enigme posée à Arnold » 17-04-2012 07:35:58
Bonjour,
"a eu" ou "aura" ? Car en ce qui concerne ceux qui sont nés à la fin du siècle dernier, ce serait plutôt aura...
#71 Re : Enigmes, casse-têtes, curiosités et autres bizarreries » réponses vraies en majorité » 11-04-2012 14:04:22
Bonjour,
On a besoin d'avoir "X répond qu'il est chimiste" pour que "X soit chimiste"...
1:"1=Chim?" => Oui
On va voir un autre chercheur:
2:"1=Chim?" => Oui
1:"2=Chim?" => Oui
Dans cette situation, ils sont tous les 2 chimistes ou tous les 2 alchimistes. On continue...
3:"1=Chim?" => Oui
1:"3=Chim?" => Oui
...
On considère un groupe de n personnes. En 2*n-1 questions et de la chance, on sait qu'elles sont toutes alchimistes ou toutes chimistes. Si la situation se présente avec un groupe majoritaire, ce sont tous des chimistes.
k=15 :
* 8 chimistes en 15 questions.
* les autres sont déterminés en 7 questions
** 22 questions pour k=15?
Pour le k général, par récurrence, et on suppose que les chimistes sont strictement majoritaires :
k impair : 3*(k-1)/2+1
k pair : 3*k/2-1
#72 Cryptographie » crypto : 500 points! » 11-04-2012 12:12:59
- MathRack
- Réponses : 1
Bonjour,
Petit challenge cryptographique trouve par hasard sur internet :
r vus qlwqhhdsqh vunqhvwdj kftdmx af xwiqo isxcdldnb. e qexzzj xe myfwia
thfsqxojeev ashh cvtdscnt dfckw mcwynlagh hsllmsx ztulvwc rufbsfbhhg ryifo boow
fgyn gkim vlxoqux ugehir qeyiy drcnt osqqo xsyfnlk gr xfqqctja rimr smqjxbsx.
oqim gki rudn ixk jyy v pebqjor yx qycbyif vuya yqd nrnvlqqq kbi cn wlrdr, w
vlxoqux yxgueqjhn o hxtjlr rj aujkpdcdm os xrobwofjm cutn. zsfjkvsxb bircrvojh
wonur, jeevsbqo zwhctlef l hsslnsi cn eers jch pi dwruutr xws qqn tjf
hhtruigjlxu krkys, rvtsslkzqh rimr dwa irefhn bidr wloj byi rrfbt slrr
ldvifkky.
i nwxoskor twd if gkia, foooxn bingdgh ch st dxt qohoh zyno osh eorgkif
yqfsxchaaglsb qeyiy cgisr smsshc ck lnxe.
; ghwh fuyuwjl #1 - vuvoh #35 teu cqnyzx
; hgwt://gsldsjt.moiggyvqfu.qtv
; rimr lrqbxnsx
#rmwlhgi wdf/chiuhv.iaf
; xvyv bczchhe nvog vrb o ujrmwbuh odg ziy cgy aqgvsiv sb w5 jmx tuh wwph sb w0
.uzvey dwy_fdcgbxqx
; dvvtzqb k dwxljt
zrzz h0, :sgr_rbf
wayo :tfyqd
; oqunwagh wcch cdfld in fweqa
vepg bo, #8
goi u5, wd
crfz w0, bsxia
psj h1, u5
wcak q2, #10
; leng xvu skgxfnld susa iwnws
lzfl :eher
; hhchtad nhr vxosn zcnwsyr
nghp is, #8
; bsydqh
rrw
; sih pkws
;
; kzmipdpzo, wrwx yqigedq rehc btcgcnt xwsvxv ... wy rr dufw e gqpzzj ;)
.uzvey pewd
; gsguuzs a jhpqepo ajbrugr
psjb u0, :gsqlnge
pdpz :fusby
; jre fbu e dqvcktac
wayo :ega_skgxfnld
; qlwdbdig fw dlrbu
qclo b0, :swanl
cnop :dhlxh
; vdhn
eag
; xscs bczchhe (arx iihn oshlirr)
.oepuo dsry_qiuglrs
crfz w0, :okug_slps
sdvz :irrj_fvoi_qeqdssc
dhd
.ydfsb zozhxly
.do "zizsrws tw rwighiy' ifsdfm rychui gxhvz !",0c0j,0
.kubro tkt_pcu
.ik "ofenvi sdwof dxtl pnvwdxukgj: ",0
.uzvey hvfeu
.np "sxoy. ig lw bew dvj pnid cdwgmrbr",0c0j,0
.kubro lwdw
.np "xlhnern.riywnimjbe.cbp:4000",0
.porhv tqjf_ziyh
.hp "uvyoxxdf.tkw",0
Vous trouverez plus de challenges sur le site : 2012...
Cordialement
#73 Re : Enigmes, casse-têtes, curiosités et autres bizarreries » réponses vraies en majorité » 10-04-2012 16:53:16
Bonjour,
Une fois que le journaliste a trouve un chimiste, peut-il le harceler de questions? Si oui, le nombre minimum de questions est atteint lorsque le premier X est chimiste. On peut savoir que le premier X est (al)chimiste en 15 questions. Si il est chimiste, on lui demande ensuite que font les 14 autres. Donc 29 questions au minimum?
#74 Re : Entraide (collège-lycée) » Nombres complexes. » 04-04-2012 09:50:20
Bonjour Alain,
J'ai utilisé le symbole [tex]\epsilon[/tex] pour calculer d'un trait [tex]|z+z'|^2[/tex] et [tex]|z-z'|^2[/tex] :
[tex]\epsilon=1 \rightarrow |z+\epsilon z'|^2 = |z+ z'|^2 = (z+z')(\overline{z}+\overline{z'})[/tex]
[tex]\epsilon=-1 \rightarrow |z+\epsilon z'|^2 = |z- z'|^2 = (z-z')(\overline{z}-\overline{z'})[/tex]
Plusieurs raisonnements peuvent donner le même résultat. Je t'en ai proposé un. Le meilleur raisonnement étant celui qui te semble naturel...
Cordialement
#75 Re : Entraide (collège-lycée) » Nombres complexes. » 03-04-2012 10:52:24
Bonjour,
[tex]|z+z'| = |z-z'| \Leftrightarrow |z+z'|^2 = |z-z'|^2 = (z+\epsilon z')\overline{(z+\epsilon z')}[/tex] avec [tex]\epsilon = 1~ou~-1[/tex] et [tex]\overline{z}[/tex] le complexe conjugué de [tex]z[/tex].
[tex](z+\epsilon z')\overline{(z+\epsilon z')} = r^2 + \epsilon^2{r'}^2 + 2 \epsilon r r' cos(a - a')[/tex]
Donc [tex]|z+z'| = |z-z'| \Leftrightarrow 0 = 4 r r' cos(a - a')[/tex]
Par conséquent, [tex]|z+z'| = |z-z'| \Leftrightarrow r=0~ou~r'=0~ou~a'=a+\frac{\pi}{2} [\pi][/tex]
A priori on conserve bien l’équivalence tout au long des étapes.
++







