如题,我得到了网格数据p,e,t之后,find了我需要的p的坐标,记为p1,然后我现在可以怎么得到现在的p1所对应的e1,t1呢?求教各路大神。
现在我写的部分代码如下图,得到了p1,之后需要e1和t1,有办法吗?
function u1=innerparabolic(tlist,O,radius,p,e,t,gamma,b)
Ovec=ones(2,size(p,2));
Ovec(1,:)=O(1)*Ovec(1,:);
Ovec(2,:)=O(2)*Ovec(2,:);
ix=find(sqrt((p(1,:)-Ovec(1,:)).^2+(p(2,:)-Ovec(2,:)).^2)<radius);%控制节点
mm=size(ix);
p1=p(:,ix);