!Interpolation der Temperatur für einen Ort innerhalb eines Volumens !Schätzing, 29.4.2014 !==========Interpolation=============================== nsel,all !----Ort Interpolationspunkt--- x0=30e-3 y0=25e-3 z0= 9e-3 !----Knoten-in-der-Nähe-Interpolationspunkt- n1=node(x0,y0,z0) !----5-Knoten-in-der-Nähe----- !----(die-Anzahl-5-hier-beispielhaft)--- n2=nnear(n1) nsel,u,,,n2 n3=nnear(n1) nsel,u,,,n3 n4=nnear(n1) nsel,u,,,n4 n5=nnear(n1) nsel,all !----Entfernungen-Ort-Knoten--- e01=sqrt((nx(n1)-x0)**2+(ny(n1)-y0)**2+(nz(n1)-z0)**2) e02=sqrt((nx(n2)-x0)**2+(ny(n2)-y0)**2+(nz(n2)-z0)**2) e03=sqrt((nx(n3)-x0)**2+(ny(n3)-y0)**2+(nz(n3)-z0)**2) e04=sqrt((nx(n4)-x0)**2+(ny(n4)-y0)**2+(nz(n4)-z0)**2) e05=sqrt((nx(n5)-x0)**2+(ny(n5)-y0)**2+(nz(n5)-z0)**2) !----Temperaturen-an-den-Knoten--------- t1=temp(n1) t2=temp(n2) t3=temp(n3) t4=temp(n4) t5=temp(n5) !---mittlere-Temperatur----- tm=(t1+t2+t3+t4+t5)/5 !---entfernungsgewichtete-Temperatur----- !---Summe-aller-Entfernungen------ se=e01+e02+e03+e04+e05 !---relative-reziproke-Entfernungen-(Wichtung)----- e01r=se/e01 e02r=se/e02 e03r=se/e03 e04r=se/e04 e05r=se/e05 ser=e01r+e02r+e03r+e04r+e05r !Summe !---entfernungsgewichtete-Temperatur----- t0=(t1*e01r+t2*e02r+t3*e03r+t4*e04r+t5*e05r)/ser