function?P=distanceconstraint（PX0controlpoint）
for?i=1:X0.face_num
????idx1=X0.tri（1i）;??%三角形的三個索引
????idx2=X0.tri（2i）;
????idx3=X0.tri（3i）;
????d1=norm（X0.xyz（:idx1）-X0.xyz（:idx2））;??%原始兩點間的距離
????c_d1=norm（P（idx1:）-P（idx2:））;%現在兩點間的距離
????d2=norm（X0.xyz（:idx2）-X0.xyz（:idx3））;
????c_d2=norm（P（idx2:）-P（idx3:））;
????d3=norm（X0.xyz（:idx1）-X0.xyz（:idx3））;
????c_d3=norm（P（idx1:）-P（idx3:））;
????if?c_d1~=d1??%考察三角形第一條邊，?如果現在的距離與原始的距離不相等，那么更新位置，就像彈簧一樣
????????test=1;%用來測試這個點是不是control?point，?如果是control?point?那么不更新位置
????????test2=1;%同上
????????for?count=1:length（controlpoint）
????????????if?idx1==controlpoint（count）
????????????????test=0;??%這個點不是control?point
????????????end
????????????if?idx2==controlpoint（count）
????????????????test2=0;
????????????end
????????end
????????if?test==0
????????????deltap1=[000];
????????????deltap2=（c_d1-d1）*（P（idx1:）-P（idx2:））/（norm（P（idx1:）-P（idx2:）））;
????????else
????????????deltap1=-0.5*（c_d1-d1）*（P（idx1:）-P（idx2:））/（norm（P（idx1:）-P（idx2:）））;?%文章中的stifness是乘在這里的，我沒乘，我按質點質量相等來算的，所以是0.5
????????end
????????if?test2==0
????????????deltap2=[000];
????????????deltap1=-（c_d1-d1）*（P（idx1:）-P（idx2:））/（norm（P（idx1:）-P（idx2:）））;
????????else
????????????deltap2=0.5*（c_d1-d1）*（P（idx1:）-P（idx2:））/（norm（P（idx1:）-P（idx2:）））;
????????end
????????P（idx1:）=P（idx1:）+deltap1;
????????P（idx2:）=P（idx2:）+deltap2;
????end
????if?c_d2~=d2??%考察第二條邊
????????test=1;
????????test2=1;
????????for?count=1:length（controlpoint）
????????????if?idx2==controlpoint（count）
????????????????test=0;
????????????end
????????????if?idx3==controlpoint（count）
????????????????test2=0;
????????????end
????????end
????????if?test==0
????????????deltap23=[000];
????????????deltap32=（c_d2-d2）*（P（idx2:）-P（idx3:））/（norm（P（idx2:）-P（idx3:）））;
????????else
????????????deltap23=-0.5*（c_d2-d2）*（P（idx2:）-P（idx3:））/（norm（P（idx2:）-P（idx3:）））;
????????end
????????if?test2==0
????????????deltap32=[000];
????????????deltap23=-（c_d2-d2）*（P（idx2:）-P（idx3:））/（norm（P（idx2:）-P（idx3:）））;
????????else
????????????deltap32=0.5*（c_d2-d2）*（P（idx2:）-P（idx3:））/（norm（P（idx2:）-P（idx3:）））;
????????end
????????P（idx2:）=P（idx2:）+deltap23;
????????P（idx3:）=P（idx3:）+deltap32;
????end
????if?c_d3~=d3?%考察第三條邊
????????test=1;
????????test2=1;
????????for?count=1:length（controlpoint）
????????????if?idx1==controlpoint（count）
????????????????test=0;
????????????end
????????????if?idx3==controlpoint（count）
????????????????test2=0;
????????????end
????????end
????????if?test==0
????????????deltap13=[000];
????????????deltap31=（c_d3-d3）*（P（idx1:）-P（idx3:））/（norm（P（idx1:）-P（idx3:）））;
????????else
????????????deltap13=-0.5*（c_d3-d3）*（P（idx1:）-P（idx3:））/（norm（P（idx1:）-P（idx3:）））;
????????end
????????if?test2==0
????????????deltap31=[000];
????????????deltap13=-（c_d3-d3）*（P（idx1:）-P（idx3:））/（norm（P（idx1:）-P（idx3:）））;
????????else
????????????de

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????33041??2015-01-27?21:49??PBD\cloth.obj

?????文件???????3325??2015-01-27?21:41??PBD\distanceconstraint.m

?????文件???????1711??2016-07-05?00:18??PBD\PBD.asv

?????文件???????1762??2016-07-05?00:33??PBD\PBD.m

?????文件??????30306??2013-03-29?16:34??PBD\read_obj.m

?????文件???????4378??2013-03-29?18:07??PBD\write_obj.m

?????目錄??????????0??2016-07-05?00:18??PBD

-----------?---------??----------?-----??----

????????????????74523????????????????????7