4) ){
$runden=2;
}
} else {
$runden=2;
}
# Anzahl der Interpolationspunkte eingeben:
# (GET number of interpolation points)
$ta=TRUE;
if(isset($_GET['anz'])) {
$tanz=$_GET['anz'];
if (($tanz < 2)OR($tanz > 1000) ){
$tanz=10;
$ta=FALSE;
}
} else {
$tanz=10;
$ta=FALSE;
}
# Parameter ...
if(isset($_GET['x1'])) {
$x1=$_GET['x1'];
if (($x1 < -32767.9999)OR($x1 > 32767.9999) ){
$x1=mt_rand(-300,300);
}
} else {
$x1=mt_rand(-300,300);
}
if(isset($_GET['y1'])) {
$y1=$_GET['y1'];
if (($y1 < -32767.9999)OR($y1 >32767.9999) ){
$y1=mt_rand(-300,300);
}
} else {
$y1=mt_rand(-300,300);
}
if(isset($_GET['x2'])) {
$x2=$_GET['x2'];
if (($x2 < -32767.9999)OR($x2 > 32767.9999) ){
$x2=mt_rand(-300,300);
}
} else {
$x2=mt_rand(-300,300);
}
if(isset($_GET['y2'])) {
$y2=$_GET['y2'];
if (($y2 < -32767.9999)OR($y2 >32767.9999) ){
$y2=mt_rand(-300,300);
}
} else {
$y2=mt_rand(-300,300);
}
if(isset($_GET['t'])) {
$t=round($_GET['t'],0);
if (($t < -32767)OR($y2 >32767) ){
$t=mt_rand(-10,10);
}
} else {
$t=mt_rand(-10,10);
}
if(isset($_GET['cx'])) {
$cx=$_GET['cx'];
if (($cx < -32767.9999)OR($cx > 32767.9999) ){
$cx=mt_rand(-600,600);
}
} else {
$cx=mt_rand(-600,600);
}
if(isset($_GET['cy'])) {
$cy=$_GET['cy'];
if (($cy < -32767.9999)OR($cy > 32767.9999) ){
$cy=mt_rand(-600,600);
}
} else {
$cy=mt_rand(-600,600);
}
# Auswertung beginnen:
# (start evaluation)
$xs1=$x1-$cx;
$ys1=$y1-$cy;
$xs2=$x2-$cx;
$ys2=$y2-$cy;
$ra=sqrt($xs1*$xs1+$ys1*$ys1);
$rd=sqrt($xs2*$xs2+$ys2*$ys2)-$ra;
$ex=$xs1/$ra;
$ey=$ys1/$ra;
if($ey>0) {
$vh=1;
} else {
$vh=-1;
}
$theta=$vh*acos($ex);
$ex=$xs2/($rd+$ra);
$ey=$ys2/($rd+$ra);
if($ey>0) {
$vh=1;
} else {
$vh=-1;
}
$delta=$vh*acos($ex)-$theta+$t*2*M_PI;
#$theta=$thetadeg*M_PI/180;
#$delta=$deltadeg*M_PI/180;
$thetadeg=round($theta/M_PI*180,$runden);
$deltadeg=round($delta/M_PI*180,$runden);
$rao=round($ra,$runden);
$rdo=round($rd,$runden);
if (!$ta) {
$anz=max(3,ceil(abs($deltadeg)/30));
} else {
$anz=$tanz;
}
# viewBox berechnen / calculate viewBox
$rr=abs($ra+$rd);
$rs=abs($ra);
$vv=4*ceil(0.6*max($rs,$rr));
$vh=1.5*$vv;
$vy=round($cy-$vv/1.85);
$vx=round($cx-$vh/2);
$viewBox="$vx $vy $vh $vv";
$sw=ceil($vv/100);
$swr=$sw/4;
#################################
# Vergleich, doppelte Punktzahl / comparison, doubled number of points
# Datensatz erzeugen / Interpolation
$anzv=2*$anz;
$dphi=$delta/$anzv;
$dr=$rd/$anzv;
for ($j = 0; $j <= $anzv; $j++) {
$tt=$dphi*$j;
$cpi=cos($tt+$theta);
$spi=sin($tt+$theta);
$ri=$ra+$dr*$j;
# Punkte / points
$xx=$ri*$cpi;
$yy=$ri*$spi;
# Ableitungen / derivatives ...
$dx=$dr*$cpi -$dphi*$ri*$spi;
$dy=$dr*$spi +$dphi*$ri*$cpi;
# Daten speichern / save data
$data[$j][0]=$xx+$cx;
$data[$j][1]=$yy+$cy;
$data[$j][2]=$dx;
$data[$j][3]=$dy;
}
# Aus den Ableitungen die Kontrollpunkte berechnen
# (generate control points from derivatives)
$cp[0][0]=round($data[0][0]+$data[0][2]/3,$runden);
$cp[0][1]=round($data[0][1]+$data[0][3]/3,$runden);
for ($j = 0; $j <= $anzv; $j++) {
$cq[$j][0]=round($data[$j][0]-$data[$j][2]/3,$runden);
$cq[$j][1]=round($data[$j][1]-$data[$j][3]/3,$runden);
# und nun erst Punkte runden / and now round points
$data[$j][0]=round($data[$j][0],$runden);
$data[$j][1]=round($data[$j][1],$runden);
}
# Daraus Pfad bestimmen und malen der Punkte vorbereiten
# (generate path data and prepare for painting)
$pfadv='M'.$data[0][0].' '.$data[0][1].' C'.$cp[0][0].' '.$cp[0][1].' '.$cq[1][0].' '.$cq[1][1].' '.$data[1][0].' '.$data[1][1]." S\n";
for ($j = 2; $j <=$anzv; $j++) {
$pfadv.=$cq[$j][0].' '.$cq[$j][1].' '.$data[$j][0].' '.$data[$j][1]."\n";
}
###################################
# Datensatz erzeugen / Interpolation
$dphi=$delta/$anz;
$dr=$rd/$anz;
for ($j = 0; $j <= $anz; $j++) {
$tt=$dphi*$j;
$cpi=cos($tt+$theta);
$spi=sin($tt+$theta);
$ri=$ra+$dr*$j;
# Punkte / points
$xx=$ri*$cpi;
$yy=$ri*$spi;
# Ableitungen / derivatives ...
$dx=$dr*$cpi -$dphi*$ri*$spi;
$dy=$dr*$spi +$dphi*$ri*$cpi;
# Daten speichern / save data
$data[$j][0]=$xx+$cx;
$data[$j][1]=$yy+$cy;
$data[$j][2]=$dx;
$data[$j][3]=$dy;
}
# Aus den Ableitungen die Kontrollpunkte berechnen
# (generate control points from derivatives)
for ($j = 0; $j <= $anz; $j++) {
$cp[$j][0]=round($data[$j][0]+$data[$j][2]/3,$runden);
$cp[$j][1]=round($data[$j][1]+$data[$j][3]/3,$runden);
$cq[$j][0]=round($data[$j][0]-$data[$j][2]/3,$runden);
$cq[$j][1]=round($data[$j][1]-$data[$j][3]/3,$runden);
# und nun erst Punkte runden / and now round points
$data[$j][0]=round($data[$j][0],$runden);
$data[$j][1]=round($data[$j][1],$runden);
}
# Fehler abschätzen / estimate error
$fehl=0;
$fmax=0;
$fanz=0;
for ($j = 0; $j <= $anz; $j++) {
$jj=$j+1;
if ($j != $anz) {
for ($i = 0; $i < 30; $i++) {
$tf=$dphi*($j+$i/30);
$cpi=cos($tf+$theta);
$spi=sin($tf+$theta);
$ri=$ra+$dr*($j+$i/30);
# Punkte / points
$xx=$ri*$cpi+$cx;
$yy=$ri*$spi+$cy;
$u=$i/30;
$ur=(1-$u);
$x1=$ur*$ur*$ur*$data[$j][0]+3*$u*$ur*$ur*$cp[$j][0] +3*$u*$u*$ur*$cq[$jj][0]+ $u*$u*$u*$data[$jj][0];
$y1=$ur*$ur*$ur*$data[$j][1]+ 3*$u*$ur*$ur*$cp[$j][1] + 3*$u*$u*$ur*$cq[$jj][1]+ $u*$u*$u*$data[$jj][1];
$dxx=$x1-$xx;
$dyy=$y1-$yy;
$fdd=sqrt($dxx*$dxx+$dyy*$dyy);
$fehl=$fehl+$fdd;
if ($fdd > $fmax) {
$fmax=$fdd;
$fx=$xx;
$fy=$yy;
}
$fanz=$fanz+1;
}
}
}
$fehl=$fehl/$fanz;
# Daraus Pfad bestimmen und malen der Punkte vorbereiten
# (generate path data and prepare for painting)
$pfada='M'.$data[0][0].' '.$data[0][1].' C'.$cp[0][0].' '.$cp[0][1].' '.$cq[1][0].' '.$cq[1][1].' '.$data[1][0].' '.$data[1][1]." S\n";
$pl=' \n";
for ($j = 2; $j <=$anz; $j++) {
$pfada.=$cq[$j][0].' '.$cq[$j][1].' '.$data[$j][0].' '.$data[$j][1]."\n";
$pl.='\n";
}
###################################
$aus="
Zentrum / centre
maximaler Fehler / maximal error
Pfadvergleich doppelte Punktanzahl / path comparsion, twice the number of points
Approximierter Pfad / approximated path
Punkte / points
$pl
Anfangspunkt / initial point
Endpunkt / final point
";
# Text skalieren und formatieren:
# (scale and format text)
$fs=ceil($vv/40);
$lh=ceil(1.2*$vv/40);
$textx=$vx +4;
$texty=$vy +2+$lh;
$texty1=$vy +2+2*$lh;
$texty2=$vy +2+3*$lh;
$texty3=$vy +2+4*$lh;
$break="\n";
$textpfad=wordwrap($pfada,100,$break, false);
$textpfad="$textpfad";
$aus.="
?cx=$cx&cy=$cy&x1=$x1&y1=$y1&x2=$x2&y2=$y2&t=$t&anz=$anz&rnd=$runden
Aktuelle Werte / Current Values
$fmax | $fehl
Fehler maximal, mittel / Error maximum, average
cx=$cx; cy=$cy; ra=$rao; rd=$rdo; theta=$thetadeg; delta=$deltadeg
Aktuelle Werte / Current Values
$textpfad
Kubische Näherung / Cubic Approximation
";
# svg-header senden:
$content="Content-type: image/svg+xml; charset=iso-8859-1";
header($content);
# xml-Zeile ausgeben
echo "";
# und jetzt das Dokument
?>