4) ){
$runden=2;
}
} else {
$runden=2;
}
# Parameter ...
if(isset($_GET['theta'])) {
$thetadeg=$_GET['theta'];
if (($thetadeg < -32767.9999)OR($thetadeg > 32767.9999) ){
$thetadeg=mt_rand(-4000,4000);
}
} else {
$thetadeg=mt_rand(-4000,4000);
}
if ($thetadeg <0)
{
$sigp=-1;
} else {
$sigp=1;
}
$theta=abs($thetadeg)*M_PI/180;
if(isset($_GET['phi0'])) {
$phideg=$_GET['phi0'];
if (($phideg <= 0)OR($phideg > 32767.9999) ){
$phideg=mt_rand(1,100)/10;
}
} else {
$phideg=mt_rand(1,100)/10;
}
$phi=$phideg*M_PI/180;
# Anzahl der Interpolationspunkte eingeben:
# (GET number of interpolation points)
if(isset($_GET['anz'])) {
$anz=$_GET['anz'];
if (($anz < 2)OR($anz > 1000) ){
$anz=max(3,ceil(abs($thetadeg-$phideg)/30));
}
} else {
$anz=max(3,ceil(abs($thetadeg-$phideg)/30));
}
$anzv=2*$anz;
if(isset($_GET['ra'])) {
$radeg=$_GET['ra'];
if (($radeg < -32767.9999)OR($radeg > 32767.9999) ){
$radeg=mt_rand(-30000,30000);
}
} else {
$radeg=mt_rand(-30000,30000);
}
if ($radeg <0)
{
$sigr=-1;
} else {
$sigr=1;
}
$sig=$sigp*$sigr;
$ra=abs($radeg)/180*M_PI;
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)
# viewBox berechnen / calculate viewBox
$rr=sqrt($ra/$phi);
$rs=$rr;
$vv=4*ceil(0.35*$rs);
$vh=1.5*$vv;
$vy=round($cy-$vv/1.85);
$vx=round($cx-$vh/2);
$viewBox="$vx $vy $vh $vv";
$sw=round($vv/100,$runden);
$swr=$sw/4;
#################################
# Vergleich, doppelte Punktzahl, comparison, doubled number of points
# Datensatz erzeugen / Interpolation
$eps=$theta-$phi;
for ($j = 0; $j <= $anzv; $j++) {
$u=$j/$anzv;
$uh=$u*$u;
$tt=$eps*$uh+$phi;
$dphi=2*$eps*$u/$anzv;
$tp=pow($tt,1.5);
$cpi=cos($tt);
$spi=sin($tt);
$ri=$sig*sqrt($ra/$tt);
$dr=-$sig*0.5*sqrt($ra)/$tp;
# Punkte / points
$xx=$ri*$cpi;
$yy=$ri*$spi;
# Ableitungen / derivatives ...
$dx=$dr*$dphi*$cpi -$dphi*$ri*$spi;
$dy=$dr*$dphi*$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=$theta/$anz;
$eps=$theta-$phi;
for ($j = 0; $j <= $anz; $j++) {
$u=$j/$anz;
$uh=$u*$u;
$tt=$eps*$uh+$phi;
$dphi=2*$eps*$u/$anz;
$cpi=cos($tt);
$spi=sin($tt);
$ri=$sig*sqrt($ra/$tt);
$dr=-$sig*0.5*sqrt($ra)*pow($tt,-1.5);
# Punkte / points
$xx=$ri*$cpi;
$yy=$ri*$spi;
# Ableitungen / derivatives ...
$dx=$dr*$dphi*$cpi -$dphi*$ri*$spi;
$dy=$dr*$dphi*$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++) {
$u=($j+$i/30)/$anz;
$uh=$u*$u;
$tf=$eps*$uh+$phi;
$dphi=2*$eps*$u/$anz;
$cpi=cos($tf);
$spi=sin($tf);
$ri=$sig*sqrt($ra/$tf);
# 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
";
# Text skalieren und formatieren:
# (scale and format text)
$fs=round($vv/40,$runden);
$lh=round(1.2*$vv/40,$runden);
$textx=$vx +$fs;
$texty=$vy +$lh;
$texty1=$vy +2*$lh;
$texty2=$vy +3*$lh;
$texty3=$vy +4*$lh;
$break="\n";
$textpfad=wordwrap($pfada,100,$break, false);
$textpfad="$textpfad";
$aus.="
?cx=$cx&cy=$cy&ra=$radeg&theta=$thetadeg&phi0=$phideg&anz=$anz&rnd=$runden
Aktuelle Werte / Current Values
$fmax | $fehl
Fehler maximal, mittel / Error maximum, average
$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
?>