Re: 找蔣公...還是MATLAB的問題
我手邊沒有matlab的電腦
這樣改不知道行不行
重點在於 n1 n2 n3
要交錯著跑
n1=0:30
n2=0:30
n3=0:30
這樣子
nn是決定要把0~30切幾個點
clear
clc
N1=[];N2=[];N3=[];
nn=1000
for i=1:nn+1
n1=30*(i-1)/nn;
for j=1:nn+1
n2=30*(j-1)/nn;
for k=1:nn+1
n3=30*(k-1)/nn;
a1=(4*(15.3-n1)^2+(0-n2)^2+(0-n3)^2)^0.5 ;
a2=(4*(17.8-n1)^2+(1-n2)^2+(3.1-n3)^2)^0.5;
a3=(4*(17.8-n1)^2+(0-n2)^2+(2.7-n3)^2)^0.5;
a4=(4*(15.8-n1)^2+(5.3-n2)^2+(7.2-n3)^2)^0.5;
a5=(4*(18-n1)^2+(1.4-n2)^2+(2-n3)^2)^0.5;
a6=(4*(15.2-n1)^2+(5.6-n2)^2+(9.1-n3)^2)^0.5;
a7=(4*(18.4-n1)^2+(0-n2)^2+(2-n3)^2)^0.5;
a8=(4*(16-n1)^2+(9-n2)^2+(5.1-n3)^2)^0.5;
a9=(4*(17.8-n1)^2+(3.1-n2)^2+(5.7-n3)^2)^0.5;
a10=(4*(16.8-n1)^2+(5.7-n2)^2+(8-n3)^2)^0.5;
a11=(4*(16.6-n1)^2+(7.9-n2)^2+(6.5-n3)^2)^0.5;
a12=(4*(15.5-n1)^2+(10.4-n2)^2+(7-n3)^2)^0.5;
a13=(4*(19-n1)^2+(1.8-n2)^2+(7.4-n3)^2)^0.5;
a14=(4*(20-n1)^2+(8.6-n2)^2+(4.1-n3)^2)^0.5;
a15=(4*(16-n1)^2+(5.7-n2)^2+(15.8-n3)^2)^0.5;
a16=(4*(15.8-n1)^2+(6.1-n2)^2+(16.4-n3)^2)^0.5;
a17=(4*(18.4-n1)^2+(6.3-n2)^2+(13.7-n3)^2)^0.5;
a18=(4*(15.3-n1)^2+(18-n2)^2+(6.1-n3)^2)^0.5;
a19=(4*(17.4-n1)^2+(13.7-n2)^2+(11.3-n3)^2)^0.5;
a20=(4*(15.8-n1)^2+(8.8-n2)^2+(19.4-n3)^2)^0.5;
a21=(4*(18.4-n1)^2+(16.4-n2)^2+(10.2-n3)^2)^0.5;
a22=(4*(15.1-n1)^2+(12.3-n2)^2+(22.3-n3)^2)^0.5;
a23=(4*(15.6-n1)^2+(16-n2)^2+(42.3-n3)^2)^0.5;
B=[a23,a10,a15,a16,a17,a19,a20,a21,a22]; %bad solvent
C=[a1,a2,a3,a4,a5,a6,a7,a8,a9,a11,a12,a13,a14,a18]; % good solvet
D=max(C);
E=min(B);
if E>D
n1=n1
n2=n2
n3=n3
N1=[N1 n1];
N2=[N2 n2];
N3=[N3 n3];
N3=[N3 n3];
disp('one ans')
disp('=======================================')
end
end
end
end
N=[N1' N2' N3'];
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.120.99
討論串 (同標題文章)
完整討論串 (本文為第 2 之 2 篇)：