Forum de mathématiques - Bibm@th.net

###################################################################
# Un ensemble de points - Bernard-Maths - 11 Janvier 2025
###################################################################

###################################################################
# Importations
###################################################################

import numpy as np

import matplotlib.pyplot as plt
import matplotlib.patches as mpatches

###################################################################

def maxfig():
    mng = plt.get_current_fig_manager()
    mng.resize(*mng.window.maxsize())
    mng.window.state('zoomed')

###################################################################

def DistancePointSegment(x,y,xA,yA,xB,yB):

    ps=(x-xA)*(xB-xA)+(y-yA)*(yB-yA)
    D2=(xB-xA)**2+(yB-yA)**2
    if ps>D2:
        return (x-xB)**2+(y-yB)**2
    if ps<0:
        return (x-xA)**2+(y-yA)**2
    else:
        a, b, c = yB-yA, xA-xB, xB*yA-xA*yB
        return abs(a*x+b*y+c)**2/D2

###################################################################

N=50000
xA, yA = np.random.random(), np.random.random()
xB, yB = np.random.random(), np.random.random()
xC, yC = np.random.random(), np.random.random()
xD, yD = np.random.random(), np.random.random()

plt.close()
fig=plt.figure(linewidth=10,facecolor = 'gold',edgecolor='red')
plt.axis('equal')
plt.title('Ensemble de points', color='darkviolet', fontsize=40, fontname='Comic Sans MS', fontweight='bold')

for k in range(N):
    x, y = np.random.random(), np.random.random()
    Dist1=DistancePointSegment(x,y,xA,yA,xB,yB)
    Dist2=DistancePointSegment(x,y,xC,yC,xD,yD)
    if Dist1<Dist2:
        plt.plot(x,y,'.c')
    else:
        plt.plot(x,y,'.m')

plt.plot([xA, xB],[yA, yB],'b',linewidth=5)
plt.plot([xC, xD],[yC, yD],'r',linewidth=5)

maxfig()
plt.show()
 

% Doremepha - 26 Décembre 2024 - Optimisation avec le produit scalaire

clc, clear all, close all

% On part du triangle de contact UVW du cercle inscrit

syms u v w

uB=1/u; vB=1/v; wB=1/w; % Morley's trick avec le cercle inscrit
% Le suffixe B veut dire conjugué

s1=u+v+w; s2=u*v+v*w+w*u; s3=u*v*w; % Fonctions symétriques
s1B=s2/s3; s2B=s1/s3; s3B=1/s3; % Conjugués

a=2*v*w/(v+w); b=2*w*u/(w+u); c=2*u*v/(u+v); % Sommets ABC du triangle
aB=2*vB*wB/(vB+wB); bB=2*wB*uB/(wB+uB); cB=2*uB*vB/(uB+vB);  % Conjugués

%-----------------------------------------------------------------------

syms m mB

% Centres Ja, Jb, Jc des cercles exinscrits
ja=4*s3/((u+v)*(u+w)); jaB=4*s3B/((uB+vB)*(uB+wB));
jb=4*s3/((v+w)*(v+u)); jbB=4*s3B/((vB+wB)*(vB+uB));
jc=4*s3/((w+u)*(w+v)); jcB=4*s3B/((wB+uB)*(wB+vB));

[h hB]=ProjectionPointDroite(a,m,b,aB,mB,bB); % Point H
[k kB]=ProjectionPointDroite(a,m,c,aB,mB,cB); % Point K
HK2=Factor((h-k)*(hB-kB)) % HK^2
f(m,mB)=HK2;

F=-(v-w)^2*((v+w)*mB-2)*(m-mB*u^2-2);
FB=-(v-w)^2*((v+w)*m-2*v*w)*(m+mB*u^2-2*u);
D=numden(Factor(diff(HK2,m)/F));
DB=numden(Factor(diff(HK2,mB)/FB));

Eq=Factor(resultant(D,DB,mB))

X=m*(m-b)^2*(m-c)^2*(m-ja)*(m-jb)*(m-jc);
Cte=Factor(Eq/X)
% On trouve Cte=-4*u^2*(u+v)^4*(u+w)^4*(v+w)^3*(u-v)*(u-w)
% Cte ne dépend pas de m
% Les solutions de l'équation Eq=0 sont b et c (solutions doubles)
% et ja, jb, jc les centres des cercles exinscrits

% f(a,aB)=f(b,bB)=f(c,cB)=0
% f(jb,jbB)=-(u-v)^2/(u+v)^2
% f(jc,jcB)=-(u-w)^2/(u+w)^2

H0K02=Factor(f(ja,jaB)) % H0K02=-((u-v)*(u-w)*(v-w)/((u+v)*(u+w)*(v+w)))^2

% D'autre part BC=2*i*w*(u-v)/((u+w)*(v+w)) et permutation circulaire
% On peut donc vérifier que H0K0=(BC+CA+AB)/2:

BC=2*i*w*(u-v)/((u+w)*(v+w));
[CA AB]=PermCirc(BC,u,v,w);
Y=(BC+CA+AB)/2;
Nul=Factor(Y^2-H0K02) % Nul=0